Next Article in Journal
Electrorheological Characteristics of Poly(diphenylamine)/magnetite Composite-Based Suspension
Next Article in Special Issue
The Effect of Surface Treatment on Structural Properties of CVD Diamond Layers with Different Grain Sizes Studied by Raman Spectroscopy
Previous Article in Journal
Hysteretic Behavior of Random Particulate Composites by the Stochastic Finite Element Method
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Ligand-Modified Boron-Doped Diamond Surface: DFT Insights into the Electronic Properties of Biofunctionalization

1
Faculty of Electronics, Telecommunications and Informatics, Gdansk University of Technology, 11/12 G. Narutowicza St., 80-233 Gdansk, Poland
2
Materials and Process Simulation Center, California Institute of Technology, 1200 East California Blvd., Pasadena, CA 91125, USA
*
Author to whom correspondence should be addressed.
Materials 2019, 12(18), 2910; https://doi.org/10.3390/ma12182910
Submission received: 2 August 2019 / Revised: 2 September 2019 / Accepted: 6 September 2019 / Published: 9 September 2019

Abstract

:
With the increasing power of computation systems, theoretical calculations provide a means for quick determination of material properties, laying out a research plan, and lowering material development costs. One of the most common is Density Functional Theory (DFT), which allows us to simulate the structure of chemical molecules or crystals and their interaction. In developing a new generation of biosensors, understanding the nature of functional linkers, antibodies, and ligands become essential. In this study, we used DFT to model a bulk boron-doped diamond slab, modified by a functional linker and a surrogate proteins ligand. DTF calculations enable the prediction of electronic transport properties in an electrochemical sensor setup, composed of a boron-doped diamond electrode functionalized by 4-amino benzoic acids and a target surrogated protein-ligand for influenza. Electron conduction pathways and other signatures associated with the detection and measurement of the target analyte are revealed.

1. Introduction

Computational power continues to increase significantly year to year, and by today’s standards, it has enabled rapid simulations, before physical experiments. This helps drive the development of high-sensitivity, high-selectivity biosensing electrochemical surfaces [1]. This involves addressing many electronic structure challenges, including surface activation for improved sensing selectivity, determination of the electrochemical windows, optimisation of surface and linker conductivities, and characterisation of analyte conductivities, among others.
Density Functional Theory (DFT) is an approximation to Schrodinger’s electronic wave equation in quantum mechanics that is used to study multi-body electronic-nuclear systems at their ground state, using functionals of the electron density. DFT is considered relatively cheap in comparison to traditional methods such as Hartree-Fock theory [2], because it directly associates a system’s ground state energy with the electron density, thereby significantly reducing the number of degrees of freedom. DFT is a powerful tool for studying nuclear fuel materials [3], electronic structure and optical properties [4], doping [5] and surface functionalization [6,7,8], among many other materials’ properties and phenomena.
Thin diamond layers are one of the most promising material for developing new electrochemical sensors [9,10,11]. These layers have high chemical resistance, biological compatibility [12] and a wide electrochemical window. These properties have enabled their broader use in electronics, after appropriate doping [13,14,15]. The most popular and well-known dopant for a diamond is boron, which leads to a p-type semiconductor or even semi-metallic behaviour [16,17]. Furthermore, the diamond surface can be easily functionalized via organic compounds to improve their selectivity and sensitivity in biosensors [18,19].
Many authors have used DFT to study surface functionalization, including Maron and Eisenstein [6,7], who investigated ligands in lanthanide chemistry and tried to answer the question about the role of f-electrons in lanthanide-ligand bonds. Boukhvalow [20] have simulated the covalent functionalization of graphene, concluding that the formation of uniform one- or two-sided functionalized graphene is an exception to the typical scenarios of graphene chemistry. Lastly, Jaramillo-Botero et al. [21] characterised non-covalent ligands on graphene nanoribbons for UDP-glucose sensors and showed the ligands do not affect the material’s transmission spectrum. Recently, Tian and Larrson reported on the adhesion energies of important biomolecules at the variously terminated (111) diamond surfaces [22].
In our recent work [23], we show utilisation of DFT to study of the refractive index of BDD resulting with good agreement with experiments carried out by spectroscopic ellipsometry [24]. The measured values were 2.55–3.37 for a wavelength in the range from 200 nm up to 1600 nm, an error less than 2.5% when compared to bulk diamond. Furthermore, we used a DFT calculation to show electrical and optical properties of allylamine [25]. Achieved results show that with proper functionals and pseudopotentials lead to accurate characterisation of bulk materials.
In this paper, we investigate the energetic stability and electronic properties of heavily boron-doped diamond (BDD) surfaces, modified by biofunctional ligands that attach to the diamond surface through aminobenzoic linkers. The former, considered here as a surrogate for a broad group of target proteins. DFT [26] with a localised basis set approach [27] was used to determine the electronic and transport properties of the BDD bio-surface, along with a generalised gradient approximation (GGA) exchange-correlation energy correction.

2. Materials and Methods

We used the QuantumATK software (version 2019.03) [28,29] to calculate the electronic and transport properties of heavily boron-doped diamond electrodes (BDDE) modified by binding of different surface ligands. We assume surface modifications on a pure (111) BDD surface, considering CVD fabricated BDDs are dominated by this crystallographic phase [30] (see Figure 1A–C). The bulk diamond and hydrogen-terminated diamond (111) surfaces were studied for reference and comparison. The diamond was doped with two boron atoms, placed within the third carbon layer from the top surface to provide an acceptor concentration of 3 × 1021 cm−3 (~2%), as described in our earlier studies [31]. This doping leads to a semiconductor material, with a bandgap of 4.39 eV.
We study the electron-donating and accepting capabilities of a surrogate protein-ligand that is part of the influenza A virus M1 antibody (BDDE-H-ABA-C3M) attached through 4-aminobenzoic acid (4-ABA) to the BDDE surface, as shown in Figure 1C. To characterise this system, we start from a diamond (111) 4 × 4 super-cell slab (Figure 1A). Then we create a hydrogenated diamond surface and remove one of hydrogens from it. Then we attach the 4-ABA linker (Figure 1B), equivalent to an amine cation radical formation of ABA grafted on the carbon-based electrode, and then we covalently attach the protein surrogate ligand to the other end of the 4-ABA (Figure 1C).
The model system is constructed as a 2-electrode configuration. The electrode region is a repetition of the left or right end of the slab. In all configurations, the left electrode (bottom) has a length of 10.09 Å. In the BDDE, the right electrode has length 10.09 Å. The right electrode length is 5.92 Å. For calculating transmission properties, we use Landauer-Buttiker Formalism (LBF) on an asymmetric electrode configuration. The proposed approach significantly reduces the total number of atoms in the model system, hence cost of calculations. This enables the simulation of a larger ligand surrogate.
Periodic boundary conditions were used in a and b directions, and the Dirichlet boundary condition was used in the transport direction of the BDDE, c.
Electronic properties simulations were done using the generalised gradient approximation (GGA) with Optimized Norm-Conserving Vanderbilt (SG15) pseudopotentials [32], a linear combination of atomic orbitals [33], a Perdew-Burke-Ernzerhof (PBE) [34] exchange-correlation correction with an energy cut-off of 185 Hartree. Calculations were carried out with k-point sampling at 3x3x101 points in Monkhorst-Pack grid [35]. For optimal comparison of electrical properties, all the slabs have the same A and B cell size.
For the convergence criteria of electronic-state calculations, the total energy change was set to less than 1.0 × 10−5 eV/atom. Calculations were carried out with k-point sampling at 3 × 3 × 101 points in Monkhorst-Pack grid [35]. For optimal comparison of electrical properties, all the slabs have the same A and B cell size.
Calculations of atomic forces used the LCAO local density approximation (LDA) with Fritz-Haber-Institute (FHI) [36] pseudopotentials and the DoubleZetaPolarized (DPZ) basis set. For optimisation of atom positions, we used the Limited Memory Broyden–Fletcher–Goldfarb–Shanno (LBFGS) optimisation algorithm [37]. The convergence of atomic forces was less than 0.01 eVÅ−1 with atomic energy change less than 1.0 × 10−5 eV/atom. Additionally, a second run with GGA and DFT-D3 dispersion correction included Van der Walls interactions, which are essential in calculating interactions between such small molecules.
The GGA-PBE algorithm was chosen to describe the electronic structure of BDD [38,39], while the delta-test described by Lejaeghre et al. [40] was used to select the individual atomic potentials, summarised in Table 1. The delta-test results compare the different implementation of pseudopotentials in Quantum ATK to state of the art all-electron pseudopotentials used in WIEN2K code [41].
The SG15 and the DZP pseudopotentials lead to significant differences in the delta-tests [42]. The SG15 medium basis set (with ×2 times higher computational costs when compared to FHI-DZP), results in improved accuracy of almost a factor of 2 in the delta-tests. Further improvement was possible using the SG15 High basis, nearly five times better delta-test (six times higher computational time cost than FHI-DZP, and significantly higher storage requirements). Consequently, we selected SG15-medium as the best cost-performance option for all the results reported in the following sections.

3. Results and Discussion

3.1. Electron Density Distribution at the BDDE-Analyte Junction

The iso-surfaces of the electron density differences obtained from a ground state DFT calculation and a constrained DFT calculation are illustrated in Figure 2. The variation in the electron density in the electrode structures is represented by the violet clouds in 2D and 3D view. The highest probability occurs in an electron close to carbon atoms. The boron doping—acceptor in diamond results in p-type semiconducting [43], which is marked in Figure 2A,B by a blue circular region.
Furthermore, the surface hydrogen termination (dangling-bond) shows a negative electron affinity which will affect band gap [44,45]. The surface functionalization by a 4-aminobenzoic linker (BDDE-H-ABA) does not show significant electron density differences (see Figure 2B) when compared to relaxed BDDE. On the contrary, adding ligand chain effects on electron-densities resulting in higher electron densities at the ligand side. We have observed that BDDE-H-ABA-C3H modification results in the higher electron densities when compared with the pristine BDDE. Thus, the ligand of the anti-M1 antibody act as n-type semiconductor due to functional groups rich in п-bonds [46,47], which make it an electrons source for the BDDE. Once polarised (approx. above 0.5 V), the entire structure works as a P-N junction created by the bio-interface (see Figure 2C) that can be switched and tuned by polarisation [48]. Usually, the electrons within the H-terminated diamond surface are mainly redistributed on the C-H bonds. Nevertheless, electron transfer from the H-terminated diamond surface to the grafted linker and then to ligand is also possible [49].
Takagi et al. [50] showed the mechanism of inducing hole carriers at the H-terminated diamond by physically adsorbed molecules. The electron transfer occurs when the energy level of unoccupied molecular orbitals in an adsorbate molecule is lower than the valence band maximum of H-terminated diamond. The concentration of doped hole carriers depends on the energy level of unoccupied molecular orbital of adsorbate molecules.
DFT computations can give specific information about the electrical properties of junction between the diamond electrode surface and organic compound chosen for its functionalization. This evidence reveals valuable data for designed biosensors based on current-voltage measurements, e.g., cyclic voltammetry in electrochemistry.

3.2. Electron Pathways Study during Device Operation

Electron pathways allow us to gauge the effect of the external potential on the structure’s properties. Electron pathways provide information about what happens during polarisation, e.g. changing character from insulating to semiconducting. To the best of our knowledge this approach is new, and there is lack of scientific literature.
Figure 3 shows a 2-terminal device structure with plotted electron pathways (attn. electrodes are not shown here for figure readability). In the first case (see Figure 3A), the electrons are transferred across bonds pathways where minimal energy is demanded to electron transfer and where are the highest electron densities. When we apply 1V of external potential across the BDDE terminals, the electronic properties shift significantly. Electrons do not travel exclusively along the bonding pathways, but they are also able to hop between neighbouring atoms due to the additional energy supplied. The applied potential spreads across the entire structure, in this case, distributed between −0.5V and +0.5V at the left and right electrodes, respectively.
BDDE surface functionalization with the amino-benzoic linker (BDDE-H-ABA) results in two specific electronic conduction pathways (see Figure 4), including one in which electrons travel around the benzene ring via the C-C bonds and another in which the exhibit a hopping-like behaviour directly from the C atoms in the benzene ring to O atoms in the COOH group, as shown in Figure 4A. The polarised structure of BDDE-H-ABA shown in Figure 4B further demonstrates evenly distributed pathways across the 4-ABA and the top BDDE surface layers (attn. electrodes are not shown here for figure readability). The 1 V external potential creates different electron conduction pathways, even from surface (BDDE) H atoms to H atoms located in the benzene ring. This provides evidence on electronic conduction through the 4-ABA linker.
Figure 5A of the BDDE-H-ABA-C3H slab shows that electrons prefer the left electrode (−). The right side (+), just minor electron pathways. This fact unambiguously shows that the ligand donates electrons to BDDE suppressing them from the right electrode. The polarised BDDE-H-ABA-C3M structure electron pathways show partial electron tunnelling occurs across the organic molecules and increased conductivity along the BDDE surface. This suggests that the nature of ligand varies thanks to electron tunnelling through structure and reveals a semiconducting behaviour. This effect is illustrated in Figure 5B. Cross-sectional electron tunnelling in Figure 5B exhibits efficient transport along protein bio-interface required for effective biosensing effect [51]. It should be also stated that DFT provides a reliable description of the intramolecular interactions of organic molecules ranging from small to relatively large. Though, calculating the vibrational modes and charge pathways in an organic molecular crystal is a much more complicated because conventional first-principles methods sometimes fail to describe weak intermolecular interactions adequately causing an unphysical expansion of the unit cell during geometry optimisations [48].

3.3. Device—Density of States

The total density of states (DOS) spectra for specific modified diamond structures are shown in Figure 6. It compares the bulk boron-doped diamond (BDD), and the pure diamond electrode surfaces modified by ad-layers BDDE-H-ABA and BDD-H-ABA-C3M. Similar DOS spectra indicating the formation of B-related acceptor impurity band were reported by Ashcheulov et al. [50] by using B dimer supercells configuration or by Wang at el. [52] applying boron substituted doping. Apparently, the Monte Carlo method was applied to estimate band gap of semiconducting boron-doped diamond [53,54]. The band structure and density of states in bulk diamond was simulated by full-band Monte Carlo code for electrons based on the empirical pseudopotential band calculation [54] revealing similar behaviour to DOS presented in Figure 6. Next, Watanabe at al. [53] computed transport properties of hot electrons in bulk diamond utilizing Monte Carlo. They demonstrated band gap energy of 5.47 eV close to experimental values making it competitive for DFT methods but just in case of a bulk diamond, not organic-inorganic interfaces. The deep boron levels could be observed in the DOS of BDDE within the gap, located in the range of 0.0–0.38 eV from the top of the valence band with peak located at 0.1 eV. This is in excellent agreement with other works [55] and the experimental results of 0.365 eV [56].
The BDDE functionalization by ABA and C3M results in the tuning of EF and vacuum level as observed in Figure 6. Accordingly, the band tail states are localised as reported by Dai et al. [45]. In case when there are more electrons at the BDDE than the adsorbate could accept then the electron transfer from BDDE to ABA-C3M occurs, inducing the holes in localized states of the band tail with a valence band bending in the surface layer as revealed by Dai et al. [45]. In such a case, band gap for BDDE-H-ABA and BDDE-H-ABA-C3M was shortened down to 4.01 eV and 2.6 eV respectively. The modification of the diamond surface by ABA and C3M shifts the vacuum level decreasing band gap, contrary to pure hydrogen termination BDDE where surface states are removed from the inside of the band gap [57].
For comparison, the bare BDD results here in the band gap of 4.39 eV. Computed band-gap of BDDE is in good agreement with previous work found in the literature [58,59], where GGA-PBE was applied. Thus, GGA-PBE could be utilized as one of the most cost-effective formalisms used for comparative boron-doped diamond interfaces band gap estimation. Therefore, we have applied here GGA-PBE for cost-effective qualitative trend comparisons, showing a band gap shift during the sequential steps of BDDE surface modification.
Nevertheless, to achieve better fitting to the experiment, you have to use semi-empirical corrections or hybrid functionals [60,61,62]. Hybrid functionals B3LYP, B3PW91, and HSE06, are computationally expensive and resource-consuming [63], and when it comes to empirical corrections they do not behave the same for each system and must be treated individually, which cannot be applied directly to the examined organic interface which have large vacuum area. Reason for that is that correction parameter is calculated with self-consistency algorithm which calculates integral over the whole unit cell with vacuum included. This vacuum is responsible for causing incorrect or even diverging results [60].
Next, Rivero et al. [57] investigated surface properties of hydrogenated diamond in the presence of adsorbates by B3LYP and PBE showing that B3LYP provides much larger band gaps than PBE. O’Donnell et al. [64] also deliberated that the DFT computed band gap of metallic surface is not convincing due to its limitations. They predicted, that bulk band gap of O-terminated diamond was 4.1 eV using GGA, while the value reaches 5.5 eV [65,66]. They proposed to apply scissor correction to GGA eigenvalues (rigid shift of unoccupied states) to achieve an agreement with both experiment and higher levels of theory in the case of diamond [67].

4. Conclusions

We show the use of DFT as a powerful tool for ab-initio predictions of electrical transport properties in the design of biofunctionalized boron-doped diamond electrodes. The simulation of BDDE and its surface functionalization via 4-amino benzoic acid, and the ligation of a protein surrogate (anti-M1 for influenza) shows that BDDE-H-ABA-C3M behaves like P-N junction, the availability of sufficient electron transport pathways and local density of states enable electronic recognition of analytes under externally applied electric fields, a band gap drop from 4.4 eV to 4.01 eV for BDDE to BDDE-ABA, and down to 2.6 eV for BDDE-H-ABA-C3M. Furthermore, by monitoring progress in functionalization, either optically or electrically (changing the band gap), we can optimise the conductivity of the sensor and search for optimal treatment. As shown here, DFT can be used to characterise, screen, and optimise the electronic transport properties of electrochemical bio-sensing devices.

Author Contributions

B.D. conducted the DFT calculations and wrote the manuscript; M.S. and R.B. contributed to the discussion and wrote the manuscript. A.J.-B. and W.A.G. contributed to the conception of the study, including DFT calculations concept.

Funding

This work was supported by The National Centre for Research and Development Techmatstrateg 347324/12/NCBR/2017. The DS funds of Faculty of Electronics, Telecommunications and Informatics of the Gdańsk University of Technology are also acknowledged.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Guo, Z.; Chen, L.; Zeng, H.; Gomez, J.A.; Plowden, J.; Fujita, T.; Katz, J.M.; Donis, R.O.; Sambhara, S. NS1 Protein of Influenza A Virus Inhibits the Function of Intracytoplasmic Pathogen Sensor, RIG-I. Am. J. Respir. Cell Mol. Biol. 2007, 36, 263–269. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Hartree, D.R. The Wave Mechanics of an Atom with a Non-Coulomb Central Field. Part II. Some Results and Discussion. Math. Proc. Camb. Philos. Soc. 1928, 24, 111–132. [Google Scholar] [CrossRef]
  3. Liu, X.-Y.; Andersson, D.A.; Uberuaga, B.P. First-principles DFT modeling of nuclear fuel materials. J. Mater. Sci. 2012, 47, 7367–7384. [Google Scholar] [CrossRef]
  4. Tian; Liu DFT Description on Electronic Structure and Optical Absorption Properties of Anionic S-Doped Anatase TiO2. J. Phys. Chem. B 2006, 110, 17866–17871. [CrossRef] [PubMed]
  5. Gurkan, Y.Y.; Kasapbasi, E.; Cinar, Z. Enhanced solar photocatalytic activity of TiO2 by selenium(IV) ion-doping: Characterization and DFT modeling of the surface. Chem. Eng. J. 2013, 214, 34–44. [Google Scholar] [CrossRef]
  6. Maron, L.; Eisenstein, O. Do f Electrons Play a Role in the Lanthanide−Ligand Bonds? A DFT Study of Ln(NR2)3; R = H, SiH3. J. Phys. Chem. A 2000, 104, 7140–7143. [Google Scholar] [CrossRef]
  7. Maron, L.; Eisenstein, O. DFT modeling of ligands in lanthanide chemistry: Is Ln[N(SiH3)2]3 a model for Ln[N(SiMe3)2]3? New J. Chem. 2001, 25, 255–258. [Google Scholar] [CrossRef]
  8. Petit, L.; Daul, C.; Adamo, C.; Maldivi, P. DFT modeling of the relative affinity of nitrogen ligands for trivalent f elements: An energetic point of view. New J. Chem. 2007, 31, 1738–1745. [Google Scholar] [CrossRef]
  9. Bogdanowicz, R. Characterization of Optical and Electrical Properties of Transparent Conductive Boron-Doped Diamond thin Films Grown on Fused Silica. Metrol. Meas. Syst. 2014, 21, 685–698. [Google Scholar] [CrossRef] [Green Version]
  10. Sobaszek, M.; Skowroński, Ł.; Bogdanowicz, R.; Siuzdak, K.; Cirocka, A.; Zięba, P.; Gnyba, M.; Naparty, M.; Gołuński, Ł.; Płotka, P. Optical and electrical properties of ultrathin transparent nanocrystalline boron-doped diamond electrodes. Opt. Mater. 2015, 42, 24–34. [Google Scholar] [CrossRef]
  11. Ficek, M.; Sobaszek, M.; Gnyba, M.; Ryl, J.; Gołuński, Ł.; Smietana, M.; Jasiński, J.; Caban, P.; Bogdanowicz, R. Optical and electrical properties of boron doped diamond thin conductive films deposited on fused silica glass substrates. Appl. Surf. Sci. 2016, 387, 846–856. [Google Scholar] [CrossRef]
  12. Tang, L.; Tsai, C.; Gerberich, W.W.; Kruckeberg, L.; Kania, D.R. Biocompatibility of chemical-vapour-deposited diamond. Biomaterials 1995, 16, 483–488. [Google Scholar] [CrossRef]
  13. Falina, S.; Kawai, S.; Oi, N.; Yamano, H.; Kageura, T.; Suaebah, E.; Inaba, M.; Shintani, Y.; Syamsul, M.; Kawarada, H. Role of Carboxyl and Amine Termination on a Boron-Doped Diamond Solution Gate Field Effect Transistor (SGFET) for pH Sensing. Sensors 2018, 18, 2178. [Google Scholar] [CrossRef] [PubMed]
  14. Karaya, R.; Baba, I.; Mori, Y.; Matsumoto, T.; Nakajima, T.; Tokuda, N.; Kawae, T. B-doped diamond field-effect transistor with ferroelectric vinylidene fluoride–trifluoroethylene gate insulator. Jpn. J. Appl. Phys. 2017, 56, 10PF06. [Google Scholar] [CrossRef]
  15. Zkria, A.; Gima, H.; Yoshitake, T. Application of nitrogen-doped ultrananocrystalline diamond/hydrogenated amorphous carbon composite films for ultraviolet detection. Appl. Phys. A 2017, 123, 167. [Google Scholar] [CrossRef]
  16. Abubakr, E.; Zkria, A.; Katamune, Y.; Ohmagari, S.; Imokawa, K.; Ikenoue, H.; Yoshitake, T. Formation of low resistivity layers on singlecrystalline diamond by excimer laser irradiation. Diam. Relat. Mater. 2019, 95, 166–173. [Google Scholar] [CrossRef]
  17. Fujimori, N.; Nakahata, H.; Imai, T. Properties of Boron-Doped Epitaxial Diamond Films. Jpn. J. Appl. Phys. 1990, 29, 824. [Google Scholar] [CrossRef]
  18. Carlisle, J.A. Precious biosensors. Nat. Mater. 2004, 3, 668. [Google Scholar] [CrossRef] [PubMed]
  19. Härtl, A.; Schmich, E.; Garrido, J.A.; Hernando, J.; Catharino, S.C.R.; Walter, S.; Feulner, P.; Kromka, A.; Steinmüller, D.; Stutzmann, M. Protein-modified nanocrystalline diamond thin films for biosensor applications. Nat. Mater. 2004, 3, 736. [Google Scholar] [CrossRef] [PubMed]
  20. Boukhvalov, D.W. DFT modeling of the covalent functionalization of graphene: From ideal to realistic models. RSC Adv. 2013, 3, 7150–7159. [Google Scholar] [CrossRef]
  21. Jaramillo-Botero, A.; Marmolejo-Tejada, J.M. All-Armchair Graphene Nanoribbon Field-Effect Uridine Diphosphate Glucose Sensor: First-Principles In-Silico Design and Characterization. IEEE Sens. J. 2019, 19, 3975–3983. [Google Scholar] [CrossRef]
  22. Tian, Y.; Larsson, K. Effect by Diamond Surface Modification on Biomolecular Adhesion. Mater. Basel Switz. 2019, 12. [Google Scholar]
  23. Dec, B.; Bogdanowicz, R. DFT studies of refractive index of boron-doped diamond. Photonics Lett. Pol. 2018, 10, 39–41. [Google Scholar] [CrossRef]
  24. Sobaszek, M.; Siuzdak, K.; Skowroński, L.; Bogdanowicz, R.; Pluciński, J. Optically transparent boron-doped nanocrystalline diamond films for spectroelectrochemical measurements on different substrates. IOP Conf. Ser. Mater. Sci. Eng. 2016, 104, 012024. [Google Scholar] [CrossRef] [Green Version]
  25. Dec, B.; Bogdanowicz, R.; Pyrchla, K. Ab-initio study of electrical and optical properties of allylamine. Photonics Lett. Pol. 2018, 10, 94–96. [Google Scholar] [CrossRef]
  26. Kohn, W.; Becke, A.D.; Parr, R.G. Density Functional Theory of Electronic Structure. J. Phys. Chem. 1996, 100, 12974–12980. [Google Scholar] [CrossRef] [Green Version]
  27. Smidstrup, S.; Stradi, D.; Wellendorff, J.; Khomyakov, P.A.; Vej-Hansen, U.G.; Lee, M.E.; Ghosh, T.; Jónsson, E.; Jónsson, H.; Stokbro, K. First-principles Green's-function method for surface calculations: A pseudopotential localized basis set approach. Phys. Rev. B 2017, 96, 195309. [Google Scholar] [CrossRef] [Green Version]
  28. Griebel, M.; Hamaekers, J. Molecular dynamics simulations of the elastic moduli of polymer–carbon nanotube composites. Comput. Methods Appl. Mech. Eng. 2004, 193, 1773–1788. [Google Scholar] [CrossRef]
  29. Griebel, M.; Knapek, S.; Zumbusch, G. Numerical Simulation in Molecular Dynamics: Numerics, Algorithms, Parallelization, Applications; Springer Science & Business Media: Heidelberg, Germany, 2007; ISBN 978-3-540-68095-6. [Google Scholar]
  30. Ekimov, E.A.; Sidorov, V.A.; Bauer, E.D.; Mel’nik, N.N.; Curro, N.J.; Thompson, J.D.; Stishov, S.M. Superconductivity in diamond. Nature 2004, 428, 542–545. [Google Scholar] [CrossRef] [Green Version]
  31. Zhao, S.; Larsson, K. Theoretical Study of the Energetic Stability and Geometry of Terminated and B-Doped Diamond (111) Surfaces. J. Phys. Chem. C 2014, 118, 1944–1957. [Google Scholar] [CrossRef]
  32. Schlipf, M.; Gygi, F. Optimization algorithm for the generation of ONCV pseudopotentials. Comput. Phys. Commun. 2015, 196, 36–44. [Google Scholar] [CrossRef] [Green Version]
  33. Chelikowsky, J.R.; Louie, S.G. First-principles linear combination of atomic orbitals method for the cohesive and structural properties of solids: Application to diamond. Phys. Rev. B 1984, 29, 3470–3481. [Google Scholar] [CrossRef]
  34. Perdew, J.P.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 1996, 77, 3865–3868. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Monkhorst, H.J.; Pack, J.D. Special points for Brillouin-zone integrations. Phys. Rev. B 1976, 13, 5188–5192. [Google Scholar] [CrossRef]
  36. Fuchs, M.; Scheffler, M. Ab initio pseudopotentials for electronic structure calculations of poly-atomic systems using density-functional theory. Comput. Phys. Commun. 1999, 119, 67–98. [Google Scholar] [CrossRef] [Green Version]
  37. Liu, D.C.; Nocedal, J. On the limited memory BFGS method for large scale optimization. Math. Program. 1989, 45, 503–528. [Google Scholar] [CrossRef] [Green Version]
  38. Turner, S.; Lu, Y.-G.; Janssens, S.D.; Da Pieve, F.; Lamoen, D.; Verbeeck, J.; Haenen, K.; Wagner, P.; Van Tendeloo, G. Local boron environment in B-doped nanocrystalline diamond films. Nanoscale 2012, 4, 5960. [Google Scholar] [CrossRef] [PubMed]
  39. Brunetta, C.D.; Karuppannan, B.; Rosmus, K.A.; Aitken, J.A. The crystal and electronic band structure of the diamond-like semiconductor Ag2ZnSiS4. J. Alloys Compd. 2012, 516, 65–72. [Google Scholar] [CrossRef]
  40. Lejaeghere, K.; Bihlmayer, G.; Björkman, T.; Blaha, P.; Blügel, S.; Blum, V.; Caliste, D.; Castelli, I.E.; Clark, S.J.; Corso, A.D.; et al. Reproducibility in density functional theory calculations of solids. Science 2016, 351, aad3000. [Google Scholar] [CrossRef]
  41. Comparing Solid State DFT Codes, Basis Sets and Potentials. Available online: https://molmod.ugent.be/deltacodesdft (accessed on 25 July 2019).
  42. Lejaeghere, K.; Speybroeck, V.; Oost, G.; Cottenier, S. Error Estimates for Solid-State Density-Functional Theory Predictions: An Overview by Means of the Ground-State Elemental Crystals. Crit. Rev. Solid State Mater. Sci. 2013, 39, 1–24. [Google Scholar] [Green Version]
  43. Bogdanowicz, R.; Fabiańska, A.; Golunski, L.; Sobaszek, M.; Gnyba, M.; Ryl, J.; Darowicki, K.; Ossowski, T.; Janssens, S.D.; Haenen, K.; et al. Influence of the boron doping level on the electrochemical oxidation of the azo dyes at Si/BDD thin film electrodes. Diam. Relat. Mater. 2013, 39, 82–88. [Google Scholar] [CrossRef]
  44. Cui, J.B.; Ristein, J.; Ley, L. Electron Affinity of the Bare and Hydrogen Covered Single Crystal Diamond (111) Surface. Phys. Rev. Lett. 1998, 81, 429–432. [Google Scholar] [CrossRef]
  45. Dai, Y.; Huang, B.; Dai, D. The role of dangling-bond, hydrogen and adsorbate in diamond surface conduction. Diam. Relat. Mater. 2003, 12, 15–19. [Google Scholar] [CrossRef]
  46. Chua, L.-L.; Zaumseil, J.; Chang, J.-F.; Ou, E.C.-W.; Ho, P.K.-H.; Sirringhaus, H.; Friend, R.H. General observation of n-type field-effect behaviour in organic semiconductors. Nature 2005, 434, 194–199. [Google Scholar] [CrossRef] [PubMed]
  47. Cornil, J.; Beljonne, D.; Calbert, J.-P.; Brédas, J.-L. Interchain Interactions in Organic π-Conjugated Materials: Impact on Electronic Structure, Optical Response, and Charge Transport. Adv. Mater. 2001, 13, 1053–1067. [Google Scholar] [CrossRef]
  48. Coropceanu, V.; Cornil, J.; da Silva Filho, D.A.; Olivier, Y.; Silbey, R.; Brédas, J.-L. Charge Transport in Organic Semiconductors. Chem. Rev. 2007, 107, 926–952. [Google Scholar] [CrossRef] [PubMed]
  49. Zhao, S.; Larsson, K. First Principle Study of the Attachment of Graphene onto Different Terminated Diamond (111) Surfaces. Available online: https://www.hindawi.com/journals/acmp/2019/9098256/ (accessed on 28 August 2019).
  50. Ashcheulov, P.; Šebera, J.; Kovalenko, A.; Petrák, V.; Fendrych, F.; Nesládek, M.; Taylor, A.; Vlčková Živcová, Z.; Frank, O.; Kavan, L.; et al. Conductivity of boron-doped polycrystalline diamond films: Influence of specific boron defects. Eur. Phys. J. B 2013, 86, 443. [Google Scholar] [CrossRef]
  51. Strange, M.; Thygesen, K.S. Towards quantitative accuracy in first-principles transport calculations: The GW method applied to alkane/gold junctions. Beilstein J. Nanotechnol. 2011, 2, 746–754. [Google Scholar] [CrossRef] [Green Version]
  52. Wang, L.; Liu, J.; Tang, T. The Influence of B, N and Si Doping on the CH3 Adsorption on the Diamond Surface Based on DFT Calculations. Crystals 2019, 9, 427. [Google Scholar] [CrossRef]
  53. Watanabe, T.; Teraji, T.; Ito, T.; Kamakura, Y.; Taniguchi, K. Monte Carlo simulations of electron transport properties of diamond in high electric fields using full band structure. J. Appl. Phys. 2004, 95, 4866–4874. [Google Scholar] [CrossRef]
  54. Kamakura, Y.; Fujita, R.; Konaga, K.; Ueoka, Y.; Mori, N.; Kotani, T. Full band Monte Carlo simulation of impact ionization in wide bandgap semiconductors based on ab initio calculation. In Proceedings of the 2016 International Conference on Simulation of Semiconductor Processes and Devices (SISPAD), Nuremberg, Germany, 6–8 September 2016; pp. 47–51. [Google Scholar]
  55. Barnard, A.S.; Russo, S.P.; Snook, I.K. Ab initio modelling of band states in doped diamond. Philos. Mag. 2003, 83, 1163–1174. [Google Scholar] [CrossRef]
  56. Collins, A.T.; Lightowlers, E.C. Photothermal Ionization and Photon-Induced Tunneling in the Acceptor Photoconductivity Spectrum of Semiconducting Diamond. Phys. Rev. 1968, 171, 843–855. [Google Scholar] [CrossRef]
  57. Rivero, P.; Shelton, W.; Meunier, V. Surface properties of hydrogenated diamond in the presence of adsorbates: A hybrid functional DFT study. Carbon 2016, 110, 469–479. [Google Scholar] [CrossRef] [Green Version]
  58. Liberman, D.A. Slater transition-state band-structure calculations. Phys. Rev. B 2000, 62, 6851–6853. [Google Scholar] [CrossRef]
  59. Ullah, M.; Ahmed, E.; Hussain, F.; Rana, A.M.; Raza, R.; Ullah, H. Electronic structure calculations of oxygen-doped diamond using DFT technique. Microelectron. Eng. 2015, 146, 26–31. [Google Scholar] [CrossRef]
  60. Tran, F.; Blaha, P. Accurate Band Gaps of Semiconductors and Insulators with a Semilocal Exchange-Correlation Potential. Phys. Rev. Lett. 2009, 102, 226401. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  61. Heyd, J.; Scuseria, G.E.; Ernzerhof, M. Hybrid functionals based on a screened Coulomb potential. J. Chem. Phys. 2003, 118, 8207–8215. [Google Scholar] [CrossRef] [Green Version]
  62. Xiao, H.; Tahir-Kheli, J.; Goddard, W.A. Accurate Band Gaps for Semiconductors from Density Functional Theory. J. Phys. Chem. Lett. 2011, 2, 212–217. [Google Scholar] [CrossRef] [Green Version]
  63. Heyd, J.; Scuseria, G.E. Efficient hybrid density functional calculations in solids: Assessment of the Heyd–Scuseria–Ernzerhof screened Coulomb hybrid functional. J. Chem. Phys. 2004, 121, 1187–1192. [Google Scholar] [CrossRef]
  64. O’Donnell, K.M.; Martin, T.L.; Allan, N.L. Light Metals on Oxygen-Terminated Diamond (100): Structure and Electronic Properties. Chem. Mater. 2015, 27, 1306–1315. [Google Scholar] [CrossRef] [Green Version]
  65. Dean, P.J.; Lightowlers, E.C.; Wight, D.R. Intrinsic and Extrinsic Recombination Radiation from Natural and Synthetic Aluminum-Doped Diamond. Phys. Rev. 1965, 140, A352–A368. [Google Scholar] [CrossRef]
  66. Ristein, J. Electronic properties of diamond surfaces — blessing or curse for devices? Diam. Relat. Mater. 2000, 9, 1129–1137. [Google Scholar] [CrossRef]
  67. Gairola, S.C. Electron and Phonon in Neutron Scattering of Semiconductor Crystals. Int. J. Theor. Phys. 2013, 12, 4256–4264. [Google Scholar] [CrossRef]
Figure 1. Optimized (4 × 4) model structures: (A)—(BDDE) relaxed boron-doped diamond (111); (B)—(BDDE-H-ABA)-H-terminated reconstructed boron-doped diamond (111) with adlayer of 4-amino-benzoic acid; (C) (BDDE-H-ABA-C3M)-H-terminated reconstructed boron-doped diamond (111) with 4-amino benzoic acid and ligand ending.
Figure 1. Optimized (4 × 4) model structures: (A)—(BDDE) relaxed boron-doped diamond (111); (B)—(BDDE-H-ABA)-H-terminated reconstructed boron-doped diamond (111) with adlayer of 4-amino-benzoic acid; (C) (BDDE-H-ABA-C3M)-H-terminated reconstructed boron-doped diamond (111) with 4-amino benzoic acid and ligand ending.
Materials 12 02910 g001
Figure 2. Electron density iso-surfaces: (A)—(BDDE) relaxed boron-doped diamond (111) slab; (B)—(BDDE-H-ABA)-H-terminated reconstructed boron-doped diamond (111) with adlayer of 4-aminobenzoic linker, (C) (BDDE-H-ABA-C3H)—BDDE modified with 4-aminobenzoic linker bonded to C3H ligand of the anti-M1 antibody protein. Azure colour defines electron density close to zero, and the purple describes values up to 1.
Figure 2. Electron density iso-surfaces: (A)—(BDDE) relaxed boron-doped diamond (111) slab; (B)—(BDDE-H-ABA)-H-terminated reconstructed boron-doped diamond (111) with adlayer of 4-aminobenzoic linker, (C) (BDDE-H-ABA-C3H)—BDDE modified with 4-aminobenzoic linker bonded to C3H ligand of the anti-M1 antibody protein. Azure colour defines electron density close to zero, and the purple describes values up to 1.
Materials 12 02910 g002
Figure 3. Electron pathways in bulk diamond: (A) structure without any external potential across the structure. (B) Structure with applied 1V potential across the structure, where on the left side is −0.5V (−), right side 0.5V (+). The 2-terminal device structure was used (attn. electrodes are not shown here for figure readability).
Figure 3. Electron pathways in bulk diamond: (A) structure without any external potential across the structure. (B) Structure with applied 1V potential across the structure, where on the left side is −0.5V (−), right side 0.5V (+). The 2-terminal device structure was used (attn. electrodes are not shown here for figure readability).
Materials 12 02910 g003
Figure 4. Electron pathways at the diamond surface modified by the attached linker (4-aminobenzoic acid): (A) structure without applied potential along C direction (from left to right), (B) structure with applied potential −1V, (−0.5V on the left (−), 0.5V on the right side (+)). The 2-terminal device structure was used (attn. electrodes are not shown here for figure readability).
Figure 4. Electron pathways at the diamond surface modified by the attached linker (4-aminobenzoic acid): (A) structure without applied potential along C direction (from left to right), (B) structure with applied potential −1V, (−0.5V on the left (−), 0.5V on the right side (+)). The 2-terminal device structure was used (attn. electrodes are not shown here for figure readability).
Materials 12 02910 g004
Figure 5. Electron pathways at the diamond surface with attached ligand through linker: (A) structure not polarised with potential, (B) structure polarised with −1V potential (−0.5V on the left (−) and 0.5V on the right side (+) of the slab). The 2-terminal device structure was used (attn. electrodes are not shown here for figure readability).
Figure 5. Electron pathways at the diamond surface with attached ligand through linker: (A) structure not polarised with potential, (B) structure polarised with −1V potential (−0.5V on the left (−) and 0.5V on the right side (+) of the slab). The 2-terminal device structure was used (attn. electrodes are not shown here for figure readability).
Materials 12 02910 g005
Figure 6. The density of states spectra for BDDE-H-ABA and BDD-H-ABA-C3M with additionally added BDD band gap for change comparison.
Figure 6. The density of states spectra for BDDE-H-ABA and BDD-H-ABA-C3M with additionally added BDD band gap for change comparison.
Materials 12 02910 g006
Table 1. Comparison of delta-test results for pseudopotentials and basis set with chosen atoms that were used in the investigated system.
Table 1. Comparison of delta-test results for pseudopotentials and basis set with chosen atoms that were used in the investigated system.
Atom TypeSG15 (Medium) [meV]SG15 (High) [meV]FHI (DPZ) [meV]
Boron1.651.883.67
Carbon0.662.805.89
Hydrogen0.960.190.54
Nitrogen6.571.042.45
Oxygen9.371.6824.28

Share and Cite

MDPI and ACS Style

Dec, B.; Sobaszek, M.; Jaramillo-Botero, A.; Goddard, W.A., III; Bogdanowicz, R. Ligand-Modified Boron-Doped Diamond Surface: DFT Insights into the Electronic Properties of Biofunctionalization. Materials 2019, 12, 2910. https://doi.org/10.3390/ma12182910

AMA Style

Dec B, Sobaszek M, Jaramillo-Botero A, Goddard WA III, Bogdanowicz R. Ligand-Modified Boron-Doped Diamond Surface: DFT Insights into the Electronic Properties of Biofunctionalization. Materials. 2019; 12(18):2910. https://doi.org/10.3390/ma12182910

Chicago/Turabian Style

Dec, Bartłomiej, Michał Sobaszek, Andrés Jaramillo-Botero, William Andrew Goddard, III, and Robert Bogdanowicz. 2019. "Ligand-Modified Boron-Doped Diamond Surface: DFT Insights into the Electronic Properties of Biofunctionalization" Materials 12, no. 18: 2910. https://doi.org/10.3390/ma12182910

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop