CN114381742B - Nonmetallic monoatomic doping BSe monolayer hydrogen evolution catalyst, design method and application - Google Patents
Nonmetallic monoatomic doping BSe monolayer hydrogen evolution catalyst, design method and application Download PDFInfo
- Publication number
- CN114381742B CN114381742B CN202210040241.3A CN202210040241A CN114381742B CN 114381742 B CN114381742 B CN 114381742B CN 202210040241 A CN202210040241 A CN 202210040241A CN 114381742 B CN114381742 B CN 114381742B
- Authority
- CN
- China
- Prior art keywords
- bse
- monolayer
- nonmetallic
- monoatomic
- configuration
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 239000002356 single layer Substances 0.000 title claims abstract description 122
- 229910052739 hydrogen Inorganic materials 0.000 title claims abstract description 77
- 239000001257 hydrogen Substances 0.000 title claims abstract description 77
- UFHFLCQGNIYNRP-UHFFFAOYSA-N Hydrogen Chemical compound [H][H] UFHFLCQGNIYNRP-UHFFFAOYSA-N 0.000 title claims abstract description 76
- 239000003054 catalyst Substances 0.000 title claims abstract description 63
- 238000000034 method Methods 0.000 title claims abstract description 31
- 238000013461 design Methods 0.000 title claims abstract description 12
- 238000006243 chemical reaction Methods 0.000 claims abstract description 30
- 230000003197 catalytic effect Effects 0.000 claims abstract description 17
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 10
- 238000004519 manufacturing process Methods 0.000 claims abstract description 5
- 229910052711 selenium Inorganic materials 0.000 claims description 28
- BUGBHKTXTAQXES-UHFFFAOYSA-N Selenium Chemical group [Se] BUGBHKTXTAQXES-UHFFFAOYSA-N 0.000 claims description 23
- 239000011669 selenium Substances 0.000 claims description 22
- ZOXJGFHDIHLPTG-UHFFFAOYSA-N Boron Chemical group [B] ZOXJGFHDIHLPTG-UHFFFAOYSA-N 0.000 claims description 17
- 229910052799 carbon Inorganic materials 0.000 claims description 16
- 125000003748 selenium group Chemical group *[Se]* 0.000 claims description 14
- 229910052757 nitrogen Inorganic materials 0.000 claims description 12
- 229910052796 boron Inorganic materials 0.000 claims description 11
- 229910052698 phosphorus Inorganic materials 0.000 claims description 11
- 239000000126 substance Substances 0.000 claims description 10
- 230000015572 biosynthetic process Effects 0.000 claims description 8
- 238000006467 substitution reaction Methods 0.000 claims description 5
- 229910052755 nonmetal Inorganic materials 0.000 claims 1
- 230000001737 promoting effect Effects 0.000 abstract 1
- 125000004429 atom Chemical group 0.000 description 20
- 238000010586 diagram Methods 0.000 description 9
- 125000004435 hydrogen atom Chemical group [H]* 0.000 description 7
- 238000011161 development Methods 0.000 description 5
- 239000002019 doping agent Substances 0.000 description 5
- 239000000463 material Substances 0.000 description 5
- 229910052710 silicon Inorganic materials 0.000 description 5
- 239000010411 electrocatalyst Substances 0.000 description 4
- 239000012535 impurity Substances 0.000 description 4
- BASFCYQUMIYNBI-UHFFFAOYSA-N platinum Chemical compound [Pt] BASFCYQUMIYNBI-UHFFFAOYSA-N 0.000 description 4
- 238000001179 sorption measurement Methods 0.000 description 4
- 229910052785 arsenic Inorganic materials 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000005868 electrolysis reaction Methods 0.000 description 3
- 229910000510 noble metal Inorganic materials 0.000 description 3
- 229910052760 oxygen Inorganic materials 0.000 description 3
- 239000004065 semiconductor Substances 0.000 description 3
- 229910052717 sulfur Inorganic materials 0.000 description 3
- CURLTUGMZLYLDI-UHFFFAOYSA-N Carbon dioxide Chemical compound O=C=O CURLTUGMZLYLDI-UHFFFAOYSA-N 0.000 description 2
- PXHVJJICTQNCMI-UHFFFAOYSA-N Nickel Chemical compound [Ni] PXHVJJICTQNCMI-UHFFFAOYSA-N 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 2
- 239000007809 chemical reaction catalyst Substances 0.000 description 2
- 230000007613 environmental effect Effects 0.000 description 2
- 239000002803 fossil fuel Substances 0.000 description 2
- 230000003993 interaction Effects 0.000 description 2
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 2
- 229910052697 platinum Inorganic materials 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000010998 test method Methods 0.000 description 2
- 229910052723 transition metal Inorganic materials 0.000 description 2
- 229910052582 BN Inorganic materials 0.000 description 1
- PZNSFCLAULLKQX-UHFFFAOYSA-N Boron nitride Chemical compound N#B PZNSFCLAULLKQX-UHFFFAOYSA-N 0.000 description 1
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 1
- UGFAIRIUMAVXCW-UHFFFAOYSA-N Carbon monoxide Chemical compound [O+]#[C-] UGFAIRIUMAVXCW-UHFFFAOYSA-N 0.000 description 1
- VGGSQFUCUMXWEO-UHFFFAOYSA-N Ethene Chemical compound C=C VGGSQFUCUMXWEO-UHFFFAOYSA-N 0.000 description 1
- 239000005977 Ethylene Substances 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 1
- 230000004888 barrier function Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 229910002092 carbon dioxide Inorganic materials 0.000 description 1
- 239000001569 carbon dioxide Substances 0.000 description 1
- 229910002091 carbon monoxide Inorganic materials 0.000 description 1
- 238000006555 catalytic reaction Methods 0.000 description 1
- 239000007806 chemical reaction intermediate Substances 0.000 description 1
- 239000003034 coal gas Substances 0.000 description 1
- 239000010941 cobalt Substances 0.000 description 1
- 229910017052 cobalt Inorganic materials 0.000 description 1
- GUTLYIVDDKVIGB-UHFFFAOYSA-N cobalt atom Chemical compound [Co] GUTLYIVDDKVIGB-UHFFFAOYSA-N 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 230000006866 deterioration Effects 0.000 description 1
- 230000009977 dual effect Effects 0.000 description 1
- 238000004146 energy storage Methods 0.000 description 1
- 238000006735 epoxidation reaction Methods 0.000 description 1
- 239000007789 gas Substances 0.000 description 1
- 229910021389 graphene Inorganic materials 0.000 description 1
- 150000002431 hydrogen Chemical class 0.000 description 1
- 238000005984 hydrogenation reaction Methods 0.000 description 1
- 230000004807 localization Effects 0.000 description 1
- 229910052751 metal Inorganic materials 0.000 description 1
- 239000002184 metal Substances 0.000 description 1
- CWQXQMHSOZUFJS-UHFFFAOYSA-N molybdenum disulfide Chemical compound S=[Mo]=S CWQXQMHSOZUFJS-UHFFFAOYSA-N 0.000 description 1
- 229910052982 molybdenum disulfide Inorganic materials 0.000 description 1
- 239000003345 natural gas Substances 0.000 description 1
- 230000007935 neutral effect Effects 0.000 description 1
- 229910052759 nickel Inorganic materials 0.000 description 1
- 239000003921 oil Substances 0.000 description 1
- 230000003647 oxidation Effects 0.000 description 1
- 238000007254 oxidation reaction Methods 0.000 description 1
- 239000001301 oxygen Substances 0.000 description 1
- 230000037361 pathway Effects 0.000 description 1
- 230000001699 photocatalysis Effects 0.000 description 1
- 231100000572 poisoning Toxicity 0.000 description 1
- 230000000607 poisoning effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000027756 respiratory electron transport chain Effects 0.000 description 1
- 239000000758 substrate Substances 0.000 description 1
- 229910052714 tellurium Inorganic materials 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- -1 transition metal carbides Chemical class 0.000 description 1
- 150000003624 transition metals Chemical class 0.000 description 1
- 238000001845 vibrational spectrum Methods 0.000 description 1
Classifications
-
- C—CHEMISTRY; METALLURGY
- C25—ELECTROLYTIC OR ELECTROPHORETIC PROCESSES; APPARATUS THEREFOR
- C25B—ELECTROLYTIC OR ELECTROPHORETIC PROCESSES FOR THE PRODUCTION OF COMPOUNDS OR NON-METALS; APPARATUS THEREFOR
- C25B1/00—Electrolytic production of inorganic compounds or non-metals
- C25B1/01—Products
- C25B1/02—Hydrogen or oxygen
- C25B1/04—Hydrogen or oxygen by electrolysis of water
-
- C—CHEMISTRY; METALLURGY
- C25—ELECTROLYTIC OR ELECTROPHORETIC PROCESSES; APPARATUS THEREFOR
- C25B—ELECTROLYTIC OR ELECTROPHORETIC PROCESSES FOR THE PRODUCTION OF COMPOUNDS OR NON-METALS; APPARATUS THEREFOR
- C25B11/00—Electrodes; Manufacture thereof not otherwise provided for
- C25B11/04—Electrodes; Manufacture thereof not otherwise provided for characterised by the material
- C25B11/051—Electrodes formed of electrocatalysts on a substrate or carrier
- C25B11/073—Electrodes formed of electrocatalysts on a substrate or carrier characterised by the electrocatalyst material
- C25B11/091—Electrodes formed of electrocatalysts on a substrate or carrier characterised by the electrocatalyst material consisting of at least one catalytic element and at least one catalytic compound; consisting of two or more catalytic elements or catalytic compounds
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
- Y02E60/30—Hydrogen technology
- Y02E60/36—Hydrogen production from non-carbon containing sources, e.g. by water electrolysis
Landscapes
- Chemical & Material Sciences (AREA)
- Engineering & Computer Science (AREA)
- Chemical Kinetics & Catalysis (AREA)
- Electrochemistry (AREA)
- Materials Engineering (AREA)
- Metallurgy (AREA)
- Organic Chemistry (AREA)
- Inorganic Chemistry (AREA)
- Catalysts (AREA)
Abstract
The design method of the nonmetallic monoatomic doping BSe monolayer hydrogen evolution catalyst comprises the following steps: (1) Selecting a monolayer BSe supercell with a set specification, and calculating the electronic configuration and the geometric structure of the supercell; (2) Constructing a nonmetallic monoatomic doping BSe monolayer configuration on the monolayer BSe supercell; (3) Determining a configuration of stable structure in the nonmetallic monoatomic doping BSe monolayer configuration; (4) The configuration for promoting hydrogen evolution reaction in the configuration with stable structure is determined to be a nonmetallic monoatomic doping BSe monolayer hydrogen evolution catalyst, and the obtained nonmetallic monoatomic doping BSe monolayer with catalytic activity is used as the hydrogen evolution catalyst, can be applied to a process for preparing hydrogen by electrolyzing water, and has important guiding significance and practical significance.
Description
Technical Field
The application belongs to the technical field of catalysts, and particularly relates to a nonmetallic monoatomic doping BSe monolayer hydrogen evolution catalyst, a design method and application.
Background
With the rapid development of the global industry, the problems of energy shortage and environmental deterioration come up from the excessive exploitation and utilization of fossil fuel oil, coal and natural gas, and the search for clean energy for sustainable development becomes a great hotspot.
Hydrogen is a high density energy carrier for zero carbon dioxide emissions, and is increasingly gaining attention because of its contribution to replacing fossil fuels and solving environmental problems. In the hydrogen production process, the electrolyzed water is used for preparing hydrogen, has the advantages of simple equipment, good product sanitary condition, high reproducibility and the like, and is a method with attractive force. Taking water decomposition reaction as an example, it includes hydrogen evolution half reaction and oxygen evolution half reaction, wherein the hydrogen evolution half reaction HER (2H++ 2e- →H) 2 ) Is a dual electron transfer pathway to a reaction intermediate comprising three possible reaction steps:
(i) Volmer step: h++ e- + →h;
(ii) Heyrovsky step: H+H++ e- →H 2 +*;
(iii)Tafel step: 2H → H 2 +*;
Wherein, represents a catalytically active site.
A great deal of theoretical and experimental researches prove that the ideal catalyst can improve the energy conversion efficiency, improve the speed of hydrogen evolution reaction, reduce the energy barrier or overpotential and provide definite active sites, so that the limiting factor of the hydrogen production by water electrolysis is the overpotential, and the overpotential is determined by different catalyst properties.
A new two-dimensional layered material monolayer BSe was successfully designed by Demirci et al in 2017. To further confirm its thermal stability, they also calculated the vibration spectrum of monolayer BSe. BSe two-dimensional material has been theoretically projected, and its unit cell structure is part of the space group D3h, similar to a monolayer InSe. Studies have shown that the monolayer BSe forms heterostructures with higher photocatalytic properties with other two-dimensional materials. At present, two-dimensional layered materials are attracting attention because of their controllable properties, low cost, large surface area, good stability, and excellent properties in energy storage, nanoelectronics, catalysis, and the like. For example, graphene, hexagonal boron nitride, molybdenum disulfide, two-dimensional transition metal carbides, and the like have been demonstrated to intercalate metal atoms to provide active sites.
The single-atom catalyst can effectively excite the activity of the two-dimensional material, is dispersed on a substrate with larger surface area, and has good catalytic performance on a series of chemical reactions such as hydrogen evolution reaction, carbon monoxide oxidation, ethylene epoxidation, water gas shift reaction, organic molecule hydrogenation and the like. The monoatomic catalysts currently used in large numbers are transition metal-based catalysts, such as platinum catalysts, cobalt catalysts and nickel catalysts. Noble metal platinum has higher catalytic activity, catalytic efficiency and catalytic index, and is an ideal hydrogen evolution reaction electrocatalyst at present. Although noble metal monoatomic catalysts can be used as water electrolysis catalysts, the noble metal monoatomic catalysts have the problems of high price, low stability, easy poisoning and the like.
The advantages of the nonmetallic catalyst such as easy acquisition of the dopant, low cost and high catalytic performance become a new development direction in the catalytic field. In recent years, various nonmetallic (C, si, N, and O) doped catalysts have been designed and demonstrated for high catalytic performance for hydrogen evolution reactions. Therefore, the development of nonmetallic doped catalysts has become a necessary trend for future development.
Disclosure of Invention
In view of this, in one aspect, some embodiments disclose a method of designing a nonmetallic monoatomic doping BSe monolayer hydrogen evolution catalyst, the method comprising the steps of:
(1) Selecting a monolayer BSe supercell with a set specification, and calculating the electronic configuration and the geometric structure of the supercell;
(2) Constructing a nonmetallic monoatomic doping BSe monolayer configuration on the monolayer BSe supercell;
(3) Determining a configuration of stable structure in the nonmetallic monoatomic doping BSe monolayer configuration;
(4) The configuration that promotes the hydrogen evolution reaction in the configuration that determines the structural stability is a nonmetallic monoatomic doping BSe monolayer hydrogen evolution catalyst.
Some embodiments disclose a method for designing a nonmetallic monoatomic doping BSe monolayer hydrogen evolution catalyst, wherein the monolayer BSe supercell selected in step (1) is a 3×3×1 monolayer BSe supercell comprising 36 atoms, and the lattice constant isThe thickness of the single layer is->B-Se bond length is->B-B bond length is->
Some embodiments disclose methods of designing nonmetallic monoatomic doping BSe monolayer hydrogen evolution catalysts, in step (2), the nonmetallic monoatomic doping configuration comprises a configuration in which nonmetallic monoatoms adsorb to the surface sites of the BSe monolayer and a configuration in which nonmetallic monoatoms replace the B or Se atoms in the BSe monolayer.
The design method of the nonmetallic monoatomic doping BSe monolayer hydrogen evolution catalyst disclosed in some embodiments is used to form a catalyst that can evaluate the structural stability of nonmetallic doping configurations, which can be calculated from the following formula:
E form =E tot [BSe * +X]-E tot [BSe]+xμ B +yμ Se -μ X
wherein E is tot [BSe*+X]For the total energy after doping with nonmetallic monoatoms, E tot [BSe]Mu, the total energy of undoped nonmetallic monoatoms X For each nonmetallic atom chemical potential, the coefficients x and y are the number of boron atoms and selenium atoms, respectively, μ B 、μ Se The chemical potentials of boron atoms and selenium atoms respectively, and X represents a substituted nonmetallic element.
Some embodiments disclose a method of designing a nonmetallic monoatomic doping BSe monolayer hydrogen evolution catalyst, the nonmetallic monoatomic doping configuration comprising: the configuration of selenium atoms in the substituted BSe monolayer, the configuration of boron atoms in the substituted BSe monolayer, the configuration of sites above boron atoms adsorbed on the BSe monolayer, the configuration of sites above selenium atoms adsorbed on the BSe monolayer and the configuration of sites above the center of the hexagonal ring of the BSe monolayer.
Some embodiments disclose methods of designing nonmetallic monoatomic doped BSe monolayer hydrogen evolution catalysts in which the selenium atom site configuration in the substituted BSe monolayer is the most stable configuration, wherein the nonmetallic monoatomic and its mass concentration is 0.258wt% c, 0.300wt% n, 0.343wt% o, 0.407wt% f, 0.600wt% si, 0.6612 wt% p, 0.685wt% s, 0.757wt% cl, 1.586wt% as, 1.689wt% br, 2.671wt% te, or 2.656wt% i. The molar ratio of B and Se atoms to nonmetallic atoms in the selenium atom site configuration in the substituted BSe monolayer is 35:1.
Some embodiments disclose a design method of a nonmetallic monoatomic doping BSe monolayer hydrogen evolution catalyst, wherein the configuration of selenium atoms in the nonmetallic monoatomic C, N, si, P, as substituted BSe monolayer has hydrogen evolution catalytic activity.
Some embodiments disclose a design method of a nonmetallic monoatomic doping BSe monolayer hydrogen evolution catalyst, wherein the nonmetallic monoatomic C, si replaces the selenium atom configuration in the BSe monolayer and has the highest catalytic activity.
In another aspect, some embodiments disclose a nonmetallic monoatomic doped BSe monolayer hydrogen evolution catalyst having a configuration of nonmetallic monoatoms, including C, N, si, P, as, substituted for selenium atoms in BSe monolayers.
In yet another aspect, some embodiments disclose the use of a nonmetallic monoatomic doped BSe monolayer hydrogen evolution catalyst, in particular, a nonmetallic monoatomic doped BSe monolayer hydrogen evolution catalyst for use in a water electrolysis hydrogen production reaction.
According to the design method of the nonmetallic monoatomic doping BSe monolayer hydrogen evolution catalyst disclosed by the embodiment of the application, a monolayer BSe supercell with a set specification is selected, a nonmetallic monoatomic doping BSe monolayer configuration is built on the monolayer BSe supercell, and the nonmetallic monoatomic doping BSe monolayer hydrogen evolution catalyst with hydrogen evolution catalytic activity is obtained by determining a structure stable in the nonmetallic monoatomic doping BSe monolayer configuration, so that the method has important guiding significance and practical significance in the process of preparing hydrogen by electrolyzing water.
Drawings
FIG. 1 is a schematic diagram of the structure and sites of a nonmetallic monoatomic doped monolayer BSe;
FIG. 2 is a schematic diagram of a monolayer BSe unit cell structure;
FIG. 3 is a graph showing the formation energy of a single layer of non-metallic monoatomic doping BSe in different configurations;
FIG. 4 is a band diagram of the configuration of selenium atoms in different nonmetallic monoatomic substituted BSe monolayers;
FIG. 5 shows Gibbs free energy diagrams of hydrogen evolution reactions with configuration of selenium atoms in different nonmetallic monoatomic substituted BSe monolayers;
FIG. 6 is a volcanic plot of hydrogen evolution reaction gibbs free energy versus exchange current for a selenium atom configuration in a different nonmetallic monoatomic substituted BSe monolayer;
FIG. 7 Gibbs free energy vs. p for selenium atom configuration in different nonmetallic monoatomic substituted BSe monolayers z A linear relationship graph with a center;
fig. 8 shows a fractional density plot and a fractional charge localization plot of selenium atom configuration in different nonmetallic monoatomic substituted BSe monolayers.
Detailed Description
The word "embodiment" as used herein does not necessarily mean that any embodiment described as "exemplary" is preferred or advantageous over other embodiments. The performance index test in the examples herein, unless otherwise specified, was performed using conventional test methods in the art. It is to be understood that the terminology used herein is for the purpose of describing particular embodiments only and is not intended to be limiting of the disclosure.
Unless otherwise defined, technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this application belongs; other test methods and techniques not specifically mentioned in this application are all those commonly used by those of ordinary skill in the art.
The terms "substantially" and "about" are used herein to describe small fluctuations. For example, they may refer to less than or equal to ±5%, such as less than or equal to ±2%, such as less than or equal to ±1%, such as less than or equal to ±0.5%, such as less than or equal to ±0.2%, such as less than or equal to ±0.1%, such as less than or equal to ±0.05%. Numerical data presented or represented herein in a range format is used only for convenience and brevity and should therefore be interpreted flexibly to include not only the numerical values explicitly recited as the limits of the range, but also to include all the individual numerical values or sub-ranges encompassed within that range. For example, a numerical range of "1 to 5%" should be interpreted to include not only the explicitly recited values of 1% to 5%, but also include individual values and sub-ranges within the indicated range. Thus, individual values, such as 2%, 3.5% and 4%, and subranges, such as 1% to 3%, 2% to 4% and 3% to 5%, etc., are included in this numerical range. The same principle applies to ranges reciting only one numerical value. Moreover, such an interpretation applies regardless of the breadth of the range or the characteristics being described.
In this document, including the claims, conjunctions such as "comprising," including, "" carrying, "" having, "" containing, "" involving, "" containing, "and the like are to be construed as open-ended, i.e., to mean" including, but not limited to. Only the conjunctions "consisting of … …" and "consisting of … …" are closed conjunctions.
Numerous specific details are set forth in the following examples in order to provide a better understanding of the present application. It will be understood by those skilled in the art that the present application may be practiced without some of these specific details. In the examples, some methods, means, instruments, devices, etc. well known to those skilled in the art are not described in detail in order to highlight the gist of the present application.
On the premise of no conflict, the technical features disclosed in the embodiments of the present application may be combined arbitrarily, and the obtained technical solution belongs to the disclosure of the embodiments of the present application.
In some embodiments, the design method of the nonmetallic monoatomic doping BSe monolayer hydrogen evolution catalyst comprises the following steps: (1) Selecting a monolayer BSe supercell with a set specification, and calculating the electronic configuration and the geometric structure of the supercell; (2) Constructing a nonmetallic monoatomic doping BSe monolayer configuration on the monolayer BSe supercell; (3) Determining a configuration of stable structure in the nonmetallic monoatomic doping BSe monolayer configuration; (4) The configuration that promotes the hydrogen evolution reaction in the configuration that determines the structural stability is a nonmetallic monoatomic doping BSe monolayer hydrogen evolution catalyst.
As an alternative embodiment, the monolayer BSe supercell selected in step (1) is a 3×3×1 monolayer BSe supercell comprising 36 atoms having a lattice constant ofThe thickness of the single layer is->B-Se bond length is->The bond length of B-B is->The supercell structure is shown in figure 1.
FIG. 2 shows a schematic structure of a BSe monolayer of 1X 1 unit cell, the top view is a top view, the side view is shown in side view, and the unit cell structure is expanded in the XY plane to obtain a 3X 1 monolayer BSe supercell.
As an alternative embodiment, in step (2), the nonmetallic monoatomic doping configuration includes a configuration in which the nonmetallic monoatoms adsorb to the surface sites of the BSe monolayer and a configuration in which the nonmetallic monoatoms replace the B or Se atoms in the BSe monolayer.
As an alternative embodiment, as shown in the schematic diagram of the nonmetallic monoatomic doping BSe site structure in fig. 1 (a), the upper diagram is a top view, the lower diagram is a side view, and the nonmetallic monoatomic doping configuration includes: substitution BSe of the selenium atom site configuration in the monolayer, which is marked as B-rich in the figure; substitution BSe of the configuration of the boron atom site in the monolayer, denoted Se-rich in the figure; the conformation of the site above BSe monolayer boron atoms is marked as B-top in the figure; the configuration of the site above the BSe monolayer selenium atom is adsorbed, and is marked as Se-top in the figure; and a configuration of a site above the center of the BSe monolayer hexagonal ring, labeled H-top in the figure. Fig. 1 (b) is a diagram of a BSe monolayer energy band structure without doped nonmetallic atoms.
As an alternative example, the formation of a doped nonmetallic monoatomic configuration can be calculated from the following formula, with the formation being able to evaluate the structural stability of the nonmetallic doping configuration:
E form =E tot [BSe * +X]-E tot [BSe]+xμ B +yμ Se -μ X
wherein E is tot [BSe*+X]For the total energy after doping with nonmetallic monoatoms, E tot [BSe]Mu, the total energy of undoped nonmetallic monoatoms X For each nonmetallic atom chemical potential, the coefficients x and y are the number of boron atoms and selenium atoms, respectively, μ B 、μ Se The chemical potentials of boron atoms and selenium atoms respectively, and X represents a substituted nonmetallic element.
BSe monolayer formation can be calculated by the following equation:
E form [BSe]=μ BSe -μ B 0 -μ Se 0
wherein mu BSe Mu, the total chemical potential of the block BSe B 0 、μ Se 0 Also of the chemical formula of each atom in the bulk B, se, respectively.
When doping nonmetallic monoatoms, the element proportion is changed, and two cases of boron-rich environment (B-rich) or selenium-rich environment (Se-rich) are mainly considered. The original selenium atom is replaced by a boron-rich environment, and the chemical potential of the boron atom is mu B B-rich =μ B 0 The method comprises the steps of carrying out a first treatment on the surface of the The selenium atom is replaced by original boron atom with chemical formula of mu Se Se-rich =μ Se 0 。
From thermodynamic equilibrium, one can assume: mu (mu) BSe =μ B +μ Se Thus, mu for boron-rich environments B B-rich 、 μ Se B-rich Calculated by the following formula:
mu for selenium-rich environment B Se-rich 、μ Se Se-rich Calculated by the following formula:
formation energy E form The more negative, the more readily the nonmetallic atoms are bonded to the monolayer BSe. As shown in fig. 3, the formation energy calculation results of doping different sites with different nonmetallic monoatoms show that the formation energy of the site configuration of the selenium atom in the substituted BSe monolayer is lower than that of other doping site configurations, which indicates that the nonmetallic monoatoms are more prone to doping at the site substituted for the selenium atom.
As an alternative example, in the non-metallic monoatomic doping configuration, the configuration of the selenium atom site in the substituted BSe monolayer is the most stable configuration, wherein the non-metallic monoatomic and its mass concentration is 0.258wt% C, 0.300wt% N, 0.343wt% O, 0.407wt% F, 0.600wt% Si, 0.662wt% P, 0.685wt% S, 0.757wt% Cl, 1.586wt% As, 1.689wt% Br, 2.671wt% Te or 2.656wt% I. As shown in fig. 4, the energy band analysis result shows that the doping configuration of the selenium atomic site in the C, si, N, P and As monoatomic substituted BSe single layer is equivalent to that of introducing a hole into the system, and the hole is a p-type semiconductor with an impurity state, wherein the p-type semiconductor is a fermi level, namely, the position of 0eV in the figure is close to the top of a valence band, and the impurity state is mainly donated by a dopant; F. the doping configuration of the selenium atom site in the Cl, br and I monoatomic substitution BSe monolayer is equivalent to introducing electrons into the system, resulting in the fermi level moving up into the rewind, and the energy band structure of the configuration doped with the single atom of the main group element O, S, te with the substituted Se atom is almost consistent with that of the undoped structure, which indicates that the single atom doping with the main group element does not effectively regulate the surface activity of the single atom.
As an alternative embodiment, the nonmetallic monoatoms C, N, si, P, as replace selenium atom configuration in BSe single layer, have hydrogen evolution catalytic activity and can be used as a hydrogen evolution reaction catalyst. As a more preferred embodiment, the nonmetallic monoatoms C, si replace the selenium atom configuration in BSe monolayers, have the highest catalytic activity and can be used as an excellent catalyst for hydrogen evolution reaction.
The hydrogen evolution reaction takes place by the Heyrovsky step [ H ] ads +(H + +e - )→H 2 ]Atomic to hydrogen combination to achieve optimal hydrogen evolution reaction performance Ji Pusi free energy Δg H The value will approach an ideal value deltag H =0 to reconcile the reaction disorder. Ji Pusi free energy ΔG H The definition is as follows:
in the above, E (H*) 、E (*) And E is (H2) Respectively, H (i.e., monolayer BSe adsorbs one hydrogen atom), x (i.e., monolayer BSe), and free H 2 Is set, is a function of the total energy of (1); e (E) ZPE And S respectively represents zero energy sum of hydrogen adsorption of different doping structuresEntropy, which can be calculated from NIST-JANAF thermochemical database and vibration frequency of adsorbed hydrogen atom configuration (H).
According to sabatier theory, the absolute value of the free energy of jeep is within 0.2, which is an ideal electrocatalyst, the corresponding overpotential is obtained by calculating the frequency of a doping system, and the result is shown in figure 5, wherein the Gibbs free energy chart shows that the monoatomic catalyst introduced with holes promotes hydrogen evolution reaction, and the monoatomic catalyst introduced with electrons inhibits hydrogen evolution reaction.
In some embodiments, the computational simulation is performed in neutral solution with absolute value of jeep free energy |Δg H* I calculates the switching current i 0 The kinetics of the hydrogen evolution reaction are described as:
wherein k is 0 、k B T is the reaction rate constant, the Boltzmann constant and the temperature under zero overpotential respectively, and represents the structure of a nonmetallic atom doping BSe; for ease of illustration we will refer to k 0 Let 1 be the value.
The activity of the hydrogen evolution reaction can be determined by the exchange current i relative to the volcanic peak 0 And ΔG H* The position of the value is measured. As shown in fig. 6, the electrocatalysts located on the left and right sides of the volcanic chart have negative and positive jeep free energy values, respectively, and the electrocatalyst closer to the standard value is located at the top of the volcanic peak. According to Sabatier theory, N, P and As monoatomic doped BSe single-layer catalyst is positioned on the left side of the volcanic diagram, the free energy value of Ji Pusi is negative, and the interaction between the system and adsorbed hydrogen atoms is strong; o, S, te, F, cl, br and I monatomic doped BSe monolayer catalysts and undoped BSe monolayer catalysts are located at the bottom right of the curve, with Ji Pusi free energy values being positive, indicating weak adsorption of hydrogen atoms to these catalysts. Ji Pusi free energy values of the C and Si monoatomic doped catalyst are close to standard values, and i is achieved 0 The maximum of the velocities, they are located near the peaks in the volcanic plot. The result shows that Si monoatomic doping BSe monolayer catalysisThe overpotential of the catalyst is 0.15eV and is close to the standard value, and the catalyst is an efficient hydrogen evolution reaction catalyst. C in FIG. 5 1 Represents a C monoatomic, si 1 Represents Si monoatoms, and the representation of the elemental atoms has similar meaning.
To further achieve a strong bonding between the monolayer system of BSe doped with different nonmetallic monoatoms and H, the nonmetallic atoms are p z Band center ε (p) z ) The definition is as follows:
wherein D (E) is a nonmetallic atom p z The DOS value for a given energy E.
As shown in FIG. 7, εp z Exhibits F, cl, br, I single-atom doped structure with negative epsilon p z Values, away from the fermi level, result in weaker adsorption to hydrogen atoms; C. si, N, P, as epsilon p of monoatomic doping structure z Near the fermi level, the binding to the hydrogen atom is stronger. Epsilon p of single-atom doped BSe monolayer catalyst with best hydrogen evolution reaction C and Si z Then it is moderate.
As shown in FIG. 8, for p-type semiconductor systems, such as C, si, N, P, as single-atom doped BSe monolayers, split-density PDOS and split-charge density PARCHG analysis, the impurity states introduced by doping are mainly determined by the p-type dopant z Track contribution. The s-orbitals of the isolated hydrogen atoms are located at the fermi level, the closer the impurity state is to the fermi level, the stronger the interaction with hydrogen. It is thus available that the N, P, as dopant is p due to its being located at the fermi level z The orbitals overlap with the s-orbitals of hydrogen so that they strongly interact, whereas the C and Si dopants p z The orbitals are located above the fermi level so that their adsorption strength with hydrogen is moderate. Therefore, the C and Si monoatomic doping BSe monolayer catalyst has high catalytic activity in hydrogen evolution reaction.
According to the design method of the nonmetallic monoatomic doping BSe monolayer hydrogen evolution catalyst, a monolayer BSe supercell with a set specification is selected, a nonmetallic monoatomic doping BSe monolayer structure is built on the monolayer BSe supercell, and the nonmetallic monoatomic doping BSe monolayer hydrogen evolution catalyst with hydrogen evolution catalytic activity is obtained by determining a structure stable in the nonmetallic monoatomic doping BSe monolayer structure, so that the method has important guiding significance and practical significance in the process of preparing hydrogen by electrolyzing water.
Technical details disclosed in the technical schemes and embodiments disclosed in the application are only illustrative of the inventive concepts of the application and are not limiting of the technical schemes of the application, and all conventional changes, substitutions or combinations of technical details disclosed in the application have the same inventive concepts as the application and are within the scope of protection of the claims of the application.
Claims (8)
1. The design method of the nonmetallic monoatomic doping BSe monolayer hydrogen evolution catalyst is characterized by comprising the following steps:
(1) Selecting a monolayer BSe supercell with a set specification, and calculating the electronic configuration and the geometric structure of the supercell; wherein the selected monolayer BSe supercell is a 3×3 monolayer BSe supercell containing 36 atoms and having a lattice constant ofThe thickness of the single layer is->B-Se bond length is->The bond length of B-B is->
(2) Constructing a nonmetallic monoatomic doping configuration on the monolayer BSe supercell;
(3) Determining a configuration of the nonmetallic monoatomic doping configuration, wherein the configuration is stable in structure; wherein the structural stability of the non-metal doped configuration is assessed by the formation energy calculated from the formula:
wherein E is tot [BSe*+X]For the total energy after doping with nonmetallic monoatoms, E tot [BSe]Mu, the total energy of undoped nonmetallic monoatoms X For each nonmetallic atom chemical potential, the coefficients x and y are the number of boron atoms and selenium atoms, respectively, μ B 、μ Se Chemical potentials of boron atoms and selenium atoms respectively, wherein X represents a substituted nonmetallic element;
(4) The configuration that promotes the hydrogen evolution reaction in the configuration that determines the structural stability is a nonmetallic monoatomic doping BSe monolayer hydrogen evolution catalyst.
2. The method for designing a nonmetallic monoatomic doped BSe monolayer hydrogen evolution catalyst according to claim 1, wherein in the step (2), the nonmetallic monoatomic doping configuration includes a configuration in which nonmetallic monoatoms adsorb to a surface site of the BSe monolayer and a configuration in which nonmetallic monoatoms replace B or Se atoms in the BSe monolayer.
3. The method for designing a nonmetallic monoatomic doping BSe monolayer hydrogen evolution catalyst according to claim 2, wherein the nonmetallic monoatomic doping configuration comprises: the configuration of selenium atoms in the substituted BSe monolayer, the configuration of boron atoms in the substituted BSe monolayer, the configuration of sites above boron atoms adsorbed on the BSe monolayer, the configuration of sites above selenium atoms adsorbed on the BSe monolayer and the configuration of sites above the center of the hexagonal ring of the BSe monolayer.
4. The method of designing a nonmetallic monoatomic doped BSe monolayer hydrogen evolution catalyst according to claim 3, wherein in the nonmetallic monoatomic doped configuration, the configuration of selenium atom sites in the substituted BSe monolayer is the most stable configuration, wherein the nonmetallic monoatomic and its mass concentration is 0.258wt% c, 0.300wt% n, 0.343wt% o, 0.407wt% f, 0.600wt% si, 0.662wt% p, 0.685wt% s, 0.757wt% cl, 1.586wt% as, 1.689wt% br, 2.671wt% te, or 2.656wt% i.
5. The method for designing a nonmetallic monoatomic doping BSe monolayer hydrogen evolution catalyst according to claim 4, wherein the nonmetallic monoatomic C, N, si, P, as replaces the selenium atom configuration in the BSe monolayer with hydrogen evolution catalytic activity.
6. The method for designing a nonmetallic monoatomic doping BSe monolayer hydrogen evolution catalyst according to claim 4, wherein the nonmetallic monoatomic C, si replaces the selenium atom configuration in the BSe monolayer with the highest catalytic activity.
7. The nonmetallic monoatomic doping BSe monolayer hydrogen evolution catalyst obtained by the design method according to any one of claims 1 to 6, wherein the nonmetallic monoatomic doping BSe monolayer hydrogen evolution catalyst has a selenium atom configuration in a nonmetallic monoatomic substitution BSe monolayer, and the nonmetallic monoatoms comprise C, N, si, P, as.
8. The use of the nonmetallic monoatomic doping BSe monolayer hydrogen evolution catalyst according to claim 7, wherein the nonmetallic monoatomic doping BSe monolayer hydrogen evolution catalyst is used for the electrolytic water hydrogen production reaction.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210040241.3A CN114381742B (en) | 2022-01-14 | 2022-01-14 | Nonmetallic monoatomic doping BSe monolayer hydrogen evolution catalyst, design method and application |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210040241.3A CN114381742B (en) | 2022-01-14 | 2022-01-14 | Nonmetallic monoatomic doping BSe monolayer hydrogen evolution catalyst, design method and application |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114381742A CN114381742A (en) | 2022-04-22 |
CN114381742B true CN114381742B (en) | 2024-03-12 |
Family
ID=81201957
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210040241.3A Active CN114381742B (en) | 2022-01-14 | 2022-01-14 | Nonmetallic monoatomic doping BSe monolayer hydrogen evolution catalyst, design method and application |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114381742B (en) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108573124A (en) * | 2018-05-15 | 2018-09-25 | 中国石油大学(华东) | A kind of method of quantitative analysis metal insertion carbon nano-tube catalyst hydrogen evolution activity |
CN109086564A (en) * | 2018-06-21 | 2018-12-25 | 太原理工大学 | A method of improving molybdenum disulfide catalytic hydrogen evolution performance |
KR20200040687A (en) * | 2018-10-10 | 2020-04-20 | 연세대학교 산학협력단 | A catalyst structure for hydrogen evolution and a method for designing the same |
KR20210053394A (en) * | 2019-11-01 | 2021-05-12 | 광주과학기술원 | Method for preparing monoatomic platinum catalyst having enhanced activity toward hydrogen evolution reaction(her) and monoatomic platinum catalyst prepared thereby |
CN113674812A (en) * | 2021-07-30 | 2021-11-19 | 哈尔滨工业大学 | Design method of metal-nitrogen-carbon nitrogen fixation electrocatalyst |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20200168300A1 (en) * | 2018-11-27 | 2020-05-28 | California Institute Of Technology | Screening methods and related catalysts, materials, compositions, methods and systems |
-
2022
- 2022-01-14 CN CN202210040241.3A patent/CN114381742B/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108573124A (en) * | 2018-05-15 | 2018-09-25 | 中国石油大学(华东) | A kind of method of quantitative analysis metal insertion carbon nano-tube catalyst hydrogen evolution activity |
CN109086564A (en) * | 2018-06-21 | 2018-12-25 | 太原理工大学 | A method of improving molybdenum disulfide catalytic hydrogen evolution performance |
KR20200040687A (en) * | 2018-10-10 | 2020-04-20 | 연세대학교 산학협력단 | A catalyst structure for hydrogen evolution and a method for designing the same |
KR20210053394A (en) * | 2019-11-01 | 2021-05-12 | 광주과학기술원 | Method for preparing monoatomic platinum catalyst having enhanced activity toward hydrogen evolution reaction(her) and monoatomic platinum catalyst prepared thereby |
CN113674812A (en) * | 2021-07-30 | 2021-11-19 | 哈尔滨工业大学 | Design method of metal-nitrogen-carbon nitrogen fixation electrocatalyst |
Non-Patent Citations (1)
Title |
---|
First-Principles Study on Transition-Metal Dichalcogenide/BSe van der Waals Heterostructures: A Promising Water-Splitting Photocatalyst;Yi Luo et al.;《J. Phys. Chem. C》;第123卷;第22742−22751页 * |
Also Published As
Publication number | Publication date |
---|---|
CN114381742A (en) | 2022-04-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Ji et al. | Graphene/MoS2/FeCoNi (OH) x and Graphene/MoS2/FeCoNiPx multilayer-stacked vertical nanosheets on carbon fibers for highly efficient overall water splitting | |
Xu et al. | Design of the synergistic rectifying interfaces in mott–schottky catalysts | |
Sun et al. | Catalysis of carbon dioxide photoreduction on nanosheets: fundamentals and challenges | |
Wang et al. | DFT study on sulfur-doped g-C3N4 nanosheets as a photocatalyst for CO2 reduction reaction | |
Huang et al. | Defect engineering of cobalt-based materials for electrocatalytic water splitting | |
Dilpazir et al. | Cobalt single atoms immobilized N-doped carbon nanotubes for enhanced bifunctional catalysis toward oxygen reduction and oxygen evolution reactions | |
Wang et al. | In situ preparation of Mo2C nanoparticles embedded in ketjenblack carbon as highly efficient electrocatalysts for hydrogen evolution | |
Pandey et al. | Interface engineering of an RGO/MoS2/Pd 2D heterostructure for electrocatalytic overall water splitting in alkaline medium | |
He et al. | Boron-doped C3N monolayer as a promising metal-free oxygen reduction reaction catalyst: a theoretical insight | |
Xiong et al. | Rational Bottom-Up Engineering of Electrocatalysts by Atomic Layer Deposition: A Case Study of Fe x Co1–x S y-Based Catalysts for Electrochemical Hydrogen Evolution | |
Pan | Metal dichalcogenides monolayers: novel catalysts for electrochemical hydrogen production | |
Zhao et al. | Carbon dots decorated hierarchical NiCo2S4/Ni3S2 composite for efficient water splitting | |
Liu et al. | Composition-tunable synthesis of “clean” syngas via a one-step synthesis of metal-free pyridinic-N-enriched self-supported CNTs: the synergy of electrocatalyst pyrolysis temperature and potential | |
Gong et al. | Thermally stable hierarchical nanostructures of ultrathin MoS2 nanosheet-coated CeO2 hollow spheres as catalyst for ammonia decomposition | |
Ma et al. | Molybdenum‐Based Co‐catalysts in Photocatalytic Hydrogen Production: Categories, Structures, and Roles | |
Garg et al. | Hexagonal CuCl monolayer for water splitting: a DFT study | |
Rhatigan et al. | Hydrogen evolution on non-metal oxide catalysts | |
Liang et al. | Design, synthesis and high HER performances of 3D Ni/Mo sulfide on Ni foam | |
Yan et al. | 2D Nanomaterial Supported Single‐Metal Atoms for Heterogeneous Photo/Electrocatalysis | |
Yu et al. | Uncovering the surface and phase effect of molybdenum carbides on hydrogen evolution: a first-principles study | |
Fang et al. | First-principles calculations on semiconducting ε-GeS and ε-SnS monolayer nanosheets with photocatalytic activity for sunlight-driven water splitting | |
Kwon et al. | Composition-tuned (MoWV) Se2 ternary alloy nanosheets as excellent hydrogen evolution reaction electrocatalysts | |
N Dhandapani et al. | Modulating the surface electronic structure of active Ni sites by engineering hierarchical NiFe-LDH/CuS over Cu foam as an efficient electrocatalyst for water splitting | |
Song et al. | New Way to Synthesize Robust and Porous Ni1–x Fe x Layered Double Hydroxide for Efficient Electrocatalytic Oxygen Evolution | |
Mahala et al. | Development of copper cobalt sulfide with Cu: Co ratio variation on carbon cloth as an efficient electrode material for the oxygen evolution reaction |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |