CN107798213B - Rapid simulation calculation method for palladium metal catalytic performance change - Google Patents

Rapid simulation calculation method for palladium metal catalytic performance change Download PDF

Info

Publication number
CN107798213B
CN107798213B CN201711000044.4A CN201711000044A CN107798213B CN 107798213 B CN107798213 B CN 107798213B CN 201711000044 A CN201711000044 A CN 201711000044A CN 107798213 B CN107798213 B CN 107798213B
Authority
CN
China
Prior art keywords
particle
palladium metal
electrons
metal atoms
lithium iron
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.)
Expired - Fee Related
Application number
CN201711000044.4A
Other languages
Chinese (zh)
Other versions
CN107798213A (en
Inventor
刘悦
庄伟彬
付凯
朱浩楠
孙跃军
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Liaoning Technical University
Original Assignee
Liaoning Technical University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Liaoning Technical University filed Critical Liaoning Technical University
Priority to CN201711000044.4A priority Critical patent/CN107798213B/en
Publication of CN107798213A publication Critical patent/CN107798213A/en
Application granted granted Critical
Publication of CN107798213B publication Critical patent/CN107798213B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16CCOMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
    • G16C20/00Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
    • G16C20/80Data visualisation
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16CCOMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
    • G16C20/00Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
    • G16C20/10Analysis or design of chemical reactions, syntheses or processes

Landscapes

  • Engineering & Computer Science (AREA)
  • Chemical & Material Sciences (AREA)
  • Crystallography & Structural Chemistry (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Computing Systems (AREA)
  • Theoretical Computer Science (AREA)
  • Data Mining & Analysis (AREA)
  • Analytical Chemistry (AREA)
  • Chemical Kinetics & Catalysis (AREA)
  • Catalysts (AREA)

Abstract

The invention provides a rapid simulation calculation method for palladium metal catalytic performance change, which comprises the following steps: constructing a ball stick model of palladium metal atoms and lithium iron niobate electrons; calculating the total energy value of the ball stick model of palladium metal atoms and lithium iron niobate electrons; determining interaction force, interaction speed and interaction position among the particles by adopting a dissipative particle dynamics method; setting the number of each particle and the deposition temperature; performing two-dimensional treatment on lithium iron niobate electrons, namely constructing a two-dimensional grid of N x N, and distributing each particle of the lithium iron niobate electrons on the intersection points of the two-dimensional grid; determining the covering positions of palladium metal atoms and two-dimensional lithium niobate ferroelectrons at different deposition temperatures according to the interaction force among the particles; and calculating the coverage rate value of the palladium metal atoms covered on the lithium iron niobate electrons at different temperatures. The surface area change trend of palladium metal atoms on different ferroelectric surfaces at different temperatures can be effectively predicted, so that the change trend of the catalytic performance of the palladium metal atoms along with the temperature can be predicted.

Description

Rapid simulation calculation method for palladium metal catalytic performance change
Technical Field
The invention belongs to the technical field of material science numerical simulation, and particularly relates to a rapid simulation calculation method for palladium metal catalytic performance change.
Background
Carrier research of supported metal catalysts is a hot research topic in the current scientific community, and some commonly used catalyst carriers researched at present comprise graphene, carbon nanofibers, gamma-Al 2O3, CeO2 nanotubes, high-molecular carriers and the like, and ferroelectrics are widely applied to the field of electronic devices, but the shape prediction of gold atom aggregation clusters formed on the surface of the ferroelectrics is few. Less research has been conducted on the use of ferroelectrics as supports, and we have conducted a simulation calculation study on the effect of using ferroelectrics as supports on the catalytic activity of catalytic metals.
Disclosure of Invention
Aiming at the defects of the prior art, the invention provides a rapid simulation calculation method for the catalytic performance change of palladium metal.
A rapid simulation calculation method for palladium metal catalytic performance change comprises the following steps:
step 1: constructing a ball stick model of palladium metal atoms and lithium iron niobate electrons;
step 2: calculating the total energy value of the ball stick model of palladium metal atoms and lithium iron niobate electrons;
and step 3: setting coordination number values of palladium metal atoms and lithium iron niobate electrons, and determining interaction force, interaction speed and interaction position among particles by adopting a dissipative particle dynamics method;
and 4, step 4: setting the number of each particle and the deposition temperature;
and 5: performing two-dimensional treatment on lithium iron niobate electrons, namely constructing a two-dimensional grid of N x N, and distributing each particle of the lithium iron niobate electrons on the intersection points of the two-dimensional grid;
step 6: determining the covering positions of palladium metal atoms and two-dimensional lithium niobate ferroelectrons at different deposition temperatures according to the interaction force among the particles;
and 7: and calculating the coverage rate value of the palladium metal atoms covered on the lithium iron niobate electrons at different temperatures.
The method for measuring the total energy value of the ball stick model of palladium metal atoms and lithium iron niobate electrons comprises the following steps: and (4) calculating the total energy value of the ball-rod model of the palladium metal atoms and the lithium iron niobate electrons by using a density functional theory calculation method.
The calculation formula for determining the interaction force, the interaction speed and the interaction position among the particles by using the dissipative particle dynamics method is as follows:
Figure BDA0001443154680000011
Figure BDA0001443154680000012
Figure BDA0001443154680000021
wherein,
Figure BDA0001443154680000022
aij(T)=aii+1.451χijkBT;
aii=15kBT;
Figure BDA0001443154680000023
Figure BDA0001443154680000024
Fijis the interaction force between particles i, j, vijIs the mutual velocity between particles i, j, rijIs the mutual position between the particles i, j,
Figure BDA0001443154680000025
is the dissipative force between the particles i, j,
Figure BDA0001443154680000026
is a random force between particles i, j, aijIs the maximum repulsive force between particles i, j, T is the temperature, kBR is a constant, ZijThe coordination number value of particle i relative to particle j, Eij(T) is the energy value of particle i relative to particle j, ZjiIs the coordination number value of particle j relative to particle i, Eji(T) is the energy value of particle j relative to particle i, ZiiIs the coordination number value of particle i, Eii(T) is the energy value of the particle i, ZjjIs the coordination number value of particle j, Ejj(T) is the energy value of particle j.
The invention has the beneficial effects that:
the invention provides a rapid simulation calculation method for the change of the catalytic performance of palladium metal, which can obviously observe that the coverage rate of palladium atoms increases along with the increase of temperature, and observe coverage rate curves obtained by each data, wherein the coverage rate curves obtained in two polarization directions have a tendency of folding along with the increase of temperature, so that the surface area change tendency of the palladium metal atoms at different temperatures can be effectively predicted, the surface area change tendency of the palladium metal atoms on different ferroelectric surfaces can be predicted, and the change tendency of the catalytic performance of the palladium metal atoms along with the temperature can be predicted.
Drawings
Fig. 1 is a flowchart of a rapid simulation calculation method for palladium metal catalytic performance change according to an embodiment of the present invention.
FIG. 2 is a two-dimensional view of palladium metal atoms at 500 and 600 ℃ with a two-dimensional lithium niobate ferroelectric coated location in accordance with an embodiment of the present invention;
FIG. 3 is a graph of palladium metal atoms and two-dimensional lithium iron niobate electron coverage at a temperature of 500 deg.C-1000 deg.C with 100 palladium metal atoms in accordance with an embodiment of the present invention;
FIG. 4 is a graph of palladium metal atoms and two-dimensional lithium iron niobate electron coverage at a temperature of 500 deg.C-1000 deg.C with 200 palladium metal atoms in accordance with an embodiment of the present invention;
FIG. 5 is a plot of palladium metal atoms at 500 deg.C-1000 deg.C and two-dimensional lithium iron niobate electron coverage for a palladium metal atom of 300 deg.C in accordance with an embodiment of the present invention;
FIG. 6 is a plot of palladium metal atoms at 400 deg.C and at 500 deg.C-1000 deg.C versus two-dimensional lithium iron niobate electron coverage for a specific embodiment of the present invention;
FIG. 7 is a graph of palladium metal atoms and bidimensionalized lithium iron niobate electron coverage at 500 deg.C-1000 deg.C with 500 palladium metal atoms in accordance with an embodiment of the present invention.
Detailed Description
The following detailed description of embodiments of the invention refers to the accompanying drawings.
A rapid simulation calculation method for palladium metal catalytic performance change is shown in figure 1, and comprises the following steps:
step 1: and constructing a ball stick model of palladium metal atoms and lithium iron niobate electrons.
In this embodiment, a polarized molecular structure of lithium iron niobate electrons (L innbo 3) is determined, a molecular structure diagram thereof is drawn by using Gauss software, and atoms are moved according to molecular polarization to obtain desired polarized molecules, thereby forming a ferroelectric substrate.
Step 2: and (4) calculating the total energy value of the ball stick model of the palladium metal atoms and the lithium iron niobate electrons.
In this embodiment, the total energy value of the ball-and-stick model of palladium metal atoms and lithium iron niobate electrons is calculated by using a study result that palladium metal atoms approach a ferroelectric substrate at a distance of 0.5 to 1000, calculating energy using Gauss energy calculation software, and calculating a Density Functional Theory (DFT) calculation method.
And step 3: and setting the coordination number value of palladium metal atoms and lithium iron niobate electrons, and determining the interaction force, the interaction speed and the interaction position among the particles by adopting a dissipative particle dynamics method.
In this embodiment, a Dissipative Particle Dynamics (DPD) is a molecular dynamics simulation algorithm for simulating the behavior of a complex fluid. The DPD method meets the condition of thermal equilibrium and is a non-lattice model method for simulating the movement of particles in continuous space and discontinuous time. The DPD method is one in which a single particle represents the entire molecule or a region of fluid containing multiple molecules, or a fragment of a macromolecule, rather than a single atom, and does not consider the behavioral details of an atom, which is considered process independent.
In this embodiment, the coordination number of the palladium metal atom and the lithium iron niobate electron is 1 to 10.
The calculation formula for determining the interaction force, the interaction speed and the interaction position among the particles by adopting a dissipative particle dynamics method is shown as formulas (1), (2) and (3):
Figure BDA0001443154680000031
Figure BDA0001443154680000032
Figure BDA0001443154680000041
wherein,
Figure BDA0001443154680000042
aij(T)=aii+1.451χijkBT;
aii=15kBT;
Figure BDA0001443154680000043
Figure BDA0001443154680000044
Fijis the interaction force between particles i, j, vijIs the mutual velocity between particles i, j, rijIs the mutual position between the particles i, j,
Figure BDA0001443154680000045
is the dissipative force between the particles i, j,
Figure BDA0001443154680000046
is a random force between particles i, j, aijIs the maximum repulsive force between particles i, j, T is the temperature, kBR is a constant, ZijThe coordination number value of particle i relative to particle j, Eij(T) is the energy value of particle i relative to particle j, ZjiIs the coordination number value of particle j relative to particle i, Eji(T) is the energy value of particle j relative to particle i, ZiiIs the coordination number value of the particle i,Eii(T) is the energy value of the particle i, ZjjIs the coordination number value of particle j, Ejj(T) is the energy value of particle j.
And 4, step 4: the number of particles and the deposition temperature were set.
In the present embodiment, the number of particles is graded from 100 to 1000 with 100 as a reference. Temperature selection: the span is 500-1000 deg.C, and the temperature gradient is 50 deg.C.
And 5: and (3) performing two-dimensional treatment on the lithium iron niobate electrons, namely constructing a two-dimensional grid of N x N, and distributing each particle of the lithium iron niobate electrons on the intersection points of the two-dimensional grid.
In this embodiment, a 100 × 100 two-dimensional grid is constructed, and each particle of lithium niobate electrons is distributed at an intersection of the two-dimensional grid.
Step 6: and determining the covering positions of palladium metal atoms and two-dimensional lithium niobate ferroelectrons at different deposition temperatures according to the interaction force among the particles.
In this embodiment, a Mathematica software is used to program and simulate a surface coverage view, a series of two-dimensional views are respectively obtained by taking the deposition temperature and the number of candidate particles as gradients, then the palladium metal atom coverage rate in the two-dimensional views is observed to represent the dispersion degree of palladium atoms in L iNbO3 polarized molecules as the substrate surface, and the catalytic activity of palladium atom catalytic metals is represented by the dispersion degree.
The two-dimensional view of the palladium metal atom and the two-dimensional lithium niobate ferroelectric electron coverage position at a temperature of 600 ℃ with 500 palladium metal atoms is shown in fig. 2.
And 7: and calculating the coverage rate value of the palladium metal atoms covered on the lithium iron niobate electrons at different temperatures.
In this embodiment, fig. 3-7 are graphs showing the coverage rate of palladium metal atoms and two-dimensional lithium iron niobate electrons at 500-1000 ℃ when the palladium metal atoms are 100-500-.

Claims (2)

1. A rapid simulation calculation method for palladium metal catalytic performance change is characterized by comprising the following steps:
step 1: constructing a ball stick model of palladium metal atoms and lithium iron niobate electrons;
step 2: calculating the total energy value of the ball stick model of palladium metal atoms and lithium iron niobate electrons;
and step 3: setting coordination number values of palladium metal atoms and lithium iron niobate electrons, and determining interaction force, interaction speed and interaction position among particles by adopting a dissipative particle dynamics method;
the calculation formula for determining the interaction force, the interaction speed and the interaction position among the particles by using the dissipative particle dynamics method is as follows:
Figure FDA0002119397190000011
Figure FDA0002119397190000012
Figure FDA0002119397190000013
wherein,
Figure FDA0002119397190000014
aij(T)=aii+1.451χijkBT;
aii=15kBT;
Figure FDA0002119397190000015
Figure FDA0002119397190000016
Fijis the interaction force between particles i, j, vijIs the mutual velocity between particles i, j, rijIs the mutual position between the particles i, j,
Figure FDA0002119397190000017
is the dissipative force between the particles i, j,
Figure FDA0002119397190000018
is a random force between particles i, j, aijIs the maximum repulsive force between particles i, j, T is the temperature, kBR is a constant, ZijThe coordination number value of particle i relative to particle j, Eij(T) is the energy value of particle i relative to particle j, ZjiIs the coordination number value of particle j relative to particle i, Eji(T) is the energy value of particle j relative to particle i, ZiiIs the coordination number value of particle i, Eii(T) is the energy value of the particle i, ZjjIs the coordination number value of particle j, Ejj(T) is the energy value of particle j; and 4, step 4: setting the number of each particle and the deposition temperature;
and 5: performing two-dimensional treatment on lithium iron niobate electrons, namely constructing a two-dimensional grid of N x N, and distributing each particle of the lithium iron niobate electrons on the intersection points of the two-dimensional grid;
step 6: determining the covering positions of palladium metal atoms and two-dimensional lithium niobate ferroelectrons at different deposition temperatures according to the interaction force among the particles;
and 7: and calculating the coverage rate value of the palladium metal atoms covered on the lithium iron niobate electrons at different temperatures.
2. The method for rapidly simulating and calculating the catalytic performance change of palladium metal according to claim 1, wherein the method for calculating the total energy value of the ball-and-stick model of palladium metal atoms and lithium iron niobate electrons specifically comprises the following steps:
and (4) calculating the total energy value of the ball-rod model of the palladium metal atoms and the lithium iron niobate electrons by using a density functional theory calculation method.
CN201711000044.4A 2017-10-24 2017-10-24 Rapid simulation calculation method for palladium metal catalytic performance change Expired - Fee Related CN107798213B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711000044.4A CN107798213B (en) 2017-10-24 2017-10-24 Rapid simulation calculation method for palladium metal catalytic performance change

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711000044.4A CN107798213B (en) 2017-10-24 2017-10-24 Rapid simulation calculation method for palladium metal catalytic performance change

Publications (2)

Publication Number Publication Date
CN107798213A CN107798213A (en) 2018-03-13
CN107798213B true CN107798213B (en) 2020-07-28

Family

ID=61534465

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711000044.4A Expired - Fee Related CN107798213B (en) 2017-10-24 2017-10-24 Rapid simulation calculation method for palladium metal catalytic performance change

Country Status (1)

Country Link
CN (1) CN107798213B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109002681B (en) * 2018-05-25 2021-08-27 胡培君 Method for simulating/predicting catalytic activity of catalyst in surface catalytic reaction and application thereof

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103678875A (en) * 2013-11-04 2014-03-26 北京理工大学 Prediction method of dispersion behaviors of flame retardant and flame retardant synergist in polymer matrix
CN104268405A (en) * 2014-09-26 2015-01-07 安徽大学 Monte carlo molecular simulation research method for kinetic process of polymerization reaction
CN105138744A (en) * 2015-08-07 2015-12-09 长春理工大学 Dissipative particle dynamics simulation method for grinding liquid particle characteristics

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102614860B (en) * 2012-02-16 2013-10-16 北京工业大学 Low-temperature plasma catalyst, and preparation method and application thereof
CN103785390B (en) * 2014-02-18 2016-03-02 北京工业大学 A kind of ferroelectric-tourmaline plasma catalyst, preparations and applicatio

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103678875A (en) * 2013-11-04 2014-03-26 北京理工大学 Prediction method of dispersion behaviors of flame retardant and flame retardant synergist in polymer matrix
CN104268405A (en) * 2014-09-26 2015-01-07 安徽大学 Monte carlo molecular simulation research method for kinetic process of polymerization reaction
CN105138744A (en) * 2015-08-07 2015-12-09 长春理工大学 Dissipative particle dynamics simulation method for grinding liquid particle characteristics

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Molecular dynamic simulations of surface morphology and pulsed laser deposition growth of lithium niobate thin films on silicon substrate;Liu Yue et al;《Chin Phys B》;20150531;第24卷(第5期);第1-4页 *
钯金属纳米颗粒催化剂在铁电表面的分散聚集模拟;刘悦 等;《第九届计算纳米科学与新能源材料国际研讨会》;20160630;第158页 *

Also Published As

Publication number Publication date
CN107798213A (en) 2018-03-13

Similar Documents

Publication Publication Date Title
Lu et al. Intervalley scattering and localization behaviors of spin-valley coupled dirac fermions
Christofides et al. Control and optimization of multiscale process systems
Mortazavi et al. Multiscale modeling of thermal conductivity of polycrystalline graphene sheets
Wu et al. Lattice mismatch induced nonlinear growth of graphene
Zhang et al. Controllable fabrication of patterned ZnO nanorod arrays: investigations into the impacts on their morphology
Jabar et al. Magnetic properties of graphene structure: a Monte Carlo simulation
Khosravian et al. Effects of a grain boundary loop on the thermal conductivity of graphene: A molecular dynamics study
CN107798213B (en) Rapid simulation calculation method for palladium metal catalytic performance change
Garcia-Lekue et al. Substrate-induced stabilization and reconstruction of zigzag edges in graphene nanoislands on Ni (111)
Ji et al. A multiscale insight into the growth of h-BN: effect of the enclosure
Jiang et al. Topology-induced thermal rectification in carbon nanodevice
Pundlik et al. First-principles investigation of the atomic and electronic structure and magnetic moments in gold nanoclusters
Annevelink et al. Grain boundary structure and migration in graphene via the displacement shift complete lattice
Piñeros et al. Design of two-dimensional particle assemblies using isotropic pair interactions with an attractive well
CN107633150B (en) Simulation calculation method for regulating and controlling catalytic performance of gold atom by using ferroelectric
Opalka et al. Relativistic Jahn-Teller effects in the photoelectron spectra of tetrahedral P4, As4, Sb4, and Bi4
Gomez-Gualdron et al. Characterization of metal nanocatalyst state and morphology during simulated single-walled carbon nanotube growth
Sutrakar et al. Universal stability and temperature dependent phase transformation in group VIIIB–IB transition metal FCC nanowires
Hu et al. Stochastic modeling and simultaneous regulation of surface roughness and porosity in thin film deposition
Zhang et al. Surface phonon localization and heat flux regulation in nanophononic metamaterials
Balcerzak et al. Hubbard pair cluster with elastic interactions. Studies of thermal expansion, magnetostriction and electrostriction
Sepehri et al. Homogeneous and heterogeneous nucleation of water vapor: a comparison using molecular dynamics simulation
Li et al. First-principles study on the structural and electronic properties of clean and atomic oxygen adsorbed ZrC (0 0 1) surface
Zhang et al. Effect of quenching temperature and size on atom movement and local structural change for small copper clusters containing 51–54 atoms during quenching processes
Lu et al. Phase-field crystal study of topological evolution during grain annihilation process in polycrystalline materials

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
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20200728