Optical Absorption and Reflectivity of a Molecular Cluster of Lithium Niobate Adsorbed on a Graphene Layer

We used density functional theory to study the adsorption of a molecular cluster of lithium niobate (LiNbO3) on a graphene layer. The cluster size is about 1.2 nm, and it has 11 molecules. We optimized the cluster, and then we calculated its interaction with the optimized graphene layer. We found that the cluster is adsorbed on the graphene layer with an adsorption energy of −2.2213 eV (−0.068 eV/atom). Afterwards, we calculated the reflectivity, and the optical absorption coefficients of the system cluster-graphene and of a graphene layer alone, to make a comparison. We found large differences in the values of these properties of the new system, with respect to the corresponding ones of graphene. We performed our calculations using the general gradient approximation (GGA), and the GGA modified for van der Waals interactions in order to take into account the long-range correlations.


Introduction
In crystalline form, LiNbO 3 has many interesting properties.Among these, we can mention that it is an anisotropic, semiconductor, birefringent, photorefractive, colorless, photo elastic, ferroelectric, and piezoelectric material.Because of its properties, this compound is very appealing for many applications [1][2][3].The photorefractive materials are used in many technological applications-such as holographic memories, phase conjugators, and dynamic holographic interferometry, among many more.Besides, in photorefractive crystals, an illumination configuration is converted into a refractive index pattern [1].
On the other hand, graphene is a two-dimensional form of carbon, which is a one atom thick, and arranged in a hexagonal distribution of atoms.It is formally a semiconductor, but it conducts electricity very well [4][5][6].The optical properties of graphene are very different from those of LiNbO 3 , in particular, absorption and reflectivity as functions of frequency [7,8].
Here, we are interested in the decoration of a graphene layer with clusters of lithium niobate to combine, in two dimensions, the properties of these materials.Thus, we would have a new material in a two-dimensional geometry which includes a cluster of a photorefractive compound.As far as we know, this is the first attempt to predict the optical properties of this kind of a two-dimensional material.
Therefore, we have approached the problem of adsorbing a molecular cluster of lithium niobate on a graphene surface.Clearly, the properties of a molecular cluster of LiNbO 3 will tend to those of crystalline LiNbO 3 as the cluster increases its dimensions.Thus, the cluster size we choose should not be very small, and in this way, its properties can resemble some of those of the lithium niobate crystal.Furthermore, a nanometer-size cluster is experimentally much easier to obtain than smaller sizes.Besides, notice that as we increase the number of atoms within the cluster, the computing time for the numerical simulation is larger.Taking into account these considerations, we included 11 molecules of LiNbO 3 in the cluster, with a size of about 1.2 nm.
We have calculated the resulting optical properties like absorption and reflectivity of the new system.We found that the cluster is, indeed, adsorbed on the pristine graphene surface and that for the new system, the absorption is quite large in the infrared and visible regions with respect to the corresponding value for graphene.In the ultraviolet region, there are important changes too, especially in the position of the peaks.Something very similar happens for the reflectivity.

Method
We started by optimizing the cluster of lithium niobate.We had 11 molecules in the cluster which had a size of around 1.2 nm.Then, we explored the possibility of adsorbing this cluster on a graphene layer.In this way, we started numerical simulations of the interaction of the cluster and the graphene surface.We utilized the DFT, with the GGA.We considered Vdw-DF2 to include long-range correlation effects as van der Waals interactions, with the Quantum Espresso code [9].The norm conserving Troullier-Martins, and the Goedecker-Hartwigsen-Hutter-Teter pseudo potentials [10], were employed.We utilized the Perdew-Burke-Ernzerhof (PBE) approximation for GGA [11] We considered as valence electronic states for carbon 2s 2 p 2 , for lithium: 2s 1 , for niobium: 4d 4 5s 2 , and for oxygen: 2s 2 2p 4 .We performed non-relativistic, and non-spin polarized calculations.We considered for the plane waves cut-off energy 1000 eV.We took 50 k points within the Monkhorst-Pack scheme [12].We took 10 −7 eV for the threshold energy convergence.In order to verify the pseudopotentials, we have selected, and calculated by energy minimization, the value of the C-C bond length 1.42 Å in the graphene layer (the experimental value is 1.42 Å).In this manner, we calculated the lattice parameters for the lithium niobate crystal with rhombohedral unit cell.We obtained the following, a = 5.4940 Å, c = 13.8625Å, and α = 55.867• (the corresponding experimental values: a = 5.494 Å, c = 13.825Å, and α = 55.867• , respectively [13]).
The graphene unit cell for the calculation has 98 carbon atoms (see Figure 1).The lithium niobate cluster has 11 molecules, i.e., 11 lithium atoms, 11 niobium atoms, and 33 oxygen atoms.know, this is the first attempt to predict the optical properties of this kind of a two-dimensional material.Therefore, we have approached the problem of adsorbing a molecular cluster of lithium niobate on a graphene surface.Clearly, the properties of a molecular cluster of LiNbO3 will tend to those of crystalline LiNbO3 as the cluster increases its dimensions.Thus, the cluster size we choose should not be very small, and in this way, its properties can resemble some of those of the lithium niobate crystal.Furthermore, a nanometer-size cluster is experimentally much easier to obtain than smaller sizes.Besides, notice that as we increase the number of atoms within the cluster, the computing time for the numerical simulation is larger.Taking into account these considerations, we included 11 molecules of LiNbO3 in the cluster, with a size of about 1.2 nm.
We have calculated the resulting optical properties like absorption and reflectivity of the new system.We found that the cluster is, indeed, adsorbed on the pristine graphene surface and that for the new system, the absorption is quite large in the infrared and visible regions with respect to the corresponding value for graphene.In the ultraviolet region, there are important changes too, especially in the position of the peaks.Something very similar happens for the reflectivity.

Method
We started by optimizing the cluster of lithium niobate.We had 11 molecules in the cluster which had a size of around 1.2 nm.Then, we explored the possibility of adsorbing this cluster on a graphene layer.In this way, we started numerical simulations of the interaction of the cluster and the graphene surface.We utilized the DFT, with the GGA.We considered Vdw-DF2 to include long-range correlation effects as van der Waals interactions, with the Quantum Espresso code [9].The norm conserving Troullier-Martins, and the Goedecker-Hartwigsen-Hutter-Teter pseudo potentials [10], were employed.We utilized the Perdew-Burke-Ernzerhof (PBE) approximation for GGA [11] We considered as valence electronic states for carbon 2s 2 p 2 , for lithium: 2s 1 , for niobium: 4d 4 5s 2 , and for oxygen: 2s 2 2p 4 .We performed non-relativistic, and non-spin polarized calculations.We considered for the plane waves cut-off energy 1000 eV.We took 50 k points within the Monkhorst-Pack scheme [12].We took 10 −7 eV for the threshold energy convergence.In order to verify the pseudopotentials, we have selected, and calculated by energy minimization, the value of the C-C bond length 1.42 Å in the graphene layer (the experimental value is 1.42 Å).In this manner, we calculated the lattice parameters for the lithium niobate crystal with rhombohedral unit cell.We obtained the following, a = 5.4940 Å, c = 13.8625Å, and α = 55.867°(the corresponding experimental values: a = 5.494 Å, c = 13.825Å, and α = 55.867°,respectively [13]).
The graphene unit cell for the calculation has 98 carbon atoms (see Figure 1).The lithium niobate cluster has 11 molecules, i.e., 11 lithium atoms, 11 niobium atoms, and 33 oxygen atoms.The system graphene-cluster of LiNbO 3 was represented using a hexagonal cell with a = 17.344Å and c = 26.04Å with periodic conditions.Note that in the z direction, we took a large enough separation to ensure that there is no interaction between the adjacent cells.
To calculate the adsorption energy of the cluster on the graphene layer, we used the expression Here, the first term E (graphene + cluster) is the energy of the new system made by the cluster of lithium niobate, plus the graphene layer; the second, [E (graphene) + E (cluster)] is the energy of the graphene alone, plus the energy of the isolated cluster of lithium niobate.
We used the method of the modified Blöch tetrahedron for integration in the first Brillouin zone [14], to calculate the band energy structure.Afterwards, we calculated the dielectric tensor, the absorption and the reflectivity.These calculations were performed for graphene and for the new system made of the cluster of lithium niobate adsorbed on graphene.

Results and Discussion
From the numerical simulations of the interaction of the cluster and the pristine graphene surface, we obtained that the cluster of lithium niobate is adsorbed on the graphene layer, with an adsorption energy of −2.2213 eV (−0.068 eV/atom), as shown in Figure 1.It was not necessary to consider defects on the graphene layer to achieve the adsorption of the cluster.The adsorption energy is large enough to have stability on the resulting system.For an adsorption energy below 10 meV, we could have instability in the new system.Anyhow, we used molecular dynamics at 300 K and atmospheric pressure to test the stability, and the system did not change its configuration under these conditions.The size of the graphene unit cell (98 carbon atoms) is large enough to avoid interactions of adjacent lithium clusters.
We should mention that the reproducibility of the results is granted by the high degree of accuracy in our numerical simulations.Furthermore, we have the validation of our pseudo potentials.This was achieved by energy minimization.We could correctly reproduce the experimental values of the lattice parameters of graphene and of LiNbO 3 .
On the other hand, it is convenient to mention that DFT theory and molecular dynamics have been used successfully to understand the optical properties of crystalline lithium niobate.This was the case in understanding the variation of these properties with the deviation from the stoichiometric composition [15].
To obtain the optical absorption and reflectivity, we calculated the energy band structure of the new system and of graphene.The calculation of the complex dielectric tensor allows to obtain optical spectra.This allows us to know how the electromagnetic radiation propagates through the material.
Then, we proceed to calculate the imaginary component of the dielectric tensor in the limit of the linear optics, and within the random phase approximation (RPA) [16] where ω is the photon frequency, n is the electron density, m is the effective mass, the indices α and β are the Cartesian components, p is the vector defining the polarization of the incident electric field, c and v represent the conduction band and valence band states, respectively.On the other hand, from the Equation (2) and using the Kramers-Kronig rule [16,17], we can obtain the real component of the dielectric tensor where P denotes the principal value.We performed the calculation over a range of 0 to 30 eV (or from 0 to rom 0 to 242,000 cm −1 ).We considered the light propagation perpendicular to the graphene layer, and calculated the reflectivity from the Fresnel formula where n and k are the real and imaginary parts of the complex refractive index (n and k are the refractive index and extinction coefficient respectively) We calculated the absorption coefficient In Figure 2, we can see the calculated absorption for the system cluster of lithium niobate adsorbed on the graphene layer (dashed curve) for the infrared region.In this figure, we present the absorption of graphene too (dotted line).Notice that the new system shows a much larger absorption than the one corresponding to graphene.The absorption in this system has two large peaks, one around 9.4 × 10 3 cm −1 , and the other around 11.25 × 10 3 cm −1 .The absorption for graphene in this region is almost zero.
where n and k are the real and imaginary parts of the complex refractive index (n and k are the refractive index and extinction coefficient respectively) We calculated the absorption coefficient In Figure 2, we can see the calculated absorption for the system cluster of lithium niobate adsorbed on the graphene layer (dashed curve) for the infrared region.In this figure, we present the absorption of graphene too (dotted line).Notice that the new system shows a much larger absorption than the one corresponding to graphene.The absorption in this system has two large peaks, one around 9.4 × 10 3 cm −1 , and the other around 11.25 × 10 3 cm −1 .The absorption for graphene in this region is almost zero.In Figure 3, we show the absorption for the new system (dashed curve) and for graphene (dotted line) in the region of visible light.The new system shows the largest absorption, which increases with increasing frequency.We notice three large peaks.The first one is around 1.2 × 10 4 cm −1 , the second is located around 2.04 × 10 4 cm −1 , and the third peak is at around 2.5 × 10 4 cm −1 .In Figure 3, we show the absorption for the new system (dashed curve) and for graphene (dotted line) in the region of visible light.The new system shows the largest absorption, which increases with increasing frequency.We notice three large peaks.The first one is around 1.2 × 10 4 cm −1 , the second is located around 2.04 × 10 4 cm −1 , and the third peak is at around 2.5 × 10 4 cm −1 .In Figure 4, we present again the absorption for the same systems (dashed curve for the new system), but for the ultraviolet region of the spectrum.We can notice that the graphene layer presents two peaks.The first one is located around 0.9 × 10 5 cm −1 , and the second is at around 1.15 × 10 5 cm −1 .
The new system has the largest absorption between 0.25 × 10 5 cm −1 and 0.75 × 10 5 cm −1 .The second peak of the new system is located at the same position than the second peak corresponding to the graphene layer.Besides, the absorption has a very similar value to the one corresponding with the graphene layer.Notice that there is a small third peak for the absorption of the new system located within a region where the corresponding absorption of the graphene layer is near zero.It is located around 1.6 × 10 5 cm −1 .Finally, in Figure 5, we show the reflectivity for the same systems for the infrared, visible, and ultraviolet regions of the spectrum.The new system (dashed curve) has the largest reflectivity from zero to around 0.8 × 10 5 cm −1 .There is a second small peak for the absorption of the new system, located at almost the same position of the second peak of the graphene layer reflectivity (around 0.9 × 10 5 cm −1 ).The second peak of the reflectivity of the new system is located at the same position than the corresponding one to the graphene layer (around 1.2 × 10 5 cm −1 ).Its value is very similar to the one corresponding to the graphene layer.In Figure 4, we present again the absorption for the same systems (dashed curve for the new system), but for the ultraviolet region of the spectrum.We can notice that the graphene layer presents two peaks.The first one is located around 0.9 × 10 5 cm −1 , and the second is at around 1.15 × 10 5 cm −1 .
The new system has the largest absorption between 0.25 × 10 5 cm −1 and 0.75 × 10 5 cm −1 .The second peak of the new system is located at the same position than the second peak corresponding to the graphene layer.Besides, the absorption has a very similar value to the one corresponding with the graphene layer.Notice that there is a small third peak for the absorption of the new system located within a region where the corresponding absorption of the graphene layer is near zero.It is located around 1.6 × 10 5 cm −1 .Figure 3. Absorption for the new system (dashed curve), and for graphene (dotted line) in the region of visible light.
In Figure 4, we present again the absorption for the same systems (dashed curve for the new system), but for the ultraviolet region of the spectrum.We can notice that the graphene layer presents two peaks.The first one is located around 0.9 × 10 5 cm −1 , and the second is at around 1.15 × 10 5 cm −1 .
The new system has the largest absorption between 0.25 × 10 5 cm −1 and 0.75 × 10 5 cm −1 .The second peak of the new system is located at the same position than the second peak corresponding to the graphene layer.Besides, the absorption has a very similar value to the one corresponding with the graphene layer.Notice that there is a small third peak for the absorption of the new system located within a region where the corresponding absorption of the graphene layer is near zero.It is located around 1.6 × 10 5 cm −1 .Finally, in Figure 5, we show the reflectivity for the same systems for the infrared, visible, and ultraviolet regions of the spectrum.The new system (dashed curve) has the largest reflectivity from zero to around 0.8 × 10 5 cm −1 .There is a second small peak for the absorption of the new system, located at almost the same position of the second peak of the graphene layer reflectivity (around 0.9 × 10 5 cm −1 ).The second peak of the reflectivity of the new system is located at the same position than the corresponding one to the graphene layer (around 1.2 × 10 5 cm −1 ).Its value is very similar to the one corresponding to the graphene layer.Finally, in Figure 5, we show the reflectivity for the same systems for the infrared, visible, and ultraviolet regions of the spectrum.The new system (dashed curve) has the largest reflectivity from zero to around 0.8 × 10 5 cm −1 .There is a second small peak for the absorption of the new system, located at almost the same position of the second peak of the graphene layer reflectivity (around 0.9 × 10 5 cm −1 ).The second peak of the reflectivity of the new system is located at the same position than the corresponding one to the graphene layer (around 1.2 × 10 5 cm −1 ).Its value is very similar to the one corresponding to the graphene layer.We must mention that the energy band structure for the pristine graphene layer corresponds to an atypical semiconductor.The size of the energy gap for the graphene layer is zero.This characteristic, and its two-dimensional configuration, are the origin of the unusual electronic properties of this material [4].Thus, graphene has an optical absorption near zero in the regions of the infrared and the visible radiation.It has two large peaks in the ultraviolet region, as we have presented in the Figures 2-4.

Conclusions
We found that the cluster of lithium niobate can be adsorbed on a pristine graphene layer, with an adsorption energy of −2.2213 eV (−0.068 eV/atom).The cluster is made of 11 molecules containing a total of 11 lithium atoms, 11 niobium atoms, and 33 oxygen atoms.
The resulting system, with the lithium niobate cluster adsorbed on the graphene surface has different electronic properties.Thus, new possibilities for the adsorption of photons appear, and others disappear.In this way, the optical absorption of the new system eliminates a peak from the absorption of the graphene layer in the ultraviolet region.It is some kind of a filter in this region of the electromagnetic spectrum.Besides, the new system absorbs much more energy in the infrared and visible regions of the spectrum.
The reflectivity of the new system changes the corresponding one to the graphene layer.It eliminates one of the two peaks shown by the graphene layer in the ultraviolet region of the spectrum.Furthermore, the reflectivity in the infrared and visible regions are much increased for the new system.
Finally, we mention that we have combined the properties of pristine graphene and a molecular cluster of lithium niobate, which is a photorefractive compound.We have calculated the optical absorption and reflectivity of this possible, novel, two-dimensional material.We must mention that the energy band structure for the pristine graphene layer corresponds to an atypical semiconductor.The size of the energy gap for the graphene layer is zero.This characteristic, and its two-dimensional configuration, are the origin of the unusual electronic properties of this material [4].Thus, graphene has an optical absorption near zero in the regions of the infrared and the visible radiation.It has two large peaks in the ultraviolet region, as we have presented in the Figures 2-4.

Conclusions
We found that the cluster of lithium niobate can be adsorbed on a pristine graphene layer, with an adsorption energy of −2.2213 eV (−0.068 eV/atom).The cluster is made of 11 molecules containing a total of 11 lithium atoms, 11 niobium atoms, and 33 oxygen atoms.
The resulting system, with the lithium niobate cluster adsorbed on the graphene surface has different electronic properties.Thus, new possibilities for the adsorption of photons appear, and others disappear.In this way, the optical absorption of the new system eliminates a peak from the absorption of the graphene layer in the ultraviolet region.It is some kind of a filter in this region of the electromagnetic spectrum.Besides, the new system absorbs much more energy in the infrared and visible regions of the spectrum.
The reflectivity of the new system changes the corresponding one to the graphene layer.It eliminates one of the two peaks shown by the graphene layer in the ultraviolet region of the spectrum.Furthermore, the reflectivity in the infrared and visible regions are much increased for the new system.
Finally, we mention that we have combined the properties of pristine graphene and a molecular cluster of lithium niobate, which is a photorefractive compound.We have calculated the optical absorption and reflectivity of this possible, novel, two-dimensional material.

Figure 1 .
Figure 1.We show the unit cell with the final configuration of the lithium cluster adsorbed on the graphene layer.The cluster has 11 molecules.This means 11 lithium atoms, 11 niobium atoms, and 33 oxygen atoms.The unit cell of graphene has 98 carbon atoms.

Figure 1 .
Figure 1.We show the unit cell with the final configuration of the lithium cluster adsorbed on the graphene layer.The cluster has 11 molecules.This means 11 lithium atoms, 11 niobium atoms, and 33 oxygen atoms.The unit cell of graphene has 98 carbon atoms.

Figure 2 .
Figure 2. Optical absorption for the system cluster of lithium niobate adsorbed on the graphene layer (dashed curve) for the infrared region.Optical absorption of graphene (dotted line).

Figure 2 .
Figure 2. Optical absorption for the system cluster of lithium niobate adsorbed on the graphene layer (dashed curve) for the infrared region.Optical absorption of graphene (dotted line).

Figure 3 .
Figure 3. Absorption for the new system (dashed curve), and for graphene (dotted line) in the region of visible light.

Figure 4 .
Figure 4. Absorption for the new system (dashed curve), and for the graphene layer (dotted line), for the ultraviolet region of the spectrum.

Figure 3 .
Figure 3. Absorption for the new system (dashed curve), and for graphene (dotted line) in the region of visible light.

Figure 4 .
Figure 4. Absorption for the new system (dashed curve), and for the graphene layer (dotted line), for the ultraviolet region of the spectrum.

Figure 4 .
Figure 4. Absorption for the new system (dashed curve), and for the graphene layer (dotted line), for the ultraviolet region of the spectrum.

Figure 5 .
Figure 5. Reflectivity for the three regions of the spectrum.For the new system (dashed line), and for the graphene layer (dotted line).

Figure 5 .
Figure 5. Reflectivity for the three regions of the spectrum.For the new system (dashed line), and for the graphene layer (dotted line).