Disclosure of Invention
The invention aims to provide a method for constructing a two-dimensional Schottky heterojunction model based on germanium alkene.
In order to achieve the purpose, the technical scheme adopted by the invention is as follows:
the invention provides a method for constructing a two-dimensional Schottky heterojunction model based on germanium alkene, which comprises the following steps:
the model establishment is completed on Materials Studio 8.0 software, which is simulation software developed by Accelrys corporation specially for the field of material science and capable of running on a PC, and integrates various advanced methods and excellent modeling and visualization in the current molecular simulation field.
The first step is as follows: setting space group P6/mmm, and lattice constant as follows:
gamma 120 deg., establish WSe
2Bulk Material model, WSe
2The crystal has a layered structure;
the second step is that: the space group is set to P6/mmm,the lattice constant is:
b is a, height of corrugation
Gamma is 120 degrees, a Germanene model is established, and the Germanene crystal has a hexagonal honeycomb structure;
the third step: the WSe prepared in the first step2Uniformly cutting (001) crystal face on the bulk material model to obtain single-layer WSe2Uniformly cutting a (001) crystal face of the Germanene model prepared in the second step to obtain a Germanene thin layer;
the fourth step: establishing a 4X 4 or 3X 3Germanene supercell according to the Germanene thin layer prepared in the third step, and preparing a single layer WSe according to the third step2Establishing a 5 × 5 or 4 × 4WSe2The supercell of (4X 4) Germanene is vertically superposed on 5X 5WSe2Or 3X 3Germanene supercells are vertically superposed on 4X 4WSe2On the supercell to obtain Germanene/WSe2A heterojunction model.
The Germanene/WSe2The heterojunction model is in a single layer WSe2Germanene is in a single-layer fluctuant germanium honeycomb structure, Ge atoms are combined by covalent bonds, and a single layer WSe is formed2Is a sandwich structure consisting of Se-W-Se, and a covalent bond is formed between a W atom and a Se atom.
The Germanene/WSe
2The lattice mismatch ratio of the heterojunction model is defined as (a)
1-a
2)/a
1Wherein a is
1And a
2Are respectively WSe
2And the lattice constant of the Germanene optimized supercell. The optimization method comprises the following steps: the structure is optimized by adopting a Broyden-Fletcher-Goldfarb-Shanno (BFGS) scheme, the plane wave function cutoff energy is set to be 400eV, and the energy convergence standard is 5.0 x 10 in the process of structure relaxation
-7eV/atom, force per atom less than
The component of atomic displacement being less than
The stress is less than 0.02 GPa. Wherein the superlattice constant a (a ═ b) of 4X 4 or 3X 3Germane is
And
5X 5 or 4X 4WSe
2Has a superlattice constant a (a ═ b) of
And
because the thickness of the vacuum layer of the slice model is selected as the thickness in the calculation
The heterojunction superlattice lattice constant c is not explained.
The Germanene/WSe2The formation energy, the binding energy and the mismatch energy of the heterojunction model are respectively defined as:
wherein E is
formTo form energy, E
cohAs binding energy, Δ E
mismatchIn order to be able to mismatch the energy,
is Germanene/WSe
2Lattice of heterojunction modelNumber is
The total energy of (2) was-64842.67 eV. E
(Germanene)Is the Germanene lattice constant
The energy at time was-3419.12 eV.
Is WSe
2Lattice constant
The energy at time was-61417.82 eV. E
(Germanene)aIs the Germanene lattice constant
The energy at time was-3419.00 eV.
Is WSe
2Lattice constant
The energy at time was-61417.80 eV. A is Germanene/WSe
2Area of the heterojunction model of
Germanene/WSe from formula (1), formula (2) and formula (3)
2Formation energy of heterojunction model E
formIs composed of
Binding energy E
cohIs composed of
Mismatch energy Δ E
mismatchIs composed of
The Germanene/WSe2The magnitude of the interlayer force of the heterojunction model is determined by lattice mismatch energy and binding energy, and is defined as:
ΔEvdW=|ΔEmismach|+|Ecoh| (4)
ΔE
vdwthe energy caused by the interlayer acting force is calculated to obtain Delta E
vdwIs composed of
Due to the interlayer forces of most van der Waals force heterojunctions
Within the interval, it can therefore be concluded that Germanene/WSe
2The heterojunction model belongs to van der waals heterojunction, and the type of bonding between two layers is physical adsorption.
The Germanene/WSe
2The calculation of the electronic performance of the heterojunction model is completed by CAStep code based on the first principle of density functional DFT, the calculation is realized in Materials Studio software, the electronic exchange related interaction is processed by adopting generalized gradient approximation GGA and functional PBE, DFT-D approximation is realized by adopting a Tkatchenko and Scheffler (TS) scheme, the interaction of ions and electrons can be described by ultra-soft pseudopotentials, the structure is optimized by adopting a Broyden-Fletcher-Goldfarb-Shanno (BFBFGS) scheme, the plane wave function truncation energy is set to be 400eV, and the energy convergence standard is 5.0 multiplied by 10 eV in the process of structure relaxation
-7eV/atom, force per atom less than
The component of atomic displacement being less than
The stress is less than 0.02 GPa. The thickness of the vacuum layer of the lamella model is selected as
The Germanene/WSe2Two-dimensional way of using GGA-PBE form for energy band of heterojunction modelDiameter G (000) -M (00.50) -K (-0.333-0.6670) -G (000), single layer WSe2Is a direct band gap with a forbidden band width of 1.532eV, see fig. 2 a; the VBM of the top of the valence band and the CBM of the bottom of the conduction band are positioned at a point K; Germanene/WSe2The heterojunction model is a typical schottky heterojunction, forming a p-type schottky contact, where germane has a band gap at point K that is nearly zero small, see fig. 2 d. This is mainly due to Germanene and WSe2The interaction forces break the band gap that occurs with the sublattice symmetry of germane. Germanene/WSe2The band structures of the heterojunctions appear to be Germanene and single layer WSe2Is an important feature of two-dimensional vdW heterostructures, so that the electronic structure of the heterojunction largely retains germane and WSe2The layers are each independent of the electronic structure.
Germanene/WSe2Germanene and WSe in heterojunctions2The band dispersion relationship of (a) is essentially unchanged and reference is made to fig. 2a, b and c, which illustrate that the heterojunction maintains its respective independent electronic structure while generating new electronic properties.
The Germanene/WSe
2The interface spacing of the heterojunction model is
The band gap opening value is 10-45 meV; germanene and WSe
2The heterojunction is formed with high carrier mobility.
The Germanene/WSe2Metallic Germanene and semiconductor monolayer WSe in heterojunction model2Germanene/WSe at equilibrium interlayer spacing2Electrostatic potential of heterojunction in Z-direction, single-layer WSe2Have deeper potentials than Germanene, Germanene and WSe2A potential difference of electrostatic (potential drop) between up to 12.3 eV; Germanene/WSe2The heterojunction model preserves germanne's high carrier mobility.
The Germanene/WSe2In the heterojunction model, the effective mass is defined as:
where k is the wave vector,
is a simplified planck constant, m is the effective mass, and E is the energy. Under the action of different interface distances, the effective masses of electrons and holes of the heterojunction are shown in figure 3b, the visible interface distance and the electric field have certain regulation effect on the carrier mobility, and in an equilibrium state, the effective mass of the electrons is 0.57m
0Effective mass of holes is-0.56 m
0When the interface spacing is
The effective mass of the electron and hole is 0.70m, respectively
0And-0.73 m
0And thus has a higher carrier mobility.
Due to the adoption of the technical scheme, the invention has the following advantages and beneficial effects:
according to the method for constructing the two-dimensional Schottky heterojunction model based on the germanium alkene, the two-dimensional Schottky heterojunction model based on the germanium alkene is suitable for a Schottky heterojunction system, especially a two-dimensional Schottky heterojunction system; the theoretical model is composed of two-dimensional materials, and is formed by combining WSe2Cutting (001) crystal face of the Germanene material to obtain Germanene/WSe2A heterojunction model; the theoretical model considers lattice mismatch of a heterojunction model, establishes a supercell and analyzes the formation stability of a heterojunction; the theoretical model adopts a first principle based on a Density Functional (DFT) to disclose the influence mechanism of an energy band structure and a layer spacing on the effective quality of a heterojunction.
The invention vertically stacks germanium alkene (Germanene) on a single layer of tungsten diselenide (WSe) for the first time2) And a density functional theory based on a first sexual principle is adopted for Germanene/WSe2The heterojunction is researched, lattice parameters, energy band structures, thermodynamic stability and the like of the supercell are calculated, and meanwhile, the influence mechanism of the interface spacing on the effective quality of the heterojunction is provided.
Detailed Description
In order to more clearly illustrate the invention, the invention is further described below in connection with preferred embodiments. It is to be understood by persons skilled in the art that the following detailed description is illustrative and not restrictive, and is not to be taken as limiting the scope of the invention.
Example 1
A construction method of a two-dimensional Schottky heterojunction model based on germanium alkene comprises the following steps:
the model establishment is completed on Materials Studio 8.0 software, which is simulation software developed by Accelrys corporation specially for the field of material science and capable of running on a PC, and integrates various advanced methods and excellent modeling and visualization in the current molecular simulation field.
The first step is as follows: setting space group P6/mmm, and lattice constant as follows:
b=a,
gamma 120 deg., establish WSe
2Bulk Material model, WSe
2The crystal has a layered structure;
the second step is that: setting space group as P6/mmm, lattice constant as:
b is a, height of corrugation
Gamma is 120 degrees, a Germanene model is established, and the Germanene crystal has a hexagonal honeycomb structure;
the third step: the WSe prepared in the first step2Uniformly cutting (001) crystal face on the bulk material model to obtain single-layer WSe2Uniformly cutting a (001) crystal face of the Germanene model prepared in the second step to obtain a Germanene thin layer;
the fourth step: establishing a 4X 4 or 3X 3Germanene supercell according to the Germanene thin layer prepared in the third step, and preparing a single layer WSe according to the third step2Establishing a 5 × 5 or 4 × 4WSe2The supercell of (4X 4) Germanene is vertically superposed on 5X 5WSe2Or 3X 3Germanene supercells are vertically superposed on 4X 4WSe2On the supercell to obtain Germanene/WSe2A heterojunction model.
The Germanene/WSe2The heterojunction model is in a single layer WSe2Germanene is in a single-layer fluctuant germanium honeycomb structure, Ge atoms are combined by covalent bonds, and a single layer WSe is formed2Is a sandwich structure consisting of Se-W-Se, and a covalent bond is formed between a W atom and a Se atom. As shown in FIG. 1, FIG. 1 shows Germanene/WSe prepared according to an embodiment of the present invention2Structural schematic diagram of heterojunction modelIn the figure, a is a supercell of 3X 3Germane vertically superposed on 4X 4WSe2Germanene/WSe obtained from Supercell of (1)2A heterojunction model schematic diagram; supercells with b 4X 4Germanene vertically superimposed on 5X 5WSe2Germanene/WSe obtained from Supercell of (1)2A schematic top view of a heterojunction model; supercell with c being Germanene vertically superposed on 5X 5WSe2Germanene/WSe obtained from Supercell of (1)2The heterojunction model is schematically shown from the side.
The lattice constants a (a ═ b) of the heterojunction supercells in a and b in FIG. 1 are respectively
And
the combination of heterojunctions with a lattice mismatch of less than 5%, the Germanene/WSe, is commonly referred to as lattice matching
2The lattice mismatch ratio of the heterojunction model is defined as (a)
1-a
2)/a
1Wherein a is
1And a
2Are respectively WSe
2And the lattice constant of the Germanene optimized supercell. The optimization method comprises the following steps: the structure is optimized by adopting a Broyden-Fletcher-Goldfarb-Shanno (BFGS) scheme, the plane wave function cutoff energy is set to be 400eV, and the energy convergence standard is 5.0 x 10 in the process of structure relaxation
-7eV/atom, force per atom less than
The component of atomic displacement being less than
The stress is less than 0.02 GPa. Wherein the superlattice constant a (a ═ b) of 4X 4 or 3X 3Germane is
And
5X 5 or 4X 4WSe
2Has a superlattice constant a (a ═ b) of
And
because the thickness of the vacuum layer of the slice model is selected as the thickness in the calculation
The heterojunction superlattice lattice constant c is not explained. The lattice mismatch ratios were calculated to be 7.2% and 1.1%, respectively. Thus, the supercells of 4X 4Germanene are vertically stacked at 5X 5WSe
2Germanene/WSe obtained from Supercell of (1)
2The heterojunction model can form a better match and the following analysis is based on the model of fig. 1 b.
In order to examine the difficulty and stability of heterojunction formation in the experiment, Germanene/WSe was calculated2The formation energy, the binding energy and the mismatch energy of the heterojunction model are respectively defined as:
wherein E is
formTo form energy, E
cohAs binding energy, Δ E
mismatchIn order to be able to mismatch the energy,
is Germanene/WSe
2The lattice constant of the heterojunction model is
The total energy of (2) was-64842.67 eV. E
(Germanene)Is a Germanene latticeConstant number
The energy at time was-3419.12 eV.
Is WSe
2Lattice constant
The energy at time was-61417.82 eV. E
(Germanene)aIs the Germanene lattice constant
The energy at time was-3419.00 eV.
Is WSe
2Lattice constant
The energy at time was-61417.80 eV. A is Germanene/WSe
2Area of the heterojunction model of
Germanene/WSe from formula (1), formula (2) and formula (3)
2Formation energy of heterojunction model E
formIs composed of
Binding energy E
cohIs composed of
Mismatch energy Δ E
mismatchIs composed of
Negative formation and binding energies indicate Germanene/WSe
2The heterojunction model is an easily formed and stable heterojunction. This result is similar in magnitude to that of two-dimensional heterostructures, such as graphene/MoS
2Graphene/phosphene and graphene/g-GaN.
The Germanene/WSe2The magnitude of the interlayer force of the heterojunction model is determined by lattice mismatch energy and binding energy, and is defined as:
ΔEvdW=|ΔEmismach|+|Ecoh| (4)
ΔE
vdwthe energy caused by the interlayer acting force is calculated to obtain Delta E
vdwIs composed of
Due to the interlayer forces of most van der Waals force heterojunctions
Within the interval, it can therefore be concluded that Germanene/WSe
2The heterojunction model belongs to van der waals heterojunction, and the type of bonding between two layers is physical adsorption.
The Germanene/WSe
2The calculation of the electronic performance of the heterojunction model is completed by CAStep code based on the first principle of density functional DFT, is realized in Materials Studio software, adopts generalized gradient approximation GGA and functional PBE to process the electronic exchange related interaction, and because Germanene and WSe
2The strong bonding interaction is lacked, and the weak vdW interaction is expected to play an important role. In order to obtain the accurate structure of the heterostructure in calculation, DFT-D approximation is realized by adopting a Tkatchenko and Scheffler (TS) scheme, the interaction of ions and electrons can be described by ultra-soft pseudopotential, the structure is optimized by adopting a Broyden-Fletcher-Goldfarb-Shanno (BFGS) scheme, the plane wave function truncation energy is set to be 400eV, and the energy convergence standard is 5.0 x 10 in the process of structure relaxation
-7eV/atom, force per atom less than
The component of atomic displacement being less than
The stress is less than 0.02 GPa. The thickness of the vacuum layer of the lamella model is selected as
As shown in FIG. 2, FIG. 2 shows Germanene/WSe prepared according to an embodiment of the present invention2Energy band diagrams of the heterojunction model, a, b and c are single-layer WSe2Germanene and Germanene/WSe2Band diagram of heterojunction model, d is Germane/WSe2The heterojunction model band is magnified at the dirac point. The Germanene/WSe2The energy band of the heterojunction model uses the two-dimensional pathway of GGA-PBE form G (000) -M (00.50) -K (-0.333-0.6670) -G (000), and it can be seen that the WSe of a single layer is2Is a direct band gap with a forbidden band width of 1.532eV, see fig. 2 a; the VBM of the top of the valence band and the CBM of the bottom of the conduction band are positioned at a point K; the pi and pi-bands of germanne cross at dirac points and cross the fermi level with a linear dirac cone, thus having semi-metallic properties, similar to graphene. As compared to Germanene alone, it can be seen that Germanene/WSe2The heterojunction model is a typical schottky heterojunction, forming a p-type schottky contact, where the dirac cone of germane is well preserved. At the same time, germane was found to have a band gap at the K point which is small, almost zero, as shown in the graph d. This is mainly due to Germanene and WSe2The interaction forces break the band gap that occurs with the sublattice symmetry of germane. Germanene/WSe2The band structures of the heterojunctions appear to be Germanene and single layer WSe2Is an important feature of two-dimensional vdW heterostructures, so that the electronic structure of the heterojunction largely retains germane and WSe2The layers are each independent of the electronic structure. Germanene/WSe2Germanene and WSe in heterojunctions2The band dispersion relationship of (a) is essentially unchanged and reference is made to fig. 2a, b and c, which illustrate that the heterojunction maintains its respective independent electronic structure while generating new electronic properties.
The Germanene/WSe
2The interface spacing of the heterojunction model is
The band gap opening value is 10-45 meV; germanene and WSe
2The heterojunction is formed with high carrier mobility. As shown in FIG. 3, FIG. 3 shows Germanene/WSe prepared according to an embodiment of the present invention
2The band gap and electron hole effective mass curve chart of the heterojunction model, a is Germanene/WSe under different interface pitches
2Band gap of heterojunction, b is Germanene/WSe under different interface pitches
2Electron hole effective mass of the heterojunction. From a, the energy dispersion relation of the conduction band bottom and valence band top near the Fermi level Dirac point can be seen at different interface pitches. The position of the pi band of Germanene is unchanged, and the Dirac cone angle formed by the pi band and the pi band is gradually reduced. Variation in the position of the pi bands of Germanene, spacing at the interface
When the band gap of Germanene is the farthest from the Fermi level, the band gap opening value is 45meV, and when the interface distance is
At this time, the band gap opening value decreases to a minimum value of 10meV, and increases from 10meV to 27meV (interface spacing: 10 meV) with an increase in interface spacing
) Then reduced to 17meV (interfacial spacing of
). Accordingly, Germanene/WSe
2The band gap opening value of the heterojunction at the dirac point is very sensitive to the lattice symmetry of the heterojunction and is influenced by the interface spacing. From b, it can be seen that when the interface spacing is
The effective mass of electrons and holes is at their maximum, 0.70m each
0And-0.73 m
0Is increased to
The effective mass of electrons and holes decreases and remains substantially unchanged. When in flatAt the equilibrium layer spacing, the effective masses of electrons and holes at the K Dirac point were calculated to be 0.53 and 0.51m
0Specification of Germanene and WSe
2The heterojunction has high carrier mobility when formed, so that the heterojunction becomes a material suitable for high-speed nano electronic and optoelectronic devices.
The Germanene/WSe2In the heterojunction model, the effective mass is defined as:
where k is the wave vector,
is a simplified planck constant. Under the action of different interface distances, the effective masses of electrons and holes of the heterojunction are shown in figure 3b, the visible interface distance and the electric field have certain regulation effect on the carrier mobility, and in an equilibrium state, the effective mass of the electrons is 0.57m
0Effective mass of holes is-0.56 m
0When the interface spacing is
The effective mass of the electron and hole is 0.70m, respectively
0And-0.73 m
0And thus has a higher carrier mobility.
FIG. 4 shows Germanene/WSe2Electrostatic potential of the heterojunction in the Z-direction is schematically shown. The Germanene/WSe2Metallic Germanene and semiconductor monolayer WSe in heterojunction model2Germanene/WSe at equilibrium interlayer spacing2Electrostatic potential of heterojunction in Z-direction, single-layer WSe2Has a deeper potential than Germanene, and in addition, due to Germanene and WSe2The electrostatic potential difference (potential drop) between them is as high as 12.3eV, so that a strong electrostatic field exists across the van der waals heterojunction, which has a large influence on the injection of charges when the germane layer is used as an electrode. Germanene/WSe2The heterojunction maintains the high carrier mobility of Germanene.
Although the present invention has been described with reference to a preferred embodiment, it should be understood that various changes, substitutions and alterations can be made herein without departing from the spirit and scope of the invention as defined by the appended claims.