Anti-Symmetric Electromagnetic Interactions’ Response in Electron Circular Dichroism and Chiral Origin of Periodic, Complementary Twisted Angle in Twisted Bilayer Graphene

Many novel physical properties of twisted bilayer graphene have been discovered and studied successively, but the physical mechanism of the chiral modulation of BLG by a twisted angle lacks theoretical research. In this work, the density functional theory, the wavefunction analysis of the excited state, and the quantum theory of atoms in molecules are used to calculate and analyze the anti-symmetric chiral characteristics of zigzag-edge twisted bilayer graphene quantum dots based on periodic complementary twisted angles. The analysis of the partial density of states shows that Moiré superlattices can effectively adjust the contribution of the atomic basis function of the fragment to the transition dipole moment. The topological analysis of electron density indicates that the Moiré superlattices structure can enhance the localization of the system, increasing the electron density of the Moiré central ring, reducing the electron surge capacity in general and inducing the reversed helical properties of the top and underlying graphene, which can be used as the origin of the chiral discrimination; it also reveals the mole in the superlattice chiral physical mechanism. On this basis, we will also study the nonlinear optical properties of twisted bilayer graphene based on a twisted angle.


Introduction
Twisted bilayer graphene (TwBLG) is a robust correlation system obtained by stacking coupling and interlayer rotation between one layer of graphene nanosheets and another layer of graphene nanosheets. This interlayer interaction will cause the symmetry of the electron to break [1]. For example, most Mott insulators exhibit an antiferromagnetic arrangement of electron spin [2,3]. In 2018, Pablo Jarrilo-Herrero et al. experimentally observed for the first time that magnetic angle graphene could achieve insulation and low-temperature superconductivity under strong electric field conditions slightly above absolute zero, which aroused considerable interest among researchers [4]. Defective bilayer graphene (BLG) with a small twisted angle will generate the early experimental preparation process [5]. The distortion angle causes the slight dislocation of the cellular lattice to produce a long period Moiré superlattice. This induces a periodic potential field and a microgap at the edge of the Briviine region of the Moiré superlattice (MS) [6]. As the interaction increases, the Quantum Anomalous Hall Effect (QAHE) and chiral edge states will occur near the insulating state [7]. The MS can be produced in a variety of ways, the latest being by cold-pressing graphite sheets [8]. The Moiré superlattice can effectively regulate the distribution of the magnetic transition dipole moment to control the chirality [9]. The MS structure has a circular dichroism due to its differential response to polarized light [10]. This optical activity originates from the relative rotation of the top and underlying layers, leading to the rotation of the helicity of Dirac Fermions [11,12]. The

Results
The UV-vis and ECD spectra are calculated to study the optical absorption characteristics of quantum dot system. As shown in Figure 1c, the spectra of Z-TwBLG-QDs-27.8° and Z-TwBLG-QDs-32.2° overlap. The main absorption peak is near 772.3 nm, which is contributed to by S15, and the second absorption peak is contributed to by two excited states (S49 and S53). The main absorption peak of Z-TwBLG-QDs-32.2° is near 772.4 nm,

Results
The UV-vis and ECD spectra are calculated to study the optical absorption characteristics of quantum dot system. As shown in Figure 1c, the spectra of Z-TwBLG-QDs-27.8 • and Z-TwBLG-QDs-32.2 • overlap. The main absorption peak is near 772.3 nm, which is contributed to by S 15 , and the second absorption peak is contributed to by two excited states (S 49 and S 53 ). The main absorption peak of Z-TwBLG-QDs-32.2 • is near 772.4 nm, which is contributed to by two excited states (S 15 and S 16 ), and the secondary peak is contributed to by two excited states (S 49 and S 53 ). As shown in Figure 1d, ECD shows two regions with high rotatory strength in the visible light region. When the twisted angle is 27.8 • , the corresponding ECD excited states are S 51 /S 53 and S 47 /S 49 . When the twisted angle is 32.2 • , the corresponding excited states are S 51 /S 53 and S 48 /S 49 . The excited states corresponding to the absorption peaks in the near-infrared region are S 13 and S 14 . In addition, the rotatory intensity of the excited states of the two systems at the same wavelength is axisymmetric, and the absolute value is almost the same (see Table S2 in supplementary). Therefore, the electromagnetic interaction intensity of their incident light-quantum dots is very close. The symmetrical rotatory strength makes the molar absorption coefficients of Z-TwBLG-QDs-27.8 • and Z-TwBLG-QDs-32.2 • satisfy an axisymmetric relationship, indicating that the molecular system exhibits symmetric circular dichroism in response to polarized light.
Circular dichroism is the asymmetric response of the electric and magnetic transition dipole moments when the system interacts with electromagnetic waves, which can reflect the chiral properties of the system. The TEDM and TMDM are plotted to reveal the physical mechanism of the axial symmetry of the ECD spectrum. The transition dipole moment density is the product of the transition density and the coordinate variables in the three directions x, y, and z of the cartesian coordinates [21]: The transition dipole moment density here is connected to the transition electric dipole moment, specifically. The contribution of the x, y, and z components to TEDM is obtained. The component of the transition magnetic dipole moment can be defined as The component of the transition magnetic dipole moment m l (r) is defined to visualize the transition magnetic dipole moment. The explicit expression is Its full-space integral is the transition magnetic dipole moment as follows: To make the distribution characteristics of TEDM and TMDM more precise and more intuitive, the smoothing electric/magnetic moments (TEDM S and TMDM S ) of the system on the right side of TEDM and TMDM in  13 , and their TEDM and TMDM on the x and y components are perpendicular to each other. ECD shows that the rotatory intensities of the two excited states contribute negatively to ECD. When the twisted angle changes to 32.2 • , the absorption peak at the same wavelength is contributed to by S 13 . Compared with S 14 , the TEDM and TMDM on the x, y, and z components are anti-symmetric, so S 13 is responsible for ECD absorption. S 14 is responsible for chirality.
Comparing the ECD excited state of the two systems at 844.5 nm, the distribution range of TEDM S in S 13 is more extensive in the x component, and the TEDM of S 14 is more extensive in the y component, see Figure 2d  Comparing the ECD excited state of the two systems at 844.5 nm, the distribution range of TEDMS in S13 is more extensive in the x component, and the TEDM of S14 is more extensive in the y component, see Figure 2d,g,e,h. On the z component, the TEDM and TMDM distributions are also anti-symmetric. The TEDMS and TMDMS show that the positive and negative isosurface face have opposite directions, see Figure 2f,i, which indicates  The ECD excited state of the two systems at 612.8 nm is contributed to by S49. When the twisted angle is 27.8°, the TEDMs on the x and y components are distributed axisym- component, the positive and negative isosurfaces also separate to a greater extent and the directions are opposite, see Figure 4j-l. Therefore, the axisymmetric ECD of the two systems at 612.7 nm is mainly caused by the anti-symmetric distribution of the positive and negative isosurfaces on the z component.  The above analysis discusses the electromagnetic interaction mechanism of Z-TwBLG-QDs based on the periodic complementary twisted angles. It is necessary to investigate the DOS analysis of the area with anti-symmetric distribution of TEDM and TMDM to reveal the physical mechanism of the anti-symmetric distribution of TEDM and TMDM, corresponding to fragments 1-5, see Figure 5a,b. Fragment definition adds the basis functions on the atom to the fragment, which is marked with different colors. As shown in Figure 5c,d, the normalized integral curve of PDOS represents the contribution of each fragment's atomic orbital, the relative positions of the blue curve, and the yellow curve alternate near the energy of −2.77 eV. According to Mulliken division, the contribution of the r-th basis function to the transition dipole moment is as follows: where P tran is the transition density matrix between certain two states and χ r |−r|χ s is the electric dipole moment integral between the r − th and s − th basis functions. The contribution of the atoms to the transition dipole moment is to sum the basis function contributions of all atoms since the transition electric dipole moment is the transition dipole moment from the ground state to the excited state. Therefore, the alternating phenomenon of the contribution of the fragment atomic orbits will cause significant changes in the response mode of the transition electric dipole during the light-quantum dot interaction of the system. Combining Formulas (1) and (2), alternating atomic orbital contributions will change the orbital wavefunctions ϕ a and ϕ b , thereby causing the changes in the transition electric dipole moment and the transition magnetic dipole moment. These changes will indicate the anti-symmetric distribution of TEDM and TMDM for the two systems. The MS structure that is based on the periodic complementary twisted angles on the atomic orbitals of the edge region are the cause of this situation. As shown in Tables S4-S6 in supplementary, the real space function of the RCP of the Coronene and the two Z-TwBLG-QDs are calculated through the topological analysis of the electron density. The of the two systems are very close. All RCP are numbered to analyze the correlation between the and the chirality of the system. The red serial numbers correspond to the RCP of the Moiré superlattice in the top graphene, and the blue serial numbers correspond to the RCP of the Moiré superlattice in the underlying graphene. Figure 6a indicates the changing trend of the normalized . The polarization of the Moiré superlattice increases as the electron density increases after removing an RCP with the smallest polarization. The red circular arrow on the top layer of the graphene is the reduced order of the electron density, see Figure 6c. The rotation direction is clockwise, so the top layer graphene is a right-handed spiral, while the underlying graphene is a lefthanded spiral, see Figure 6b,d. The Laplace function of the electron density (∇ ) is to calculate the gradient of : The above discussion qualitatively reveals the physical mechanism of the axisymmetric ECD of the quantum dot system. To investigate the origin of the symmetric circular dichroism quantitatively, it is necessary to calculate the helical properties of chiral interactions in molecular systems. These can be investigated by the polarization (the sensitivity of electron density to deformation in an electric field) of atomic fragments in chiral molecules according to the theory of molecular polarizabilities [21]. The electron density is defined as where η i is the occupied number of orbital i, ϕ is the orbital wavefunction, χ is the basis function, and C is the coefficient matrix. The element in the i − th row and j − th column correspond to the expansion coefficient of orbit j relative to the basis function i. When the twisted angle is 27.8 • , the MS structure contains an axisymmetrically distributed Moiré superlattice along the central benzene ring of the graphene layer. The RCP of the Moiré superlattice is found in the chemical bond formed by the benzene ring on the top and underlying graphene. As shown in Tables S4-S6 in supplementary, the real space function of the RCP of the Coronene and the two Z-TwBLG-QDs are calculated through the topological analysis of the electron density. The ρ RCP of the two systems are very close. All RCP are numbered to analyze the correlation between the ρ RCP and the chirality of the system. The red serial numbers correspond to the RCP of the Moiré superlattice in the top graphene, and the blue serial numbers correspond to the RCP of the Moiré superlattice in the underlying graphene. Figure 6a indicates the changing trend of the normalized ρ RCP . The polarization of the Moiré superlattice increases as the electron density increases after removing an RCP with the smallest polarization. The red circular arrow on the top layer of the graphene is the reduced order of the electron density, see Figure 6c. The rotation direction is clockwise, so the top layer graphene is a right-handed spiral, while the underlying graphene is a left-handed spiral, see Figure 6b,d. The Laplace function of the electron density (∇ 2 ρ RCP ) is to calculate the gradient of ρ RCP : Molecules 2022, 27, 6525 10 of 14 The normalized ∇ is the same as the curve of normalized because ∇ is the second derivative of , see Figure S3(a and b). The numerical value indicates that ∇ is positive, indicating that the interlayer interaction originates from the closed-shell interaction of the π-π bond in the Moiré superlattice. The interaction that is different from the covalent bond comes from the shared electron pair, so there is no electron accumulation in the bonding area. The increase in the degree of aggregation will reduce the surging ability of electrons and weaken the role of electrons in the electrostatic interaction of the region of MS, making the nuclear electrostatic potential of Z-TwBLG-QDs relative to Coronene increase, see Figure 6e,f.
The values of between the Coronene and the two Z-TwBLG-QDs reflect that the MS will enhance the localization. It is necessary to analyze the localized orbital locator (LOL) of the Coronene and the two systems of Z-TwBLG-QDs to further reveal the influence of the MS on the localization of the system. Tables S4-S6 in supplementary show that The normalized ∇ 2 ρ RCP is the same as the curve of normalized ρ RCP because ∇ 2 ρ RCP is the second derivative of ρ RCP , see Figure S3a,b. The numerical value indicates that ∇ 2 ρ RCP is positive, indicating that the interlayer interaction originates from the closed-shell interaction of the π-π bond in the Moiré superlattice. The interaction that is different from the covalent bond comes from the shared electron pair, so there is no electron accumulation in the bonding area. The increase in the degree of aggregation will reduce the surging ability of electrons and weaken the role of electrons in the electrostatic interaction of the region of MS, making the nuclear electrostatic potential of Z-TwBLG-QDs relative to Coronene increase, see Figure 6e,f.
The values of ρ RCP between the Coronene and the two Z-TwBLG-QDs reflect that the MS will enhance the localization. It is necessary to analyze the localized orbital locator (LOL) of the Coronene and the two systems of Z-TwBLG-QDs to further reveal the influence of the MS on the localization of the system. Tables S4-S6 in supplementary show that the value of LOL at the central RCP of the two systems is significantly greater than that of Coronene. The color-filled map of LOL-pi at 1.2 Bohr above the plane of the molecular system was drawn to visually analyze this difference. As shown in Figure 7c, the LOL-pi coloring diagram of Coronene shows that the difference from RCP is that the outer-shell electrons of the carbon atom of the central six-membered ring have no significant changes due to the stable conjugated structure. On the contrary, the π electrons delocalize between adjacent carbon atoms at the six protruding points of the zigzag edge. After adding the Moiré, the delocalization path of the π electron does not include the six-cornered outer benzene ring. Additionally, the delocalization gradually weakens from the outside to the inside because the gradient of orbital wavefunction of the system at the boundary is larger than the inside, see Figure 7d,e. Another tool for analyzing the delocalization of electrons is the electron localization function (ELF), which has a similar physical meaning to LOL. The changing trend of the ELF value of the three systems' RCP is the same as that of LOL, see Figure 7b. The previous discussion indicates that the electrons at the edge of Z-TwBLG-QDs are more active than the internal ones under the edge effect. The ELF_pi isosurface can intuitively reflect the delocalization characteristics of the molecular system contributed to by π electrons, see Figure 7f. When the pi isovalue of the six-membered ring in the Coronene center is broken, On the contrary, the π electrons delocalize between adjacent carbon atoms at the six protruding points of the zigzag edge. After adding the Moiré, the delocalization path of the π electron does not include the six-cornered outer benzene ring. Additionally, the delocalization gradually weakens from the outside to the inside because the gradient of orbital wavefunction of the system at the boundary is larger than the inside, see Figure 7d,e.
Another tool for analyzing the delocalization of electrons is the electron localization function (ELF), which has a similar physical meaning to LOL. The changing trend of the ELF value of the three systems' RCP is the same as that of LOL, see Figure 7b. The previous discussion indicates that the electrons at the edge of Z-TwBLG-QDs are more active than the internal ones under the edge effect. The ELF_pi isosurface can intuitively reflect the delocalization characteristics of the molecular system contributed to by π electrons, see Figure 7f. When the pi isovalue of the six-membered ring in the Coronene center is broken, the isovalue is 0.495. When the pi isovalue of the six-member ring at the edge is broken, the isovalue is 0.518. Therefore, the conjugation degree of the central six-membered ring is lower than the edge, see Figure S4a,b. When the isovalue is set to 0.7, it can be clearly seen that the π isosurface at the edge is more significant than the internal one, and this law also applies to the more extensive system such as Z-TwBLG-QDs. As shown in Figure 7g,h, the π isosurface of Z-TwBLG-QDs-27.8 gradually decreases from the edge to the internal, and the degree of π conjugation decreases. This situation is the same when the twisted angle is 32.2 • .

Methods
The quantitative calculation of the two molecular systems consists of two steps. The first step is using Gaussian16 [22] combined with GFN2-xTB [23] to optimize the ground state geometry, and the GFN2 XTB method is similar to the semi empirical DFT of DFTB, with good accuracy and universality. The further optimization of the ground state structure uses the DFT method [24], B3LYP functional [25], the 6-31G(d) basis set [26], and DFT-D3 dispersion-correction [27]. The step of excited state calculation is to calculate the Boltzmann distribution of the model with internal energy, and the final structure is calculated by time-dependent density functional (TDDFT). The output of this step is used to calculate and visualize the electron-hole pair density, the transition density matrix (TDM), and all the configuration coefficients of the TMDM and TEDM. The wavefunction of excited state generates the electron-hole pair density through the linear combination of the single excitation configuration function [28]. TEDM and TMDM are determined by the tangential vector of the transition density space [29,30]. The molecular orbitals of the PDOS are calculated by the Hartree-Fock equation theoretical method, which is based on a single-electron approximation [31]. The topological analysis of electron density is based on the QTAIM [32], and the real space function (a function with three-dimensional space coordinates as variables) of the ring critical point (RCP) of the Moiré superlattice was obtained.

Conclusions
In this work, various first-principle calculation methods are used to analyze the chiral physical mechanism of Z-TwBLG-QDs. The analysis of partial density of states indicates that the MS based on the periodic complementary twisted angle will change the contribution of the basis functions of the edge atoms of the system to the transition dipole moment. That will cause the orbital wavefunction to change and induce an anti-symmetric transition electric dipole moment and a transition magnetic dipole moment, resulting in axisymmetric circular dichroism in ECD. The topological analysis of electron density shows that the chirality is caused by the opposite helical characteristics of top graphene nanosheets and the underlying graphene nanosheets of the Moiré superlattices system. The structure of the Moiré superlattices weakens the conjugation degree of the edge π electrons, which has a positive effect on the localization of the molecular system. This research is of great significance in revealing the chirality generation and regulation mechanism of layered two-dimensional chiral materials and elucidating the structure-function relationship of chiral molecules.

Supplementary Materials:
The following supporting information can be downloaded at: https: //www.mdpi.com/article/10.3390/molecules27196525/s1, Figure S1: The normalized of Z-TwBLG-QDs-27.8 • and Z-TwBLG-QDs-32.2 • (a, b), Figure S2: The isovalue of the central six-membered ring of Coronene when the connected region of the ELF_pi isosurface was just fractured (0.495). The isovalue of the edge of Coronene when the connected region of the ELF_pi isosurface was just fractured (0.518), Table S1: The wavelength and oscillator intensity of main one-photon absorption excited state in Z-TwBLG-QDs-27.8 • and Z-TwBLG-QDs-32.2 • , Table S2: The wavelength and rotatory intensity of main ECD excited states in Z-TwBLG-QDs-27.8 • and Z-TwBLG-QDs-32.2 • , Table S3: The real space function values of the RCP in Coronene, Table S4