Electronic structure of 1,4-Phenylenediacrylic acid on graphene and bilayer graphite: from experiments to DFT and ab initio molecular dynamics simulations
{"title":"Electronic structure of 1,4-Phenylenediacrylic acid on graphene and bilayer graphite: from experiments to DFT and ab initio molecular dynamics simulations","authors":"Elaheh Mohebbi , Eleonora Pavoni , Pierluigi Stipa , Marina Petroselli , Cristina Minnelli , Luca Pierantoni , Davide Mencarelli , Martino Aldrigo , Emiliano Laudadio , Mir Masoud Seyyed Fakhrabadi","doi":"10.1016/j.chphi.2025.100947","DOIUrl":null,"url":null,"abstract":"<div><div>In this work, density functional theory (DFT) calculations and ab initio molecular dynamics (AIMD) simulations were implemented to expand the knowledge about the interaction of graphene and bilayer graphite surface with 1,4-Phenylenediacrylic acid (C<sub>12</sub>H<sub>10</sub>O<sub>4</sub>). The DFT calculations demonstrated that C<sub>12</sub>H<sub>10</sub>O<sub>4</sub> molecule has an opening band gap of 0.0062 eV at the top position over the graphene sheet higher than the cross and bridge sites with lower band gaps of 0.0050 eV and 0.0046 eV, respectively. The HOMO-LUMO splitting calculations confirmed more mixture of LUMO states of the C<sub>12</sub>H<sub>10</sub>O<sub>4</sub> and graphene in the carbon-carbon double bond in vinyl segment and the COOH functional group in the C<sub>12</sub>H<sub>10</sub>O<sub>4</sub>@Graphene (top) adsorption site. Then, the increasing of the molecule units on the graphene substrate resulted in a higher electronic band gap of 0.0068 eV and LUMO energy level of 0.9528 than 0.9383 eV for the monomer ones. The AIMD calculations were used to mimic the self-assembly process of the C<sub>12</sub>H<sub>10</sub>O<sub>4</sub> molecules on the graphene layer at room temperature, remarking high adsorption capabilities of the latter one. The imaginary and real parts of dielectric constant have been evaluated and for all cases the maximum intensity of the main first peak has been found at 2.43 THz. The results of the static part of dielectric constant showed high <em>Re</em>(<em>ω</em>) for the adsorption of C<sub>12</sub>H<sub>10</sub>O<sub>4</sub> monomers on the graphene surface, while by increasing the number of C<sub>12</sub>H<sub>10</sub>O<sub>4</sub> units <em>Re</em>(<em>ω</em>) resulted remarkably reduced. The maximum value predicted is 7817 in C<sub>12</sub>H<sub>10</sub>O<sub>4</sub>@Graphene (cross) along the in-plane xx polarization and 2747 for 4C<sub>12</sub>H<sub>10</sub>O<sub>4</sub>@Graphene along the in-plane yy directions. Finally, the adsorption of C<sub>12</sub>H<sub>10</sub>O<sub>4</sub> layer on the AB stacking bilayer graphite has been considered to simulate the experimental scanning tunnelling microscopy (STM) image of self-assembled C<sub>12</sub>H<sub>10</sub>O<sub>4</sub> on highly oriented pyrolytic graphite (HOPG) surface. The zero-band gap has been predicted since the electronic structure of graphene near the K point varies by increasing its thickness.</div></div>","PeriodicalId":9758,"journal":{"name":"Chemical Physics Impact","volume":"11 ","pages":"Article 100947"},"PeriodicalIF":4.3000,"publicationDate":"2025-09-26","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"Chemical Physics Impact","FirstCategoryId":"1085","ListUrlMain":"https://www.sciencedirect.com/science/article/pii/S2667022425001331","RegionNum":0,"RegionCategory":null,"ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"Q2","JCRName":"CHEMISTRY, PHYSICAL","Score":null,"Total":0}
引用次数: 0
Abstract
In this work, density functional theory (DFT) calculations and ab initio molecular dynamics (AIMD) simulations were implemented to expand the knowledge about the interaction of graphene and bilayer graphite surface with 1,4-Phenylenediacrylic acid (C12H10O4). The DFT calculations demonstrated that C12H10O4 molecule has an opening band gap of 0.0062 eV at the top position over the graphene sheet higher than the cross and bridge sites with lower band gaps of 0.0050 eV and 0.0046 eV, respectively. The HOMO-LUMO splitting calculations confirmed more mixture of LUMO states of the C12H10O4 and graphene in the carbon-carbon double bond in vinyl segment and the COOH functional group in the C12H10O4@Graphene (top) adsorption site. Then, the increasing of the molecule units on the graphene substrate resulted in a higher electronic band gap of 0.0068 eV and LUMO energy level of 0.9528 than 0.9383 eV for the monomer ones. The AIMD calculations were used to mimic the self-assembly process of the C12H10O4 molecules on the graphene layer at room temperature, remarking high adsorption capabilities of the latter one. The imaginary and real parts of dielectric constant have been evaluated and for all cases the maximum intensity of the main first peak has been found at 2.43 THz. The results of the static part of dielectric constant showed high Re(ω) for the adsorption of C12H10O4 monomers on the graphene surface, while by increasing the number of C12H10O4 units Re(ω) resulted remarkably reduced. The maximum value predicted is 7817 in C12H10O4@Graphene (cross) along the in-plane xx polarization and 2747 for 4C12H10O4@Graphene along the in-plane yy directions. Finally, the adsorption of C12H10O4 layer on the AB stacking bilayer graphite has been considered to simulate the experimental scanning tunnelling microscopy (STM) image of self-assembled C12H10O4 on highly oriented pyrolytic graphite (HOPG) surface. The zero-band gap has been predicted since the electronic structure of graphene near the K point varies by increasing its thickness.