Next Article in Journal
Effective Removal of Different Heavy Metals Ion (Cu, Pb, and Cd) from Aqueous Solutions by Various Molecular Weight and Salt Types of Poly-γ-Glutamic Acid
Next Article in Special Issue
Preparation and Support Effect of Graphdiyne Nanotubes with Abundant Cu Quantum Dots
Previous Article in Journal
Unravelling the Flotation Performance of 1-Hydroxy-2-naphthyl hydroxamic Acid and Styrene Phosphonic Acid Collectors on Monazite Using Experiments and DFT Calculations
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Physical Mechanism of Linear and Nonlinear Optical Properties of Nanographene-Induced Chiral Inversion

1
College of Science, Liaoning Petrochemical University, Fushun 113001, China
2
Liaoning Provincial Key Laboratory of Novel Micro-Nano Functional Materials, Fushun 113001, China
*
Authors to whom correspondence should be addressed.
Molecules 2024, 29(5), 1053; https://doi.org/10.3390/molecules29051053
Submission received: 3 January 2024 / Revised: 26 February 2024 / Accepted: 26 February 2024 / Published: 28 February 2024
(This article belongs to the Special Issue Advanced Carbon Nanomaterials and Their Applications)

Abstract

:
Based on density functional theory (DFT) and wave function analysis, the ultraviolet and visible spectrophotometry (UV-Vis) spectra and Raman spectra of 1-meso and 1-rac obtained by the chiral separation of chiral nanographenes are theoretically investigated. The electron excitation properties of 1-meso and 1-rac are studied by means of transition density matrix (TDM) and charge density difference (CDD) diagrams. The intermolecular interaction is discussed based on an independent gradient model based on Hirshfeld partition (IGMH). The interaction of 1-meso and 1-rac with the external environment is studied using the electrostatic potential (ESP), and the electron delocalization degree of 1-meso and 1-rac is studied based on the magnetically induced current under the external magnetic field. Through the chiral separation of 1-rac, two enantiomers, 1-(P, P) and 1-(M, M), were obtained. The electrical–magnetic interaction of the molecule is revealed by analyzing the electron circular dichroism (ECD) spectra of 1-meso, 1-(P, P) and 1-(M, M), the transition electric dipole moment (TEDM) and the transition magnetic dipole moment (TMDM). It is found that 1-(P, P) and 1-(M, M) have opposite chiral properties due to the inversion of the structure.

1. Introduction

Graphene is a two-dimensional carbon nanomaterial composed of carbon atoms with sp2 hybrid orbitals in a hexagonal honeycomb lattice [1]. It is ultra thin, ultra light and super strong, and has many interesting properties, thus being known as the “king of new materials” [2,3]. Due to the hexagonal lattice structure formed by carbon atoms, graphene has two typical types of edges, namely, armchair edges and zigzag edges [4]. Graphene has attracted increasing attention due to its unique physical properties (such as strong mechanical strength, large surface area and good thermal, optical and electrical conductivity). The discovery of graphene has significantly accelerated research on the manufacture of low-cost electrode materials [5,6,7,8,9,10]. Graphene and its derivatives, including graphene oxide (GO) and reduced graphene oxide (r GO), are emerging as an important class of nanomaterials in the field of biosensors [11,12,13,14,15].
Since the discovery of graphene manufacturing methods by K. S. Novoselov et al. in 2003, the study of nanographene has become a major trend as the size of graphene devices reduces [16,17,18,19,20]. Nanographene is a type of graphene fragment with dimensions in the nanometer range. In recent years, with the development of nanographene chemistry, it has been increasingly applied in electronic devices, energy storage, sensing and other fields [21,22,23,24].
Chirality is widespread in nature. If an object is different from its mirror image, it is referred to as “chiral”, and its mirror image cannot coincide with the original object [25,26]. Chiral objects and their mirror images are called enantiomers. Chiral nanographenes with a high fluorescence quantum yield ( ϕ F ) and large asymmetry factor ( g lum ) are essential for the development of circularly polarized luminescence (CPL) materials. Wang at el. proposed a new strategy to extend the helical π of luminous groups with a high diameter ϕ F to obtain chiral nanographenes with excellent luminous properties while maintaining the distribution of frontier molecular orbitals (FMOs). Using perylene as the parent core, the authors synthesized chiral nanographenes with an ϕ F content of 93% and separated them by chiral HPLC to obtain 1-meso and 1-rac [27].
Besides experimental studies, there is a lack of theoretical work on 1-meso and 1-rac, particularly in the analysis and interpretation of their UV-Vis spectra, charge transfer properties and chiral mechanism, as well as other physical mechanisms. In this paper, the two chiral structures are introduced. Based on quantum chemical calculations and wave function analysis, the structure is theoretically studied by generating UV-Vis, ECD and Raman spectra, and the physical mechanism of light absorption is explained and combined with the visualization of the charge transfer process. This study provides a theoretical basis for the application of chiral nanographene structures in optoelectronics, linear optics and nonlinear optics.

2. Results

By taking perylene as the core, 1-meso and 1-rac were designed and synthesized as chiral nanographenes. The two enantiomers of chiral nanographenes were obtained through chiral separation. Their structures are presented in Figure 1a,b. The directions and angles of molecular torsion on both sides are different, as shown in Figure 1c,d. Figure S1 draws the IGMH diagram between the upper and lower segments of 1-meso and 1-rac at the torsion angle position. There are green isosurfaces in the close distance between the two intersecting benzene rings, suggesting the presence of interlayer interaction between the upper and lower layers. However, the electron density in the interaction region is very low. In this section, the optical properties, chirality mechanism and vibration modes of 1-meso and 1-rac will be investigated.

2.1. HOMO-LUMO Molecular Orbital Analysis

Figure 2a and Figure 2b, respectively, present the HOMO-LUMO molecular orbitals of 1-meso and 1-rac, where the positive and negative phases of the highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO) are shown in red and blue, respectively. The HOMO and LUMO energies of both 1-meso and 1-rac are depicted in Figure 2. It is worth noting that the HOMO-LUMO gap of 1-meso is measured at 2.535 eV, while that of 1-rac is 2.551 eV. The HOMO and LUMO energies indicate the ease with which a molecule can gain or lose electrons, and the HOMO-LUMO gap reflects the ability of electrons to transition from occupied to empty orbitals. A higher gap value indicates a lower probability of electron transition. Comparatively, 1-meso exhibits a greater tendency for excitation than 1-rac. By examining the molecular orbital diagrams, both the HOMO and LUMO diagrams demonstrate a similar distribution of electrons and holes, mainly concentrated on the central perylene. In addition, the transition from S0 to S1 represents a local excitation process.

2.2. The Physical Mechanism Analysis of UV-Vis Spectra

We have plotted the UV-Vis spectra of 1-meso and 1-rac, as shown in Figure S2a,b. Figure 3 presents the combined UV-Vis spectra of 1-meso and 1-rac. The compounds 1-rac and 1-meso exhibit similar absorption spectra, with a slight redshift in the absorption spectra of 1-meso. This is consistent with the results of the HOMO-LUMO gap analysis. In the experiment, the maximum absorption values of the lowest energies of the compounds 1-rac and 1-meso are 538 nm and 539 nm, respectively, which is relatively small compared to the maximum absorption values of 545 nm and 546 nm of the lowest energies of 1-rac and 1-meso in the calculation of TD-DFT, reflecting the accuracy of the calculation [27]. The UV-Vis spectra show that the absorption spectra bands of 1-meso and 1-rac are mainly within 220–700 nm. They both have five absorption peaks. Due to their similar electronic structure, the photophysical properties in the ground state are not affected by the conformation.

2.3. Electronic Transition Properties

To further analyze the excitation properties of electron spectra, TDM and CDD maps of excited states where 1-meso and 1-rac contribute significantly to the absorption peaks are drawn. In the CDD diagram, the red isosurfaces represent electrons, and the blue isosurfaces represent holes. The expression for electrons and holes is as follows:
ρ hole ( r ) = ρ ( loc ) hole ( r ) + ρ ( cross ) hole ( r )
ρ ( loc ) hole ( r ) = i a ( w i a ) 2 φ i φ i i a ( w i a ) 2 φ i φ i
ρ ( cross ) hole ( r ) = i a j i a w i a w j a φ i φ j i a j i a w i a w j a φ i φ j
ρ ele ( r ) = ρ ( loc ) ele ( r ) + ρ ( cross ) ele ( r )
ρ ( loc ) ele ( r ) = i a ( w i a ) 2 φ a φ a i a ( w i a ) 2 φ a φ a
ρ ( cross ) ele ( r ) = i a i b a w i a w i b φ a φ b i a i b a w i a w i b φ a φ b
where r is the coordinate vector, φ is the orbital wave function, i or j is the occupied orbital label, and a or b is the empty orbital label. i a represents each excitation configuration of the cycle, and i a represents each de-excitation configuration of the cycle.
Figure 4 shows the TDM and CDD diagrams of 1-meso. The absorption peak at 545 nm is contributed by S1, as shown in Figure S2a. The TDM and CDD of S1 are shown in Figure 4a. The transition density is mainly distributed in the lower left corner, and the electrons and holes are evenly distributed in the perylene in the middle of 1-meso. The absorption peak at 368 nm primarily originates from S9, which is characterized by the transition density being predominantly distributed on the left and lower sides. Local excitation occurs on the perylene molecule located in the center of 1-meso, with electrons from both sides transferring to the perylene center. The electronic transition mode observed under S9 involves local excitation accompanied by charge transfer, as depicted in Figure 4b. The absorption peak at 291 nm is primarily attributed to transitions occurring in S21 and S30, with the transition density primarily concentrated in the lower left corner of the diagonal, albeit with a small amount dispersed in other positions. Electrons and holes are evenly distributed across the perylene molecule, with a small presence in other regions. These characteristics of local excitation are depicted in Figure 4c,d. Similarly, the absorption peak at 273 nm is ascribed to S47, while the absorption peak at 241 nm is primarily contributed by S100. The electron transition mode in these two excited states primarily involves charge transfer excitation. Specifically, electrons in S47 transition from the middle of 1-meso to both the upper and lower sides, while electrons from both sides of 1-meso transfer to perylene in S100.
Figure S3 presents the TDM and CDD of 1-rac in various excited states. The excited states corresponding to the UV-Vis absorption peaks of 1-rac are highly similar to those of 1-meso molecules, except for the S31 and S49 excited states. A comparison of the TDM and CDD diagrams reveals nearly identical electronic transition properties, which can be attributed to their similar electronic structures.

2.4. ESP Analysis

In the ESP diagram of the molecular surface, the ESP value of different surface regions is represented by different colors, making the ESP distribution on the molecular surface clear. The value of the ESP reflects the ability of the molecule to interact with other molecules. Figure 5, respectively, shows the ESP of the molecular surface for 1-meso and 1-rac. Red indicates a positive ESP value, while blue indicates a negative ESP value. The yellow dots represent the maximum, and the blue dots represent the minimum. It can be observed that the ESP in the region where the C atom is located is negative, with the minimum value distributed near the C atom; in contrast, the ESP in the region where the H atom is located is positive, with the maximum value distributed above the H atom. The maximum value of 1-meso is significantly higher than that of 1-rac, while there is almost no difference in the minimum value. Consequently, the H atom corresponding to the maximum ESP of 1-meso is more prone to undergoing nucleophilic reactions.

2.5. Response to External Magnetic Field and Magnetic Induced Current Density

If the electrons in a molecular system have strong delocalization, a significant induced ring current will be generated when a magnetic field is applied. Therefore, we studied the induced currents of 1-meso and 1-rac induced by an external magnetic field, which is perpendicular to the molecular plane from inside to outside, as shown in Figure 6, Figures S4 and S5 (in the Supplementary Materials). The induced currents of 1-meso and 1-rac are clearly circular, and the red circular arrow indicates the direction of the system’s circular current. It can be seen from the figure that the current isosurface in the yellow region is mainly distributed on the inner side of the benzene ring. The ring current in each benzene ring is clockwise, and a clockwise ring current is also formed at the edge of the molecule. This indicates that 1-meso and 1-rac have strong electron dissociation.

2.6. Raman Spectral Analysis

The Raman spectra are a type of scattering spectra, where the horizontal axis represents the frequency of the scattered light in relation to the incident light, and the vertical axis indicates the intensity of the scattered light. The Raman activity is a unique property of each vibration mode. The Raman spectrum can reflect the characteristic information of the molecular structure of substances, and is a unique property in the study of molecular chemical structure and chemical composition analysis.
In the UV-Vis spectra, the most prominent absorption peaks for the 1-meso and 1-rac conformers within the ultraviolet wavelength range are observed at 291 nm, while the most intense absorption peak within the visible light band occurs at 545 nm. Therefore, we select two laser wavelengths (325 nm and 532 nm) that closely match these peaks to irradiate the respective structures, facilitating the acquisition of their resonance Raman spectra. Figure 7 shows the Raman and resonance Raman spectra of 1-meso (a) and 1-rac (b). It can be seen from Figure 7 that the resonance Raman spectra of 1-meso and 1-rac after laser irradiation have the same Raman peak wave number as the static Raman spectra, and the strongest Raman peak in both the static Raman spectra and the resonance Raman spectra is located near 1600 cm−1. However, the Raman intensity is significantly enhanced, indicating that both wavelengths of laser can significantly improve the Raman intensity of 1-meso and 1-rac by about two orders of magnitude.
Then, the vibration attribution of 1-meso and 1-rac Raman characteristic peaks are analyzed and explained, respectively. With different frequencies, the vibration modes of the molecule are also different. In the vibration mode diagram, the direction of the vector arrow indicates the main direction of the atomic vibration, and the length of the arrow indicates the intensity of the vibration, which is the energy involved in the vibration. Figure 8 shows the vibration modes corresponding to 1-meso in different bands. The vibration mode at 1578.29 cm−1 corresponds with the contraction of the central benzene ring and the vibration of the surrounding hydrogen atoms. On the contrary, the vibration mode at 1595.75 cm−1 represents the lateral vibration of hydrogen atoms, as shown in Figure 8a,b. Figure 8c,d depict the vibration modes observed at the prominent peak position of 1-meso in the 325 nm resonance Raman spectra, where the vibration mode at 1600.64 cm−1 signifies the vibration of hydrogen atoms and the stretching vibration of the carbon atoms in the benzene ring. Additionally, the vibration mode at 1629.95 cm−1 indicates the lateral oscillation of H atoms. Figure 8e,f display the vibration modes of 1-meso at the position of the prominent peak in the 532 nm resonance Raman spectra. The corresponding vibration patterns at the higher peaks are attributed to the vibration of the H atoms. Notably, the vibration mode at 1600.64 cm−1 corresponds to the same position under 325 nm laser irradiation.
Figure 9a,b depict the vibration mode of 1-rac observed at the peak position in the static Raman spectra. Both the peaks at 1581.80 cm−1 and 1601.47 cm−1 exhibit similar vibration patterns, which are attributed to the movement of C and H atoms. Figure 9c,d illustrate the vibration mode of 1-rac corresponding to the prominent peak in the 325nm resonance Raman spectra. Specifically, the vibration observed at 1514.38 cm−1 is related to the motion of C atoms and H atoms on both sides, excluding those associated with the intermediate perylene. In Figure 9e,f, the vibration mode of 1-rac is observed at the position of the strong peak in the 532 nm resonance Raman spectra. Here, the vibration observed at 1373.28 cm−1 is associated with the stretching of the central benzene ring and the lateral vibrations of additional H atoms. Notably, the vibration pattern observed at 1601.47 cm−1 remains consistent under the irradiation of three different lasers. Due to the certain symmetry of the molecules, their vibration patterns are also symmetrically stretched.

2.7. ECD Spectra and Transition Electric\Magnetic Dipole Moment

Ji-Kun Li et al. experimentally obtained two enantiomers of 1-(P, P) and 1-(M, M) by the chiral separation of 1-rac at room temperature [27], as shown in Figure S6. Therefore, we judge the chirality of 1-meso and 1-rac enantiomers by calculating their electron circular dichroism, as shown in Figure 10.
ECD can effectively characterize the chirality of chromophores in molecules. Theoretically, the strength of ECD can be defined as
φ j | μ e | φ i φ j | μ m | φ i B
μ e represents the TEDM and μ m represents the TMDM. Here, the first term in the absolute value is light absorption and the second term is circular dichroism. The dominant ECD spectra are the tensor products of the TEDM and TMDM. Figure 10a shows the ECD spectra of 1-meso, while Figure 10b,c present the ECD spectra of two enantiomers of 1-rac. Figure 10d is their combined diagram. It can be observed from the figure that the chirality of 1-meso is weak, the chirality of the two enantiomers of 1-rac are strong, and the images of 1-(P, P) and 1-(M, M) are symmetric and have opposite Cotton effects.
Next, the properties of these three structural chirality were analyzed by studying the density of the TEDM\TMDM in different directions. As shown in Figure 10a, 1-meso has a positive rotatory strength in S10, S46 and S100, and a negative rotatory strength in S1, S8 and S43. The TEDM density of the six excited states successively along the X, Y and Z directions (Figure 11a and Figures S7–S11 in the Supplementary Materials). The TMDM density increases in all three directions, exhibiting an alternating distribution of positive and negative regions. Moreover, both the positive and negative isosurfaces of the TEDM and TMDM are predominantly distributed on the left and right sides of the molecule. Notably, the distributions of the TEDM and TMDM for states S47 and S100 are similar, with the positive and negative TMDM densities in the Z direction clearly separated. This indicates that the chiral effects of states S47 and S100 have a similar nature. In structure 1-(P, P), the rotatory strength of states S10 and S94 is positive, whereas that of states S1 and S100 is negative. Conversely, in 1-(M, M), the rotatory strength of S10 and S94 is negative, while that of S1 and S100 is positive. It is worth nothing that within the same excited state, the positive and negative isosurfaces of 1-(M, M) for both the TEDM and TMDM along the X, Y and Z directions are nearly opposite to the positive and negative isosurfaces of 1-(P, P); refer to Figure 11b,c and Figures S12–S17 in the Supplementary Materials. This discrepancy is due to the inherent symmetry between the chirality of 1-(P, P) and 1-(M, M). It is remarkable that the TEDM and TMDM densities of the three structures are mainly distributed on the left and right sides of the structure, rather than on perylene, indicating that the chirality is mainly attributed to the asymmetry on both sides of the structure.
Table 1 shows the TEDM and TMEM components and their eigenvalues of the excited states at X, Y and Z of 1-meso, 1-(P, P) and 1-(M, M), which make the highest contribution to the ECD spectra. The absolute value of the molecular tensor product represents the light excitation intensity, and the greater the absolute value, the stronger the excitation intensity. From the table, it can be seen that when the absolute value of the 1-meso tensor product is low, 1-(P, P) and 1-(M, M) are relatively high. Thus, the two structures obtained from 1-rac by chiral separation have strong chirality. The intensity and direction of each ECD spectra in the table were completely matched. The results show that the TEDM\TMDM is complete and self-consistent in inferring the chiral mechanism.

3. Calculation Method

In this research, we constructed the 1-meso and 1-rac models using GaussView 6.0.16 software [28] and Gaussian 16 (A.03) program [29]. The optimization of the structures was carried out based on the density functional theory (DFT) [30] calculation method, utilizing the B3LYP [31] functional and 6–31g(d) [32] basis set, along with DFT-D3 [33] correction. Subsequently, the B3LYP [34,35] function set and 6–31g(d) base set were used for the electronic excitation calculation using the TD-DFT calculation method, followed by the visualization and analysis of the excitation results using Multiwfn [36]. The UV-Vis spectra, ECD spectra, TMD and Raman spectra of the system were generated using Origin 2022 software [37]. Furthermore, employing the VMD program [38], we plotted the CDD diagram, TEDM, TMDM, IGMH and ESP of the system [39,40,41]. Additionally, the vibration modes of 1-meso and 1-rac were depicted using GaussView, while the magnetic induced current density was illustrated using the AICD program [42].

4. Conclusions

In summary, we employ quantum chemistry calculations and wave function analysis to study the HOMO and LUMO molecular orbitals, UV-Vis and Raman spectra, and the ESP and magnetic induced current density of 1-meso and 1-rac. 1-(M, M) and 1-(P, P) are obtained through the chiral separation of 1-rac, and their chiral mechanisms are investigated. The HOMO and LUMO diagrams both exhibit a similar distribution of electrons and holes, primarily concentrated around the central perylene, where S0→S1 transitions represent local excitations. The TDM and CDD plots in the spectral analysis reflect the linear optical properties of the 1-meso and 1-rac. First, the UV-Vis spectra shows that 1-meso and 1-rac have similar electronic properties, resulting in consistent electronic excitation properties. The form of electron excitation is valence shell excitation, comprising local excitation and charge transfer excitation. The IGMH analysis reveals a weak interaction between the benzene rings. The ESP maxima of 1-meso and 1-rac are evidently different, with the maximum of 1-meso being significantly higher than that of 1-rac. The magnetic induced current density demonstrates strong electron delocalization in both 1-meso and 1-rac. Raman spectra analysis shows that lasers can significantly enhance the Raman intensity of 1-meso and 1-rac. Furthermore, chiral separation is employed to divide 1-rac into 1-(M, M) and 1-(P, P) structures. These resulting structures are mirror images of each other, exhibiting a stronger level of chirality compared to 1-meso. The analysis of the TEDM and TMDM reveals that the positive and negative phases of 1-(M, M) and 1-(P, P) are opposite in the X, Y and Z directions, indicating the presence of chiral symmetry between the two. This research provides theoretical insights into the photophysical properties of other chiral nanographene materials, potentially guiding future studies in this field.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/molecules29051053/s1, Figure S1: IGMH isosurfaces of 1-meso and 1-rac (a, b), and their side views (c, d). Figure S2: UV-vis spectra of 1-meso (a) and 1-rac (b). Figure S3: TDM and CDD diagrams of 1-rac under S1 (a), S9 (b), S21 (c), S31 (d), S49 (e), S100 (f). In CDD diagram, blue represents holes and red represents electrons. The isovalue is 0.0005. Figure S4: Magnetic induced current density of 1-meso. Figure S5: Magnetic induced current density of 1-rac. Figure S6: Structure diagram of 1- (M, M) and 1- (P, P). Figure S7: TEDM and TMDM of 1-meso in S1 excited state. Figure S8: TEDM and TMDM of 1-meso in S10 excited state. Figure S9: TEDM and TMDM of 1-meso in S43 excited state. Figure S10: TEDM and TMDM of 1-meso in S47 excited state. Figure S11: TEDM and TMDM of 1-meso in S100 excited state. Figure S12: TEDM and TMDM of 1- (P, P) in S1 excited state. Figure S13: TEDM and TMDM of 1- (P, P) in S10 excited state. Figure S14: TEDM and TMDM of 1- (P, P) in S94 excited state. Figure S15: TEDM and TMDM of 1- (M, M) in S1 excited state. Figure S16: TEDM and TMDM of 1- (M, M) in S10 excited state. Figure S17: TEDM and TMDM of 1- (M, M) in S94 excited state.

Author Contributions

Conceptualization, Z.Y.; methodology, X.G.; software, X.G.; formal analysis, Y.Z.; investigation, Z.Y.; writing—original draft preparation, Y.Z. and Y.J.; visualization, Y.J. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Scientific Research Fund Project of Education Department of Liaoning province (No. L2019028, L2017LQN037, JYTMS20231431, JYTMS20231446).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are contained within the article and Supplementary Materials.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Nakada, K.; Fujita, M.; Dresselhaus, G.; Dresselhaus, M.S. Edge state in graphene ribbons: Nanometer size effect and edge shape dependence. Phys. Rev. B 1996, 54, 17954. [Google Scholar] [CrossRef]
  2. Liang, T.; Jiang, H.; Yang, Z.; Yang, Z.; Bai, H.; Shangguan, Y.; Wang, H. Research progress of graphene used in national defense. In Proceedings of the Eighth International Conference on Energy Materials and Electrical Engineering (ICEMEE 2022), Guangzhou, China, 25 April 2023; SPIE: Bellingham, WA, USA, 2023; Volume 12598, pp. 81–86. [Google Scholar]
  3. Zhang, J.; Li, J. Status and Prospect of Application about the New Materials Graphene. In MATEC Web of Conferences, Proceedings of the 13th Global Congress on Manufacturing and Management (GCMM 2016), Zhengzhou, China, 28–30 November 2016; EDP Sciences: Les Ulis, France, 2017; Volume 100, p. 04029. [Google Scholar]
  4. Kim, J.; Lee, N.; Min, Y.H.; Noh, S.; Kim, N.K.; Jung, S.; Joo, M.; Yamada, Y. Distinguishing zigzag and armchair edges on graphene nanoribbons by X-ray photoelectron and Raman spectroscopies. ACS Omega 2018, 3, 17789–17796. [Google Scholar] [CrossRef]
  5. Geim, A.K. Graphene: Status and prospects. Science 2009, 324, 1530–1534. [Google Scholar] [CrossRef]
  6. Li, X.; Yu, J.; Wageh, S.; Al-Ghamdi, A.A.; Xie, J. Graphene in photocatalysis: A review. Small 2016, 12, 6640–6696. [Google Scholar] [CrossRef] [PubMed]
  7. Grigorenko, A.N.; Polini, M.; Novoselov, K.S. Graphene plasmonics. Nat. Photonics 2012, 6, 749–758. [Google Scholar] [CrossRef]
  8. Wu, J.; Lin, H.; Moss, D.J.; Loh, K.P.; Jia, B. Graphene oxide for photonics, electronics and optoelectronics. Nat. Rev. Chem. 2023, 7, 162–183. [Google Scholar] [CrossRef] [PubMed]
  9. Shu, J.C.; Cao, M.S.; Zhang, Y.L.; Wang, Y.Z.; Zhao, Q.L.; Fang, X.Y.; Yang, S.H.; Qin, Y.; Yuan, J. Atomic-molecular engineering tailoring graphene microlaminates to tune multifunctional antennas. Adv. Funct. Mater. 2023, 2212379. [Google Scholar] [CrossRef]
  10. Razaq, A.; Bibi, F.; Zheng, X.; Papadakis, R.; Jafri SH, M.; Li, H. Review on graphene-, graphene oxide-, reduced graphene oxide-based flexible composites: From fabrication to applications. Materials 2022, 15, 1012. [Google Scholar] [CrossRef]
  11. Alam, S.N.; Sharma, N.; Kumar, L. Synthesis of graphene oxide (GO) by modified hummers method and its thermal reduction to obtain reduced graphene oxide (rGO). Graphene 2017, 6, 1–18. [Google Scholar] [CrossRef]
  12. Habte, A.T.; Ayele, D.W. Synthesis and characterization of reduced graphene oxide (rGO) started from graphene oxide (GO) using the tour method with different parameters. Adv. Mater. Sci. Eng. 2019, 2019, 5058163. [Google Scholar] [CrossRef]
  13. Konios, D.; Stylianakis, M.M.; Stratakis, E.; Kymakis, E. Dispersion behaviour of graphene oxide and reduced graphene oxide. J. Colloid Interface Sci. 2014, 430, 108–112. [Google Scholar] [CrossRef]
  14. Zhang, S.; Pei, Y.; Zhao, Z.; Guan, C.; Wu, G. Simultaneous manipulation of polarization relaxation and conductivity toward self-repairing reduced graphene oxide based ternary hybrids for efficient electromagnetic wave absorption. J. Colloid Interface Sci. 2023, 630, 453–464. [Google Scholar] [CrossRef]
  15. Rajapaksha, P.; Orrell-Trigg, R.; Shah, D.; Cheeseman, S.; Vu, K.B.; Ngo, S.T.; Murdoch, B.J.; Choudhury, N.R.; Yin, H.; Cozzolino, D.; et al. Broad spectrum antibacterial zinc oxide-reduced graphene oxide nanocomposite for water depollution. Mater. Today Chem. 2023, 27, 101242. [Google Scholar] [CrossRef]
  16. Novoselov, K.S.; Geim, A.K.; Morozov, S.V.; Jiang, D.; Zhang, Y.; Dubonos, S.V.; Grigorieva, I.V.; Firsov, A.A. Electric field effect in atomically thin carbon films. Science 2004, 306, 666–669. [Google Scholar] [CrossRef]
  17. Novoselov, K.S.; Geim, A.K.; Morozov, S.V.; Jiang, D.; Katsnelson, M.I.; Grigorieva, I.V.; Dubonos, S.V.; Firsov, A.A. Two-dimensional gas of massless Dirac fermions in graphene. Nature 2005, 438, 197–200. [Google Scholar] [CrossRef]
  18. Neto, A.C.; Guinea, F.; Peres, N.M.; Novoselov, K.S.; Geim, A.K. The electronic properties of graphene. Rev. Mod. Phys. 2009, 81, 109. [Google Scholar] [CrossRef]
  19. Ares, P.; Novoselov, K.S. Recent advances in graphene and other 2D materials. Nano Mater. Sci. 2022, 4, 3–9. [Google Scholar] [CrossRef]
  20. Kumar, N.; Salehiyan, R.; Chauke, V.; Botlhoko, O.J.; Setshedi, K.; Scriba, M.; Masukume, M.; Ray, S.S. Top-down synthesis of graphene: A comprehensive review. FlatChem 2021, 27, 100224. [Google Scholar] [CrossRef]
  21. Tarcan, R.; Todor-Boer, O.; Petrovai, I.; Leordean, C.; Astilean, S.; Botiz, I. Reduced graphene oxide today. J. Mater. Chem. C 2020, 8, 1198–1224. [Google Scholar] [CrossRef]
  22. Balandin, A.A. Phononics of graphene and related materials. ACS Nano 2020, 14, 5170–5178. [Google Scholar] [CrossRef] [PubMed]
  23. Olabi, A.G.; Abdelkareem, M.A.; Wilberforce, T.; Sayed, E.T. Application of graphene in energy storage device—A review. Renew. Sustain. Energy Rev. 2021, 135, 110026. [Google Scholar] [CrossRef]
  24. Syama, S.; Mohanan, P.V. Comprehensive application of graphene: Emphasis on biomedical concerns. Nano-Micro Lett. 2019, 11, 6. [Google Scholar] [CrossRef] [PubMed]
  25. Bonner, W.A. The origin and amplification of biomolecular chirality. Orig. Life Evol. Biosph. 1991, 21, 59–111. [Google Scholar] [CrossRef]
  26. McNaught, A.D. Compendium of Chemical Terminology; Blackwell Science: Oxford, UK, 1997. [Google Scholar]
  27. Li, J.K.; Chen, X.Y.; Zhao, W.L.; Guo, Y.; Zhang, Y.; Wang, X.; Sue, A.C.; Cao, X.; Li, M.; Chen, C.; et al. Synthesis of Highly Luminescent Chiral Nanographene. Angew. Chem. 2023, 135, e202215367. [Google Scholar] [CrossRef]
  28. Baughman, R.H.; Eckhardt, H.; Kertesz, M. Structure-property predictions for new planar forms of carbon: Layered phases containing sp2 and sp atoms. J. Chem. Phys. 1987, 87, 6687–6699. [Google Scholar] [CrossRef]
  29. Frisch, M.J.; Trucks, G.W.; Schlegel, H.B.; Scuseria, G.E.; Robb, M.A.; Cheeseman, J.R.; Scalmani, G.; Barone, V.; Petersson, G.A.; Nakatsuji, H.; et al. Gaussian 16, Revision A.03; Gaussian, Inc.: Wallingford, CT, USA, 2016. [Google Scholar]
  30. Kohn, W.; Sham, L.J. Self-consistent equations including exchange and correlation effects. Phys. Rev. 1965, 140, A1133. [Google Scholar] [CrossRef]
  31. Becke, A.D. Density-functional thermochemistry. IV. A new dynamical correlation functional and implications for exact-exchange mixing. J. Chem. Phys. 1996, 104, 1040–1046. [Google Scholar] [CrossRef]
  32. Rassolov, V.A.; Pople, J.A.; Ratner, M.A.; Windus, T.L. 6-31G* basis set for atoms K through Zn. J. Chem. Phys. 1998, 109, 1223–1229. [Google Scholar] [CrossRef]
  33. Grimme, S. Density functional theory with London dispersion corrections. Wiley Interdiscip. Rev. Comput. Mol. Sci. 2011, 1, 211–228. [Google Scholar] [CrossRef]
  34. Macernis, M.; Mourokh, L.; Vinciunas, L.; Valkunas, L. Mechanism of proton transfer in bacteriorhodopsin. Lith. J. Phys. 2022, 62, 179–193. [Google Scholar] [CrossRef]
  35. Mačernis, M. DFT functional analysis for the modeling Raman bands and absorption spectra of the lycopene structure. Open Access J. Chem. 2018, 2, 1–10. [Google Scholar] [CrossRef]
  36. Lu, T.; Chen, F. Multiwfn: A multifunctional wavefunction analyzer. J. Comput. Chem. 2012, 33, 580–592. [Google Scholar] [CrossRef]
  37. Liu, Z.; Lu, T.; Chen, Q. Vibrational Spectra and Molecular Vibrational Behaviors of All-Carboatomic Rings, cyclo [18] carbon and Its Analogues. Chem.–Asian J. 2021, 16, 56–63. [Google Scholar] [CrossRef]
  38. Humphrey, W.; Dalke, A.; Schulten, K. VMD: Visual molecular dynamics. J. Mol. Graph. 1996, 14, 33–38. [Google Scholar] [CrossRef]
  39. Gai, X.; Sheng, H.; Wang, J. Physical mechanism on the linear spectrum and nonlinear spectrum in a twist bilayer graphdiyne nanodisk. Phys. Chem. Chem. Phys. 2023, 25, 20049–20065. [Google Scholar] [CrossRef] [PubMed]
  40. Li, N.; Zhang, L.; Wang, J. Modulation of chiral spectral deflection by van der Waals force-induced molecular electropolarization in catenane oligomers. RSC Adv. 2023, 13, 11055–11061. [Google Scholar] [CrossRef] [PubMed]
  41. Lu, C.; Chen, P.; Sheng, H.; Li, C.; Wang, J. Physical mechanism on linear spectrum and nonlinear spectrum in double helical carbon nanomolecule–infinitene. Spectrochim. Acta Part A Mol. Biomol. Spectrosc. 2022, 282, 121674. [Google Scholar] [CrossRef]
  42. Geuenich, D.; Hess, K.; Köhler, F.; Herges, R. Anisotropy of the induced current density (ACID), a general method to quantify and visualize electronic delocalization. Chem. Rev. 2005, 105, 3758–3772. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Structure diagram of 1-meso (a); structure diagram of 1-rac (b); structural side view of 1-meso (c); structural side view of 1-rac (d).
Figure 1. Structure diagram of 1-meso (a); structure diagram of 1-rac (b); structural side view of 1-meso (c); structural side view of 1-rac (d).
Molecules 29 01053 g001
Figure 2. HOMO-LUMO molecular orbitals of 1meso (a) and 1-rac (b).
Figure 2. HOMO-LUMO molecular orbitals of 1meso (a) and 1-rac (b).
Molecules 29 01053 g002
Figure 3. Combined absorption spectra of 1-meso and 1-rac.
Figure 3. Combined absorption spectra of 1-meso and 1-rac.
Molecules 29 01053 g003
Figure 4. TDM and CDD diagrams of 1-meso under S1 (a), S9 (b), S21 (c), S30 (d), S47 (e), S100 (f). In CDD diagram, blue represents holes and red represents electrons. The isovalue is 0.0005.
Figure 4. TDM and CDD diagrams of 1-meso under S1 (a), S9 (b), S21 (c), S30 (d), S47 (e), S100 (f). In CDD diagram, blue represents holes and red represents electrons. The isovalue is 0.0005.
Molecules 29 01053 g004
Figure 5. ESP isosurface diagram of 1-meso (a) and 1-rac (b).
Figure 5. ESP isosurface diagram of 1-meso (a) and 1-rac (b).
Molecules 29 01053 g005
Figure 6. Magnetic induced current densities of 1-meso (a) and 1-rac (b).
Figure 6. Magnetic induced current densities of 1-meso (a) and 1-rac (b).
Molecules 29 01053 g006
Figure 7. Raman and resonance Raman spectra of 1-meso (a) and 1-rac (b).
Figure 7. Raman and resonance Raman spectra of 1-meso (a) and 1-rac (b).
Molecules 29 01053 g007
Figure 8. (a) The Vibration mode at position 1578.29 cm−1. (b) The vibration mode at position 1595.75 cm−1. (c) The vibration mode at 1600.64 cm−1 positions under laser irradiation at 325 nm. (d) The vibration mode at 1629.95 cm−1 positions under 325 nm laser irradiation. (e) The vibration mode at 1451.01 cm−1 positions under 532 nm laser irradiation. (f) The vibration mode at 1600.64 cm−1 positions under 532 nm laser irradiation.
Figure 8. (a) The Vibration mode at position 1578.29 cm−1. (b) The vibration mode at position 1595.75 cm−1. (c) The vibration mode at 1600.64 cm−1 positions under laser irradiation at 325 nm. (d) The vibration mode at 1629.95 cm−1 positions under 325 nm laser irradiation. (e) The vibration mode at 1451.01 cm−1 positions under 532 nm laser irradiation. (f) The vibration mode at 1600.64 cm−1 positions under 532 nm laser irradiation.
Molecules 29 01053 g008
Figure 9. (a) The Vibration mode at position 1581.80 cm−1. (b) The vibration mode at position 1601.47 cm−1. (c) The vibration mode at 1514.38 cm−1 positions under laser irradiation at 325 nm. (d) The vibration mode at 1601.47 cm−1 positions under 325 nm laser irradiation. (e) The vibration mode at 1373.28 cm−1 positions under 532 nm laser irradiation. (f) The vibration mode at 1601.47 cm−1 positions under 532 nm laser irradiation.
Figure 9. (a) The Vibration mode at position 1581.80 cm−1. (b) The vibration mode at position 1601.47 cm−1. (c) The vibration mode at 1514.38 cm−1 positions under laser irradiation at 325 nm. (d) The vibration mode at 1601.47 cm−1 positions under 325 nm laser irradiation. (e) The vibration mode at 1373.28 cm−1 positions under 532 nm laser irradiation. (f) The vibration mode at 1601.47 cm−1 positions under 532 nm laser irradiation.
Molecules 29 01053 g009
Figure 10. ECD spectra of 1-meso (a), 1-(P, P) (b), 1-(M, M) (c) and their combination diagram (d).
Figure 10. ECD spectra of 1-meso (a), 1-(P, P) (b), 1-(M, M) (c) and their combination diagram (d).
Molecules 29 01053 g010
Figure 11. 1-meso and S8 (a), 1-(P, P) and S100 (b) and 1-(M, M) and S100 (c) are TEDMs\TMDMs in different directions. The green (orange) isosurfaces represent the positive (negative) TEDM, and yellow (purple) isosurfaces represent the positive (negative) TMDM. X, Y, Z represent the direction.
Figure 11. 1-meso and S8 (a), 1-(P, P) and S100 (b) and 1-(M, M) and S100 (c) are TEDMs\TMDMs in different directions. The green (orange) isosurfaces represent the positive (negative) TEDM, and yellow (purple) isosurfaces represent the positive (negative) TMDM. X, Y, Z represent the direction.
Molecules 29 01053 g011
Table 1. The value of the transition electric/magnetic dipole moment and the eigenvalue of their tensor product.
Table 1. The value of the transition electric/magnetic dipole moment and the eigenvalue of their tensor product.
XYZEigenvalue
1-mesoTEDM0.0000−0.44170.1305−0.6414
S8TMDM0.0000−0.60162.8786
1-(P,P)TEDM−0.00032.4479−0.0000−3.6456
S100TMDM−0.00011.48930.0001
1-(M,M)TEDM−0.0005−2.44780.00013.6450
S100TMDM0.00021.48910.0002
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Yang, Z.; Gai, X.; Zou, Y.; Jiang, Y. The Physical Mechanism of Linear and Nonlinear Optical Properties of Nanographene-Induced Chiral Inversion. Molecules 2024, 29, 1053. https://doi.org/10.3390/molecules29051053

AMA Style

Yang Z, Gai X, Zou Y, Jiang Y. The Physical Mechanism of Linear and Nonlinear Optical Properties of Nanographene-Induced Chiral Inversion. Molecules. 2024; 29(5):1053. https://doi.org/10.3390/molecules29051053

Chicago/Turabian Style

Yang, Zhiyuan, Xinwen Gai, Yi Zou, and Yongjian Jiang. 2024. "The Physical Mechanism of Linear and Nonlinear Optical Properties of Nanographene-Induced Chiral Inversion" Molecules 29, no. 5: 1053. https://doi.org/10.3390/molecules29051053

APA Style

Yang, Z., Gai, X., Zou, Y., & Jiang, Y. (2024). The Physical Mechanism of Linear and Nonlinear Optical Properties of Nanographene-Induced Chiral Inversion. Molecules, 29(5), 1053. https://doi.org/10.3390/molecules29051053

Article Metrics

Back to TopTop