CN110600084A - Dissipative particle dynamics method for researching interfacial polymerization reaction mechanism of composite film - Google Patents
Dissipative particle dynamics method for researching interfacial polymerization reaction mechanism of composite film Download PDFInfo
- Publication number
- CN110600084A CN110600084A CN201910858591.9A CN201910858591A CN110600084A CN 110600084 A CN110600084 A CN 110600084A CN 201910858591 A CN201910858591 A CN 201910858591A CN 110600084 A CN110600084 A CN 110600084A
- Authority
- CN
- China
- Prior art keywords
- interfacial polymerization
- dpd
- water
- chloride
- polymerization reaction
- 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
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
Landscapes
- Engineering & Computer Science (AREA)
- Computing Systems (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention belongs to the field of high-performance membrane materials, and particularly relates to a dissipative particle dynamics method for researching a composite membrane interfacial polymerization reaction mechanism by utilizing a computer simulation technology. The method comprises the following steps: (1) determining the composition of a solvent and a monomer in a water phase and an oil phase of an interfacial polymerization system; (2) respectively constructing a DPD model of a solvent and a water-soluble monomer in a water phase and a DPD model of an organic solvent and an oil-soluble monomer in an oil phase; (3) establishing an interface polymerization reaction system DPD model consisting of a water phase and an oil phase; (4) calculating interaction parameters among all DPD beads, namely conservative force parameters; (5) carrying out DPD simulation by using Materials Studio software, and obtaining a motion trail file and a related calculation file of each DPD bead after system balance; (6) and observing the structural characteristics of the separation layer of the composite membrane generated by the interfacial polymerization reaction according to the simulation result, and researching influence factors determining the performance of the separation layer in the interfacial polymerization process by combining with a calculation file. The invention provides a theoretical basis for systematically researching the main control factors of the interfacial polymerization reaction of the separation layer of the composite membrane, and has important guiding significance for preparing the high-performance separation membrane by the experimental process for improving the interfacial polymerization reaction.
Description
Technical Field
The invention belongs to the field of high-performance membrane materials, and particularly relates to a dissipative particle dynamics method for researching a composite membrane interfacial polymerization reaction mechanism by utilizing a computer simulation technology.
Background
With the occurrence of the problems of water resource shortage, increasingly serious water pollution and the like, the membrane separation technology is used as one of economic and efficient technologies for sewage treatment, seawater desalination and brackish water desalination, and has wide market application prospect. The membrane material, as the core of the membrane separation technology, will directly affect the membrane separation performance and the application of the membrane technology. The preparation of high-performance membrane materials is a hot spot of continuous development and research in the industry and academia. Currently commercialized reverse osmosis membranes, nanofiltration membranes and organic solvent resistant nanofiltration composite membranes are generally prepared by interfacial polymerization of a water-soluble monomer in a water phase and a polybasic acid chloride monomer in an oil phase on the surface of a base membrane to form a polyamide selective separation layer. In the interfacial polymerization process, the monomer concentration, the reaction time and the structure of the base membrane are key factors influencing the performance of the finally prepared polyamide composite membrane. Although great progress has been made in the synthesis of novel membrane materials and in the modification of membrane materials, there is still insufficient research and explanation on the microstructure properties and mechanism of interfacial polymerization during the preparation of high-performance membrane materials. Therefore, the method is very important for researching the microstructure characteristics and mechanism in the process of preparing the high-performance composite membrane by the interfacial polymerization reaction.
At present, widely applied experimental characterization methods and detection methods (such as SEM, TEM, AFM and the like) are difficult to meet the requirements of quantitatively analyzing the surface microscopic characteristics and the dynamic change process of the water-soluble monomer in the interface polymerization process at the atomic molecular level, and are also difficult to explain the interface polymerization mechanism at the molecular atomic level. Therefore, it remains a great challenge to experimentally precisely control the structure and performance of a polymeric composite membrane. Compared with experimental means, computer simulation can help researchers understand possible influences of interfacial polymerization reaction processes, clear monomer/solvent concentrations, structures, chemical properties and the like on polymerization reactions from a microscopic view angle, and is a powerful means for researching interfacial polymerization reaction mechanisms of separation membranes. The Dissipative Particle Dynamics (DPD) simulation technology can make up the defects of experimental means, provides a new thought for deeply researching thermodynamic and kinetic influence factors of interfacial polymerization reaction of a composite membrane separation layer, and provides a reliable theoretical basis for exploring interfacial polymerization reaction mechanisms so as to further improve the interfacial polymerization process and obtain high-performance membrane materials.
Disclosure of Invention
Aiming at the defects of the prior art, the invention provides a method for researching and researching thermodynamic and kinetic influence factors of interfacial polymerization reaction of a composite membrane separation layer from a microscopic angle by using Materials Studio calculation software and DPD simulation so as to develop the understanding of the interfacial polymerization reaction mechanism of a high-performance separation membrane.
The purpose of the invention is realized by the following technical scheme:
a dissipative particle dynamics method for researching a composite film interfacial polymerization reaction mechanism comprises the following steps:
determining various substance components and chemical structures thereof in a water phase and an oil phase of an interfacial polymerization system, wherein the water phase comprises water-soluble water and water-soluble monomers, and the oil phase comprises an organic solvent and oil-soluble monomers;
step two, respectively constructing a DPD model of the water phase solvent and the water-soluble monomer molecules and a DPD model of the organic solvent and the oil-soluble monomer molecules in the oil phase;
step three, establishing a DPD structural model of an interfacial polymerization reaction system consisting of a water phase and an oil phase;
step four, calculating interaction parameters among all DPD beads, namely conservative force parameters;
fifthly, carrying out DPD simulation by using Materials Studio software, and obtaining a motion track file and a related calculation file of each DPD bead after the system is balanced;
and step six, observing the structural characteristics of the separation layer of the composite membrane generated by the interfacial polymerization reaction according to the simulation result, and researching influence factors determining the performance of the separation layer in the interfacial polymerization process by combining the calculation file.
Further, the water-soluble monomer is one or more of piperazine, 2-methylpiperazine, 2, 5-dimethylpiperazine, 4-aminomethylpiperazine, 2, 5-diethylpiperazine, alpha-cyclodextrin, beta-cyclodextrin, gamma-cyclodextrin, delta-cyclodextrin, p-phenylenediamine, m-phenylenediamine, trimesamine, diaminotoluene, ethylenediamine, propylenediamine, xylylenediamine, 1, 3-diaminocyclohexane or 1, 4-diaminocyclohexane, and the concentration of the water-soluble monomer is 0.01 to 8.0 wt%.
Further, the water-soluble monomer is preferably piperazine, m-phenylenediamine or cyclodextrin.
Further, the organic solvent in the oil phase is one or more of n-hexane, cyclohexane, heptane, octane, naphtha, Isopar-E, Isopar-G, Isopar-L or mineral oil.
Further, the oil-soluble monomer in the oil phase is a polyacyl chloride monomer, and comprises one or more of trimesoyl chloride, terephthaloyl chloride, isophthaloyl chloride, terephthaloyl chloride, benzenetrisulfonyl chloride, tricaprylyl chloride, butanetriacyl chloride, pentatriacyl chloride, glutaroyl chloride, adipoyl chloride, maleoyl chloride, cyclopropanetriacyl chloride, cyclobutanetriacyl chloride, cyclobutane tetracoyl chloride, cyclopentanedioyl chloride, cyclopentane triacyl chloride, cyclopentane tetracoyl chloride, cyclohexane diacyl chloride, cyclohexane triacyl chloride or cyclohexane tetracoyl chloride, and the concentration of the oil-soluble monomer is 0.01-4.0 wt%.
Further, the specific steps of the second step are as follows:
(1) coarse granulating each substance in the system according to the chemical structures of the solvent and the monomer in the water phase and the oil phase, and defining DPD beads of different types;
(2) the types of the beads are set by using a Materials Visualizer module of Materials Studio software, and corresponding DPD beads are used for constructing DPD molecular models of water molecules, water-soluble monomer molecules, organic solvent molecules and oil-soluble monomer molecules.
Further, the third step comprises the following specific steps:
(1) building a cubic box by using Materials Studio software, averagely dividing the box into an upper layer and a lower layer, setting the upper layer as an organic phase in the interfacial polymerization reaction process, placing an organic solvent and an oil-soluble monomer, setting the lower layer as a water phase in the interfacial polymerization reaction process, and placing water molecules and a water-soluble monomer;
(2) the number of water molecules, water-soluble monomer molecules, organic solvent molecules, and oil-soluble monomer molecules is determined by the monomer concentration required for the interfacial polymerization reaction.
Further, the fourth step is as follows:
the Flory-Huggins parameters of the beads are obtained by molecular dynamics simulation or through a reference document, and then the interaction force parameters among the coarse-grained beads are calculated according to the DPD theory.
Further, the step five is specifically as follows:
(1) selecting a Geometry Optimization task in a Mesocite module to carry out structural Optimization on the constructed interface aggregation system;
(2) carrying out DPD simulation on the interface polymerization system DPD model obtained by the construction method in the third step and the conservative force parameter obtained in the fourth step by utilizing a Mesocite module of Materials Studio software to obtain a stable equilibrium state interface polymerization structure;
(3) and outputting and storing a motion track file and a related calculation file of each bead in the DPD simulation, wherein the related files refer to an interaction energy file, a concentration file, a density file, a radial distribution function file, a mean square displacement file and a mutual distance file.
Further, the sixth step is as follows:
(1) outputting the structure of the interface polymerization system DPD model obtained in the fifth step when the model reaches a stable equilibrium state, and observing the motion tracks of all DPD beads;
(2) according to a motion track file and a related calculation file of beads in which water-soluble monomers and oil-soluble monomers react with each other, making a snapshot to reflect a track evolution diagram of a structure of a polymer layer formed by the interfacial polymerization reaction of the two monomers at an interface along with time, calculating an evolution diagram of the concentrations of the water-soluble monomers and the oil-soluble monomers in the interface polymerization layer along with time, and inspecting the concentration distribution of the two monomers near the interface polymerization layer at different times;
(3) through the analysis, the factors determining the structure and the performance of the polymer in the interfacial polymerization reaction process and the influence regulation thereof are analyzed.
Aiming at the interfacial polymerization reaction process of the composite membrane separation layer, the DPD simulation is combined with experimental research, the influence of each factor on the structure and the performance of the composite membrane separation layer is researched, the reaction mechanism of the interfacial polymerization process is explored from a microscopic angle, and a theoretical basis is laid for obtaining a high-performance membrane material by improving the interfacial polymerization process.
The DPD simulation method disclosed by the invention is used for researching the interfacial polymerization reaction mechanism of the composite film separation layer, and compared with the traditional method, the DPD simulation method has the following remarkable advantages: (1) the structure and performance of a composite membrane separation layer formed by interfacial polymerization reaction can be researched on a mesoscale, and the influence mechanism of each factor in the interfacial polymerization reaction process can be researched; (2) a dynamic visual effect diagram of the interfacial polymerization reaction can be provided, which is beneficial to further and deeply knowing the mechanism of the interfacial polymerization reaction; (3) the research result makes up the defects of experimental means, and the forming process of the composite membrane separation layer and the influence factors of the structure and the performance of the composite membrane separation layer can be visually and vividly observed from a microscopic angle; (4) the method can be applied to the fields of chemical, environmental and life sciences related to membrane material preparation, water treatment and the like.
Drawings
FIG. 1 is a schematic diagram of DPD coarse granulation model of solvent, gel film material and monomer in hydrogel phase and organic phase in interfacial polymerization process in example 2;
FIG. 2 is a graph showing the distribution of each substance in the aqueous phase and the oil phase in the initial state of interfacial polymerization in example 2.
Figure 3. example 2, the composite membrane separation layer formed by interfacial polymerization of PIP in the aqueous phase and TMC in the oil phase at equilibrium is shown.
Detailed Description
For better understanding of the present invention, the technical solution of the present invention is further described in detail with reference to the specific embodiments, but the present invention is not limited thereto, and modifications or equivalent substitutions may be made to the technical solution of the present invention without departing from the spirit and scope of the technical solution of the present invention.
Example 1
In the embodiment, the interface polymerization reaction process of the composite film separation layer is simulated and the reaction mechanism is explored on the computing server by using Materials Studio software and a DPD simulation method, and the method mainly comprises the following aspects:
determining various substance components and chemical structures thereof in a water phase and an oil phase of an interfacial polymerization system, wherein the water phase comprises water and water-soluble monomers, and the oil phase comprises an organic solvent and oil-soluble monomers;
secondly, respectively constructing a DPD model of the solvent and the water-soluble monomer molecules in the water phase and a DPD model of the organic solvent and the oil-soluble monomer molecules in the oil phase;
establishing a DPD structural model of an interfacial polymerization reaction system consisting of a water phase and an oil phase;
fourthly, calculating interaction parameter force parameters among all DPD beads;
fifthly, carrying out DPD simulation by using Materials Studio software, and obtaining a motion track file and a related calculation file of each DPD bead after the system is balanced;
and sixthly, observing the structural characteristics of the separation layer of the composite membrane generated by the interfacial polymerization reaction according to the simulation result, and researching influence factors determining the performance of the separation layer in the interfacial polymerization process by combining with a calculation file.
The method comprises the following specific steps:
(1) determining various substance components and chemical structures thereof in a water phase and an oil phase of an interfacial polymerization system, wherein the water phase comprises water and water-soluble monomers, the oil phase comprises an organic solvent and oil-soluble monomers, different interfacial polymerization monomers can form composite membrane separation layers with different structures and properties, and proper solvents and monomers can be selected according to actual conditions;
(2) coarse granulating each substance in the system according to the chemical structures of the solvent and the monomer in the water phase and the oil phase, and defining DPD beads of different types;
(3) setting the types of beads by using a Materials Visualizer module of Materials Studio software, and constructing a DPD molecular model of water molecules, water-soluble monomer molecules, organic solvent molecules and oil-soluble monomer molecules by using corresponding DPD beads;
(4) building a cube box with a proper size by using Materials Studio software, averagely dividing the box into an upper layer and a lower layer, setting the upper layer as an organic phase in the interfacial polymerization reaction process, placing an organic solvent and an oil-soluble monomer, setting the lower layer as a water phase in the interfacial polymerization reaction process, and placing water molecules and a water-soluble monomer; the number of water molecules, water-soluble monomer molecules, organic solvent molecules and oil-soluble monomer molecules is determined by the monomer concentration required for the interfacial polymerization reaction;
(5) obtaining Flory-Huggins parameters of the beads by molecular dynamics simulation or through a reference, further calculating interaction force parameters among the coarse-grained beads according to a DPD theory, wherein the relationship between the Flory-Huggins parameters and the DPD conservative force parameters is as shown in formulas (1) and (2):
wherein, aiiRepresents the interaction force parameter between identical DPD beads;
aijrepresenting the interaction force parameters between different DPD beads.
In the formula, NmRepresents the level of coarse grain in the DPD simulation, namely: the number of water molecules contained in one DPD bead;
kBt represents an energy unit in the DPD simulation;
ρ represents the density of the DPD analog system, and ρ is 3 in this embodiment;
χijindicating Flory between different DPD beadsThe Huggins parameter, which can be obtained by molecular dynamics simulation or from the authentic scientific literature.
(6) Selecting a Geometry Optimization task in a Mesocite module, and carrying out structural Optimization on the constructed interface polymerization system, wherein the Optimization conditions are set as follows: selecting a self-calculated DPD force field, optimizing the mass and selecting the Customized method, selecting the Ewald method by electrostatic action and selecting the Beadbased method by the Van der Waals force, wherein the truncation distance is
(7) Carrying out DPD simulation on the interface polymerization system DPD model obtained by the construction method in the step (4) and the conservative force parameter obtained in the step (5) by utilizing a Mesocite module of Materials Studio software to obtain a stable equilibrium state interface polymerization structure; the force fields selected in the DPD simulation process are consistent with those in the step (6) and are interaction force parameters obtained by calculation according to formulas (1) and (2), and the simulation time and the step length can be adjusted according to the actual condition and the size of a simulation system;
(8) outputting and storing a motion track file and a related calculation file of each bead in the DPD simulation, wherein the related files refer to an interaction energy file, a concentration file, a density file, a radial distribution function file, a mean square displacement file and a mutual distance file;
(9) outputting the structure of the interface polymerization system DPD model obtained in the step (7) when the model reaches a stable equilibrium state, and observing the motion tracks of all DPD beads;
(10) according to a motion track file and a related calculation file of beads in which water-soluble monomers and oil-soluble monomers react with each other, making a fast graph which reflects a track evolution graph of a structure of a polymer layer formed by the interface polymerization reaction of the two monomers at an interface along with time, calculating an evolution graph of the concentrations of the water-soluble monomers and the oil-soluble monomers in the interface polymerization layer along with time, and inspecting the concentration distribution of the two monomers near the interface polymerization layer at different times;
(11) through the analysis, the factors determining the structure and the performance of the polymer and the influence rule thereof in the interfacial polymerization reaction process are analyzed.
Example 2
In the embodiment, the method takes the example that piperazine monomer in aqueous solution and trimesoyl chloride monomer in normal hexane are subjected to interfacial polymerization to form a composite membrane separation layer, and comprises the following specific steps:
(1) determining the components and chemical structures of the materials in the water phase and the oil phase of the interfacial polymerization system, wherein the solvent and the water-soluble monomer in the water phase are respectively water (H)2O) and piperazine (PIP), the organic solvent and the oil-soluble monomer in the oil phase are n-Hexane (Hexane) and trimesoyl chloride (TMC), respectively;
(2) according to the molecular structures of water, PIP, Hexane and TMC, the substances in the system are coarsely granulated, and DPD beads of different types are defined, as shown in the Bead type column in FIG. 1:
(3) the method comprises the steps of setting the types of beads by using a Materials Visualizer module of Materials Studio software, and constructing DPD molecular models of water molecules, PIP molecules, Hexane molecules and TMC molecules by using corresponding DPD beads, as shown in the column of Coarse-grained molecules in FIG. 1;
(4) the volume is constructed by using Materials Studio softwareThe square box is divided into an upper layer and a lower layer (the volume is respectively)) The upper layer is set as an organic phase in the interfacial polymerization reaction process, Hexane and TMC molecules are placed in the organic phase, the lower layer is set as a water phase in the interfacial polymerization reaction process, water molecules and PIP molecules are placed in the organic phase, and the constructed interfacial polymerization system is shown in figure 2, wherein: the composition of the oil phase is Hexane/TMC 0.95:0.05, the composition of the water phase is water/PIP 0.9:0.1, and for observation, water molecules in the water phase and Hexane molecules in the oil phase are set to be invisible;
(5) the Flory-Huggins parameters of the beads are obtained by molecular dynamics simulation or through a reference, the interaction force parameters of the coarse-grained beads are further calculated according to a DPD theory, and the relationship between the Flory-Huggins parameters and the DPD conservative force parameters is as shown in formulas (1) and (2), and the method is the same as the embodiment 1:
(6) selecting a Geometry Optimization task in a Mesocite module, and carrying out structural Optimization on the constructed interface polymerization system, wherein the Optimization conditions are set as follows: selecting a self-calculated DPD force field, optimizing the mass and selecting 'Customized', selecting 'EWald' method by electrostatic interaction and selecting 'Bead based' by Van der Waals force, wherein the truncation distance is
(7) Carrying out DPD simulation of 500000 steps (about 75 ns) on a system with an optimized structure by utilizing a Mesocite module of Materials Studio software to obtain a stable equilibrium state interface aggregation structure, wherein force fields selected in the DPD simulation process are consistent with those in the step (6) and are conservative force parameters obtained by calculation according to formulas 1 and 2, and motion tracks of all DPD beads of one frame are stored in every 500 steps;
(8) after the equilibrium state is reached, outputting and storing a motion trail file of PIP molecules and TMC molecules in DPD simulation and a related calculation file, wherein the related file refers to an interaction energy file, a concentration file, a density file, a radial distribution function file, a mean square displacement file and a mutual distance file;
(9) outputting the structure of the interface polymerization system DPD model obtained in the step (7) when the model reaches a stable equilibrium state, and observing the motion tracks of all DPD beads;
(10) according to a motion trail file and a related calculation file of DPD molecules of PIP and TMC monomers which generate interface polymerization reaction, making a snapshot map reflecting a time-dependent locus evolution diagram of structures of polymer layers formed by the interface polymerization reaction of the two monomers at an interface, calculating a time-dependent evolution diagram of concentrations of PIP and TMC in the interface polymerization layer, and inspecting concentration distribution of the two monomers near the interface polymerization layer at different times;
(11) the separation layer of the composite membrane formed by the interfacial polymerization of PIP in the aqueous phase and TMC in the oil phase at equilibrium is shown in fig. 3. Through the analysis, the D beads in the PIP molecules and the C beads in the TMC molecules are determined to generate the polyamide polymer film through the interfacial polymerization reaction, and the concentration of the PIP and TMC monomers has important influence on the thickness, the pore diameter and the porosity of the formed polymer film.
Although the preferred embodiments of the present patent have been described in detail, the present patent is not limited to the above embodiments, and various changes and modifications can be made without departing from the spirit and scope of the present patent within the knowledge of those skilled in the art.
Claims (10)
1. A dissipative particle dynamics method for researching a composite membrane interfacial polymerization reaction mechanism is characterized by comprising the following steps:
determining various substance components and chemical structures thereof in a water phase and an oil phase of an interfacial polymerization system, wherein the water phase comprises water and water-soluble monomers, and the oil phase comprises an organic solvent and oil-soluble monomers;
step two, respectively constructing a DPD model of the water phase solvent and the water-soluble monomer molecules and a DPD model of the organic solvent and the oil-soluble monomer molecules in the oil phase;
step three, establishing a DPD structural model of an interfacial polymerization reaction system consisting of a water phase and an oil phase;
step four, calculating interaction force parameters among all DPD beads;
fifthly, carrying out DPD simulation by using Materials Studio software, and after the system is balanced, obtaining a motion trail file and a related calculation file of each DPD bead;
and step six, observing the structural characteristics of the separation layer of the composite membrane generated by the interfacial polymerization reaction according to the simulation result, and researching influence factors determining the performance of the separation layer in the interfacial polymerization process by combining with a calculation file.
2. The dissipative particle dynamics method for studying interfacial polymerization mechanism of composite membrane according to claim 1, wherein: the water-soluble monomer is one or more of piperazine, 2-methyl piperazine, 2, 5-dimethyl piperazine, 4-amino methyl piperazine, 2, 5-diethyl piperazine, alpha-cyclodextrin, beta-cyclodextrin, gamma-cyclodextrin, delta-cyclodextrin, p-phenylenediamine, m-phenylenediamine, sym-benzenetriamine, diaminotoluene, ethylene diamine, propylene diamine, phenyl dimethyl diamine, 1, 3-diaminocyclohexane or 1, 4-diaminocyclohexane, and the concentration of the water-soluble monomer is 0.01-8.0 wt%.
3. The dissipative particle dynamics method for studying interfacial polymerization mechanism of composite membrane according to claim 2, wherein: the water-soluble monomer is preferably piperazine, m-phenylenediamine or cyclodextrin.
4. The dissipative particle dynamics method for studying interfacial polymerization mechanism of composite membrane according to claim 1, wherein: the organic solvent in the oil phase is one or more of n-hexane, cyclohexane, heptane, octane, naphtha, Isopar-E, Isopar-G, Isopar-L or mineral oil.
5. The dissipative particle dynamics method for studying interfacial polymerization mechanism of composite membrane according to claim 1, wherein: the oil-soluble monomer in the oil phase is a polybasic acyl chloride monomer and comprises one or more of trimesoyl chloride, paraphthaloyl chloride, isophthaloyl chloride, paraphthaloyl chloride, benzene trisulfonyl chloride, tricaprylyl chloride, butanetriacyl chloride, pentatricolyl chloride, glutaryl chloride, adipoyl chloride, maleoyl chloride, cyclopropanetriacyl chloride, cyclobutane trisulfonyl chloride, cyclobutane tetrasulfonyl chloride, cyclopentane diacyl chloride, cyclopentane triacyl chloride, cyclopentane tetraacyl chloride, cyclohexane diacyl chloride, cyclohexane triacyl chloride or cyclohexane tetraacyl chloride, and the concentration of the oil-soluble monomer is 0.01-4.0 wt%.
6. The dissipative particle dynamics method for studying interfacial polymerization reaction mechanism of composite membrane according to claim 1, wherein the second step comprises the following steps:
(1) coarse granulating each substance in the system according to the chemical structures of the solvent and the monomer in the water phase and the oil phase, and defining DPD beads of different types;
(2) the types of the beads are set by using a Materials Visualizer module of Materials Studio software, and corresponding DPD beads are used for constructing DPD molecular models of water molecules, water-soluble monomer molecules, organic solvent molecules and oil-soluble monomer molecules.
7. The dissipative particle dynamics method for studying interfacial polymerization reaction mechanism of composite membrane according to claim 1, wherein the third step comprises the following steps:
(1) building a cubic box by using Materials Studio software, averagely dividing the box into an upper layer and a lower layer, setting the upper layer as an organic phase in the interfacial polymerization reaction process, placing an organic solvent and an oil-soluble monomer, setting the lower layer as a water phase in the interfacial polymerization reaction process, and placing water molecules and a water-soluble monomer;
(2) the number of water molecules, water-soluble monomer molecules, organic solvent molecules, and oil-soluble monomer molecules is determined by the monomer concentration required for the interfacial polymerization reaction.
8. The dissipative particle dynamics method for studying interfacial polymerization reaction mechanism of composite membrane according to claim 1, wherein the fourth step is as follows:
the Flory-Huggins parameters of the beads are obtained by molecular dynamics simulation or through a reference document, and then the interaction force parameters among the coarse-grained beads are calculated according to a DPD theory.
9. The dissipative particle dynamics method for studying interfacial polymerization reaction mechanism of composite membrane according to claim 1, wherein the fifth step is as follows:
(1) selecting a Geometry Optimization task in a Mesocite module to carry out structural Optimization on the constructed interface aggregation system;
(2) carrying out DPD simulation on the interface polymerization system DPD model obtained by the construction method in the third step and the acting force parameter obtained in the fourth step by utilizing a Mesocite module of Materials Studio software to obtain a stable equilibrium state interface polymerization structure;
(3) and outputting and storing a motion track file and a related calculation file of each bead in the DPD simulation, wherein the related files refer to an interaction energy file, a concentration file, a density file, a radial distribution function file, a mean square displacement file and a mutual distance file.
10. The dissipative particle dynamics method for studying interfacial polymerization reaction mechanism of composite membrane according to claim 1, wherein the sixth step is as follows:
(1) outputting the structure of the interface polymerization system DPD model obtained in the fifth step when the model reaches a stable equilibrium state, and observing the motion tracks of all DPD beads;
(2) according to a motion track file and a related calculation file of beads in which water-soluble monomers and oil-soluble monomers react with each other, making a fast graph which reflects a track evolution graph of a structure of a polymer layer formed by the interface polymerization reaction of the two monomers at an interface along with time, calculating an evolution graph of the concentrations of the water-soluble monomers and the oil-soluble monomers in the interface polymerization layer along with time, and inspecting the concentration distribution of the two monomers near the interface polymerization layer at different times;
(3) through the analysis, the factors determining the structure and the performance of the polymer and the influence rule thereof in the interfacial polymerization reaction process are analyzed.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910858591.9A CN110600084B (en) | 2019-09-11 | 2019-09-11 | Dissipative particle dynamics method for researching interfacial polymerization reaction mechanism of composite film |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910858591.9A CN110600084B (en) | 2019-09-11 | 2019-09-11 | Dissipative particle dynamics method for researching interfacial polymerization reaction mechanism of composite film |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110600084A true CN110600084A (en) | 2019-12-20 |
| CN110600084B CN110600084B (en) | 2022-02-15 |
Family
ID=68858808
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910858591.9A Active CN110600084B (en) | 2019-09-11 | 2019-09-11 | Dissipative particle dynamics method for researching interfacial polymerization reaction mechanism of composite film |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110600084B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111739589A (en) * | 2020-06-29 | 2020-10-02 | 青岛科技大学 | Mesoscale-based simulation method for dissolving lignin by eutectic solvent |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105138744A (en) * | 2015-08-07 | 2015-12-09 | 长春理工大学 | Dissipative particle dynamics simulation method for grinding liquid particle characteristics |
| CN110021380A (en) * | 2019-04-25 | 2019-07-16 | 济南大学 | A method of each atom scattering nature in glass system is probed into based on molecular dynamics simulation |
| CN110064312A (en) * | 2019-04-29 | 2019-07-30 | 袁书珊 | A kind of high throughput solvent resistant interfacial polymerization composite membrane and preparation method thereof |
-
2019
- 2019-09-11 CN CN201910858591.9A patent/CN110600084B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105138744A (en) * | 2015-08-07 | 2015-12-09 | 长春理工大学 | Dissipative particle dynamics simulation method for grinding liquid particle characteristics |
| CN110021380A (en) * | 2019-04-25 | 2019-07-16 | 济南大学 | A method of each atom scattering nature in glass system is probed into based on molecular dynamics simulation |
| CN110064312A (en) * | 2019-04-29 | 2019-07-30 | 袁书珊 | A kind of high throughput solvent resistant interfacial polymerization composite membrane and preparation method thereof |
Non-Patent Citations (4)
| Title |
|---|
| JIAHUI HU ET AL.: "Fabrication of a high-flux sulfonated polyamide nanofiltration membrane: Experimental and dissipative particle dynamics studies", 《JOURNAL OF MEMBRANE SCIENCE》 * |
| ZHAOHUAN MAI ET AL.: "Parameterization of a mesoscopic model for the self-assembly of linear sodium alkyl sulfates", 《THE JOURNAL OF CHEMICAL PHYSICS》 * |
| 孙德林 等: "耗散粒子动力学模拟Nafion膜和PVA/Nafion共混膜的介观结构", 《物理化学学报》 * |
| 孙红 等: "水在 Nafion 膜中扩散行为介观分析", 《沈阳建筑大学学报（自然科学版）》 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111739589A (en) * | 2020-06-29 | 2020-10-02 | 青岛科技大学 | Mesoscale-based simulation method for dissolving lignin by eutectic solvent |
| CN111739589B (en) * | 2020-06-29 | 2022-07-05 | 青岛科技大学 | Mesoscale-based simulation method for dissolving lignin by eutectic solvent |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110600084B (en) | 2022-02-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110517735B (en) | Dissipative particle dynamics method for simulating gel film interface polymerization reaction process | |
| Yan et al. | Screening and design of covalent organic framework membranes for CO2/CH4 separation | |
| Muscatello et al. | Multiscale molecular simulations of the formation and structure of polyamide membranes created by interfacial polymerization | |
| Hu et al. | A microscopic characterization of wettability in shale kerogen with varying maturity levels | |
| Phillip et al. | Diffusion and flow across nanoporous polydicyclopentadiene-based membranes | |
| Wang et al. | Molecular simulation of gas adsorption in shale nanopores: A critical review | |
| Zhang et al. | Molecular insights into the composition–structure–property relationships of polyamide thin films for reverse osmosis desalination | |
| Mollahosseini et al. | Molecular dynamics simulation for membrane separation and porous materials: A current state of art review | |
| Keskin et al. | Assessment of a metal− organic framework membrane for gas separations using atomically detailed calculations: CO2, CH4, N2, H2 mixtures in MOF-5 | |
| Rahaman et al. | Dynamics of water adsorption onto a calcite surface as a function of relative humidity | |
| CN110556168B (en) | Simulation method for analyzing diffusion property of water-soluble monomer in gel film | |
| Priya et al. | Machine learning assisted screening of two-dimensional materials for water desalination | |
| Lin et al. | Relative importance of geometrical and intrinsic water transport properties of active layers in the water permeability of polyamide thin-film composite membranes | |
| Liu et al. | Interfacial resistance and length-dependent transport diffusivities in carbon nanotubes | |
| Kong et al. | An atomistic simulation study on POC/PIM mixed-matrix membranes for gas separation | |
| CN110600084B (en) | Dissipative particle dynamics method for researching interfacial polymerization reaction mechanism of composite film | |
| He et al. | Molecular insights into the structure-property relationships of 3D printed polyamide reverse-osmosis membrane for desalination | |
| Liyana-Arachchi et al. | Ultrathin molecular-layer-by-layer polyamide membranes: insights from atomistic molecular simulations | |
| Liu et al. | A highly rigid and conjugated microporous polymer membrane for solvent permeation and biofuel purification: a molecular simulation study | |
| Tang et al. | Crystallization kinetics and mechanisms of calcium sulfate dihydrate: experimental investigation and theoretical analysis | |
| Lin et al. | Effects of void morphology on fracturing characteristics of porous rock through a finite-discrete element method | |
| Qiu et al. | Large-scale metal–organic framework nanoparticle monolayers with controlled orientation for selective transport of rare-earth elements | |
| Bushuev et al. | Subnanometer topological tuning of the liquid intrusion/extrusion characteristics of hydrophobic micropores | |
| Zhang et al. | Unveiling the molecular mechanisms of thickness-dependent water dynamics in an ultrathin free-standing polyamide membrane | |
| Dutta et al. | System size-dependent transport properties in materials of nanoscale dimension |
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 |