CN112164425A - Simulation method of silicon dioxide coated calcium fluoride based on Materials Studio - Google Patents
Simulation method of silicon dioxide coated calcium fluoride based on Materials Studio Download PDFInfo
- Publication number
- CN112164425A CN112164425A CN202011031502.2A CN202011031502A CN112164425A CN 112164425 A CN112164425 A CN 112164425A CN 202011031502 A CN202011031502 A CN 202011031502A CN 112164425 A CN112164425 A CN 112164425A
- Authority
- CN
- China
- Prior art keywords
- model
- calcium fluoride
- materials studio
- silicon dioxide
- force field
- 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.)
- Granted
Links
Images
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C10/00—Computational theoretical chemistry, i.e. ICT specially adapted for theoretical aspects of quantum chemistry, molecular mechanics, molecular dynamics or the like
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C20/00—Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
- G16C20/10—Analysis or design of chemical reactions, syntheses or processes
Landscapes
- Theoretical Computer Science (AREA)
- Computing Systems (AREA)
- Engineering & Computer Science (AREA)
- Chemical & Material Sciences (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Health & Medical Sciences (AREA)
- Physics & Mathematics (AREA)
- Analytical Chemistry (AREA)
- Chemical Kinetics & Catalysis (AREA)
- Crystallography & Structural Chemistry (AREA)
- Silicon Compounds (AREA)
Abstract
The invention particularly relates to a method for simulating silicon dioxide coated calcium fluoride by using Materials Studio software, which comprises the following steps: 1) constructing an initial model to be calculated by using software to obtain a silicon dioxide-calcium fluoride molecular structure data file; 2) performing structural optimization on the initial model by using a Forcite module in software; 3) performing molecular dynamics operation on the model by using a Forcite module in software; 4) performing kinetic analysis on the model subjected to the kinetic treatment; 5) calculating the interface distance and energy of the silicon dioxide-calcium fluoride molecular structure to obtain distance and energy parameters; 6) the adsorption energy between the silicon dioxide and the calcium fluoride is obtained by calculation. According to the method, the silicon dioxide coated calcium fluoride with complex procedures in the experimental process is simulated by a computer, and detailed analysis is carried out from an atomic level, so that an experimenter is helped to analyze the simulated coating result, the experimental period is effectively shortened, the experimental cost is saved, a large number of experiments are avoided, and the safety of the experimenter is ensured.
Description
Technical Field
The invention relates to the technical field of computer simulation, in particular to a simulation method of silicon dioxide coated calcium fluoride based on Materials Studio.
Background
The information in this background section is only for enhancement of understanding of the general background of the invention and is not necessarily to be construed as an admission or any form of suggestion that this information forms the prior art that is already known to a person of ordinary skill in the art.
Compared with the traditional ceramic material, the strength and the fracture of the nano composite ceramic material are greatly improved, and the wear resistance, the hardness and the high-temperature performance of the nano composite ceramic material are greatly improved. The ceramic material has a series of excellent special properties such as high temperature resistance, abrasion resistance, corrosion resistance and the like, but the brittleness of the ceramic material limits the application.
Calcium fluoride is a good solid lubricant. The addition of the solid lubricant can ensure that the cutter material has lubricating and antifriction wear-resisting properties, and the use of lubricating oil is avoided. However, the poor mechanical properties of the calcium fluoride solid lubricant lead to the reduced reliability of the self-lubricating ceramic material. The development of the surface modification technology enables the solid lubricant to have good lubricating and antifriction effects and good mechanical properties. The solid lubricant of calcium fluoride coated by silicon dioxide has both antifriction effect and toughening effect.
However, the prior art has shortcomings in studying the mechanism of the nanocomposite ceramic material. The coating only stays in the experimental stage, and the research on the molecular level is still poor. The computer simulation technology is utilized to carry out multi-level and multi-scale coupling simulation on the microstructure and the mechanical property of the nano composite ceramic material, the relation between the organization structure and the mechanical property of the material is known from multiple levels, so that the understanding on the evolution mechanism of the microstructure can be deepened, the development of the existing theory and the proposal of a new theory can be promoted, the microstructure of the material can be obviously improved, the preparation process can be optimized, the mechanical property of the nano composite ceramic material can be improved, the experimental resources can be saved, the research and development cost can be reduced, and the like.
Disclosure of Invention
Aiming at the defects of the prior art, the invention aims to provide a simulation method of silicon dioxide coated calcium fluoride based on Materials Studio, which simulates a coating process through a computer and explains the coating process from a mechanism aspect.
Specifically, the technical scheme of the invention is as follows:
in a first aspect of the present invention, the present invention provides a simulation method of a Materials Studio-based silica coated calcium fluoride, comprising the steps of:
s1: in Materials Studio, build SiO2And CaF2And cutting the surface of the material with a tool;
s2: in Materials Studio, the SiO created is2And CaF2The crystal model establishes lamellar connection to form three-layer conformation of upper, middle and lower parts;
s3: selecting a force field for the model constructed at S2 in Materials Studio;
s4: in Materials Studio, carrying out structural optimization on the model constructed by S2;
s5: in Materials Studio, performing molecular dynamics calculation on the model after structure optimization in S4;
s6: in Materials Studio, the model after molecular dynamics calculation in S5 was analyzed in terms of conformation, adsorption energy, radial distribution function.
Further, in S1, by cutting the crystal face, SiO is obtained2And CaF2Crystal model to obtain SiO with different surface properties2And CaF2(ii) a The crystal model is enlarged, and the size of the crystal model is not less than twice of the truncation radius.
Further, in S2, since the end atomic species are different, a plurality of different interface models can be created, and the SiO after cutting can be modeled2And CaF2The crystal model establishes layered connection to form a multi-layer molecular structure model from top to bottomThe layers react with each other.
Further, in S3, SiO is selected to cover the inorganic crystal model2And CaF2The type of force field of (d); preferably, the type of the force field is a COMPASS force field, the COMPASS force field can be applied to an inorganic crystal, and SiO2And CaF2All in the range covered by the force field, so that each atom and bond can be reasonably distributed with reasonable data.
Further, in S4, when performing structural optimization on the model constructed in S2, the algorithm adopted is selected from: smart, Steepest device, concrete gradient, Quasi-Newton, ABNR; the convergence accuracy is selected from: coarse, Medium, Fine and Ultra-Fine; the maximum iteration times can be set by self;
the force field distribution method is selected from automatic distribution and manual distribution; the Charge distribution method is selected from the group consisting of Use current, Charge using QEq, Charge using gater, and Forcefield associated; the superposition method is divided into EWald, Atom based, Group based and PPPM.
The static term is an EWald summation method, the precision is 0.001, the van der Waals term adopts an Atom-based summation method, a local server is selected, the number of parallel operation is selected to be 8, and after parameters are set, the structure of the model is optimized.
Preferably, the selection algorithm is Smart, the convergence precision is Fine, and the maximum iteration number is 500, so that the optimized model can meet the convergence requirement, and the model is prevented from not converging after reaching the maximum iteration number;
preferably, the force field allocation method is selected to be automatic allocation and the charge allocation method is Forcefield assigned.
Further, in S5, the parameters set by the molecular dynamics calculation are as follows:
the ensemble is selected from NVE, NVT, NPT and NPH, preferably NVT and NVE;
the initial speed allocation is selected from Current, Random, preferably Random;
the pressure is 0, the temperature is 298k, the time step is 1fs, the simulation time is 100ps, and the energy deviation is 500000000;
the temperature control method is selected from the group consisting of Velocity Scale, Nose, Andersen, Berendsen and NHL, and the temperature control method is preferably Nose;
selecting a force field distribution method as automatic force field distribution, a charge distribution method as Forcefield assigned, an electrostatic term as an EWald summation method with the precision of 0.001, and an Atom-based summation method for van der Waals; selecting a local server, wherein the number of the local servers running in parallel is 8; and after the parameters are set, performing molecular dynamics calculation on the model after structure optimization.
Further, in S6, it is observed whether the model after kinetic run meets the criteria for expected coating, and if the criteria for coating are met, molecular dynamics analysis is performed.
Wherein, the standard judgment basis of the coating is as follows: if the distance between the final model layers is not more than 4A, the binding energy is negative, and the conformation conforms to the concept of coating, the model meets the standard of coating.
The specific embodiment of the invention has the following beneficial effects:
the process of coating calcium fluoride by silicon dioxide is simulated by a computer, so that an experimenter is helped to analyze the coated result, the experimental resources are effectively saved, and the research and development cost is reduced;
the process of coating is more clear and clear for a user by displaying the crystal surface, the coating environment and the simulation result.
Drawings
The accompanying drawings, which are incorporated in and constitute a part of this specification, are included to provide a further understanding of the invention, and are incorporated in and constitute a part of this specification, illustrate exemplary embodiments of the invention and together with the description serve to explain the invention and not to limit the invention.
FIG. 1 is a block flow diagram of a method for investigating a silica-coated calcium fluoride according to the present invention;
FIG. 2 is a Model of the initial molecular structure of Model 1 of silica-calcium fluoride according to example 1 of the present invention;
FIG. 3 is a Model of molecular structure at the interface after the Model 1 structure optimization treatment of silica-calcium fluoride in example 1 of the present invention;
FIG. 4 is a Model of the molecular structure at the interface after Model 1 kinetic simulation of silica-calcium fluoride in example 1 of the present invention.
Detailed Description
The invention will be further illustrated with reference to the following specific examples. It should be understood that these examples are for illustrative purposes only and are not intended to limit the scope of the present invention. Unless defined otherwise, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art. Any methods and materials similar or equivalent to those described herein can be used in the methods of the present invention. The preferred embodiments and materials described herein are intended to be exemplary only.
Example 1
As shown in fig. 1, a simulation method of silica-coated calcium fluoride based on Materials Studio, which takes silica and calcium fluoride as research objects in this embodiment, specifically includes the following steps:
(1) building an initial model of silicon dioxide and calcium fluoride;
building a silicon dioxide and calcium fluoride molecular structure model by using a Visualizer interface in Materials Studio software, and cutting a crystal face of the built model, wherein a layer1 and a layer 3 of the model are silicon dioxide, and two terminal surfaces can be obtained after cutting: o terminal and Si-O, taking the O terminal surface as layer1 and layer 3 of Model 1, taking the Si-O terminal surface as layer1 and layer 3 of Model 2, taking the upper surface atoms of layer 2 as F, Ca mixed surface, wherein the numbers of F and Ca are both 100, the lower surface atoms are F, and the number is 100; and (3) establishing layered connection by using a Build layer tool in the Build, and arranging a vacuum layer between interfaces to obtain a silicon dioxide-calcium fluoride molecular structure data file to be calculated. In the initial structure Model, Model 1: a, b, 41A, c, α, β, 90, γ, 120 °; model 2: a, b, 41A, c, α, β, 90, γ, 120. Each of the initial structural models contained Ca 900, F1800, Si 768, and O1536, and the entire model showed electric neutrality.
(2) Carrying out structural optimization on the molecular structure of the silicon dioxide-calcium fluoride;
setting parameters for structure optimization, adopting a Geometry optimization task in a Forcit module when optimizing a silicon dioxide-calcium fluoride molecular structure model, adopting a Smart algorithm, wherein Energy is 0.001, Force is 0.5, Displacement is 0.015A, the number of calculation steps is 500, a COMPASS Force field is adopted in the Force field, charges are Forcefield assembled, charge groups are automatically distributed in the Force field, an electrostatic term is an EWALD (equal-sum-of-weight) method is adopted in the accuracy of 0.001, an Atom-based sum-of-weight method is adopted in a van der Waals term, a truncation method adopts a Cubic spline, the truncation radius is 12.5A, the bond width is 1A, and Long range correction is adopted in the calculation process. And after the calculation is finished, obtaining a silicon dioxide-calcium fluoride molecular structure data file after energy minimization treatment.
(3) Performing kinetic calculation on the molecular structure of the adsorbent-adsorbent subjected to optimization treatment;
setting parameters of kinetic calculation, adopting a Dynamics task in a Forcite module when optimizing a silicon dioxide-calcium fluoride molecular structure model, firstly adopting a regular ensemble, setting the temperature to be 297K, adopting a Nose temperature control method, setting the time step to be 1fs, the time duration to be 100ps, adopting a COMPASSASS force field for a force field, automatically distributing charges for a charge group for the force field, adopting an EWALD (electro-static summation) method for an Ewald summation) and the precision to be 0.001, adopting an Atom-based summation method for a Van der Waals, adopting a Cubic spline for a truncation method, setting the truncation radius to be 12.5A and the bond width to be 1A, and adopting a Long range correction in the calculation process; then adopting a micro-canonical ensemble, setting the temperature to be 297K, adopting a Nose temperature control method, setting the time step to be 1fs, setting the time duration to be 100ps, adopting a COMPASS force field for the force field, setting the charges to be Forcefield assigned, automatically distributing the charge groups to the force field, adopting an EWald summation method for the static term, setting the precision to be 0.001, adopting an Atom-based summation method for the Van der Waals term, adopting a Cubic spline for the truncation method, setting the truncation radius to be 12.5A and setting the bond width to be 1A, and adopting a Long range iteration in the calculation process; and finally, adopting a regular ensemble, setting the temperature to be 297K, adopting a Nose temperature control method, setting the time step to be 1fs, setting the time duration to be 1000ps, adopting a COMPASS force field for a force field, setting the charges to be Forcefield assigned, automatically distributing the charge groups to the force field, adopting an EWALD (equal-wavelet-based summation) for the static term, setting the precision to be 0.001, adopting an Atom-based summation for the Van der Waals term, adopting a Cubic spline for a truncation method, setting the truncation radius to be 12.5A and setting the bond width to be 1A, and adopting Long range rectification in the calculation process.
(4) Performing kinetic analysis on the model subjected to the kinetic treatment;
and extracting each index parameter of dynamics by using Analysis in a Forcite module in Materials Studio software to obtain angular distribution, angular evolution, concentration distribution, density and density field and state parameter information of each output frame.
(5) Calculating the interface distance and energy of the silicon dioxide-calcium fluoride molecular structure to obtain distance and energy parameters;
the Distance between the interfaces was measured using the Distance of SD viewer in Materials Studio software.
Calculating Energy by utilizing an Energy task in a Forcite module in Materials Studio software, firstly calculating the Energy of a silicon dioxide-calcium fluoride integral system, outputting a data file containing the Energy, deleting the calcium fluoride in the silicon dioxide-calcium fluoride system obtained in the step (3), calculating the Energy after removing the calcium fluoride, outputting the data file containing the Energy, finally deleting the silicon dioxide in the silicon dioxide-calcium fluoride system obtained in the step (3), calculating the Energy after removing the silicon dioxide, and outputting the data file containing the Energy.
(6) The adsorption energy between the silicon dioxide and the calcium fluoride is obtained by calculation.
TABLE 1 energy and adsorption energy of silica-calcium fluoride molecular structural system
By measuring the interlayer Distance between interfaces by using the Distance of SD viewer in Materials Studio software, the interface Distance on Model 1 is less than 0A, the migration phenomenon of each atom in the system is obvious, and the surface atomic structure Model is fully relaxed. The Ca atom activity of the surface layer is far greater than that of the F atom, so that the surface layer is more likely to react with the outside. In the upper layer of the layer 2, a fault phenomenon also occurs, only a small amount of Ca atoms with positive charges move to the surface of negative charges of the layer1 in the lower layer, and the atom migration phenomenon of the upper layer of the layer 2 is more violent than that of the lower layer, which indicates that the upper layer is more likely to act; model 2 Upper and lower boundsThe surface distance is 2A, the atom position of the interface layer vertical to the interface direction is changed, and the longitudinal relaxation phenomenon of atoms is generated. CaF2Overall stability, and SiO2A large degree of positional change occurs. SiO 22The whole shrinks towards the interface. The atomic layer closer to the interface has a larger relaxation amount, and the atomic layer farther from the interface has a smaller relaxation amount. In Model 1, the upper surface of layer 2 is a positive charge surface, the lower surface of layer 3 is a negative charge surface, the interaction between the interface atoms promotes the layers 2 and 3 to approach each other, and on the contrary, the layers 1 and 2 are far away from each other.
In Model 2, the charge theory does not work. This is probably because the interfacial oxygen atom in Model 1 is in an unbound, dangling state, has high activity and significant charge effect, and is liable to react with CaF2An effect occurs. The interface in Model 2 is in an O-Si bond and state, and the interface is relatively stable.
And as shown in table 1, the adsorption energy is negative, indicating that the adsorption is an exothermic process, and the adsorption is stable. Furthermore, the adsorption energy of Model 1 was lower than that of Model 2, indicating that Model 1 was more stable. Therefore, adsorption phenomenon occurred in both Model 1 and Model 2.
Example 2
As shown in fig. 1, a simulation method of silica-coated calcium fluoride based on Materials Studio, which takes silica and calcium fluoride as research objects in this embodiment, specifically includes the following steps:
(1) building an initial model of silicon dioxide and calcium fluoride;
building a silicon dioxide and calcium fluoride molecular structure model by using a Visualizer interface in Materials Studio software, and cutting a crystal face of the built model, wherein a layer1 and a layer 3 of the model are silicon dioxide, and two terminal surfaces can be obtained after cutting: o terminal and Si-O, the O terminal surface is taken as layer1 and layer 3 of Model 3, the Si-O terminal surface is taken as layer1 and layer 3 of Model 4, the upper surface atoms of layer 2 are F, the number of which is 200, the lower surface atoms are Ca, and the number of which is 100; and (3) establishing layered connection by using a Build layer tool in the Build, and arranging a vacuum layer between interfaces to obtain a silicon dioxide-calcium fluoride molecular structure data file to be calculated. In the initial structure Model, Model 3: a, b, 41A, c, α, β, 90, γ, 120 °; model 4: a, b, 41A, c, α, β, 90, γ, 120. Each of the initial structural models contained Ca 900, F1800, Si 768, and O1536, and the entire model showed electric neutrality.
(2) Carrying out structural optimization on the molecular structure of the silicon dioxide-calcium fluoride;
setting parameters for structure optimization, adopting a Geometry optimization task in a Forcit module when optimizing a silicon dioxide-calcium fluoride molecular structure model, adopting a Smart algorithm, wherein Energy is 0.001, Force is 0.5, Displacement is 0.015A, the number of calculation steps is 500, a COMPASS Force field is adopted in the Force field, charges are Forcefield assembled, charge groups are automatically distributed in the Force field, an electrostatic term is an EWALD (equal-sum-of-weight) method is adopted in the accuracy of 0.001, an Atom-based sum-of-weight method is adopted in a van der Waals term, a truncation method adopts a Cubic spline, the truncation radius is 12.5A, the bond width is 1A, and Long range correction is adopted in the calculation process. And after the calculation is finished, obtaining a silicon dioxide-calcium fluoride molecular structure data file after energy minimization treatment.
(3) Performing kinetic calculation on the molecular structure of the adsorbent-adsorbent subjected to optimization treatment;
setting parameters of kinetic calculation, adopting a Dynamics task in a Forcite module when optimizing a silicon dioxide-calcium fluoride molecular structure model, firstly adopting a regular ensemble, setting the temperature to be 297K, adopting a Nose temperature control method, setting the time step to be 1fs, the time duration to be 100ps, adopting a COMPASSASS force field for a force field, automatically distributing charges for a charge group for the force field, adopting an EWALD (electro-static summation) method for an Ewald summation) and the precision to be 0.001, adopting an Atom-based summation method for a Van der Waals, adopting a Cubic spline for a truncation method, setting the truncation radius to be 12.5A and the bond width to be 1A, and adopting a Long range correction in the calculation process; then adopting a micro-canonical ensemble, setting the temperature to be 297K, adopting a Nose temperature control method, setting the time step to be 1fs, setting the time duration to be 100ps, adopting a COMPASS force field for the force field, setting the charges to be Forcefield assigned, automatically distributing the charge groups to the force field, adopting an EWald summation method for the static term, setting the precision to be 0.001, adopting an Atom-based summation method for the Van der Waals term, adopting a Cubic spline for the truncation method, setting the truncation radius to be 12.5A and setting the bond width to be 1A, and adopting a Long range iteration in the calculation process; and finally, adopting a regular ensemble, setting the temperature to be 297K, adopting a Nose temperature control method, setting the time step to be 1fs, setting the time duration to be 1000ps, adopting a COMPASS force field for a force field, setting the charges to be Forcefield assigned, automatically distributing the charge groups to the force field, adopting an EWALD (equal-wavelet-based summation) for the static term, setting the precision to be 0.001, adopting an Atom-based summation for the Van der Waals term, adopting a Cubic spline for a truncation method, setting the truncation radius to be 12.5A and setting the bond width to be 1A, and adopting Long range rectification in the calculation process.
(4) Performing kinetic analysis on the model subjected to the kinetic treatment;
and extracting each index parameter of dynamics by using Analysis in a Forcite module in Materials Studio software to obtain angular distribution, angular evolution, concentration distribution, density and density field and state parameter information of each output frame.
(5) Calculating the interface distance and energy of the silicon dioxide-calcium fluoride molecular structure to obtain distance and energy parameters;
the Distance between the interfaces was measured using the Distance of SD viewer in Materials Studio software.
Calculating Energy by utilizing an Energy task in a Forcite module in Materials Studio software, firstly calculating the Energy of a silicon dioxide-calcium fluoride integral system, outputting a data file containing the Energy, deleting the calcium fluoride in the silicon dioxide-calcium fluoride system obtained in the step (3), calculating the Energy after removing the calcium fluoride, outputting the data file containing the Energy, finally deleting the silicon dioxide in the silicon dioxide-calcium fluoride system obtained in the step (3), calculating the Energy after removing the silicon dioxide, and outputting the data file containing the Energy.
(6) The adsorption energy between the silicon dioxide and the calcium fluoride is obtained by calculation.
TABLE 2 energy and adsorption energy of silica-calcium fluoride molecular structural system
By measuring the interlayer Distance between interfaces by using the Distance of SD viewer in Materials Studio software, the Distance between interfaces under Model 3 is less than 0A, and the surface atomic structure Model is fully relaxed after geometric structure optimization and dynamics. The migration phenomenon of the lower layer of the layer 2 is obvious, and Ca atoms are easy to react with the upper layer of the layer 1. Only the surface atoms of the upper layer move significantly. Only a few atoms move to the lower layer of layer 3, where there are many positively charged Ca atoms. Layer 3 migrates away from the interface. (ii) a The distance between the upper interface and the lower interface of the Model 4 is 2A, and after geometric structure optimization and dynamics are carried out, the final configuration is similar to that of the Model 2. Layer 2 is overall more stable, and both layers 1, 3 shrink towards the interface. In Model 3, the upper layer surface of layer 2 is a negatively charged surface, the lower layer surface of layer 3 is a negatively charged surface, and the interaction between the interface atoms promotes the separation of layers 2 and 3, and conversely, layers 1 and 2 are close to each other. Model 4, the transition to the interface.
And as shown in table 2, the adsorption energy is negative, indicating that the adsorption is an exothermic process, and the adsorption is stable. And the adsorption energy of Model 3 is lower than that of Model 4, indicating that Model 3 is more stable. Therefore, adsorption phenomenon occurred in both Model 3 and Model 4.
Example 3
As shown in fig. 1, a simulation method of silica-coated calcium fluoride based on Materials Studio, which takes silica and calcium fluoride as research objects in this embodiment, specifically includes the following steps:
(1) building an initial model of silicon dioxide and calcium fluoride;
building a silicon dioxide and calcium fluoride molecular structure model by using a Visualizer interface in Materials Studio software, and cutting a crystal face of the built model, wherein a layer1 and a layer 3 of the model are silicon dioxide, and two terminal surfaces can be obtained after cutting: o terminal and Si-O, taking the O terminal surface as layer1 and layer 3 of Model 1, taking the Si-O terminal surface as layer1 and layer 3 of Model 2, taking the upper surface layer atoms of layer 2 as F, Ca mixed surface, taking the numbers of F and Ca as 100, taking the lower surface layer atoms as F, and taking the numbers as 100; and (3) establishing layered connection by using a Build layer tool in the Build, and arranging a vacuum layer between interfaces to obtain a silicon dioxide-calcium fluoride molecular structure data file to be calculated. In the initial structure Model, Model 5: a, b, 41A, c, α, β, 90, γ, 120 °; model 6: a, b, 41A, c, α, β, 90, γ, 120. Each of the initial structural models contained Ca 900, F1800, Si 768, and O1536, and the entire model showed electric neutrality.
(2) Carrying out structural optimization on the molecular structure of the silicon dioxide-calcium fluoride;
setting parameters for structure optimization, adopting a Geometry optimization task in a Forcit module when optimizing a silicon dioxide-calcium fluoride molecular structure model, adopting a Smart algorithm, wherein Energy is 0.001, Force is 0.5, Displacement is 0.015A, the number of calculation steps is 500, a COMPASS Force field is adopted in the Force field, charges are Forcefield assembled, charge groups are automatically distributed in the Force field, an electrostatic term is an EWALD (equal-sum-of-weight) method is adopted in the accuracy of 0.001, an Atom-based sum-of-weight method is adopted in a van der Waals term, a truncation method adopts a Cubic spline, the truncation radius is 12.5A, the bond width is 1A, and Long range correction is adopted in the calculation process. And after the calculation is finished, obtaining a silicon dioxide-calcium fluoride molecular structure data file after energy minimization treatment.
(3) Performing kinetic calculation on the molecular structure of the adsorbent-adsorbent subjected to optimization treatment;
setting parameters of kinetic calculation, adopting a Dynamics task in a Forcite module when optimizing a silicon dioxide-calcium fluoride molecular structure model, firstly adopting a regular ensemble, setting the temperature to be 297K, adopting a Nose temperature control method, setting the time step to be 1fs, the time duration to be 100ps, adopting a COMPASSASS force field for a force field, automatically distributing charges for a charge group for the force field, adopting an EWALD (electro-static summation) method for an Ewald summation) and the precision to be 0.001, adopting an Atom-based summation method for a Van der Waals, adopting a Cubic spline for a truncation method, setting the truncation radius to be 12.5A and the bond width to be 1A, and adopting a Long range correction in the calculation process; then adopting a micro-canonical ensemble, setting the temperature to be 297K, adopting a Nose temperature control method, setting the time step to be 1fs, setting the time duration to be 100ps, adopting a COMPASS force field for the force field, setting the charges to be Forcefield assigned, automatically distributing the charge groups to the force field, adopting an EWald summation method for the static term, setting the precision to be 0.001, adopting an Atom-based summation method for the Van der Waals term, adopting a Cubic spline for the truncation method, setting the truncation radius to be 12.5A and setting the bond width to be 1A, and adopting a Long range iteration in the calculation process; and finally, adopting a regular ensemble, setting the temperature to be 297K, adopting a Nose temperature control method, setting the time step to be 1fs, setting the time duration to be 1000ps, adopting a COMPASS force field for a force field, setting the charges to be Forcefield assigned, automatically distributing the charge groups to the force field, adopting an EWALD (equal-wavelet-based summation) for the static term, setting the precision to be 0.001, adopting an Atom-based summation for the Van der Waals term, adopting a Cubic spline for a truncation method, setting the truncation radius to be 12.5A and setting the bond width to be 1A, and adopting Long range rectification in the calculation process.
(4) Performing kinetic analysis on the model subjected to the kinetic treatment;
and extracting each index parameter of dynamics by using Analysis in a Forcite module in Materials Studio software to obtain angular distribution, angular evolution, concentration distribution, density and density field and state parameter information of each output frame.
(5) Calculating the interface distance and energy of the silicon dioxide-calcium fluoride molecular structure to obtain distance and energy parameters;
the Distance between the interfaces was measured using the Distance of SD viewer in Materials Studio software.
Calculating Energy by utilizing an Energy task in a Forcite module in Materials Studio software, firstly calculating the Energy of a silicon dioxide-calcium fluoride integral system, outputting a data file containing the Energy, deleting the calcium fluoride in the silicon dioxide-calcium fluoride system obtained in the step (3), calculating the Energy after removing the calcium fluoride, outputting the data file containing the Energy, finally deleting the silicon dioxide in the silicon dioxide-calcium fluoride system obtained in the step (3), calculating the Energy after removing the silicon dioxide, and outputting the data file containing the Energy.
(6) The adsorption energy between the silicon dioxide and the calcium fluoride is obtained by calculation.
TABLE 3 energy and adsorption energy of silica-calcium fluoride molecular structural system
By measuring the interlayer Distance between interfaces by using the Distance of SD viewer in Materials Studio software, the interface Distance on the Model 5 is less than 0A, the migration phenomenon of each atom in the system is obvious, and the surface atomic structure Model is fully relaxed. The lower layer only has a small amount of Ca atoms with positive charges moving to the surface of the negative charges of the layer1, and the phenomenon of atom migration of the upper layer of the layer 2 is more violent than that of the lower layer, which shows that the upper layer is more likely to act; the distance between the upper interface and the lower interface of the Model 6 is 2A, and the atomic position of the interface layer vertical to the interface direction is changed, so that the longitudinal relaxation phenomenon of atoms is generated. CaF2Overall stability, and SiO2A large degree of positional change occurs. SiO 22The whole shrinks towards the interface. The atomic layer closer to the interface has a larger relaxation amount, and the atomic layer farther from the interface has a smaller relaxation amount.
And as shown in table 3, the adsorption energy is negative, indicating that the adsorption is an exothermic process, and the adsorption is stable. Further, the adsorption energy of Model 5 was lower than that of Model 6, indicating that Model 5 was more stable. Therefore, adsorption phenomenon occurred in both Model 5 and Model 6.
The above description is only a preferred embodiment of the present invention and is not intended to limit the present invention, and various modifications and changes may be made by those skilled in the art. Any modification, equivalent replacement, or improvement made within the spirit and principle of the present invention should be included in the protection scope of the present invention.
Claims (10)
1. A simulation method of silicon dioxide coated calcium fluoride based on Materials Studio is characterized by comprising the following steps:
s1: in Materials Studio, build SiO2And CaF2And cutting the surface of the material with a tool;
s2: in Materials Studio, the SiO created is2And CaF2The crystal model establishes lamellar connection to form three-layer conformation of upper, middle and lower parts;
s3: selecting a force field for the model constructed at S2 in Materials Studio;
s4: in Materials Studio, carrying out structural optimization on the model constructed by S2;
s5: in Materials Studio, performing molecular dynamics calculation on the model after structure optimization in S4;
s6: in Materials Studio, the model after molecular dynamics calculation in S5 was analyzed in terms of conformation, adsorption energy, radial distribution function.
2. The simulation method of Materials Studio-based silica coated calcium fluoride as set forth in claim 1, wherein SiO is obtained by cutting crystal planes in S12And CaF2Crystal model to obtain SiO with different surface properties2And CaF2(ii) a The crystal model is enlarged, and the size of the crystal model is not less than twice of the truncation radius.
3. The simulation method of Materials Studio-based silica-coated calcium fluoride as set forth in claim 1, wherein in S2, since the terminal atoms are different in kind, a plurality of different interfaces can be modeled, and the cut SiO is2And CaF2The crystal model establishes layered connection to form a multi-layer molecular structure model of upper, middle and lower parts for interlayer mutual reaction.
4. The simulation method of Materials Studio based silica coated calcium fluoride as set forth in claim 1, wherein in S3, SiO covering the inorganic crystal model is selected2And CaF2The type of force field of (d); preferably, the force field type is a COMPASS force field.
5. The Materials Studio-based simulation method of silica-coated calcium fluoride as set forth in claim 1, wherein in S4, a suitable algorithm, convergence accuracy, maximum number of iterations are selected; selecting a proper force field distribution method, a proper charge distribution method and a proper superposition method; selecting the number of servers and running in parallel, and carrying out structural optimization on the model after setting parameters;
wherein the algorithm is selected from: smart, Steepest device, concrete gradient, Quasi-Newton, ABNR; the convergence accuracy is selected from: coarse, Medium, Fine and Ultra-Fine; the maximum number of iterations is 500;
the force field distribution method is selected from automatic distribution and manual distribution; the Charge distribution method is selected from the group consisting of Use current, Charge using QEq, Charge using gater, and Forcefield associated; the superposition method comprises the steps of EWald, Atom based, Group based and PPPM;
the static term is an EWald summation method, the precision is 0.001, the van der Waals term adopts an Atom-based summation method, and the number of the selected local servers and the selected parallel operation is 8.
6. The Materials Studio-based simulation method of silica-coated calcium fluoride as set forth in claim 5, wherein the algorithm is selected to be Smart, and the convergence accuracy is selected to be Fine; the force field distribution method is selected to be automatic distribution, and the charge distribution method is Forcefield assigned.
7. The Materials Studio-based simulation of silica-coated calcium fluoride of claim 1, wherein in S5, the appropriate ensemble, initial velocity distribution, pressure, temperature, time step, simulation time, energy bias are selected; selecting a proper temperature control method and a proper pressure control method; selecting a proper force field distribution method, a proper charge distribution method and a proper superposition method; and selecting the number of the servers and the parallel operation, setting parameters, and then carrying out molecular dynamics calculation on the model with the optimized structure.
8. The simulation method of Materials Studio based silica coated calcium fluoride according to claim 1,
the ensemble is selected from NVE, NVT, NPT and NPH, preferably NVT and NVE;
the initial speed allocation is selected from Current, Random, preferably Random;
the pressure is 0, the temperature is 298k, the time step is 1fs, the simulation time is 100ps, and the energy deviation is 500000000;
the temperature control method is selected from the group consisting of Velocity Scale, Nose, Andersen, Berendsen and NHL, and the temperature control method is preferably Nose;
selecting a force field distribution method as automatic force field distribution, a charge distribution method as Forcefield assigned, an electrostatic term as an EWald summation method with the precision of 0.001, and an Atom-based summation method for van der Waals; the number of the local servers selected and operated in parallel is 8.
9. The Materials Studio-based simulation method of calcium fluoride coating with silica according to claim 1, wherein in S6, it is observed whether the model after kinetic run meets the criteria for expected coating, and if so, molecular dynamics analysis is performed.
10. The material Studio-based simulation method of calcium fluoride coating with silica according to claim 6, wherein the standard criteria for coating are: if the distance between the final model layers is not more than 4A, the binding energy is negative, and the conformation conforms to the concept of coating, the model meets the standard of coating.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011031502.2A CN112164425B (en) | 2020-09-27 | 2020-09-27 | Simulation method of silicon dioxide coated calcium fluoride based on Materials Studio |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011031502.2A CN112164425B (en) | 2020-09-27 | 2020-09-27 | Simulation method of silicon dioxide coated calcium fluoride based on Materials Studio |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112164425A true CN112164425A (en) | 2021-01-01 |
| CN112164425B CN112164425B (en) | 2023-03-28 |
Family
ID=73864142
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011031502.2A Active CN112164425B (en) | 2020-09-27 | 2020-09-27 | Simulation method of silicon dioxide coated calcium fluoride based on Materials Studio |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112164425B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113515854A (en) * | 2021-06-10 | 2021-10-19 | 太原理工大学 | Molecular dynamics method for evaluating charged activated water mist atomization dust fall performance |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20100023473A1 (en) * | 2006-12-11 | 2010-01-28 | Neumann Marcus A | Tailor-made force fields for crystal structure prediction |
| CN102446235A (en) * | 2010-10-11 | 2012-05-09 | 中国石油化工股份有限公司 | Method for simulating and calculating interaction parameters among chemical components by using computer |
| CN106096258A (en) * | 2016-06-07 | 2016-11-09 | 上海理工大学 | A kind of analysis method of fluorochemical urethane surface hydrophobic |
| CN108846253A (en) * | 2018-06-21 | 2018-11-20 | 上海理工大学 | A method of utilizing computer simulation building graphene oxide and polycarbodiimide cross-linked structure model |
| CN110436942A (en) * | 2019-08-28 | 2019-11-12 | 齐鲁工业大学 | The preparation method of silicon dioxide coated nano sheet calcirm-fluoride composite granule |
| CN111540413A (en) * | 2020-04-16 | 2020-08-14 | 南京航空航天大学 | Novel method for predicting interface thermal resistance of silicon dioxide/epoxy resin composite material |
-
2020
- 2020-09-27 CN CN202011031502.2A patent/CN112164425B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20100023473A1 (en) * | 2006-12-11 | 2010-01-28 | Neumann Marcus A | Tailor-made force fields for crystal structure prediction |
| CN102446235A (en) * | 2010-10-11 | 2012-05-09 | 中国石油化工股份有限公司 | Method for simulating and calculating interaction parameters among chemical components by using computer |
| CN106096258A (en) * | 2016-06-07 | 2016-11-09 | 上海理工大学 | A kind of analysis method of fluorochemical urethane surface hydrophobic |
| CN108846253A (en) * | 2018-06-21 | 2018-11-20 | 上海理工大学 | A method of utilizing computer simulation building graphene oxide and polycarbodiimide cross-linked structure model |
| CN110436942A (en) * | 2019-08-28 | 2019-11-12 | 齐鲁工业大学 | The preparation method of silicon dioxide coated nano sheet calcirm-fluoride composite granule |
| CN111540413A (en) * | 2020-04-16 | 2020-08-14 | 南京航空航天大学 | Novel method for predicting interface thermal resistance of silicon dioxide/epoxy resin composite material |
Non-Patent Citations (3)
| Title |
|---|
| ZHAOQIANG CHEN 等: "Synthesis and Simulation of CaF2@Al(OH)3 Core-Shell Coated Solid Lubricant Composite Powder", 《MDPI》 * |
| 衣明东 等: "纳米CaF２对自润滑陶瓷材料力学性能的影响", 《硅酸盐学报》 * |
| 陈辉 等: "纳米SiO2包覆片状微米h-BN复合粉体的制备与表征", 《人工晶体学报》 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113515854A (en) * | 2021-06-10 | 2021-10-19 | 太原理工大学 | Molecular dynamics method for evaluating charged activated water mist atomization dust fall performance |
| CN113515854B (en) * | 2021-06-10 | 2022-09-27 | 太原理工大学 | Molecular dynamics method for evaluating charged activated water mist atomization dust-settling performance |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112164425B (en) | 2023-03-28 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Li et al. | Digital twin in aerospace industry: A gentle introduction | |
| Lemos | Discrete element modeling of masonry structures | |
| Balaguer et al. | Robotics and automation in construction | |
| Zhong et al. | Geographic information system-based visual simulation methodology and its application in concrete dam construction processes | |
| CN112164425B (en) | Simulation method of silicon dioxide coated calcium fluoride based on Materials Studio | |
| CN103097657A (en) | Machine, computer program product and method to generate unstructured grids and carry out parallel reservoir simulation | |
| CN102141084A (en) | Numerical simulation method for temperature and thickness relation of static thrust bearing gap oil film based on film thickness variation | |
| CN111858821B (en) | Processing method and device for space analysis task, computer equipment and storage medium | |
| Chen | Computational geomechanics and hydraulic structures | |
| Gong et al. | Modelling rock failure with a novel continuous to discontinuous method | |
| CN107909287A (en) | A kind of method of the infrastructure management based on BIM | |
| Douglas et al. | BIM, digital twin and cyber-physical systems: crossing and blurring boundaries | |
| Du et al. | Construction simulation for a core rockfill dam based on optimal construction stages and zones: case study | |
| Zheng et al. | A robust potential-based contact force solution approach for discontinuous deformation analysis of irregular convex polygonal block/particle systems | |
| CN113033060B (en) | Optimization method for predicting complex coal seam mining structure | |
| Wu et al. | An application of space syntax to critical working space analysis: The case of building construction | |
| Zhao et al. | Construction and application of a high precision 3D simulation model for geomechanics of the complex coal seam | |
| Luque-Sala et al. | A Virtual reconstruction of Gaudi’s skyscraper hotel attraction using physics-based simulation | |
| Karhu | Product model based design of precast facades | |
| CN105447228A (en) | Earth disturbing gravity assignment model applicability evaluating method | |
| Lim et al. | Four cases of production-installation simulation for free-form concrete panels | |
| Liu | Research on stability of hydraulic system based on nonlinear PID control | |
| Budkowski | Multidimensional cadastre as an element of participation in modern space management | |
| Zhang et al. | Dynamic Modeling and CAE Cosimulation Method for Heavy‐Duty Concrete Spreader | |
| Li et al. | Addictive Manufacturing, Prosumption and Construction Safety |
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 |