Molecular properties controlling chirality transfer to halide perovskite: computational insights †

Interfacing chiral molecules with achiral semiconductor nanomaterials has been shown to induce chiroptical ( i.e. , preferential absorption of circularly polarized light) signatures. This chirality transfer phenomena is attractive for enhancing the functionality of nanomaterials by making them sensitive to the circular polarization states of photons, and attractive for selective optical sensing applications. Chirality transfer via interfacing chiral molecules with nanomaterial surfaces generally gives small sensitivities as determined by anisotropy factors, the ratio of polarized and linear absorption. Here we use atomistic time-dependent density functional theory simulations to investigate the molecular properties which influence the chiroptical signatures of a lead-halide perovskite cluster. First, we find that for chiral molecules that contain aryl groups, such as ( R / S )-methylbenzylamine and ( R / S )-methylphenyl acetic acid, modulating their molecular dipole via meso functionalization with strong withdrawing groups, such as NO 2 , can enhance intensities of chiroptical signatures. Second, using the lactic acid series (( S )-lactic, ( S )-malic, and ( S,S ) tartaric acid) we average over conformations of these molecules on the cluster surface. We find that limited conformational flexibility of lactic acid provides the largest chiroptical intensities while signal from tartaric acid quenches after conformational average. These results demonstrate an importance of chemical functionalization leading to a polarized chiral system. Altogether, either through acid–base equilibrium or strategic functionalization, and limited conformationally degrees of freedom are important for obtaining high-intensity nanomaterial chiroptical signatures through chirality transfer.


Introduction
Chiroptical signatures have fundamental and applied applications. This is realized for small molecules where spectroscopy methods, such as electronic and vibronic circular dichroism, can be used to assign absolute molecular configurations. Crystalline materials show chiroptical signatures when crystallizing in symmetry breaking crystal groups 1,2 or when the chirality is induced (i.e., chirality transfer) from the semiconductors environment. 3,4 Chiroptical signatures can also emerge when an achiral material is placed within an applied magnetic field 5 or under applied stress. 6 For applications related to optical sensing an important device parameter is the so-called anisotropy factor which is the ratio of circularly polarized absorption and linear absorption. The anisotropy factor sets the upper-bound on the ability for a material to reliably distinguish between circular polarization states of photons in linear spectroscopy.
An attractive strategy for inducing chiroptical signatures into semiconductor nanomaterials is using chiral organic molecules as surface passivation ligands on semiconductor nanomaterial surfaces. 3,[7][8][9][10] The field of surface chemistry, specifically regarding colloidal stability [11][12][13] and surface composition 14,15 has undergone significant advancements and is now highly developed. Using chiral molecules would allow new materials to be synthesized and processed using established methods, but they are expected to provide enhanced optical functionality. An example of this is using the nanomaterials sensitivity of photon circular polarization to construct a quantum dot based spin bit device. 16 Many studies focused primarily on binding chiral amino acids to semiconductor nanomaterials, such as (X)-penicillamine, 7 (X)-cystine 8,[17][18][19] and (X)-glutathione, 3 or other types of amines or ammoniums, such as (X,X), diaminocyclohexane, 20 (X)-methylbenzylammonium, 10 and (X)-octylbutylammonium where X = R,S denoting molecular chirality. 10 Chiral carboxylic acids have been also been explored with quantum dots that exhibit excess metal on the surface that their conjugate bases, such as lactate, malate, and tartartate, can bind to. 21 For chirality induction through chiral molecules on 0D quantum dot/nanocrystal surfaces the reported anisotropy factors are relatively low, on the order of 10 À5 -10 À3 . 21 This motivates the exploration of design principles for how to enhance anisotropy spectra in semiconductor quantum dot systems.
There have been three conjectures as to the mechanism of chirality transfer between a semiconductor nanomaterial and a chiral molecule. One possible explanation is that templating a surface with chiral molecules introduces forces which coherently alter the structures. 22,23 Charge transfer excitations have also been observed to result in emergent CD spectra between molecular dimers, so that hybridized semiconductor -chiral molecule excitations are also possible. 24 Dipole interactions between semiconductor and chiral molecules has also been shown to result in emergent CD spectra with the orientation of the dipole field dictating preference of polarization in optical absorption. 25 Our previous work has shown the chirality transfer mechanism is primarily due to the dipole interaction 26 which has also been used to explain emergent circular dichroism in polymer blends. 27 This motivates the investigation of structure-property relationships for these chiral molecular dipole induced chiroptical signatures. The contribution and manipulation of surface chemistry in relation to the impact of chirality-transfer induced circular dichroism is an aspect that has not been thoroughly explored.
Chemical modifications of the molecular ligand can enhance the optical properties of the quantum dot. Functionalizing with withdrawing groups is expected to increase the molecular dipole which can couple to the quantum dot and enhance optical features. 28 Another aspect that has not appeared to be explored is the role of molecular conformations on the surface of the nanomaterial. For isolated molecules, such as in dilute solutions, it is well established that the circular dichroism is sensitive to changes in conformations and an averaging procedure is needed to accurately predict circular dichroism spectra from ab initio methods. 29 Previous studies have noticed a binding configuration effect where there was a noticeable distinction between semiconductor CD when the molecule is in a different absolute configuration. 17,26 Using atomistic time-dependent density functional theory (TD-DFT) simulations, here we set out to explore two scenarios: (i) how chemical modifications and (ii) molecular conformations of the chiral molecule on the surface of a nanomaterial impact the resultant circular dichroism and anisotropy spectra from a halide perovskite nanomaterial (i.e. chirality transfer). For (i) chemical modification we initially explore a chiral methyl-aryl tetrahedrally coordinated carbon with either amine, ammonium and carboxylic acid functional groups. Here we provide a comment on chemical notation. To highlight the charge of the chemical species bound to the perovskite cluster we explicitly place a charged symbol on cation and anion moieties. When discussing the role of aryl functionalization at the meso location on increasing the polarization of the chiral molecule we will refer to as an added character to the end of the common molecular name, so instead of using the cumbersome (R)-4-chloro-methylbenylammonium naming convention we use (R)-MBAmm + -Cl. In Fig. 1 we show (a) (R)-methylbenzyl amine (MBA-H), (b) (R)-methylbenzyl ammonium (MBAmm + -H) cation, and (c) (R)-methylphenyl acetate (MPAc À -H) anion, to bind to the perovskite surface. We then test how substituting the meso position of the aryl groups with electron donating and withdrawing groups, such as hydroxy (OH) and nitro (NO 2 ) groups, alter the chirality induction effect with the expectation that the stronger ground-state dipole of electron withdrawing groups induces higher intensity CD peaks. We find that when comparing the chirality transfer induced circular dichroism of (R)-MBA-H, (R)-MBAmm + -H, and (R)-MPAc À -H to a CsPbI 3 cluster that (R)-MPAc À -H gives the largest CD peaks and anisotropy intensities. Modifying the meso position of the aryl groups for this set of chiral molecules with electron donating OH or electron withdrawing Cl ions, the circular dichroism remains unchanged. In contrast, functionalizing with NO 2 results in larger CD intensities for all models, particularly most pronounced for (R)-MBA-NO 2 . We interpret this as a trade-off between charged molecules and their polarizability. The charged anionic (R)-MPAc À -H and (R)-MBAmm + -H surface ligand provides greater polarization of the perovskite cluster compared to the neutral (R)-MBA-H. However, the (R)-MBA-H molecule is more polarizable than (R)-MPAc À -Hand (R)-MBAmm + -H and, consequently, shows greater enhancement in CD peak intensities with aryl substitution of withdrawing NO 2 groups.
To address (ii) exploring the influence of conformational degrees of freedom of chiral molecules on the surface, we use the molecules from the lactic acid series Fig. 1(d) (S)-lactic, (e) (S)-malic, and (f) (S,S) tartaric acids. Experimentally, it is observed that these molecules show relatively large anisotropies for their lowest optical band at approximately 10 À2 . 30 An increase of molecular size lead to a concomitant increase of molecular conformational degrees of freedom. As such, (S)-malic and (S,S)-tartaric ligands adapt three principle conformations of trans, gaunche(À), and gaunche(+). When averaging over conformations of these molecules on the surface, we find (L)-lactic acid provides comparable circular dichroism and anisotropy intensities as (R)-MPAc À , while the effect drastically decreases when going to (S)-malic and (S,S) tartaric acids. We attribute this to chiroptically 'dark' conformers which are stabilized by the surface. Overall, this work provides modeling insight into molecular properties and degrees of freedom which impact the intensity of emergent CD from the perovskite quantum dots through chirality transfer. Namely, charged and polarizable molecular ligands provide enhanced chirality transfer and conformationally simple molecules are likely to provide more stable CD spectra due to the limited influence of chiroptically 'dark' conformations, which can be stabilized by the surface.  Fig. 1(g) illustrates the binding site of chiral molecules to the perovskite surface. The chiral molecules are bound to the surface as a chemical moiety that has overall +1 charge, as either neutral (R)-MBA-H bound to Cs + , a ammonium cation (R)-MBAmm + -H, or a Pb 2+ (COOR) À compound which corresponds to a deprotonated (R)-MPAc À , (S)-lactate, (S)malate, or (S,S)-tartate. We explore excess Pb 2+ on the surface of LHP nanocrystals as it has been shown to be a prevalent defect for as-synthesized 31 LHP nanocrystals which introduces shallow-defect states. 32,33 When binding (R)-MBA-H, (R)-MBAmm + -H, and (R)-MPAc À -H we neglect conformer averaging due to the conformation simplicity of these molecules, noting that (R)-MBAmm + -H shows no change in CD spectra when rotated on the cluster surface (as discussed later in text). From the optimized structures we substitute electron withdrawing groups, such as NO 2 and Cl, and donating groups OH at the meso position of the aryl group and keep the nuclear coordinates frozen, but allowing the donating/withdrawing groups to optimize.

(i) Considered structures
When binding the lactic acid series to the perovskite cluster surface we account for a distribution of initial conditions based on the trans, gaunche(À), and gaunche(+) conformers for (S)-malic and (S,S)-tartaric acid. For (S)-malic acid we consider six distinct initial conditions while for (S,S)-tartaric acid we model 3 distinct initial conditions. For the conformationally simpler (S)-lactic acid we consider 3 various initial conditions based on rotation about the chiral center. The optimized structures do not necessarily reflect the initial conditions and are characterized by the dihedral angle between the carboxylic acid groups.

(ii) Chemical modifications
We first explore the chirality transfer properties based on binding different functional groups to a CsPbI 3 halide perovskite cluster. Halide perovskite nanocrystals have a diverse surface chemistry due to the binary cationic-anionic surface termination of (100) crystal planes, as illustrated in Fig. 1(g). This allows for a variety of functional groups to bind to the surface, such as alkylamines, 34  Even though the chiral centers of the molecules are the same, we can see that the functional group and binding mode to the surface have a large impact on the CD spectrum with (R)-MPAc À -H providing an order of magnitude larger response compared to that of (R)-MBA and multiples larger than (R)-MBAmm + -H. This trend is also demonstrated in the anisotropy spectra, Fig. 2(e), with the two lowest energy CD peaks providing the largest anisotropy maxima.
From the optimized structures, we further explore functionalizing the meso positions of the aryl group with electron donating (withdrawing) OH (Cl,NO 2 ) groups. This is done to modulate the molecular dipole of the chiral molecule, which we hypothesize can enhance the intensity of chiroptical signatures due to their dependence on the magnitude of transition dipoles, see eqn (4) in Methods. To illustrate the effect of chemical functionalization on tuning the molecular dipole we first inspect the isolated molecules. Specifically, we examine the change in molecular polarization by functionalizing the meso position of the aryl group for the neutral molecules  Table S1 (ESI †). While for the ionic complex (R)-MBAmm-I we observe that the dipole magnitude increases significantly, the meso functionalization provides a narrower range of enhancement from 11.55 to 14.01 Debye.
When the functionalized chiral molecules are bound to the halide perovskite cluster surface, shown in ESI, † Fig. S2, the OH and Cl functionalization's provide a small change to the overall dipole moment, but NO 2 can provide 1.5-2 times increased dipole moment, as shown in Table S1 (ESI †). For (R)-MBA-H, we observe a significant enhancement of the CD peak intensities with NO 2 functionalization to (R)-MBA-NO 2 , while OH and Cl functionalizations (R)-MBA-OH and (R)-MBA-Cl provide almost identical CD spectra to the pristine models, as shown in Fig. 3. We note that the same trend emerges for functionalized (R)-MBAmm + -H and (R)-MPAc À -H, shown in ESI, † Fig. S3 and S4 (ESI †), respectively, but with NO 2 functionalized (R)-MBAmm + -NO 2 and (R)-MPAc À -NO 2 providing less pronounced increase in CD peak intensities. This is attributed to the trend observed for the chiral molecules that (R)-MBA-H as a neutral species is more polarizable while the cation (R)-MBAmm + -H and anion (R)-MPAc À -H are less polarizable species.
To confirm that the observed chiroptical signatures of the molecule-cluster system originate from the halide perovskite cluster, we overlay the CD and absorption spectra of the isolated cluster, isolated molecule and the combined system for (R)-MBA-H, (R)-MBAmm + -H, and for (R)-MPAc À -H in ESI, † Fig. S5-S7, respectively. These plots show that the absorption spectra of the combined system are identical to the isolated cluster and that CD signals emerge below the optical gap of the isolated chiral molecule. For the NO 2 functionalized models, we do detect a spectral feature below the optical gap of the isolated cluster, as shown in Table S1 (ESI †) the optical gap of the combined system shows a smaller optical gap compared to the unfunctionalized models. This feature is attributed to a partial intra-molecular charge transfer transition due to the strong electron-withdrawing nature of the NO 2 group, see ESI, † Fig. S8.
An interesting observation is that for the neutral (R)-MBA-H molecule introducing an electron-withdrawing group NO 2 through aryl substitution of the meso-position giving (R)-MBA-NO 2 , increases the maximum |g anisotropy | by a factor of 5. One possible explanation for this is due to the increased polarization at the bonding interface. This would consequently polarize the molecular orbitals and increase the magnetic dipole response. To investigate this trend, we analyze the charge distribution using Bader charge analysis and compare the charge associated to the  cluster-ligand interface, illustrated in ESI, † Fig. S9. From comparing the Bader charges, shown in ESI, † Table S4, we find there is a negligible change in charge associated with atoms associated with the molecular-cluster bond. Thus, bond polarization does not contribute to the changes in the chiroptical response with a polar group functionalization.
To further investigate the origin of CD intensity enhancement for a functionalized (R)-MBA-H system, we inspect the CD spectra calculated with a narrow linewidth of 100 meV. This allows us to examine contributions of specific electronic transitions. In Fig. 4(a) we overlay the CD spectra broadened by 350 meV with the narrower broadening of 100 meV for (R)-MBA-H and (R)-MBA-NO 2 . From this comparison, we observe that the increase in the |g anisotropy | is attributable to a significant increase in intensity between 270 nm and 290 nm, as seen in Fig. 4(a), which is identified as peak (i) and peak (ii) for (R)-MBA-H and (R)-MBA-NO 2 , respectively, in Fig. 4(b). For this CD peak, we further single out the most intense rotatory strengths and analyze the NTOs, shown in Fig. 4(c) and (d). We observe that the transitions are identical in terms of orbital density distribution. Upon further investigation we find that the enhancement is attributable to an increase of the transition electric dipole, as evident by the increased oscillator strength in ESI, † Table S3, while the transition magnetic dipole remains largely unchanged. This is an illustration of the concept of tuning the electronic polarization of the molecule to enhance the magnitude of chirality transfer driven by dipole coupling.
Another interesting observation is that the anionic (R)-MPAc À -H bound to the surface gives larger CD and |g anisotropy | intensities compared to that of cationic (R)-MBAmm + -H. Computing the root-mean-square for the rotatory strengths up to 5 eV for (R)-MBAmm + -H (N = 68) and (R)-MPAc À -H (N = 68) gives values of 440 and 480, respectively. This indicates that, on average, it would be expected the CD and |g anisotropy | maxima should be 10% larger for (R)-MPAc À -H compared to (R)-MBAmm + -H. The values in Table 1 show the first band and maximum band intensities increase by 910% and 62%, respectively. We find that the main differences in the CD spectra arise from the change in density of excited-states (DOES), as shown in Fig. 5. To directly compare the CD spectrum of these models we plot the spectrum normalized to the optical gap. Fig. 5 shows the bandgap normalized CD spectrum overlayed with their DOES for (a) (R)-MBAmm + -H and (b) (R)-MPAc À -H. Their EDOS are overlayed in Fig. 5(c). Comparing their DOES shows greater splitting of the lower energy transitions (below 0.25 eV) for (R)-MPAc À -H than that for (R)-MBAmm + -H. Thus the greater splitting between excited-states underpins the larger intensity CD spectra in (R)-MPAc À -H than (R)-MBAmm + -H due to reduced overlap between positive and negative contributions to the thermally broadened CD spectrum. This illustrates that not only the magnitude of chiroptical signatures matter, but the distribution of excited states, which contribute to chiroptical activity, have to be considered for trying to maximize the anisotropy signal.
We note that for these chiral aryl molecules we neglect conformational averaging on the surface. This is due to their inherent lack of conformational degrees of freedom aside from rotation of the functional groups around the chiral center. As a test, we explored various rotations of the (R)-MBAmm + -H molecule on the CsPbI 3 cluster surface and find that the CD signal remains virtually unchanged, as shown in ESI, † Fig. S10.

(iii) Molecular conformations
Motivated by the observation that (R)-MPAc À -H provides the highest intensity CD peaks and most resolvable anisotropy peaks, see Fig. 2, we next explore the lactic acid series, i.e., (S)-lactic, (S)-malic, and (S,S) tartaric acids. The chiroptical signatures of these molecules have been explored recently, which shows that in aqueous solutions tartaric acid provides 4 times the CD intensity enhancements compared to Lactic and Malic Acids along with each of them showing |g anisotropy | maxima greater than 0.01 which is large for organic molecules. 42 We hypothesize that these desirable chiroptical features can be transferred to the perovskite cluster. Similar to (R)-MPAc À -H these molecules are bound to the cluster in the de- When geometry is optimized, as shown in ESI, † Fig. S11-S13, the trans and gaunche+ conformers remain in these conformers, while the gauncheÀ conformers are optimized into a structure in between trans and gauncheÀ with dihedrals ranging between 24 and 28 degrees, as seen in ESI, † Table S4.
We further compute Boltzmann weights and thermal average for each system with the values shown in Table S5 (ESI †). These thermal averages are then used to calculate the conformationally averaged CD spectra for models, which are shown in Fig. 6(a). For (S,S) Tartaric Acid, the gaunche(+) conformers are the most energetically favorable, in contrast to the solution phase where trans structure is typically observed. 43,44 This is attributed to the ability for gaunche(+) to form additional hydrogen bonds between the carboxylic group and surface  halides as well as hydroxyl groups forming coordination with surface Cs + cations. In Fig. 6(a) we observe that with (S)-lactate (red, solid) bound to the surface, the induced CD spectra from the perovskite cluster has much more intensity compared to the CD spectra induced by (S)-malate (green, solid) and (S,S)-tartartate (blue, solid). Fig. 4(b) shows computed the anisotropy factors for the conformationally averaged systems. Here we find a similar trend: Pb(lactate) (red, solid) bound to the surface has higher maximum intensity peaks than the induced by (S)-malate (green, solid) and (S,S)-tartartate (blue, solid) spectra. Overall, the molecules which have conformational flexibility, effectively 'quench' the CD signal from the CsPbI 3 perovskite cluster. This trend could have two possible origins. On one hand, each conformer could have 'mirrored' CD spectra and similar Boltzman weights. Hence, the signal would be quenched after averaging. On the other hand, there could be conformers that provide 'bright' and 'dark' chirality transfer with the nanomaterial surface stabilizing a 'dark' conformer. From the Boltzmann weights shown in ESI, † Table S5, we find that the gaunche(+) conformers of (S)-Malate and (S,S) Tartarate provide greater than 90% weight. Further, when examining the computed CD for 1, gaunche(+) for both molecules, as shown in ESI, † Fig. S14, it can be seen that these conformers induce 'dark' CD into the CsPbI 3 perovskite cluster. In contrast, trans and gaunche(À) conformers of (S)-malate and (S,S)tartartate on the surface can induce 'bright' chiroptical signatures comparable to intensity similar to (S)-lactate. Thus, we observe that simply binding chiral molecules to a nanomaterial surface does not guarantee 'bright' induced chirality transfer. Conformations of the chiral molecule greatly impact the induced nanomaterial's CD. Conformational effects have also shown to be important in determining tthe metal cluster circular dichroism spectrum as well. 45 To explain the conformational dependence of chirality transfer and induced CD, we borrow ideas from literature. In a recent paper a curious correlation was found between the intensity of circular dichroism of the material and degree of spin-polarization of charge-carriers through the chirality induced spin selectivity (CISS). 46 In the same spirit, we try to identify a mechanism for the observed quenching of circular dichroism spectra for perovskite clusters with gaunche(+) conformers on their surface based on CD spectra from the isolated molecule (i.e. correlation between isolated molecule CD intensity and CD induced by chirality transfer). ESI, † Table S6 shows the rotatory strength for the lowest energy S 1 excited-state for the anionic molecule, the Root-Mean-Square (RMS) of the rotatory strength for all computed excited-states (N = 20), and the RMS for all computed excited-states for the joint cluster/ molecule system (N = 300). In short, there is no clear correlation between isolated molecule conformation and the chirality transfer strength into the perovskite cluster. One observation is that the gaunche(+) conformers (5-M, 6-M, 3-T) and other chiroptically 'dark' models (4-M) have negative R S1 values. In contrast, for chiroptically active models, the isolated molecules have positive R S1 value. The exact coupling mechanism describing conformer's dependence is difficult to determine. Analysis such as non-degenerate coupled oscillator 19 or machine learning of molecular descriptors 47 are tools to be used to gain more insight.
Finally, we remark on a trend related to molecules giving the largest chiroptical intensities. The largest computed anisotropies are summarized in Table 1. Within the lactic acid series, we find that (L)-Lactate acid provides the largest anisotropy intensities. The other molecular systems, which gave consistently large anisotropy intensities include (R)-MPAc À -H and (R)-MBAmm + -H with (R)-MPAc-NO 2 providing the largest anisotropy across all models investigated. All these molecules are conformationally limited, as they have one chiral center and are composed of a single central carbon atom. This indicates that conformational rigidity is important for inducing large chiroptical signatures via chirality transfer.

Conclusion
Using chiral molecules to induce chiroptical signatures into semiconductor nanomaterials (via chirality transfer) is an attractive prospect as it allows to leverage developments in colloidal surface chemistry to introduce enhanced functionality for opto-electronic applications, such as polarized photodetectors. Such applications require high fidelity to distinguish between binary circularly polarized states, as characterized by anisotropy spectra. Typically for colloidal semiconductor quantum dots the maximum anisotropy value is in the range of 10 À5 -10 À3 . Here we use time-dependent density functional theory calculations to investigate molecular properties, such as charge state and conformational flexibility, which influence the chirality transfer between a chiral molecule (ligand) bound to an achiral lead halide perovskite cluster.
First, we find that for chiral molecules that contain aryl groups, such as (R/S)-methylbenzylamine and (R/S)-methylphenyl acetic acid, modulating their molecule dipole via meso functionalization with strong withdrawing groups, such as NO 2 , can significantly enhance intensities of chiroptical signatures. Along the same line of thought, modulating the acid-base equilibrium of the chiral molecular to charged species, which binds to the nanomaterial surface, also increases the chirality transfer effect. Second, using the lactic acid series ((S)-lactic, (S)-malic, and (S,S) tartaric acids) we explore conformational ensemble average of these molecules on the cluster surface. With increasing molecular size and functionalization, there are also increasing number of chiral centers, which could contribute to enhancement of chiroptical signal. We find that limited conformational flexibility of lactic acid provides the largest chiroptical intensities while (S,S) tartaric acid effectively quenches them. Overall, these results demonstrate that chemical modification emphasizing a polarized chiral molecule, either through acid-base equilibrium or strategic functionalization, and limited conformationally degrees of freedom, are important for obtaining highintensity nanomaterial chiroptical signatures featuring enhanced chirality transfer.

Methods
Ground-state geometry and electronic structure are found from the self-consistent DFT equations which minimizes total energies E k for each conformer k. Excited-electronic states are found using time-dependent density functional theory (TDDFT) by solving the Casida equations. Briefly, the solutions of the Casida equation produce eigenvectors X n and Y n , which represent transition densities, with eigenvalues O n which represent the transition frequencies. For chiroptical activity the optical transitions of interest are the electronic transition dipole, eqn (1) and magnetic transition dipole, eqn (2). These transition dipoles are represented as where c j E and c i j i are Kohn-Sham wavefunctions of occupied and virtual electronic states, respectively, e represents the fundamental charge, m e is electron mass, and c is the speed of light. Within TDDFT framework, these dipole moments are found from expectation values of the transition densities 48 m elec 0n;k ¼ X ijm elec ij;k X ij;n;k þ Y ij;n;k À Á (3) m mag 0n;k ¼ X ijm mag ij;k X ij;n;k À Y ij;n;k À Á where 0 and n represent ground and excited-states, respectively. From the transition dipoles, the corresponding oscillator strength, eqn (5), and rotatory strength, eqn (6), for an electronic transition of frequency O n can be determined where they describe the probability of absorbing linearly polarized light and circularly polarized light The extinction coefficient e and differential extinction coefficient De, describing circular dichroism, for each conformer k are then generated with eqn (7) and (8) e k e À E k À E min kT (10) De k e À E k À E min kT (11) where E min and E k are the conformer with the lowest total energy and for a specific conformer, respectively. To compare the relative intensities of the linear and polarized absorption, we further compute the anisotropy factor 4,49 from the computed extinction and differential extinction coefficients as All calculations are done using Gaussian 16 package. 50 For chiral molecules the structures are optimized with dispersion corrected ob97xd functional and a 6-311++g** basis set. Both nuclear coordinates of the perovskite cluster and chiral molecule are allowed to relax during optimization. For the pristine perovskite cluster the stoichiometry Cs 19 Pb 8 I 36 was chosen giving a total charge of À1. This structure was chosen to minimize any asymmetry induced circular dichroism signal. When adding molecules to the cluster surface the charge is kept at À1. TDDFT calculations for 20 excited states are then done using the same model chemistry combined with a polarized continuum model (PCM) to account for solvent effects, namely a conductor-like PCM (CPCM) formalism 51,52 is used with solvent dielectric constant corresponding to toluene (e = 2.37). For geometry optimization of the combined perovskite cluster/ chiral molecule models ob97xd functional and a split basis set is used with heavy elements {[Cs, Pb, I]} being treated with LANL2DZ effective core potentials 53 while light elements [C,H,O,N,Cl] being treated with 6-311++g(d,p) basis set. 54 300 electronic excited-states are then computed for each structure using TDDFT formalism and C-PCM solvent model corresponding to toluene dielectric constant. To visualize the distribution of the computed transition density we used the concept of natural transition orbitals 55 (NTOs). When computing spectra for the isolated molecules, as shown in ESI, † 20 states are calculated. When computing the perovskite cluster with chiral molecules on the surface we compute 300 states for the shown circular dichroism spectra.
To ensure that the anisotropy factor 'behaves nicely' near the absorption onset transition energy, we add a constant number to the absorption spectra to avoid divergences for the singular denominator in the regions of band absorption onsets. For consistency the constant is chosen to be 10% of the first absorption peak intensity for each conformer averaged system. Bader charge analysis is done with the charge densities computed from the pseudo-potential 56,57 plane-wave based Vienna ab Initio Software Program (VASP). 58 A plane-wave basis set with a cut-off of 500 eV is used in a simulation cell with 11 A of vacuum in each cardinal direction. The Bader charges 59 where determined using the Bader Charge Analysis program. 60

Author contributions
All authors have given approval to the final version of the manuscript.

Conflicts of interest
There are no conflicts to declare.