Main

Organic crystals that demonstrate an intermolecular distance at the equilibrium point (3–5 Å; Fig. 1(i)) are one of the most stable aggregated phases. In this equilibrium state, organic molecules generally exhibit strong intermolecular interactions (for example, ππ stacking, energy transfer, charge transfer and self-absorption), which result in a notable drop in their photoluminescence quantum yield (PLQY)11. This phenomenon is widely known as aggregation-caused quenching (ACQ), which is a substantial restriction for the practical use of aggregated emitters. An effective strategy to alleviate ACQ was to weaken or eliminate intermolecular interactions by spatially disengaging the emitters far away from each other (>10 Å; Fig. 1(ii)), such as by dissolving them in solvents, diluting them in a foreigner matrix12,13,14, embedding them into porous hosts15,16,17 or growing macrocycle-hosts-based supramolecular single crystals18,19. These well-isolated single molecules have allowed efficient emission and set the foundation for organic light-emitting diodes2,20, lasers and many other optoelectronic applications. Nonetheless, scarce reports exist on organic emitters with an intermolecular distance between the equilibrium and dilution states and their behaviours have received limited research attention.

Fig. 1: Schematic illustration for different scenarios of molecular emitters with varying intermolecular distances.
figure 1

(i) Equilibrium state. The intermolecular distance in organic crystals is typically within 3–5 Å, which enables the lowest potential energy. This can be considered as the equilibrium state with aggregate behaviour. (ii) Dilution state. The molecules are separated far away from each other (>10 Å), leading to a dilution state with single-molecule behaviour. (iii) Near-equilibrium state. Organic emitters, confined in organic–inorganic 2D SLs with an intermolecular distance near the equilibrium state, behave as single molecules but also deliver some aggregate properties, thus being termed SMAs.

Two-dimensional (2D) layered hybrid perovskites are emerging solution-processable semiconducting materials, combining the advantages of both organic and inorganic components21,22,23,24,25. These unique 2D organic–inorganic superlattices (SLs) offer a platform to investigate the behaviours of organic emitters near the equilibrium state because the inorganic sheets intrinsically demonstrate a square lattice with an approximately 6-Å pitch to house organic molecules (Fig. 1(iii)), which is potentially sufficient to modulate or suppress intermolecular interactions. Numerous studies have been focusing on how the incorporation of organic moieties improves the luminescent efficiency26, charge-transport capability27 and stability22,27,28 of inorganic slabs, which have led to many breakthroughs in high-performance perovskite electronic and optoelectronic devices29,30,31,32,33. However, making use of the inorganic sublattice to tune the intermolecular interactions, molecular packing and emission properties of organic molecules remains largely unexplored. Starting from the late 1990s, several groups have reported the formation of organic semiconductor–perovskite SLs and confirmed that the emitting species could be organic chromophores22,34,35,36,37. However, the range of organic molecular emitters that can be incorporated into layered perovskites is rather limited and their PLQYs are often low (typically below 10%).

In this study, we demonstrate a new phase of molecular aggregates, which we refer to as a SMA, that was achieved near the equilibrium state by combining a 2D inorganic sublattice with suitably tailored organic chromophores. In this hybrid SL, the behaviour of organic emitters closely resembles that of individual single molecules, as evidenced by their similar emission wavelengths and lifetimes, as well as near-unity PLQYs. Theoretical and experimental investigations highlight the pivotal role of backbone dihedrals of the organic emitter in maintaining this single-molecule behaviour. Notably, despite exhibiting single-molecule-like properties, the strong alignment and close packing of organic molecules within the 2D SLs resulted in intense directional emission, ultrafast radiative recombination and efficient lasing, which are linked to the traits of ordered molecular ensembles or aggregates.

SMA materials design

We designed and synthesized two new organic molecular emitters, 2-(5-(7-(9,9-dimethyl-9H-fluoren-2-yl)benzo[c][1,2,5]thiadiazol-4-yl)thiophen-2-yl)ethan-1-ammonium bromide (referred to as FBTT; Fig. 2a) and 2-(4-(7-(9,9-dimethyl-9H-fluoren-2-yl)benzo[c][1,2,5]thiadiazol-4-yl)phenyl)ethan-1-ammonium bromide (referred to as FBTP; Fig. 2e), and incorporated them into the 2D perovskite lattices. The synthetic procedures of the organic molecules are detailed in the Supplementary Information (Supplementary Figs. 14). The chromophore core was adapted from the repeating unit of a bright-green-emitting polymer, poly(9,9-dioctylfluorene-alt-benzothiadiazole) (F8BT)38. The alkyl chain length was cut to one carbon to ensure that they can be accommodated in the perovskite lattice. Linkers, such as thienyl and phenyl groups, were further introduced to connect the chromophore core with the ammonium tail and then ionically bonded with the lead halide matrix.

Fig. 2: Single-molecule behaviour of molecular emitters in perovskite 2D SLs.
figure 2

a, Molecular structure of FBTT. b, GIWAXS pattern of FBTT 2D SL thin film. c,d, PL spectra and PLQY (c) and FLIM mapping (d) of FBTT monomer in PMMA, FBTT 2D SLs and FBTT aggregate films. e, Molecular structure of FBTP. f, GIWAXS pattern of FBTP 2D SL thin film. g,h, PL spectra and PLQY (g) and FLIM mapping (h) of FBTP monomer in PMMA, FBTP 2D SLs and FBTP aggregate films. Note that the colour scales in b and f are adjusted differently to effectively highlight the scattering along the qxy direction; the dashed white rings highlight the scattering pattern at a typical q near 0.5 Å1 with their intensity shown in Extended Data Fig. 2; the solid white lines are blank spaces derived from the raw GIWAXS data, which are collected by splicing the detector images; and the white arrows indicate the (100) diffraction and confirm the 2D perovskite structure. The dashed violet curves in c and g are the PL spectra of (PEA)2PbBr4 (PEA is phenylethylammonium), which implies the emission feature of the 2D lead bromide lattice, and the grey dash-dot lines indicate the PL peaks of the 2D SLs compared with those of the monomers and aggregates. The insets are PL images of the corresponding thin-film samples and all the image scales are 10 μm × 10 μm. The scales of all FLIM images in d and h are 50 μm × 50 μm.

Source Data

Lead bromide was selected as the inorganic framework to construct a reversed type-I band alignment for organic compounds to emit efficiently39, which is confirmed by density functional theory (DFT) and time-dependent density functional theory (TDDFT) simulations (Extended Data Fig. 1). We then carried out thin-film studies on the aggregates (neat organic emitters), 2D SLs (organic-emitter-incorporated perovskites) and isolated monomers (2 wt% organic emitters doped poly(methyl methacrylate) (PMMA) films). Grazing-incidence wide-angle X-ray scattering (GIWAXS) and X-ray diffraction (XRD) characterizations helped us verify that both FBTT and FBTP were successfully incorporated, forming a layered perovskite structure (Fig. 2b,f, Extended Data Fig. 2 and Supplementary Fig. 5). The sharp ultraviolet–visible (UV–vis) peak around 400 nm and the nearby shoulder peak from 2D SL films (Supplementary Fig. 5) can be indexed to the excitonic peak of perovskites and organic emitters, respectively, further supporting the formation of 2D perovskites.

The PL spectrum of FBTT SLs (peaked at about 594 nm; Fig. 2c) is similar to that of its aggregates and redshifted compared with its monomer state (PL peak at 559 nm). This redshift is usually ascribed to the intermolecular interaction, especially ππ stacking, induced bandgap narrowing. Notably, by contrast, the PL spectrum of FBTP SLs (peaked at about 540 nm; Fig. 2g) is similar to the FBTP monomers and clearly blueshifted compared with the aggregates (peaked at 553 nm). The inset PL images in Fig. 2c,g further support the above results. The absence of perovskite emission (violet profiles in Fig. 2c,g) in both FBTT and FBTP SLs substantiates the fact that both 2D SLs are subject to a reversed type-I band alignment. Moreover, the FBTT SLs and aggregates seem to have a similar PLQY of about 15%, which is much lower than its monomer (81.2%). This can be explained by the notorious ACQ effect. Surprisingly, the ACQ effect is absent in the FBTP SLs sample, which features a much higher PLQY (92.8%) than that of its aggregates (42.4%). Such a high PLQY is very close to that of its monomers (PLQY = 95.7%). To the best of our knowledge, this is the highest PLQY reported so far for a 2D perovskite material22,26,35. These results indicate that the FBTT molecules in perovskites behave more like aggregates, whereas the FBTP molecules in perovskites behave like single molecules.

Fluorescence lifetime imaging microscopy (FLIM) measurements show that FBTT SLs manifest similar PL lifetime as its aggregates (Fig. 2d), which is much shorter than that of the FBTT monomer, again validating the aggregation behaviour of FBTT in the perovskite SLs. The PL lifetime of FBTP SLs matches well with its monomer case (Fig. 2h). It should be noted that the pitch size of the inorganic square lattice can be fine-tuned by introducing different halides, such as Cl or I, which would enable different behaviours of organic emitters. After mixing with I, the inorganic sublattice will slightly expand and the inter-emitter distance will increase towards the dilution state. Consequently, the single-molecule behaviour is maintained, as evidenced by an identical PL emission peak (solid blue curve in Extended Data Fig. 3) to that of single molecules and the pure bromide case, whereas the incorporation of Cl will constrain the inorganic sublattice and the organic emitter will be brought even closer towards the equilibrium state, thus leading to a similar PL emission as that of aggregates (solid brown and red curves in Extended Data Fig. 3). Furthermore, the FBTT and FBTP yielded 2D SL thin films exhibiting similar surface morphology and thickness (Extended Data Fig. 4, Supplementary Fig. 6 and Methods). Also, the emission spectra of single-crystalline 2D perovskite crystals derived from FBTT and FBTP resemble those of their corresponding polycrystalline thin films (Extended Data Fig. 5). These observations eliminate film quality and crystallinity as confounding factors governing the differential emission properties.

Molecular insights

To understand the relationship between the emitter structures and their behaviours in the lattice, we performed equilibrium molecular dynamics (MD) simulations on individual FBTT-based and FBTP-based layered perovskites, following equilibration and statistical sampling over finite-temperature equilibrium structures (see Methods for simulation details). As shown in Fig. 3a,b, FBTT adopts a more ordered packing style than FBTP (side view in Supplementary Fig. 7). Further investigations on the configuration of each ligand reveal that FBTT ligands are more planar than FBTP, characterized by the fact that the benzothiadiazole and thienyl rings in FBTT ligands are parallel (Fig. 3c), whereas the benzothiadiazole and phenyl rings in FBTP are not (Fig. 3d). Co-planar molecules usually experience stronger intermolecular interactions, such as ππ stacking, that promote emission quenching and wavelength redshifting, agreeing well with the experimental observations in FBTT.

Fig. 3: Molecular insights into the rotational and vibrational behaviours of organic molecules confined in 2D SLs.
figure 3

a,b, Side view of thienyl-based (a) and phenyl-based (b) 2D SLs during MD simulations. c,d, Top view of thienyl-based (c) and phenyl-based (d) 2D SLs during MD simulations. Fluorene units are omitted for clarity. e, Probability distribution function (PDF) of each ligand’s backbone dihedral torsion between benzothiadiazole and the respective aromatic ring (Ar) in the perovskite lattice. f, The energy PDF of each ligand in the perovskite lattice. g, PDF of each ligand’s molecular planarity parameter (MPP) in the perovskite lattice. h, Free energy for neighbouring ligand rotation in the perovskite lattice. i, Zoomed-in top-view images of FBTT (left) and FBTP (right) molecules extracted from their corresponding 2D SL single-crystal structures. Hydrogen atoms are omitted for clarity. In ad,i, blue lines denote benzothiadiazole (BT); green lines denote thienyl (Th) rings and phenyl (Ph) rings; purple denotes fluorene (Fl) units; blue spheres denote nitrogen, yellow spheres sulfur, pink and purple spheres bromine, orange spheres lead, and grey spheres carbon. j,k, FTIR of FBTT aggregates and 2D SL films (j) and FBTP aggregates and 2D SL films (k), in which the peak between 1,450 cm−1 and 1,370 cm−1 can be assigned to the bending mode of the aliphatic CH2 group (α(CH2)) and the multipeaks in the range 1,620–1,540 cm−1 could come from the C=C stretching of aromatic linkers.

Source Data

To quantitatively understand the configuration and intermolecular interactions of the ligands, the distribution of the inter-ring dihedral, the site energy (calculated as a sum over all intermolecular terms involving the ligands) and the molecular planarity parameter (calculated as the root-mean-square deviation of atoms from the fitting plane to the ligand) of each ligand were calculated. All distribution data were collected from each frame of the MD trajectory based on different ligands. The simulations on the dihedral between fluorene and benzothiadiazole give similar distribution profiles for both FBTT and FBTP (Supplementary Fig. 8), whereas the dihedral torsions between benzothiadiazole and phenyl/thienyl rings are very different. As shown in Fig. 3e, FBTT manifests a narrow dihedral distribution converging around 0°, whereas the distribution of FBTP spreads out and peaked at approximately 45° and 135°, which is consistent with the DFT energy calculation for the backbone dihedral torsion of the ligands (Supplementary Fig. 9). The energy (Fig. 3f) and molecular planarity parameter (Fig. 3g) distributions show similar results, in which FBTT has a narrower distribution, suggesting that FBTT is more planar and has a more ordered packing than FBTP (see the span of deviation from plane in Supplementary Fig. 10 for more information). These match our observations in the MD simulation in Fig. 3c,d. To examine the packing behaviour of each ligand, we conducted a free-energy analysis with respect to ligand rotation in the lattice (Fig. 3h). The reaction coordinate of the free energy is defined as the packing angle between the ligand being rotated and the reference ligand (Extended Data Fig. 6). It clearly shows that FBTT requires more energy than FBTP to break from parallel stacking (angle = 0°). Furthermore, the FBTP case has a plateau around 50–70°, indicating that FBTP ligands are prone to reorganize themselves to a non-parallel packing pattern in the lattice (Fig. 3h). All of the above analyses imply that FBTT is more planar and tends to have stronger intermolecular interactions, whereas FBTP is less planar and tends to behave like a single molecule in the lattice. These molecular-scale differences reveal the subtle aspects of ligand design that conduce the formation of SMAs on designated lattice incorporation.

To experimentally examine the molecular configuration and packing behaviours of organic emitters in the inorganic lattices, we resolved their single-crystal structures (Supplementary Table 1). Supplementary Figs. 11 and 12 show the single-crystal structures of FBTT-contained and FBTP-contained hybrid materials; respectively. The FBTT molecules in the perovskite lattice adopt a herringbone packing (Extended Data Fig. 7), which is a common motif for organic semiconductors. The torsion angle between the benzothiadiazole and thienyl rings is roughly 2° (Fig. 3i, left). By strong contrast, the packing motif of organic FBTP molecules in the perovskite lattice is dominated by two different stacking styles, criss-cross and herringbone (Extended Data Fig. 7). This unique stacking style of organic FBTP molecules is radically different from the normally observed herringbone or lamellar 2D ππ stacking (Supplementary Fig. 13) in organic-semiconductor crystals40 and organic-molecule-incorporated 2D perovskite structures22,27,41. Furthermore, the torsional angle between the benzothiadiazole and phenyl rings is found to be approximately 45° (Fig. 3i, right). All of the experimental outputs from single-crystal XRD results demonstrate that FBTT is more planar and has a more ordered packing than FBTP, which is in excellent agreement with our MD simulation predictions.

We performed Fourier-transform infrared (FTIR) spectroscopy to inspect the footprint of molecular structures and intermolecular interactions. Infrared absorption from 3,100 to 2,900 cm−1 was typically assigned to the C–H stretching mode42 (Fig. 3j,k). In the FBTT case, the 2D SLs show a broad feature in this stretching region (Fig. 3j, red curve), which is identical to that of its aggregates (Fig. 3j, blue curve). Notably, we observed clear fine structures from FBTP SLs in this stretching region (Fig. 3k, red curve), which are absent in the corresponding aggregates (Fig. 3k, blue curve). The sharp peak at 2,960 cm−1 corresponds to the asymmetrical stretching of CH3 groups that are attached to the fluorene unit and the peaks around 3,050 cm−1 (Fig. 3k) seem to come from the C–H stretch of the aromatic ring in FBTP molecules. These sharp peaks and fine structures indicate that FBTP molecules possess sufficient freedom to vibrate like single molecules in the perovskite lattice.

Chemical tunability of SMAs

We further developed two new single-molecule-like emitters with different emission colours. Specifically, replacing the fluorene unit with a weaker donor (tolyl group) generates a sky-blue emitter, PBTP (Fig. 4a; 2-(4-(7-(o-tolyl)benzo[c][1,2,5]thiadiazol-4-yl)phenyl)ethan-1-ammonium bromide). Also, substituting the benzothiadiazole group with a stronger acceptor (benzo[1,2-c:4,5-c′]bis([1,2,5]thiadiazole)) produces a red emitter, BBTP (Fig. 4c; 8-bromobenzo[1,2-c:4,5-c′]bis([1,2,5]thiadiazole-4-yl)phenyl)ethan-1-ammonium bromide). TDDFT calculations confirmed reversed type I band alignment relative to [PbBr4]2− lattice (Extended Data Fig. 1). In parallel, XRD, GIWAXS and UV–vis measurements (Extended Data Fig. 8) indicate that organic emitters were successfully incorporated into the lead bromide matrix. The PLQYs of 2D SLs are also much higher than that of the corresponding aggregates (Fig. 4b,d). Both PL emission wavelengths and PLQYs of 2D SLs are similar to that of the corresponding monomers and blueshifted from the aggregates, demonstrating the generality of our strategy. It is noteworthy that our observation is conceptually different from aggregation-induced emission43, as aggregation-induced emission relies on the restriction of the intramolecular rotation of individual molecules, whereas our perovskite SLs maintain the single-molecular signature by regulating intermolecular distance at near-equilibrium state. Also, the newly developed 2D SLs exhibit greatly improved photostability when compared with the well-known (PEA)2PbBr4 (Supplementary Fig. 14).

Fig. 4: Single-molecule behaviour of other newly developed organic emitters in 2D SLs.
figure 4

a, Molecular structure of PBTP. b, PL spectra and PLQY of PBTP monomer in PMMA film, PBTP 2D SLs and PBTP aggregate films. c, Molecular structure of BBTP. d, PL spectra and PLQY of BBTP monomer in PMMA film, BBTP 2D SLs and BBTP aggregate films. The insets in b and d are PL images of the corresponding thin-film samples and all the image scales are 10 μm × 10 μm; the dashed violet curves are the PL spectra of (PEA)2PbBr4, which implies the emission feature of the 2D lead bromide lattice; and the grey dash-dot lines indicate the PL peaks of the monomers compared with those of the 2D SLs and aggregates.

Source Data

Molecular-ensemble properties

Owing to the ordered molecular arrangement in the SLs, we investigated the angle-dependent PL emissions of FBTP-based films (Fig. 5a and Supplementary Fig. 15). The angular PL of monomer and aggregate films exhibit a typical emission profile of isotropic emitters (Fig. 5b and Supplementary Fig. 16). Notably, the FBTP SL demonstrates a strong anisotropic feature, in which the directional emissions centre around 50° (Fig. 5c). This specific angle is indeed perpendicular to the transition dipole of the organic emitters in the perovskite SLs (Supplementary Figs. 12 and 17). Moreover, we unambiguously observed a short-lived PL component in the FBTP SLs film, which is absent in FBTP monomers film that shows constant longer decay (Supplementary Fig. 18). To precisely resolve the timescale of the short-lifetime component, we conducted temperature-dependent streak-camera characterizations. Similarly, the FBTP monomers show consistent slow PL decay at different temperatures (Fig. 5d,f and Extended Data Fig. 9a,b); whereas a short-lived component (below 100 ps) in FBTP SLs is clearly presented across wide temperature regions (Fig. 5e,f, Extended Data Fig. 9c,d and Supplementary Table 2). Notably, such an ultrafast decay has not been reported in a conventional organic single-molecule system, which might be ascribed to the mutual interaction of the radiation field of those well-aligned neighbouring single molecules44, linking to the behaviours of ordered molecular ensembles or aggregates45,46,47,48,49.

Fig. 5: SMA behaviour of organic emitters in 2D SLs.
figure 5

a, Schematic of the angle-dependent PL measurement setup configuration. b,c, Angle-resolved PL of FBTP monomers (b) and 2D SLs (c). Note that an s-polarized laser was used as the excitation source. d,e, Streak-camera measurement results on FBTP monomers (d) and 2D SLs (e) at different temperatures. f, Extracted PL decay from steak-camera data shown in d and e. The decay curves for the monomers all overlapped at different temperatures. g, Power-dependent PL intensity of FBTP 2D SLs at room temperature. h, Relative PLQY of FBTP SLs at different temperatures, using room-temperature PLQY as a reference. i, Schematic structure of lasing device with emission layer (that is, FBTP 2D SLs, monomers and aggregates) sandwiched between two DBRs. j, PL spectra evolution of FBTP 2D SL under different pump fluences. Inset, polar plots of emission intensity at pump fluences below (red) and above (blue) the threshold. k, Corresponding PL intensity against increasing pump fluences in log–log scale, showing a clear ‘kink’ at a threshold energy density of \({P}_{{\rm{th}}}^{{\rm{2D}}}\) = 0.80 μJ cm−2. The intensity dependence is fitted to a power law xp with p2D = 11.22 above the threshold. a.u., arbitrary units.

Source Data

Note that all of the temperature-dependent streak-camera results were measured with a fluence of 0.1 μJ cm−2, which is below the exciton–exciton annihilation regime (Fig. 5g), thus excluding the possibility of bimolecular-process-induced fast decay. Moreover, the relative PLQYs of 2D SLs at different temperatures are all greater than 90% and reach unity (100%) at around 150–200 K (Fig. 5h and Supplementary Fig. 19). This implies that the non-radiative channels play an unimportant role in this ultrafast decay. One may notice the slight PLQY drop from 150 to 100 K, which can be ascribed to the inorganic sublattice contraction-induced emission quenching50,51. This observation also indirectly supports our hypothesis that the emitters are sustained as single molecules by the room-temperature perovskite lattice. Last, the PL spectra of short-lived components of FBTP 2D SLs were slightly blueshifted compared with the longer components (Fig. 5e). On the basis of the intramolecular charge transfer nature of FBTP molecules (Supplementary Figs. 20 and 21), we believe that the short-wavelength emissions are probably coming from the locally excited states, which would be more feasible to undergo cooperatively radiative recombination.

We further conducted lasing characterizations by placing the emission layer between two high-reflectivity distributed Bragg reflectors (DBRs; Fig. 5i). As the pump power increases, the emission intensity from the FBTP 2D SL device steadily rises (Fig. 5j). When plotting the PL intensities against the pump fluences, a distinct ‘kink’ emerges at a threshold \({P}_{{\rm{th}}}^{{\rm{2D}}}\) = 0.80 μJ cm−2, which is accompanied by a sharp decline in emission linewidth (Fig. 5k), linearly polarized output (inset in Fig. 5j and Supplementary Fig. 22) and outstanding coherence (Supplementary Fig. 23), demonstrating the transition from spontaneous emission to full lasing action. Furthermore, the superlinear intensity dependence is fitted to a power law y = xp with p2D = 11.22 above the threshold for the 2D SL, which greatly exceeds that of the FBTP monomer (pmon = 1.70 with \({P}_{{\rm{th}}}^{{\rm{mon}}}\) = 2.71 μJ cm−2; Extended Data Fig. 10a–e). Notably, lasing emission was not observed from the FBTP aggregates-based device, probably attributable to disorder and the ACQ effect (Extended Data Fig. 10f,g). These findings together validate that the exceptional gain performance of FBTP 2D SLs resulted from the organized molecular arrangement at a near-equilibrium intermolecular distance (Supplementary Table 3).

Discussion

We have successfully incorporated a wide range of organic emitters into the 2D perovskite lattice with tunable emissions spanning from blue to green and red. We found that the molecular emitters with a suitable intramolecular twist in the perovskite SLs could retain the characteristics of single molecules. Surprisingly, these molecular emitters in perovskite SLs also exhibit dense packing and strong alignment resembling aggregates, which leads to unusual emission behaviours such as directional emission, enhanced radiative recombination rates and low-threshold lasing. With a vast selection of organic emitters of desirable properties, the hybrid SL defines a rich family of optoelectronic materials for solid-state lighting applications. For instance, preliminary investigations on light-emitting diode devices demonstrate more than 50-fold enhancements in external quantum efficiency when the FBTP molecules are confined in the perovskite 2D SLs compared with their true aggregates (Supplementary Fig. 24). Finally, this approach could be applicable to other inorganic motifs (Supplementary Fig. 25), such as layered metal halide–organic heterostructures52, molecule-intercalated layered 2D atomic crystal SLs53 and 1D or 0D organic–inorganic hybrid clusters54,55, which are yet to be further explored. In brief, the SMA confined in perovskite 2D SLs go beyond the present classification of organic matter, such as typical H, J or null aggregates56,57, representing a previously undiscovered phase at a near-equilibrium distance.

Methods

Chemicals and reagents

Organic solvents including anhydrous N,N-dimethylformamide (DMF), chlorobenzene (CB) and dichloromethane (DCM), acetonitrile (ACN), o-dichlorobenzene (DCB) and solid chemicals including lead bromide (PbBr2) were purchased from Sigma-Aldrich. All of the above chemicals were used as received. All chemical reagents and solvents for organic molecular emitter synthesis were purchased from Combi-Blocks and used as received. More details on the synthesis of organic molecular emitters, FBTT, FBTP, PBTP and BBTP, can be found in the Supplementary Information.

Thin-film fabrication

(1) The precursor solution for perovskites (2D SLs) was formulated with a concentration of 50 mM by dissolving organic emitters (that is, FBTT, FBTP, PBTP, BBTP, PEA) and PbBr2 in DMF with a stoichiometry ratio of 2:1 inside the nitrogen-filled glovebox. (2) The precursor solution for monomers was prepared by mixing 1 mg of organic emitters and 50 mg of PMMA with 1 ml mixed solvent (CB:DMFv/v = 1:1) and then stirring continuously at 100 °C overnight to obtain homogeneous polymer solutions. The weight percentage of organic emitters relative to PMMA is set at 2 wt% to ensure that organic emitters can be considered as isolated single molecules without experiencing obvious ACQ effect. (3) The precursor solution for aggregates was made by dissolving 10 mg of organic emitters in 1 ml DMF. All of the solutions were then saved for further use. Bare Si/SiO2 wafer, quartz, glass slides or DBRs were cleaned by ultrasonication in detergent, deionized water, acetone and isopropanol for 15 min each and then dried with dry air. (4) The substrates were treated with UV–ozone for 20 min and then transferred into a glovebox for spin coating. The above precursor solutions were spin-coated onto the pre-cleaned substrates at 2,000 rpm for 60 s, followed by thermal annealing on a hot plate at 150 °C for 10 min. These obtained films were used for further characterizations.

Bulk single-crystal growth

The FBTP-contained 2D SL single crystal was obtained through a vapour-diffusion method. Specifically, a 2:1 molar ratio of FBTP and PbBr2 was dissolved in DMF to make a 25-mM solution by heating at 70 °C for 1 h. Then, 0.1 ml of precursor solution was injected into a 4-ml small vial and placed in a 20-ml large vial containing 3 ml of a mixed solvent of CB and DCM with a volume ratio of 2:1, which was immediately sealed with a cap. The system was left undisturbed in a refrigerator (about 4 °C) for one month, yielding thin green plates. The FBTT-contained 2D SL single crystals were obtained through slow-cooling crystallization using a solution composed of approximately 1 mg of FBTT, 20 mg of lead bromide (PbBr2), 200 μl of hydrobromide acid (HBr, 48 wt% in H2O) and 300 μl of ethanol. Ethanol was added to assist the dissolution of the organic cations and crystallization. After mixing the precursors and solution, the contents of the sample vial were heated to over 100 °C by a heat gun until all of the materials were completely dissolved and the solution was clear. The vials were then moved to a Dewar flask water bath at 95 °C to cool down for 72 h until it reached room temperature. With this process, yellowish-orange bulk single crystals were obtained in the form of thin plates.

Single-crystalline nanocrystal growth

2D perovskite nanocrystals were synthesized using a modified co-solvent evaporation method58. 0.02 mmol of LBr (L = PEA, FBTT, FBTP) and 0.01 mmol of PbBr2 precursors were dissolved in a 2-ml solution of DMF and CB mixed in a 1:1 volume ratio to prepare 5-mM stock solutions. The concentrated PEA precursor solution was then diluted 120 times using a co-solvent system of CB, AN and DCB mixed in a volume ratio of 2.5:1:0.01. For (FBTT)2PbBr4 and (FBTP)2PbBr4, the stock solution was diluted 720 times by CB/AN/DCB co-solvent with a volume ratio of 7.4:1:0.01. A Si/SiO2 substrate was placed inside a 20-ml glass vial kept on a hot plate at 70 °C. Approximately 10 μl of the diluted precursor was then dropped on the Si/SiO2 substrate. The solvent evaporation is associated with the nucleation and growth of the nanocrystals on the substrate. The substrate was then removed from the hot plate in about 10 min, once all the solvent was evaporated.

DBR device fabrication

Three DBR devices were fabricated by sandwiching thin-film layers of FBTP 2D SLs, FBTP aggregates and FBTP monomers, respectively, between two highly reflective DBRs. The bottom DBR was first fabricated by an e-beam evaporator, consisting of 12.5 pairs of silicon dioxide (92.4 nm) and tantalum pentoxide (61.8 nm) capped by silicon dioxide. The emission layer was then spin-coated on the bottom DBR, for which the thickness was adjusted to approximately 140 nm by controlling the spin speed or the concentration of precursor solutions. Finally, the bottom DBR with the emission layer on top was transferred into the e-beam evaporator chamber again to complete the fabrication of a top DBR, which consists of nine pairs of silicon dioxide (91.8 nm) and tantalum pentoxide (65.6 nm).

Characterizations

NMR spectra

Nuclear magnetic resonance (NMR) spectra were acquired at room temperature using a Bruker AV 400-MHz spectrometer with CDCl3 or DMSO-d6 as the solvent and tetramethylsilane (TMS) as an internal standard. Chemical shifts of 1H NMR and 13C NMR signals were reported as values (ppm) relative to the TMS standard.

Mass spectra

High-resolution mass spectrometry was acquired in positive electrospray mode on an LTQ Orbitrap XL instrument (Thermo Fisher Scientific).

Single-crystal XRD analysis

Single crystals were analysed using a Bruker Quest diffractometer with kappa geometry, an I-μ-S microsource X-ray tube, a laterally graded multilayer Göbel mirror for single-crystal monochromatization and an area detector (Photon2 CMOS). Data collections were conducted at 150 K with Cu Kα radiation (λ = 1.54178 Å).

UV–vis absorption spectra

Thin-film absorption spectra were recorded on an Agilent UV-Vis-NIR Cary 5000 spectrometer in transmission mode.

PL spectra

Steady-state PL spectra were obtained with an Olympus microscope system (BX53) integrated with an X-CITE 120Q UV lamp. The filter cube contains a band-pass filter (330–385 nm) for excitation, a dichroic mirror (cut-off wavelength 400 nm) for light splitting and a 420-nm long-pass filter for emission collection. The collected PL signals were analysed by a spectrometer (SpectraPro HRS-300).

PLQY

The thin-film samples for PLQY measurements were deposited onto quartz substrates by following the preparation of precursor solutions and the fabrication procedures detailed in the ‘Thin-film fabrication’ subsection. The absolute PLQYs at room temperature were obtained by a three-step technique with a home-designed system, which consists of a continuous-wave laser (375 nm), an integrating sphere, optical fibre and a spectrometer. The relative PLQYs at low temperatures were estimated on the basis of the integrated emission intensity of the PL spectra at different temperatures for the 2D SLs film. By taking the PLQY of the sample at room temperature as a reference and correcting for absorption48, the relative PLQYs of the film were then calculated.

Powder XRD

Thin-film XRD was collected with Rigaku SmartLab (Cu Kα, λ = 1.54056 Å) in Bragg–Brentano mode.

Thickness measurement

The thickness of the thin-film samples was measured with a Bruker DektakXT stylus profilometer. Here the thickness of FBTT, FBTP, PBTP, BBTP and PEA 2D SLs samples were determined to be 55.3, 58.2, 24.3, 23.1 and 20.8 nm, respectively.

Atomic force microscopy

The surface morphology and roughness were obtained with a MultiMode 8-HR AFM (Bruker) in tapping mode.

Scanning electron microscopy

Thin-film scanning electron microscopy images were acquired with a high-resolution field-emission scanning electron microscope (SU8010).

FLIM

FLIM measurements were performed using a Nikon TE2000 confocal microscope with water immersion objective (60×, NA = 1.2) equipped with an Alba Fast FLIM system (ISS). Specifically, samples were excited using a 440-nm pulsed laser with modulation frequency of 10 MHz and imaged through a 506-nm long-pass filter, followed by MPD APD detectors. After image collection, biexponential fitting of FLIM images was performed using VistaVision software (ISS) to obtain the fluorescent lifetimes of each pixel.

GIWAXS

GIWAXS spectra were collected at beamline 7.3.3 at the Advanced Light Source at Lawrence Berkeley National Laboratory using an incident angle of 0.18° and wavelength of 1.24 Å (energy 10 keV). The detector used was PILATUS 2M (Dectris, Inc.) and the data were calibrated using silver behenate as a standard using the Nika Igor Pro package59.

FTIR

Attenuated total reflectance FTIR spectroscopy was conducted on a Thermo Nicolet Nexus 470 FTIR, equipped with a diamond-attenuated total reflectance crystal sampling accessory, with N2 purging.

Low-temperature PL and time-resolved PL measurements

A home-built confocal micro-PL setup was used to carry out temperature-dependent steady-state and time-resolved optical measurements. A 447-nm picosecond pulsed diode laser (LDH-P-C-450B, PicoQuant, 50 ps, 5 MHz) was used as the excitation source and focused onto the surface of the samples using a Nikon objective (40×, NA = 0.6). The emitted signal was collected by the same objective, dispersed with a monochromator (Andor Technology) and detected by a spectrometer (Andor Shamrock 3030i) and charge-coupled device (CCD; Andor Newton 920). The excitation scatter was rejected using a suitable optical filter placed before the detector. For the time-resolved PL dynamics, the signal was detected by a single-photon avalanche diode (PicoQuant, PDM Series) with a single-photon counting module (PicoQuant), with a time resolution of about 100 ps. For most of the measurements, a low excitation fluence of 0.1 μJ cm−2 was used to prevent the onset of parasitic bimolecular processes such exciton–exciton annihilation.

For the temperature-dependent PL measurements, a closed-cycle optical cryostat Cryostation s50 (Montana Instruments) was used. The sample was placed on the holder and a hard vacuum was established in the sample chamber. Temperatures in the range 100–295 K, with a temperature stability <10 mK, were then attainable using the control units.

Angle-dependent PL measurements

Angle-resolved PL measurements were performed on a home-built system. The sample is illuminated with a pulsed laser (100 fs, 400 nm) produced by the second-harmonic generation from a Ti:sapphire laser system. The laser pulse polarization is controlled by a polarizer and half-wave plate and impinges on the sample at normal incidence. The excited PL light is collected by a lens coupled to an optical fibre and finally directed into a spectrometer. The angle subtended by the lens is approximately 1°.

Temperature-dependent streak-camera measurements

The samples were excited with 440-nm-wavelength light pulses from an ORPHEUS optical parametric amplifier powered by a Pharos amplifier with a 2-kHz repetition rate and 170-fs pulse duration. The beam was focused with a lens into a 1-mm spot size on the sample that was mounted in a cryostat and held under a vacuum of <8 × 10−5 mbar. The PL was collected by an achromatic lens and guided into a spectrograph (Andor Kymera 328i), which was connected to a charge-coupled device camera (Andor iXon Life 888) and streak camera (Hamamatsu C10910) for performing time-integrated PL and time-resolved PL measurements, respectively.

Photostability test

The photostability test on the 2D perovskite thin-film samples were carried out by tracking their absorption spectra under the irradiation of a UV curing lamp in a glovebox, in which the UV lamp had an output power of 0.31 W cm−2 and was 5 cm away from the samples.

Lasing characterizations

Optically pumped lasing measurements were carried out on a home-built far-field micro-PL system in ambient conditions. The excitation pulses (400 nm, about 100 fs, 1 kHz) were generated from the second harmonic of the fundamental output of a regenerative amplifier (Solstice, Spectra-Physics, 800 nm, about 100 fs, 1 kHz), which was in turn seeded by a mode-locked Ti:sapphire laser (Mai Tai, Spectra-Physics, 800 nm, about 100 fs, 80 MHz). The DBR devices with emission layer sandwiched between two DBRs were locally excited with a laser beam focused down to about 50 µm in diameter through an objective (Nikon CFLU Plan, 5×, NA = 0.15), with input power altered by neutral density filters. After passing through a 420-nm long-pass emission filter, the collected PL signal from the DBR devices was subsequently coupled to a grating spectrometer (Acton SP-2358) and recorded with a thermoelectrically cooled CCD (Princeton Instruments, ProEm 1600B).

Coherence measurements

The spatial coherence of the lasing emission from the DBR device was evaluated by a Michelson interferometer setup. Initially, the emission from the DBR device was divided into two beams with a beam splitter, which were subsequently directed to two separate arms of the interferometer. These beams were then reflected by the interferometer mirrors and overlapped on a CCD camera. The length of one interferometer arm was precisely adjusted to ensure that both beams travelled the same distance before reaching the CCD camera. Clear interference fringes will be recorded if the initial emission exhibits coherence.

Simulations

MD simulations

Unbiased MD sampling: the modified MYP model was used for all simulations, as described in our previous work. LAMMPS and PLUMED were used to perform the MD simulations. All simulations used a 1-fs integration time step and periodic boundary conditions. Long-range electrostatics was modelled using the particle–particle–particle–mesh (PPPM) algorithm and Lennard-Jones interactions were truncated at 15 Å. The initial structure of 2D perovskites was generated by constructing representative unit cells of ideal perovskites lattice with the bulky organic cations placed at the surface. The simulation was first relaxed in the NVE ensemble with restrained atomic displacements of 0.01 Å per time step for 50 ps, followed by a 100-ps NPT equilibration with the Nosé–Hoover thermostat and barostat. The boundary of the y direction, which is parallel to the bulky organic cations, is extended by 20 Å to prevent the interaction of organic cations from different sides as a result of periodic boundary conditions. During the NPT equilibration, the barostat was only applied to the x and z directions, which are normal to the bulky organic cations. Finally, a 100-ps NVT simulation is conducted to evaluate the distribution of backbone dihedrals, the site energy of ligands and the molecular planarity parameter.

Free-energy calculation on the ligand rotation: the reaction coordinate of the free-energy calculation is defined as the angle between the fitted planes of two randomly picked neighbouring ligands. Steered MD was then used to calculate the free-energy curve. In steered MD, we use a spring constant of 1,000 kcal (mol-Å)−1 and a constant velocity of 0.0125 rad ps−1 to steer the ligand into the target-packing angle. All of the reported results are calculated on the basis of five independent simulation runs.

DFT calculations

Geometry optimizations and excited-state calculations for FBTT, FBTP, PBTP and BBTP molecules were carried out by means of DFT and TDDFT as implemented in the Gaussian 16 package60 using the B3LYP functional61 and the def2-TZVP basis set62 in vacuum. The alkylammonium tail is not included for brevity, which tends to be flexible and varies its geometry depending on the matrix. In our experience, an expensive basis set such as def2-TZVP is necessary to achieve adequate accuracy in TDDFT calculations when considering the band alignment in perovskites. Transition dipole moments were visualized with the Multiwfn software63.

Optical-field-distribution simulation

The electric-field intensity distribution in the 2D SLs-based DBR device was calculated in the electromagnetic-wave-frequency domain by using the commercial software COMSOL Multiphysics. A periodic boundary condition was applied to the modelling configuration. The transmitted spectrum was retrieved by setting the bottom side of the 2D SLs layer as a periodic excitation port.

Reporting summary

Further information on research design is available in the Nature Portfolio Reporting Summary linked to this article.