MOFs material structure performance analysis method, storage medium and electronic equipment

Abstract

The invention relates to an MOFs material structure performance analysis method, a storage medium and electronic equipment. The method comprises the steps that a first structure file of a target MOFs material model is acquired, and a first input parameter file is made; reading the first structure file and the first input parameter file by using DFTB + software for geometric optimization to obtain a second structure file and compound energy and unit cell parameters thereof; if the target MOFs material is of a newly designed structure, forming energy is calculated based on a theoretical chemical equation, simple substance energy and compound energy; if the target MOFs material is of an existing structure, forming energy is calculated based on an existing chemical equation and energy of reactants and products of the equation; making a second input parameter file; reading the second structure file and the second input parameter file by using DFTB + software to carry out state density calculation and visualization to obtain a state density map; and based on the cell parameters, forming an energy and state density map to analyze the target MOFs material.

Description

Technical Field

[0001] This application relates to the technical field of material analysis, and specifically, to a method for analyzing the structural properties of MOFs materials, a storage medium, and an electronic device. Background Art

[0002] Metal-organic frameworks (MOFs), also known as metal-organic skeletons, are a new type of porous material and have attracted much attention due to their wide application potential in the fields of gas adsorption, separation, catalysis, etc. The structure of MOFs consists of metal centers (secondary building units, SBUs) and organic ligands, and their pore size, topological structure, and functional characteristics can be precisely regulated by adjusting the types of metal centers and organic ligands. With the progress of MOFs synthesis and modification technologies, researchers can flexibly change their structures and chemical properties by means such as introducing functional groups, exchanging metals and ligands, and encapsulating nano-guests, further expanding their application fields. Currently, more than 100,000 different MOFs have been successfully synthesized.

[0003] However, with the rapid increase in the types of MOFs materials, the number of MOFs that can be synthesized in the future may reach trillions, far exceeding the scope that traditional experimental methods can handle. It is unrealistic to rely solely on experimental synthesis and testing of all potential MOFs materials. Especially under limited time and resources, the efficiency of experimental screening faces severe challenges. The research and development of MOFs materials urgently need to rely on computational simulations. Simulations can use chemical calculation software developed by commercial companies such as Material studio, etc., but these software have problems such as high fees. Summary of the Invention

[0004] The purpose of this application is to provide a method for analyzing the structural properties of MOFs materials, a storage medium, and an electronic device, aiming to analyze the structural properties of MOFs materials through simulation calculations to provide analysis data for MOF material screening, improve the efficiency of screening MOFs materials, reduce the screening cost, and this process is implemented based on the free and open-source software DFTB+, reducing the cost of simulation calculations.

[0005] To achieve the above purpose, in the first aspect of this application, a method for analyzing the structural properties of MOFs materials is provided. The method includes: Obtain the first structure file of the target MOFs material model and create the first input parameter file for geometric optimization of the target MOFs material model; Run the geometric optimization command in the DFTB+ software to read the first structure file and the first input parameter file and perform geometric optimization simulations to obtain the second structure file, compound energy, and unit cell parameters of the geometrically optimized target MOFs material model; If the target MOF material has a newly designed structure, calculate the formation energy of the target MOF material based on the theoretical chemical equation for synthesizing the target MOF material from the elemental substances of the constituent elements of the target MOF material, the energies of the elemental substances, and the energy of the compound; If the target MOF material has an existing structure, calculate the formation energy of the target MOF material based on the existing chemical equation for synthesizing the target MOF material and the energies of the reactants and products in this existing chemical equation, where the energy of the product includes the energy of the compound; Create a second input parameter file for the target MOF material model for density of states calculation; Run the density of states calculation command in the DFTB+ software to read the second structure file and the second input parameter file and perform density of states calculation, and visualize the density of states calculation results to obtain the density of states diagram of the geometrically optimized target MOF material model; Analyze the target MOF material based on the unit cell parameters, formation energy, and density of states diagram.

[0006] Optionally, the calculating the formation energy of the target MOF material based on the theoretical chemical equation for synthesizing the target MOF material from the elemental substances of the constituent elements of the target MOF material, the energies of the elemental substances, and the energy of the compound includes: Determine the elemental substances of all constituent elements of the target MOF material based on the chemical formula of the target MOF material; Determine the reaction equation for generating the target MOF material using all the elemental substances and balance it to obtain the theoretical chemical equation; Calculate the energies of all the elemental substances; Multiply the product of the energy of the compound by the stoichiometric number of the target MOF material in the theoretical chemical equation, subtract the product of the energy of each elemental substance and the stoichiometric number of the elemental substance in the theoretical chemical equation, to obtain the formation energy.

[0007] Optionally, the calculating the formation energy of the target MOF material based on the existing chemical equation for synthesizing the target MOF material and the energies of the reactants and products in this existing chemical equation includes: Obtain the existing chemical equation for synthesizing the target MOF material; Calculate the energies of all reactants and all products except the target MOF material in the existing chemical equation; Multiply the product of the energy of the compound by the stoichiometric number of the target MOF material in the existing chemical equation, add the product of the energy of each product except the target MOF material and the stoichiometric number of the product in the existing chemical equation, subtract the product of the energy of each reactant and the stoichiometric number of the reactant in the existing chemical equation, to obtain the formation energy.

[0008] Optionally, calculating the energies of all the simple substances includes: For each simple substance, obtaining the third structure file of the simple substance and creating a third input parameter file for geometric optimization of the simple substance; Running a geometric optimization command in the DFTB+ software to read the third structure file and the third input parameter file and perform geometric optimization simulation to obtain the energies of each geometrically optimized simple substance.

[0009] Optionally, calculating the energies of all reactants and all products except the target MOFs material in the existing chemical equation includes: For each reactant and each product except the target MOFs material in the existing chemical equation, obtaining the fourth structure file of the reactant or product and creating a fourth input parameter file for geometric optimization; Running a geometric optimization command in the DFTB+ software to read the fourth structure file and the fourth input parameter file and perform geometric optimization simulation to obtain the energies of each of the reactants and each of the products after geometric optimization.

[0010] Optionally, creating the first input parameter file for geometric optimization of the target MOFs material model includes: Creating a first initial file in hsd format; Writing the Geometry module in the first initial file and setting the object of geometric optimization to the first structure file in this module; Writing the Driver module in the first initial file and setting the calculation type to geometric optimization in this module, setting to enable full lattice optimization and use a rational function optimizer, setting the maximum number of optimization steps to any step between 100 and 200 steps, and setting the convergence criterion to GradElem of the gradient element being less than 1e-4; Writing the Hamiltonian module in the first initial file and setting to enable self-consistent charge calculation in this module, with the convergence value and the maximum number of iterations of this calculation being 1e-6 and 2000 times respectively, setting the Slater-Koster parameters to use the 3OB parameter set, setting the maximum angular momentum quantum number of each element of the target MOFs material model, setting the k-point sampling in the Brillouin zone to use the Monkhorst-Pack method and only selecting the Γ point as the k-point sampling point; Writing the Analysis module in the first initial file and setting the calculation properties in this module; Writing the ParserOptions module in the first initial file and setting the parser version in this module, finally obtaining the first input parameter file for geometric optimization of the target MOFs material model.

[0011] Optionally, the second input parameter file for calculating the density of states in the target MOFs material model includes: Create a second initial file in hsd format; Write the Geometry module in the second initial file, and set the object for density of states calculation as the second structure file in this module; Write the Hamiltonian module in the second initial file, and in this module, set to enable self-consistent charge calculation with a convergence value of 1e-5 for this calculation, set the maximum angular momentum quantum number for each element of the target MOFs material model, and set the k-point sampling density in the Brillouin zone as an 8×8×8 grid; Write the Analysis module in the second initial file, and in this module, set the projection region for density of states calculation of each element of the target MOFs material model; Write the ParserOptions module in the second initial file, and in this module, set the parser version, and finally obtain the second input parameter file for calculating the density of states in the target MOFs material model.

[0012] Optionally, the analysis of the target MOFs material based on the unit cell parameters, formation energy, and density of states diagram includes: Analyze the unit cell shape and size of the target MOFs material based on the unit cell parameters; When the target MOFs material is a newly designed structure, if the formation energy is greater than 0, determine that the target MOFs material cannot be formed by spontaneous reaction; When the target MOFs material is an existing structure, determine one or more existing chemical equations with the minimum formation energy for experimental verification; Conduct total density of states analysis and partial density of states analysis on the target MOFs material based on the density of states diagram.

[0013] In the second aspect of the present application, there is provided a computer-readable storage medium, on which a computer program is stored, and when the program is executed by a processor, the steps of the method described in any one of the above first aspects are implemented.

[0014] In the third aspect of the present application, there is provided an electronic device, including: A memory, on which a computer program is stored; A processor, configured to execute the computer program in the memory to implement the steps of the method described in any one of the above first aspects.

[0015] Through the above technical solution, for the target MOFs material model to be analyzed, obtain its first structure file and create its first input parameter file. Read the first structure file and the first input parameter file through DFTB+ for geometric optimization simulation to perform geometric optimization on the target MOFs material model, and obtain the compound energy and unit cell parameters. For different situations of the target MOFs material, if it is a newly designed structure, calculate the formation energy based on the theoretical chemical equation for the synthesis of elemental components, elemental energy, and compound energy; if it is an existing structure, calculate the formation energy based on the existing chemical equation, reactant, and product energies. Finally, create a second input parameter file, and use DFTB+ to read this file and the second structure file of the geometrically optimized target MOFs material model for density of states calculation and visualization to obtain the density of states diagram. After the operation is completed, the target MOFs material can be analyzed based on the calculated unit cell parameters, formation energy, and density of states diagram. Thus, the simulation calculation and analysis of MOFs materials are realized based on the open-source and free software DFTB+, which can provide more data support for screening MOFs materials, improve the screening efficiency, and reduce the time and cost of screening.

[0016] Other features and advantages of the present application will be described in detail in the subsequent specific implementation section. Brief Description of the Drawings

[0017] The drawings are used to provide a further understanding of the present application, and constitute a part of the specification. Together with the following specific implementation, they are used to explain the present application, but do not constitute a limitation to the present application. In the drawings: Figure 1 is a flowchart of a method for analyzing the structural performance of a MOFs material shown according to an exemplary embodiment.

[0018] Figure 2 is a structural model diagram of Mg-MOF-74 after geometric optimization shown according to an exemplary embodiment.

[0019] Figure 3 is a density of states diagram of Mg-MOF-74 shown according to an exemplary embodiment.

[0020] Figure 4 is a block diagram of an electronic device shown according to an exemplary embodiment. Detailed Description of the Invention

[0021] The following detailed description of the specific implementation of the present application is provided in conjunction with the drawings. It should be understood that the specific implementation described herein is only used to illustrate and explain the present application, and is not used to limit the present application.

[0022] Metal-organic framework materials (MOFs) have shown broad application prospects in many fields such as air and water capture, carbon capture (carbon dioxide capture), and adsorption and separation of volatile organic compounds (VOCs) due to their highly tunable pore structures, large specific surface areas, and excellent adsorption properties, and possess important industrial and environmental practical values. However, the number of reported MOF compounds is huge (for example, more than 100,000 MOFs are included in the Cambridge Structural Database CSD). If only relying on traditional experimental methods for performance analysis and screening, there are many inconveniences: 1. Time and cost limitations of experimental screening: Although MOF materials have broad application prospects in many fields, the large and increasing number of MOF materials makes it infeasible to rely on traditional experiments to synthesize and test one by one. Under limited time and resources, simply relying on experimental screening methods cannot cover all potential MOF materials, resulting in low efficiency and high cost in the screening process.

[0023] 2. High complexity of structure and performance: The structure and performance of MOF materials are affected by the complex interactions of multiple factors such as metal centers, ligands, pore structures, and functional groups, and it is difficult to comprehensively evaluate them by traditional experimental means. With the continuous progress of synthesis technologies, the variety of MOF materials has increased, making the systematic study of their performance more difficult.

[0024] 3. Lack of efficient screening and optimization tools: Currently, there is a lack of efficient screening and optimization tools in the research and development process of MOFs. Traditional experimental methods not only require a large number of experiments and time, but also need to be repeatedly debugged and improved, increasing the investment of time and cost. In addition, existing research methods are difficult to efficiently process the huge MOF material library, restricting the potential of MOF materials in applications.

[0025] Therefore, with the rapid increase in the variety of MOF materials and their high structural complexity, traditional experimental analysis methods are difficult to meet the requirements of rapid, efficient screening and optimization. In order to accelerate the research and development process of MOF materials, new methods based on computational simulation need to be introduced to analyze the structural performance of MOF materials, providing a basis for MOF material screening, thereby reducing the experimental workload, lowering costs, and improving the research and development efficiency.

[0026] DFTB+ (Density Functional Tight Binding Plus) is a fast quantum mechanics calculation software based on the density functional tight binding (DFTB) theory, suitable for electronic structure calculations and molecular dynamics simulations of systems such as molecules, clusters, crystals, surfaces, and nanomaterials. As an approximation method based on the density functional theory (DFT), DFTB+ greatly improves the calculation efficiency by using pre-computed DFT parameter files, especially suitable for large-scale molecular dynamics simulations.

[0027] Figure 1 is a flowchart of a method for analyzing the structural properties of a MOFs material shown according to an exemplary embodiment. Refer to Figure 1 , the method includes: S101, obtain the first structure file of the target MOFs material model, and make the first input parameter file for geometric optimization of the target MOFs material model.

[0028] S102, run the geometric optimization command in the DFTB+ software to read the first structure file and the first input parameter file and perform geometric optimization simulation to obtain the second structure file, compound energy, and unit cell parameters of the geometrically optimized target MOFs material model.

[0029] S103, if the target MOFs material is a newly designed structure, calculate the formation energy of the target MOFs material based on the theoretical chemical equation for synthesizing the target MOFs material from the elemental substances of the constituent elements of the target MOFs material, the energy of the elemental substances, and the compound energy.

[0030] S104, if the target MOFs material is an existing structure, calculate the formation energy of the target MOFs material based on the existing chemical equation for synthesizing the target MOFs material and the energies of the reactants and products of the existing chemical equation, where the energy of the products includes the compound energy.

[0031] S105, make the second input parameter file for density of states calculation of the target MOFs material model.

[0032] S106, run the density of states calculation command in the DFTB+ software to read the second structure file and the second input parameter file and perform density of states calculation, and visualize the density of states calculation results to obtain the density of states diagram of the geometrically optimized target MOFs material model.

[0033] S107, analyze the target MOFs material based on the unit cell parameters, formation energy, and density of states diagram.

[0034] In step S101, the target MOFs material model can be a computer modeling model of the MOFs material to be analyzed, which is reflected by the first structure file. The first structure file can be obtained by constructing the target MOFs material model in chemical calculation software. In addition, it is also possible to download the existing structural data of the target MOFs material model from a database and obtain the first structure file through preprocessing. For example, the crystal structure data file in cif format of the corresponding target MOFs material can be downloaded from the Cambridge Crystallographic Data Centre (CCDC), and the solvent coordinated with metal ions is removed to generate a reasonable input structure for calculation. Then, the cif format crystal structure data file after desolvation is converted into the gen format recognizable by the DFTB+ software to obtain the first structure file. For example, when the target MOFs material is Mg-MOF-74, the cif file with the CCDC number 1863524 of this material can be downloaded from the CCDC, and after removing the solvent coordinated with metal ions, it is converted into the gen format using a self-written script, and this gen format file can be used as the first structure file.

[0035] The input parameter file includes the first input parameter file, the second input parameter file, the third input parameter file, and the fourth input parameter file, which can be input parameter files in hsd format and record multiple parameter settings for simulation calculations.

[0036] In step S101, the first input parameter file is made as a preparation file for geometric optimization of the target MOFs material model in the DFTB+ software, and the calculation parameter settings required for geometric optimization are written in the first input parameter file.

[0037] Optionally, in step S101, the first input parameter file for geometric optimization of the target MOFs material model is made, including: S1011, create a first initial file in hsd format.

[0038] S1012, write the Geometry module in the first initial file, and set the object of geometric optimization as the first structure file in this module.

[0039] S1013, write the Driver module in the first initial file, and set the calculation type as geometric optimization in this module, set to enable full lattice optimization and use a rational function optimizer, set the maximum number of optimization steps as any step between 100 and 200 steps, and set the convergence criterion as the gradient element GradElem being less than 1e-4.

[0040] S1014. Write the Hamiltonian module in the first initial file, and in this module, set to enable self-consistent charge calculation, with the convergence value and the maximum number of iterations of this calculation being 1e-6 and 2000 times respectively. Set the Slater-Koster parameters to use the 3OB parameter set. Set the maximum angular momentum quantum number of each element of the target MOFs material model. Set the k-point sampling in the Brillouin zone to use the Monkhorst-Pack method and only select the Γ point as the k-point sampling point.

[0041] S1015. Write the Analysis module in the first initial file, and in this module, set the calculation properties.

[0042] S1016. Write the ParserOptions module in the first initial file, and in this module, set the parser version, and finally obtain the first input parameter file for the geometric optimization of the target MOFs material model.

[0043] In step S1011, a blank first initial file in hsd format can be created first, and then the corresponding parameter settings can be written. After step S1011 is executed, proceed to execute steps S1012 to S1016. Steps S1012 to S1016 are all to write the corresponding modules in the first initial file. The application does not limit their execution order. After all modules are written, the first input parameter file can be obtained.

[0044] In step S1012, write the Geometry module in the first initial file. Through this module, set the object of the geometric optimization simulation calculation to the first structure file. For example, the file name of the first structure file in gen format can be directly introduced to specify the geometric optimization object as the first structure file, or the structure parameters in the first structure file can be copied to the Geometry module to complete the setting of the geometric optimization object.

[0045] In step S1013, write the Driver module in the first initial file. Through this module, specify that the simulation calculation performed on the first structure file is geometric optimization, and optimize the ground state structure of the target MOFs material model by the full geometric optimization method to obtain a structure with completely relaxed atomic positions and lattice parameters, without applying any symmetry constraints.

[0046] To improve the accuracy of calculations, lattice optimization is enabled in the module by the code "LatticeOpt = Yes", and the rational function optimizer (Rational Optimizer) is set to perform full lattice optimization through the code "Optimizer = Rational{}", which can effectively handle complex systems containing a large number of atoms. By setting the maximum number of optimization steps between 100 and 200 steps in the module, for example, 150 steps, and setting the convergence criterion as the gradient element GradElem being less than 1e-4, the termination condition for geometric optimization is set. In addition, other contents can also be set in the module. For example, the file name of the file output after the geometric optimization calculation can be set, etc. This application does not make specific restrictions on this.

[0047] In step S1014, the Hamiltonian module is written into the first initial file, and the parameter settings for the DFTB operation are written through this module. Specifically, the self-consistent charge operation can be enabled in this module by the code "Scc=Yes" to obtain an accurate electronic structure, and the convergence value for the self-consistent charge calculation is set to 1e-6, indicating that the calculation stops when the charge density converges to 1e-6. The maximum number of iterations for the self-consistent charge calculation is set to 2000, indicating that at most 2000 self-consistent iterative calculations are performed.

[0048] The Slater-Koster parameters are model parameters that describe the transition of electrons between atoms, and they play a crucial role in calculating and simulating electronic structures. The Slater-Koster parameter file set contains multiple files recording the Slater-Koster parameters, and each Slater-Koster parameter file can correspond to a set of element pairs of molecular structures. Since the geometric optimization operation of the DFTB+ software requires the use of Slater-Koster parameters, for MOFs materials, the 3OB parameter set can be set in the Hamiltonian module, and this parameter set provides a good balance between calculation efficiency and accuracy.

[0049] In the Hamiltonian module, it is also possible to specify the maximum angular momentum quantum number for the constituent elements of the target MOF material, that is, to define the orbital type used for each element. Taking the aforementioned Mg-MOF-74 as an example, it can be set that the magnesium element (Mg) uses the p orbital, the hydrogen element (H) uses the s orbital, the carbon element (C) uses the p orbital, and the oxygen element (O) uses the p orbital, which can be achieved by setting through the MaxAngularMomentum sub-module of the Hamiltonian module. Since the primitive unit cells of MOF materials are generally numerous and the structures are relatively complex, for example, the total number of atoms in Mg-MOF-74 is 162, so the k-point sampling in the Brillouin zone can be set to use the Monkhorst-Pack method in the Hamiltonian module, and only the Γ point is selected as the k-point sampling point. The Monkhorst-Pack method is a commonly used method for generating k-point grids when calculating supercells. By generating a uniform k-point grid to sample the Brillouin zone, the calculation accuracy can be improved. The relatively large periodic unit of MOF materials makes sampling at the Γ point sufficient to accurately describe its electronic structure, and in this way, the calculation efficiency can also be improved. Specifically, it can be achieved in the KPointsAndWeights sub-module of the Hamiltonian module.

[0050] In step S1015, the Analysis module can be written into the first initial file. The Analysis module can be used to configure and control how to analyze calculation results and calculation properties. For example, "CalculateForces = Yes" can be written to set the calculation properties. This module provides users with a variety of options to extract important information from the simulation results and further process and analyze the calculation data. Through this module, users can select which calculation properties need to be analyzed or extracted and generate relevant result outputs.

[0051] In step S1016, the ParserOptions module can be written into the first initial file. For example, the parser version can be set by the command "ParserVersion = 14". After all the above settings are written into the first initial file, the first input parameter file for the geometric optimization of the target MOF material model is finally obtained.

[0052] After step S101 is completed, step S102 can be executed. In step S102, a geometry optimization command can be run in the DFTB+ software. The DFTB+ software runs commands to read the first structure file and the first input parameter file, and performs a geometry optimization simulation calculation on the target MOFs material model of the first structure file based on the settings of the first input parameter file. After the calculation is completed, a second structure file is obtained. The second structure file contains the geometrically optimized target MOFs material model. In addition, after the calculation is completed, the compound energy and unit cell parameters of the target MOFs material model can also be obtained. The unit cell parameters include the edge length and angles of the unit cell. For example, for Mg-MOF-74, the calculated compound energy is -8688.13 eV, and the unit cell parameters can be seen in Table 1. The geometrically optimized model can be seen in Figure 2 。

[0053] Table 1 Unit cell parameter table of Mg-MOF-74

[0054] After step S102 is completed, step S103 or step S104 is executed. If the target MOFs material is a newly designed structure and there is no current experimental synthesis route, step S103 can be entered. Optionally, in step S103, the formation energy of the target MOFs material is calculated based on the theoretical chemical equation for synthesizing the target MOFs material from the elemental substances of the composition elements of the target MOFs material, the energy of the elemental substances, and the compound energy, including: S1031, determine the elemental substances of all the composition elements of the target MOFs material based on the chemical formula of the target MOFs material.

[0055] S1032, determine the reaction formula for generating the target MOFs material using all the elemental substances and balance it to obtain the theoretical chemical equation.

[0056] S1033, calculate the energy of all the elemental substances.

[0057] S1034, multiply the product of the compound energy and the stoichiometric coefficient of the target MOFs material in the theoretical chemical equation, and subtract the product of the energy of each elemental substance and its stoichiometric coefficient in the theoretical chemical equation to obtain the formation energy.

[0058] In step S1031, first determine the elemental substances of its composition elements according to the chemical formula of the target MOFs material. Continuing with the previous example, assuming that Mg-MOF-74 is a newly designed structure and its chemical formula is C 72 H 18 O 54 Mg 18 , it can be determined that the elemental substances of its composition elements are C, H 2 , O2 and Mg.

[0059] After step S1031 is completed, step S1032 can be executed. In step S1032, determine the reaction formula for generating the target MOFs material using all the simple substances and balance it to obtain the theoretical chemical equation, that is, the reaction formula for theoretically generating the target MOFs material using simple substances. Through this theoretical chemical equation, the stoichiometric coefficients of the reactants and products can be determined, that is, the coefficients in front of the substances participating in the reaction in the chemical reaction equation, to provide data for subsequent calculations. Continuing with the previous example, using C, H 2 , O 2 and Mg to generate the theoretical chemical equation of Mg-MOF-74 can be: 18Mg + 27O 2 + 72C + 9H 2 → C 72 H 18 O 54 Mg 18 (1) After step S1032 is completed, step S1033 can be executed. In step S1033, calculate the energies of all simple substances. In step S1033, calculating the energies of all simple substances can include: S10331, for each simple substance, obtain the third structure file of the simple substance and create the third input parameter file for geometric optimization of the simple substance.

[0060] S10332, run the geometric optimization command in the DFTB+ software to read the third structure file and the third input parameter file and perform geometric optimization simulation to obtain the energies of each geometrically optimized simple substance.

[0061] In step S10331, for each simple substance, the third structure file of the simple substance can be obtained first, and the third input parameter file for geometric optimization can be created. The specific execution process can refer to step S101.

[0062] For example, a model of a simple substance can be constructed using chemical modeling software and converted into the gen format to obtain a third structure file, or the existing model structure file can be directly downloaded from a database. Since the third input parameter file is also used for geometric optimization, reference can be made to the description of the production process of the first input parameter file in step S101. By creating a new hsd format file and writing multiple settings, the third input parameter file can be obtained. The settings of its modules can also refer to the first input parameter file, using the same modules, and still using a rational function optimizer, using a 3OB parameter set, etc. The k-point sampling in the Brillouin zone can still use the Monkhorst-Pack method, and only the Γ point is selected as the k-point sampling point. The difference is that for the third input parameter file of each simple substance, the orbital type setting for the constituent elements only includes the elements of that simple substance.

[0063] After step S10331 is completed, step S10332 can be executed. In step S10332, by running the command for geometric optimization in the DFTB+ software, the software reads the third structure file and the third input parameter file of each simple substance, and performs geometric optimization simulation calculations on the simple substance model of the third structure file based on the settings of the third input parameter file. After the calculation is completed, the energy of the geometrically optimized simple substance model can be obtained. For example, for Mg-MOF-74, the energy of the simple substance after the calculation is completed can be seen in Table 2. Of course, in other possible implementation manners, other methods can also be used to calculate the energy of the simple substance, and the present application does not make specific limitations on this.

[0064] Table 2 Total energy of Mg-MOF-74 and energies of C, O 2 , H 2 and Mg

[0065] After step S1033 is completed, step S1034 is executed. In step S1034, based on the theoretical chemical equation, the energy of the product is subtracted from the energy of the reactant to obtain the formation energy, that is, the product of the compound energy of the target MOFs material and the stoichiometric coefficient of the target MOFs material in the theoretical chemical equation is subtracted from the product of the first energy of each simple substance in the theoretical chemical equation. The first energy product is the product of the energy of the simple substance and the stoichiometric coefficient of the simple substance in the theoretical chemical equation.

[0066] Continuing with the previous example, for Mg-MOF-74, referring to Table 2 and the theoretical chemical equation (1), its formation energy can be calculated by the formula E = -8688.13 - (-38.06)×72 - (-18.3)×9 - (-9.4)×18 - (-175.5)×27.

[0067] If the target MOF material has an existing structure, that is, there is currently a path for experimental synthesis, then the execution can proceed to step S104. Optionally, in step S104, if the target MOF material has an existing structure, calculate the formation energy of the target MOF material based on the existing chemical equation for synthesizing the target MOF material and the energies of the reactants and products in this existing chemical equation, including: S1041. Obtain the existing chemical equation for synthesizing the target MOF material.

[0068] S1042. Calculate the energies of all reactants and all products except the target MOF material in the existing chemical equation.

[0069] S1043. Multiply the energy of the compound by the stoichiometric coefficient of the target MOF material in the existing chemical equation, add the product of the energy of each product except the target MOF material and its stoichiometric coefficient in the existing chemical equation, and subtract the product of the energy of each reactant and its stoichiometric coefficient in the existing chemical equation to obtain the formation energy.

[0070] In step S1041, obtain the existing chemical equation for synthesizing the target MOF material, that is, the actually achievable synthesis path. Taking Mg-MOF-74 as an example, select the following existing chemical equation for it: 18MgO + 9C 8 H 6 O 6 → C 72 H 18 O 54 Mg 18 + 18H 2 O (2) After step S1041 is completed, proceed to execute step S1042. In step S1042, calculate the energies of all reactants and products in the existing chemical equation, except that the energy of the compound of the target MOF material has already been calculated. Optionally, in step S1042, calculate the energies of all reactants and all products except the target MOF material in the existing chemical equation, including: S10421. For each reactant and each product except the target MOF material in the existing chemical equation, obtain the fourth structure file of this reactant or product, and create a fourth input parameter file for geometric optimization.

[0071] S10422. Run the geometric optimization command in the DFTB+ software to read the fourth structure file and the fourth input parameter file and perform geometric optimization simulations to obtain the energies of each optimized reactant and each optimized product.

[0072] In step S10421, for each reactant in the existing chemical equation and each product other than the target MOF material, a fourth structure file can be obtained by constructing a model or downloading from an existing database, and a fourth input parameter file for its geometry optimization can be made. Specifically, reference can be made to the description of step S10331. The difference lies in that for the fourth input parameter file of each reactant or product, the orbital type setting for the constituent elements only includes the constituent elements of that reactant or product.

[0073] After step S10421 is completed, step S10422 is executed. In step S10422, by running the command for geometry optimization in the DFTB+ software, the software reads the fourth structure file and the fourth input parameter file of each reactant or product, and performs a geometry optimization simulation calculation on the model of the fourth structure file based on the settings of the fourth input parameter file. After the calculation is completed, the energies of the reactants and products after geometry optimization can be obtained. For example, for the existing chemical equation (2), the calculation of the energies of the reactants and products can be seen in Table 3. Of course, in other possible implementation manners, other methods can also be used to calculate the energies of the reactants and products, and the present application does not make specific limitations on this.

[0074] Table 3 Total energy of Mg-MOF-74 and C 8 H 6 O 6 , energy table of MgO and H2O

[0075] After step S1042 is completed, step S1043 is executed. In step S1043, the product of the compound energy of the target MOF material multiplied by its stoichiometric number in the existing chemical equation is added to the sum of the products of the energies of each product, and then the sum of the products of the energies of each reactant is subtracted from this sum to obtain the formation energy. The product of the energy of the product is the product of the energy of the product and its stoichiometric number in the existing chemical equation, and the product of the energy of the reactant is the product of the energy of the reactant and its stoichiometric number in the existing chemical equation.

[0076] Continuing with the aforementioned example, for Mg-MOF-74, referring to Table 3 and the existing chemical equation (2), its formation energy can be calculated by the formula E = -8688.13 + (-110.45)×18 - (-99.3)×18 - (-975.17)×9.

[0077] After obtaining the formation energy, step S105 can be entered. In step S105, a second input parameter file is made for the calculation of the density of states of the target MOFs material model. Optionally, in step S105, making the second input parameter file for the density of states calculation of the target MOFs material model includes: S1051, create a second initial file in hsd format.

[0078] S1052, write the Geometry module in the second initial file, and in this module, set the object for the density of states calculation to the second structure file.

[0079] S1053, write the Hamiltonian module in the second initial file, and in this module, set to turn on the self-consistent charge calculation with a convergence value of 1e-5 for this calculation, set the maximum angular momentum quantum number for each element of the target MOFs material model, and set the k-point sampling density in the Brillouin zone to an 8×8×8 grid.

[0080] S1054, write the Analysis module in the second initial file, and in this module, set the projection region for the density of states calculation of each element of the target MOFs material model.

[0081] S1055, write the ParserOptions module in the second initial file, and in this module, set the parser version, and finally obtain the second input parameter file for the density of states calculation of the target MOFs material model.

[0082] In step S1051, create a second initial file in hsd format, the initial content of which can be empty, and then corresponding parameter settings can be written. After step S1051 is executed, steps S1052 to S1055 are entered. Steps S1052 to S1055 are all to write corresponding module settings in the second initial file. The present application does not limit the execution order thereof. After all modules are written, the second input parameter file can be obtained.

[0083] In step S1052, write the Geometry module in the second initial file. By this module, set the object for the density of states calculation to the second structure file. For example, the file name of the second structure file in gen format can be directly introduced to specify that the object for the density of states calculation is the second structure file, or the structure parameters in the second structure file can be copied to the Geometry module to complete the setting of the object for the density of states calculation.

[0084] In step S1053, the Hamiltonian module is written into the second initial file. The self-consistent charge calculation can be enabled in this module, and the convergence value of this calculation is set to 1e-5. In this way, the charge density will be automatically updated in each calculation until the change in the charge density in each calculation is less than this value, and then the self-consistent calculation is considered to have converged. In a possible implementation, this module may include two sub-modules, MaxAngularMomentum and KPointsAndWeights. The maximum angular momentum quantum number of each constituent element of the target MOFs material can be set through the MaxAngularMomentum sub-module, which is achieved by defining the orbital type adopted for each element. Continuing with the previous example, for Mg-MOF-74, it can be set that the magnesium element (Mg) adopts the p orbital, the hydrogen element (H) adopts the s orbital, the carbon element (C) adopts the p orbital, and the oxygen element (O) adopts the p orbital. In addition, the 8×8×8 grid with a higher k-point density can be set as the k-point sampling density of the Brillouin zone through the KPointsAndWeights sub-module. For example, the KPointsAndWeights sub-module can be set as follows: KPointsAndWeights = SupercellFolding{ 8 0 0 0 8 0 0 0 8 0.5 0.5 0.5 } In step S1054, the Analysis module can be continuously written into the second initial file to create a projection region for each element (such as Mg, H, C, and O in Mg-MOF-74), and the density of states calculation for different elements is performed within this region. Specifically, the ProjectStates sub-module can be set in this module, and multiple Region sub-sub-modules are written in this sub-module to create projection regions for each element.

[0085] In step S1055, the ParserOptions module can be continuously written into the second initial file. For example, the parser version can be set through the command "ParserVersion = 12". After all the above settings are written into the second initial file, the second input parameter file for the density of states calculation of the target MOFs material model is finally obtained.

[0086] After step S105 is completed, step S106 can be entered. In step S106, the density of states calculation command can be run in the DFTB+ software. The DFTB+ software running command reads the second structure file and the second input parameter file, and performs the density of states calculation on the target MOFs material model of the second structure file based on the settings of the second input parameter file. After the calculation is completed, the density of states calculation results can be visualized to obtain the density of states diagram. For example, a visualization tool such as python can be used to draw the density of states diagram. See Figure 3 , taking Mg-MOF-74 as an example, the density of states diagram can show the total density of states (TDOS) and the partial density of states (PDOS) under the equilibrium structure of Mg-MOF-74. The total density of states is used to represent the distribution of all electron states of the target MOFs material such as Mg-MOF-74 in the energy space, and the partial density of states is used to represent the contribution of elements of specific atomic orbitals to the total density of states.

[0087] The density of states, or the density of electronic states (DOS), describes the distribution of electronic states in the energy space. The density of states can give information such as the interaction between local partial wave orbitals, the energy level shift and energy level dispersion of electronic states, and can be directly used to analyze doping effects, etc. In the calculation, the density of states can intuitively give the electron energy spectrum diagram obtained according to the k vector along the high-symmetry direction in the Brillouin zone, giving us various information we want. Therefore, its status is very important. The density of states (DOS) can represent the ratio of the number of electrons at the corresponding energy. Having a peak at a certain energy level means that there are more electrons distributed here. The abscissa of the density of states diagram is energy, and the ordinate represents the relative quantity, which is a relative value. It is very helpful for revealing the interaction between atomic orbitals and the energy level shift information of electrons.

[0088] After step S106 is completed, step S107 can be entered. In step S107, the target MOFs material can be analyzed in terms of corresponding material structure characteristics, stability, and electronic behavior based on the unit cell parameters, formation energy, and density of states diagram obtained from the previous steps. Optionally, in step S107, the analysis of the target MOFs material based on the unit cell parameters, formation energy, and density of states diagram includes: S1071, analyzing the unit cell shape and size of the target MOFs material based on the unit cell parameters.

[0089] S1072, when the target MOFs material is a newly designed structure, if the formation energy is greater than 0, it is determined that the target MOFs material cannot be formed by spontaneous reaction.

[0090] S1073, when the target MOFs material is an existing structure, determining one or more existing chemical equations with the minimum formation energy for experimental verification.

[0091] S1074. Perform the total density of states analysis and the partial density of states analysis on the target MOFs material based on the density of states diagram.

[0092] In step S1071, the unit cell parameters (including the unit cell edge lengths a, b, c and the unit cell angles α, β, γ) are the basic geometric parameters characterizing the crystal structure. Through them, not only can the shape, size and volume of the unit cell be accurately determined, but also the seven crystal systems (such as cubic, hexagonal, orthorhombic, etc.) to which the crystal belongs can be judged based on their numerical relationships, and further the space group symmetry can be deduced. These parameters provide the core basis for material structure analysis and can be further used for calculations such as atomic packing density.

[0093] After step S1071 is completed, step S1072 or step S1073 can be executed. If the target MOFs material is a newly designed structure, step S1072 can be entered. In step S1072, the stability of the MOFs material with the newly designed structure is analyzed according to the formation energy. If its formation energy is greater than 0, it means that its stability is poor and it cannot be formed through spontaneous reaction, and additional driving forces (such as high temperature, high pressure or catalyst) are needed, which makes its synthesis more complex and consumes more. In some cases, MOFs structures with similar performance but lower formation energy may be screened, and further research on MOFs materials with formation energy greater than 0 may be abandoned. If the formation energy is less than 0, the stability of the target MOFs material is good. If its performance in other aspects also meets the requirements, it can be included in further research, such as research in the experimental stage.

[0094] If the target MOFs material is an existing structure, step S1073 can be entered. In step S1073, there may be one or more synthesis routes for the target MOFs material. When there are multiple synthesis routes, the formation energies corresponding to multiple synthesis routes can be calculated according to the previous steps, and then one or more existing chemical equations corresponding to the synthesis routes with the minimum formation energy can be selected and further verified and analyzed in experiments. If there is only one synthesis route, the existing chemical equation corresponding to this synthesis route can be directly used for experimental verification and analysis.

[0095] After step S1072 or step S1073 is completed, step S1074 can be entered. In step S1074, the total density of states analysis and the partial density of states analysis can be performed according to the density of states diagram. See Figure 3 , whose abscissa is energy and ordinate is the density of states DOS. Taking Mg-MOF-74 as an example, its total density of states (TDOS) can be seen in Figure 3 the fluctuating curve data corresponding to DOS in, and the partial density of states (PDOS) can be obtained through Figure 3Characterize the fluctuation curve regions of the respective orbitals of C, H, Mg, and O elements corresponding to the total density of states, and show the specific contributions of different atomic orbitals to the total density of states.

[0096] See Figure 3 , in Mg-MOF-74, there is no obvious density of states (DOS) crossing the Fermi level, indicating that Mg-MOF-74 is a semiconductor. The partial density of states (PDOS) of Mg-MOF-74 can characterize the distribution of various electronic states in the valence band. The p-states and s-states of Mg, O, and C atoms play a major role in the valence band (VB), respectively. At the same time, there are three internal gaps in Mg-MOF-74: upper (0 - 6 eV), (8 - 12 eV), and lower (-6 - -12 eV). Studies have shown that Mg atoms make an important contribution to the electronic states above the higher gap (>8 eV), O atoms make a large contribution below the lower internal gap (<-6 eV), and near the Fermi level (E = 0 eV), considering the valence band, C atoms dominate the contribution to the electronic states.

[0097] Through the above technical solutions, for the target MOFs material model to be analyzed, obtain its first structure file and create its first input parameter file. Read the first structure file and the first input parameter file through DFTB+ for geometric optimization simulation to optimize the geometry of the target MOFs material model, and obtain the compound energy and unit cell parameters. For different situations of the target MOFs material, if it is a newly designed structure, calculate the formation energy based on the theoretical chemical equation for the synthesis of the constituent element simple substances, the energy of the simple substances, and the compound energy; if it is an existing structure, calculate the formation energy based on the existing chemical equation, the energies of the reactants and products. Finally, create the second input parameter file, and use DFTB+ to read this file and the second structure file of the geometrically optimized target MOFs material model for density of states calculation and visualization to obtain the density of states diagram. After the operation is completed, the target MOFs material can be analyzed based on the calculated unit cell parameters, formation energy, and density of states diagram. Thus, the simulation calculation and analysis of MOFs materials are realized based on the open-source and free software DFTB+, which can provide more data support for screening MOFs materials, improve the screening efficiency, and reduce the time and cost of screening.

[0098] Molecular simulation technology provides atomic-level thermodynamic and kinetic data for the study of MOF materials, which are difficult to obtain through experiments. Through simulation, researchers can obtain more accurate structural, energy, and performance information in highly symmetric and complex systems, thus making up for the limitations of experimental research. In addition, molecular simulation has significant advantages in large-scale material screening. It can complete a large number of computational screenings in a short time, greatly improving the research efficiency and reducing costs. Compared with traditional experimental screening methods, molecular simulation can not only verify experimental phenomena but also provide a theoretical basis for the phenomena, providing support for the performance optimization of MOFs and the design of new materials. The results of molecular simulation can be a powerful supplement to experimental data. By comparing with experimental results, it can verify experimental phenomena and provide a theoretical basis for experimental results. Especially when it is difficult to obtain some complex data through experiments, DFTB+ can provide a feasible computational approach to help explain the experimentally observed phenomena.

[0099] Computational simulation, especially DFTB+, provides powerful technical support for the structural optimization and performance regulation of MOF materials. The high crystallinity and symmetry of MOFs enable the accurate characterization of atomic positions, reducing computational uncertainty and thus improving the accuracy of simulation. With the help of advanced computational hardware and software, simulation methods can optimize structures and properties based on low-resolution experimental data and obtain more reliable parameters. This not only enhances the understanding of MOF materials but also provides a solid theoretical support for their applications in the fields of energy storage, catalysis, adsorption, etc. Since DFTB+ is an approximate method based on density functional theory, it can significantly reduce the computational cost while maintaining a relatively high computational accuracy. This advantage is particularly important for the screening and study of a large number of different MOF structures, enabling efficient computational exploration without sacrificing accuracy.

[0100] Figure 4 is a block diagram of an electronic device shown according to an exemplary embodiment. As Figure 4 shown, the electronic device 700 may include: a processor 701, a memory 702. The electronic device 700 may also include one or more of a multimedia component 703, an input / output (I / O) interface 704, and a communication component 705.

[0101] Among them, the processor 701 is used to control the overall operation of the electronic device 700 to complete all or part of the steps in the above-mentioned MOFs material structure performance analysis method. The memory 702 is used to store various types of data to support the operation of the electronic device 700. These data may include, for example, instructions for any application or method operating on the electronic device 700, as well as application-related data, such as contact data, sent and received messages, pictures, audio, video, and so on. The memory 702 can be implemented by any type of volatile or non-volatile storage device or a combination thereof, such as Static Random Access Memory (SRAM), Electrically Erasable Programmable Read-Only Memory (EEPROM), Erasable Programmable Read-Only Memory (EPROM), Programmable Read-Only Memory (PROM), Read-Only Memory (ROM), magnetic memory, flash memory, magnetic disk or optical disk. The multimedia component 703 may include a screen and an audio component. Among them, the screen may be a touch screen, and the audio component is used to output and / or input audio signals. For example, the audio component may include a microphone for receiving external audio signals. The received audio signal may be further stored in the memory 702 or sent through the communication component 705. The audio component also includes at least one speaker for outputting audio signals. The I / O interface 704 provides an interface between the processor 701 and other interface modules, and the above other interface modules may be a keyboard, a mouse, buttons, etc. These buttons may be virtual buttons or physical buttons. The communication component 705 is used for wired or wireless communication between the electronic device 700 and other devices. Wireless communication, such as Wi-Fi, Bluetooth, Near Field Communication (NFC), 2G, 3G, 4G, NB-IOT, eMTC, or other 5G, etc., or a combination of one or more of them, is not limited here. Therefore, the corresponding communication component 705 may include: a Wi-Fi module, a Bluetooth module, an NFC module, and so on.

[0102] In an exemplary embodiment, the electronic device 700 can be implemented by one or more application specific integrated circuits (ASICs), digital signal processors (DSPs), digital signal processing devices (DSPDs), programmable logic devices (PLDs), field programmable gate arrays (FPGAs), controllers, microcontrollers, microprocessors, or other electronic components, and is used to execute the above-mentioned MOFs material structure performance analysis method.

[0103] In another exemplary embodiment, a computer-readable storage medium including program instructions is further provided. When the program instructions are executed by a processor, the steps of the above-mentioned MOFs material structure performance analysis method are implemented. For example, the computer-readable storage medium can be the above-mentioned memory 702 including program instructions, and the above program instructions can be executed by the processor 701 of the electronic device 700 to complete the above-mentioned MOFs material structure performance analysis method.

[0104] The preferred embodiments of the present application have been described in detail above in conjunction with the accompanying drawings. However, the present application is not limited to the specific details in the above embodiments. Within the scope of the technical concept of the present application, various simple modifications can be made to the technical solutions of the present application, and these simple modifications all belong to the protection scope of the present application.

[0105] In addition, it should be noted that, in the above specific embodiments, the various specific technical features described can be combined in any suitable manner without conflict. To avoid unnecessary repetition, the present application does not separately describe various possible combination methods.

[0106] In addition, any combination can be made between different embodiments of the present application as long as it does not violate the idea of the present application, and it should also be regarded as the content recorded in the present application.

Claims

1. A method for analyzing the structure and performance of MOFs materials, characterized in that: The method comprises: Obtaining a first structure file of a target MOFs material model, and preparing a first input parameter file for geometry optimization of the target MOFs material model; Run the geometry optimization command in the DFTB+ software to read the first structure file and the first input parameter file and perform geometry optimization simulation to obtain the second structure file, compound energy and unit cell parameters of the target MOFs material model after geometry optimization; If the target MOFs material is a newly designed structure, the formation energy of the target MOFs material is calculated based on the theoretical chemical equation for synthesizing the target MOFs material from the constituent elements of the target MOFs material, the energy of the single substance, and the energy of the compound; If the target MOFs material is an existing structure, the formation energy of the target MOFs material is calculated based on an existing chemical equation for synthesizing the target MOFs material, and the energy of reactants and products of the existing chemical equation, wherein the energy of the product includes the energy of the compound; Prepare a second input parameter file for state density calculation of the target MOFs material model; Run the state density calculation command in the DFTB+ software to read the second structure file and the second input parameter file and perform state density calculation, and visualize the state density calculation results to obtain a state density diagram of the target MOFs material model after geometry optimization; The target MOFs materials were analyzed based on unit cell parameters, formation energy and state density diagrams.

2. The method according to claim 1, characterized in that The method of calculating the formation energy of the target MOFs material based on the theoretical chemical equation for synthesizing the target MOFs material using the constituent elements of the target MOFs material, the energy of the single substance and the energy of the compound comprises: Determine the single substance of all constituent elements of the target MOFs material based on the chemical formula of the target MOFs material; Determine the reaction formula for generating the target MOFs material using all the single substances and balance it to obtain a theoretical chemical equation; Calculate the energies of all said elements; The energy of the compound is multiplied by the product of the stoichiometric number of the target MOFs material in the theoretical chemical equation, and the product of the energy of each of the single substances and the stoichiometric number of the single substance in the theoretical chemical equation is subtracted to obtain the formation energy.

3. The method according to claim 1, characterized in that The calculation of the formation energy of the target MOFs material based on the existing chemical equation for synthesizing the target MOFs material, the energy of the reactants and the products of the existing chemical equation comprises: Obtain the existing chemical formula for synthesizing the target MOFs material; Calculate the energy of all reactants and all products except the target MOFs material in the existing chemical equation; The formation energy is obtained by multiplying the energy of the compound by the stoichiometric number of the target MOFs material in the existing chemical equation, adding the energy of each of the products except the target MOFs material and the stoichiometric number of the product in the existing chemical equation, and subtracting the energy of each reactant by the stoichiometric number of the reactant in the existing chemical equation.

4. The method according to claim 2, characterized in that: The calculating of the energy of all the single substances comprises: For each element, a third structure file of the element is obtained, and a third input parameter file for geometry optimization of the element is prepared; Run the geometry optimization command in the DFTB+ software to read the third structure file and the third input parameter file and perform geometry optimization simulation to obtain the energy of each element after geometry optimization.

5. The method according to claim 3, characterized in that: The energy of all reactants and all products except the target MOFs material in the existing chemical equation is calculated, including: For each reactant in the existing chemical equation and each product other than the target MOFs material, a fourth structure file of the reactant or product is obtained, and a fourth input parameter file for geometry optimization is prepared; The geometry optimization command is run in the DFTB+ software to read the fourth structure file and the fourth input parameter file and perform a geometry optimization simulation to obtain the energy of each of the reactants and each of the products after geometry optimization.

6. The method according to claim 1, characterized in that The first input parameter file for geometric optimization of the target MOFs material model comprises: Create the first initial file in hsd format; Writing a Geometry module into the first initial file, and setting the object of geometry optimization in the module as the first structure file; Write a Driver module in the first initial file, and set the calculation type to geometry optimization in the module, set full lattice optimization to be enabled and use rational function optimizer, set the maximum number of optimization steps to any number of steps between 100 and 200, and set the convergence criterion to the gradient element GradElem being less than 1e-4; Write a Hamiltonian module in the first initialization file, and set the self-consistent charge calculation to be enabled in the module, and the convergence value and maximum number of iterations of the calculation are 1e-6 and 2000 respectively, set the Slater-Koster parameters to use the 3OB parameter set, set the maximum angular momentum quantum number of each element of the target MOFs material model, set the Brillouin zone k-point sampling to use the Monkhorst-Pack method and select only the Γ point as the k-point sampling point; Write the Analysis module in the first initialization file and set the calculation properties in the module; The ParserOptions module is written into the first initial file, and the parser version is set in the module, and finally the first input parameter file for geometry optimization of the target MOFs material model is obtained.

7. The method according to claim 1, characterized in that The second input parameter file for calculating the state density of the target MOFs material model comprises: Create a second initial file in hsd format; Write a Geometry module into the second initial file, and set the object of state density calculation in the module to be the second structure file; Write the Hamiltonian module in the second initial file, and set the self-consistent charge calculation to be enabled in the module and the convergence value of the calculation to 1e-5, set the maximum angular momentum quantum number of each element of the target MOFs material model, and set the k-point sampling density of the Brillouin zone to a grid of 8×8×8; Writing an Analysis module into the second initial file, and setting a projection region for calculating the density of states of each element of the target MOFs material model in the module; The ParserOptions module is written into the second initial file, and the parser version is set in the module, and finally the second input parameter file for state density calculation of the target MOFs material model is obtained.

8. The method according to claim 1, characterized in that The target MOFs material is analyzed based on unit cell parameters, formation energy and state density diagram, including: Analyze the unit cell shape and size of the target MOFs material based on the unit cell parameters; When the target MOFs material is a newly designed structure, if the formation energy is greater than 0, it is determined that the target MOFs material cannot be formed by a spontaneous reaction; When the target MOFs material is an existing structure, one or more existing chemical equations with the minimum formation energy are determined for experimental verification; Based on the state density diagram, the total state density analysis and partial wave state density analysis are performed on the target MOFs material.

9. A computer-readable storage medium having a computer program stored thereon, characterized in that: When the program is executed by a processor, the steps of the method described in any one of claims 1 to 8 are implemented.

10. An electronic device, characterized in that: include: a memory having a computer program stored thereon; A processor, configured to execute the computer program in the memory to implement the steps of the method according to any one of claims 1 to 8.

Images