MXenes material performance prediction method based on graph neural network

Abstract

The invention belongs to the technical field of material performance prediction, and discloses an MXenes material performance prediction method based on a graph neural network, and the method comprises the steps: obtaining the structural feature parameters of a to-be-predicted MXenes material; respectively matching node derivative information, edge derivative information and global derivative information based on the structural feature parameters; based on the node derivative information, the edge derivative information and the global derivative information, constructing an initial feature vector set, based on edge feature vectors and based on global feature vectors; and performing hierarchical processing on each feature vector in the initial feature vector set through each layer of the graph neural network to obtain the performance index of the to-be-predicted MXenes material target. According to the scheme, the key technical problems that in the prior art, geometric information coding is insufficient, and intra-layer and inter-layer interaction cannot be effectively distinguished are effectively solved.

Description

Technical Field

[0001] This invention belongs to the field of material property prediction technology, specifically relating to a method for predicting material properties based on graph neural networks (MXenes). Background Technology

[0002] Traditional machine learning methods are only applicable to MXenes with single functional groups, and their effectiveness is limited for complex features; they cannot accurately predict MXenes with mixed surface terminations. Existing graph neural network models suffer from insufficient geometric information encoding and inability to effectively distinguish intra-layer and inter-layer interactions when dealing with the unique layered structure and complex interactions of MXenes, resulting in limited prediction accuracy and performance saturation.

[0003] In the prior art, there is a method of predicting the dynamic viscosity of MXene-soybean oil nanofluid by using long short-term memory network (LSTM) combined with Bayesian optimization technology. First, different concentrations of MXene nanofluid are prepared through experiments, and a rheometer is used to measure the dynamic viscosity in the temperature range of 25-85°C to obtain 200 sample data. Then, an LSTM deep learning model is constructed, with temperature and nanoplatelet concentration as input parameters, and the long-term dependence relationship in the sequence data is processed through the memory unit. The Bayesian optimization method is used to automatically adjust the hyperparameters of LSTM, and the network is trained through the back propagation algorithm, finally realizing the prediction of the dynamic viscosity of MXene nanofluid. However, this method is only for the prediction of a single performance parameter, and cannot simultaneously predict multiple key performances of the material, so the application range is limited. Using the traditional LSTM network structure, mainly processing time series data, cannot effectively represent the atomic-level spatial structure of MXene material and the complex interaction relationship between atoms, lacks deep understanding of the geometric structure information of the material; it is heavily dependent on a large amount of experimental data, and different ratio nanofluid samples need to be prepared and subjected to tedious rheological tests, which has high experimental cost and long cycle, and cannot realize high-throughput material screening; the input features are limited to macroscopic temperature and concentration parameters, ignoring the microstructure characteristics of MXene materials, resulting in insufficient understanding of the nature of the material by the model, and limited prediction accuracy. Further, in the prior art, there is also a method of using VICReg graph contrastive learning framework. The core technology is to generate enhanced graph pairs through 6 kinds of random graph enhancement technologies, and use VICReg loss function for contrastive learning training. After training, linear regression, random forest, LightGBM, XGBoost and other downstream models are used to predict the performance of embedded features, which can predict 7 kinds of MXenes performance such as work function, Fermi level density of states, Poisson ratio, bulk modulus, shear modulus, Young's modulus and magnetism. However, this method relies on random graph enhancement strategy to change the molecular graph structure, which may destroy the key structure information of MXenes material, especially for MXenes materials with fine layered structure, random deletion of nodes or edges will lose important intra-layer-inter-layer interaction information; it needs to be pre-trained unsupervisedly before using traditional machine learning models for downstream prediction, the training process is complex and the calculation efficiency is low; it is not specially designed for the unique two-dimensional layered structure characteristics of MXenes, and lacks the ability to distinguish different interaction modes of covalent bond and van der Waals force, which limits the further improvement of prediction accuracy.

[0004] Further, the patent CN120340696A discloses a high polymer material performance prediction method based on deep learning, which is mainly for predicting the aging performance of high polymer materials. The method first prepares a group of high polymer material samples with different variable parameter combinations, respectively performs initial tensile test and tensile test after accelerated aging test, and obtains tensile performance parameters such as breaking strength and elongation. Then two independent multilayer perceptron deep learning network models are constructed: the initial performance prediction model adopts the structure of 5 feature inputs→128 neurons in the first hidden layer→64 neurons in the second hidden layer→32 neurons in the third hidden layer→1 neuron output layer, to predict the initial tensile performance of the material; the aging performance prediction model adopts the structure of 7 feature inputs→128 neurons in the first hidden layer→64 neurons in the second hidden layer→32 neurons in the third hidden layer→1 neuron output layer, to predict the final tensile performance of the material. All hidden layers use ReLU activation function, and mean square error is used as the loss function for training. Finally, the comprehensive performance of the material is evaluated through the environmental harshness coefficient calculation and the performance evaluation index. This method is only applicable to the specific performance prediction of high polymer materials and cannot be extended to the multi-performance prediction of other types of materials such as two-dimensional materials MXenes. Using traditional multilayer perceptron network, it lacks effective representation of atomic-level structure information of materials and cannot capture the complex interactions and geometric spatial relationships between atoms. The environmental harshness coefficient calculation uses a simple linear weighting method and does not consider the nonlinear coupling effects between different environmental factors and the intrinsic characteristics of the material structure. This method is mainly based on experimental data and relies on a large number of physical test sample preparations, which is time-consuming and cannot meet the demand of high-throughput material screening and rapid performance prediction. SUMMARY

[0005] In view of the deficiencies of the prior art, the present application provides a MXenes material performance prediction method based on a graph neural network.

[0006] The method comprises: obtaining structure characteristic parameters of a MXenes material to be predicted, the structure characteristic parameters comprising atomic coordinates, atomic types, lattice parameters and edge characteristic information, the edge characteristic information comprising interatomic distances, angle information and chemical bond types; based on the structure characteristic parameters, respectively matching node-derived information, edge-derived information and global-derived information; based on the node-derived information, the edge-derived information and the global-derived information, constructing an initial feature vector set, wherein the initial feature vector set comprises a node feature vector constructed based on the node-derived information, an edge feature vector constructed based on the edge-derived information and a global feature vector constructed based on the global-derived information; The feature vectors in the initial feature vector set are processed in layers by layers of the graph neural network to obtain the target performance index of the MXenes material to be predicted.

[0007] In some embodiments, the obtaining of the structural characteristic parameter of the MXenes material to be predicted comprises: obtaining the MXenes material to be predicted from a database; performing structural file analysis on the MXenes material to be predicted by the density functional theory, calculating performance data, and performing secondary verification on the performance data measured by a model to obtain the structural characteristic parameter.

[0008] In some embodiments, the matching of the node-derived information, the edge-derived information and the global-derived information respectively based on the structural characteristic parameter comprises: matching atomic attribute information of each atom based on the atomic types to obtain the node-derived information, the atomic attribute information at least including atomic number, valence and atomic radius; encoding interaction features between atoms based on the edge feature information to obtain edge-derived information for describing atom pair interaction information; determining a target extraction dimension of the global-derived information; extracting each dimension information matching the target extraction dimension based on the atomic coordinates and the atomic types to obtain the global-derived information.

[0009] In some embodiments, the constructing of the initial node feature vector based on the node-derived information comprises: inputting the node-derived information into a multi-layer perception machine based on a mapping function to map the node-derived information into a high-dimensional hidden space vector to obtain the initial node feature vector represented based on the high-dimensional hidden space vector; wherein the mapping function is: (1); wherein, is a node feature vector of atomic node i, ELU is an exponential linear unit activation function, is an original feature vector of atomic node i, is an embedding weight matrix, is a bias term.

[0010] In some embodiments, the constructing of the initial edge feature vector based on the edge-derived information comprises: combining and encoding the edge-derived information based on a Gaussian radial basis function and a spherical harmonic function to generate an edge feature vector representing atom pair interaction, wherein the edge feature vector is represented as: (2); wherein, is the eigenvector of edge (i, j), RBF is a Gaussian radial basis function, and SBF is a spherical harmonic function, is the distance encoding, which refers to the Euclidean distance between atom node i and atom node j: (3); is the angle encoding, which is a key for the MXenes layered structure: (4).

[0011] In some embodiments, constructing a global feature vector based on the global derived information comprises: aggregating and mapping the global derived information through a multi-layer perception mechanism or a nonlinear change structure to obtain the global feature vector, wherein the global feature vector is represented as follows: (5); wherein, represents a global feature vector of the MXenes material to be predicted, is an aggregation function, is a node feature vector of the i-th atom node.

[0012] In some embodiments, the layered processing of each feature vector in the initial feature vector set through each layer of the graph neural network to obtain the target performance indicator of the MXenes material to be predicted comprises: processing the input node feature vector and the edge feature vector through a multi-head attention graph convolution layer of the graph neural network combined with an attention weight to obtain a graph convolution feature vector of the MXenes material to be predicted, and inputting the graph convolution feature vector to a global attention layer of the graph neural network; the global attention layer of the graph neural network receives the graph convolution feature vector and performs correlation analysis to generate a global attention enhanced feature, processes the global attention enhanced feature through a feature reweighting mechanism, obtains a global correlation feature vector of the MXenes material to be predicted, and inputs the global correlation feature vector to a fusion prediction layer of the graph neural network; the fusion prediction layer of the graph neural network acquires the global feature vector and receives the global correlation feature vector input by the global attention layer, performs fusion prediction on the global feature vector and the global correlation feature vector, and outputs the predicted MXenes material target performance indicator.

[0013] In some embodiments, the attention weight is determined based on the following formula: (6); wherein, is the attention weight of atom node i to atom node j in the kth attention head, , are query and key vectors, respectively, is the attention head dimension, is the attention weight matrix, is the neighbor set of atom node i.

[0014] In some embodiments, the processing of the global attention enhanced feature through the feature reweighting mechanism to obtain the global correlation feature vector of the to-be-predicted MXenes material comprises: combining the multi-head attention function, reweighting the covalent bond interaction of atoms in the same layer; combining the cross-layer attention function, reweighting the van der Waals force interaction of atoms across layers; wherein, the multi-head attention function is: (7); the cross-layer attention function is: (8); wherein, is the intra-layer attention feature of the lth layer, is the inter-layer attention feature of the lth layer, is the multi-head attention function, is the cross-layer attention function.

[0015] In some embodiments, the global feature vector and the global correlation feature vector are fused and predicted to output the target performance index of the predicted MXenes material, comprising: using attention weighted pooling and a fully connected network layer, the global correlation feature vector is nonlinearly transformed through a multi-layer MLP to obtain a global transformed feature vector; the global transformed feature vector and the global feature vector are fused and predicted to output the target performance index; wherein, the fully connected network layer is according to the following formula: (9); (10); (11); wherein, is the output of the first layer, r is the fusion feature vector, is the output of the lth layer, y is the final prediction value, and are the weight matrix and bias term of the 0th layer, respectively, and are the weight matrix and bias term of the first layer respectively, and are the weight matrix and bias term of the output layer.

[0016] The present application includes but is not limited to the following advantages: (1) The present scheme combines the graph neural network technology with the unique layered structure characteristics of MXenes materials, effectively solves the key technical problems such as insufficient geometric information coding and inability to effectively distinguish intra-layer and inter-layer interactions in the prior art through the combination of Gaussian radial basis function and spherical harmonic function coding, multi-level feature fusion mechanism and other technical means; (2) The present scheme uses phased input and hierarchical processing, first lets the graph convolution layer focus on learning the local topological features of node features and edge features, avoids the interference of global features, and then gradually integrates global features through subsequent layers, which not only retains the fineness of local features, but also supplements the integrity of global information, and finally realizes the collaborative modeling of local interaction and global structure attributes through the fusion layer; (3) The present application is specially designed for the unique two-dimensional layered structure characteristics of MXenes materials, and the geometric information of the layered structure is specially coded through the angle coding formula, which can accurately capture the differences between the intra-layer covalent bond interaction and the inter-layer van der Waals force interaction of MXenes materials, and solves the technical problem that traditional methods cannot effectively represent the layered structure of MXenes; (4) Compared with the traditional density functional theory calculation method which requires a large amount of computing resources and time cost, the present application can realize fast prediction based on the trained graph neural network model, greatly reducing the computing cost and improving the material screening efficiency. BRIEF DESCRIPTION OF DRAWINGS

[0017] In order to more clearly illustrate the technical solutions in the embodiments of the present application or the prior art, the drawings needed to be used in the embodiments or prior art description will be briefly introduced as follows.

[0018] Figure 1 is a flowchart of the MXenes material performance prediction method based on the graph neural network of the present application embodiment; Figure 2 is a flowchart of obtaining the structure characteristic parameters of the MXenes material to be predicted of the present application embodiment; Figure 3 is a flowchart of matching and deriving information based on the structure characteristic parameters of the present application embodiment; Figure 4 is a flowchart of determining the target performance index of the predicted MXenes material of the present application embodiment; Figure 5 is a prediction result verification diagram of the MXenes material performance prediction method based on the graph neural network of the present application embodiment. DETAILED DESCRIPTION

[0019] The embodiment of the present application provides a kind of based on the performance prediction method of MXenes material of graph neural network, method includes: the structural characteristic parameter of the MXenes material to be predicted is acquired, structural characteristic parameter includes atomic coordinates, atomic type, lattice parameter and edge feature information, edge feature information includes interatomic distance, angle information and chemical bond type;Based on structural characteristic parameter, respectively match node derivative information, edge derivative information and global derivative information;Based on node derivative information, edge derivative information and global derivative information, construct initial feature vector set, wherein, initial feature vector set includes the node feature vector based on node derivative information, the edge feature vector based on edge derivative information and the global feature vector based on global derivative information;Each feature vector in initial feature vector set is handled by layer by each layer of graph neural network, to obtain the target performance index of the MXenes material to be predicted.This scheme combines the graph neural network technology with the layer structure characteristics unique to MXenes material, and solves the key technical problems such as insufficient geometry information coding and inability to effectively distinguish intra-layer and inter-layer interaction in the prior art by means of combination coding of Gaussian radial basis function and spherical harmonic function, multi-level feature fusion mechanism and other technical means.

[0020] The terms "first", "second", "third", "fourth" and the like in the description and claims of the present application and in the above DETAILED DESCRIPTION, if any, are used for distinguishing between similar objects and not necessarily for describing a particular sequential or chronological order. It is to be understood that the use of these terms in the description and claims of the present application is solely for distinguishing between the comparable objects and not necessarily for describing a particular sequential or chronological order. It is to be understood that the data so used in the description and claims is interchangeable under appropriate circumstances, such that the embodiments of the application described herein are capable of accomplishing functionalities using either the described data sequences or alternative data sequences that do not necessarily accomplish exactly the same functionalities in exactly the same manner. Additionally, the term "comprising" or "having" and variations thereof as used in the detailed description and claims of the present application are used generically and are intended to cover the recited elements, steps or components, among others, without excluding or excluding the presence of additional elements, steps, components, among others.

[0021] For the sake of understanding, the specific flow of the embodiment of the present application is described below, specifically, Figure 1 The flow chart of the MXenes material performance prediction method based on graph neural network of the embodiment of the present application is shown in Figure 1 , comprising the following steps: S100, the structural characteristic parameter of the MXenes material to be predicted is acquired.

[0022] Among them, structural characteristic parameter includes atomic coordinates, atomic type, lattice parameter and edge feature information, edge feature information includes interatomic distance, angle information and chemical bond type.

[0023] Specifically,Figure 2 As another flowchart of the MXenes material performance prediction method based on the graph neural network of the embodiments of the present application, the flowchart is an exemplary illustration of step S100, see 2, comprising the following steps: S200, obtaining a to-be-predicted MXenes material from a database.

[0024] S202, performing structural file analysis on the to-be-predicted MXenes material by density functional theory, calculating performance data, and performing secondary verification with the performance data measured by the model to obtain structural characteristic parameters.

[0025] In an example, 23,857 MXenes structure samples can be obtained from the MXenes Database database (website https: / / anant.mrc.iisc.ac.in / apps / mxene), each sample containing structural parameters and performance parameters calculated by density functional theory (DFT).

[0026] The specific operation is as follows: a) Analyze the structural file to obtain atomic coordinates, atomic types and lattice parameters; b) Represent each MXenes crystal structure as G = (V, E), where node V represents an atom and edge E represents an interaction relationship between atoms with a distance less than 8 Å; c) Obtain performance parameters: magnetic moment (μB), band gap (eV) and lattice constant (Å); d) Divide the data set into a training set of 19,085, a validation set of 2,386 and a test set of 2,386 in a ratio of 8:1:1.

[0027] S102, based on the structural characteristic parameters, respectively matching node derived information, edge derived information and global derived information.

[0028] Specifically, Figure 3 As another flowchart of the MXenes material performance prediction method based on the graph neural network of the embodiments of the present application, the flowchart is an exemplary illustration of step S102, see 3, comprising the following steps: S300, based on the atomic type, matching the atomic attribute information of each atom to obtain the node derived information, and the atomic attribute information at least includes the atomic number, the valence and the atomic radius.

[0029] S302, based on the edge feature information, encoding the interaction features between atoms to obtain the edge derived information for describing the interaction information of the atom pair.

[0030] S304, determining the target extraction dimension of the global derived information.

[0031] S306. Based on atomic coordinates and atomic type, extract information of each dimension that matches the target extraction dimension to obtain global derived information.

[0032] S104. Construct an initial feature vector set based on node-derived information, edge-derived information, and global-derived information.

[0033] The initial feature vector set includes node feature vectors constructed based on node-derived information, edge feature vectors constructed based on edge-derived information, and global feature vectors constructed based on global-derived information.

[0034] Specifically, constructing the initial node feature vector based on node-derived information includes: Based on the mapping function, the node-derived information is input into the multilayer perceptron to map the node-derived information into a high-dimensional hidden space vector, thereby obtaining the initial node feature vector based on the high-dimensional hidden space vector representation. The mapping function is: (1); In the formula, Let be the node feature vector of atomic node i, and ELU be the exponential linear unit activation function. Let be the original feature vector of atomic node i. To embed the weight matrix, This is a bias term.

[0035] In one example, the 14 original features of the atom, such as atomic number, chemical valence, and atomic radius, are mapped to a 300-dimensional hidden space through a multilayer perceptron, and the node feature vector is calculated using formula (1).

[0036] Constructing the initial edge feature vector based on edge-derived information includes: Edge-derived information is combined and encoded using Gaussian radial basis functions and spherical harmonic functions to generate edge feature vectors characterizing the interactions between atomic pairs. These edge feature vectors are represented as follows: (2); In the formula, Let be the eigenvector of edge (i,j), RBF be the Gaussian radial basis function, and SBF be the spherical harmonic function. It is distance encoding, referring to the Euclidean distance between atomic node i and atomic node j: (3); It's angular encoding, which is key for the layered structure of MXenes: (4).

[0037] In an example, 64 Gaussian radial basis functions are adopted; angle coding: the interatomic angle information is coded by using spherical harmonic function, and the edge feature vector is obtained by using formulas (2), (3) and (4) for calculation, especially for the MXenes layered structure design.

[0038] The global feature vector is constructed based on the global derived information, including: The global derived information is aggregated and mapped by a multi-layer perception mechanism or a nonlinear change structure to obtain the global feature vector, wherein the global feature vector is represented as follows: (5) ; In the formula, the global feature vector of the MXenes material to be predicted, is an aggregation function, is the node feature vector of the i th atomic node.

[0039] In an example, c) 6-dimensional global features such as total number of atoms, average atomic mass, and total number of electrons of the material are extracted, and the global feature vector is obtained by using formula (5) for calculation.

[0040] S106, each feature vector in the initial feature vector set is processed by each layer of the graph neural network to obtain the target performance index of the MXenes material to be predicted.

[0041] Specifically, Figure 4 is another flowchart of the MXenes material performance prediction method based on the graph neural network of the embodiments of the present application, which is a specific description of step S106, and the specific description of step S106 is described with reference to Figure 4 , including the following steps: S400, the node feature vector and the edge feature vector input are processed by the multi-head attention graph convolution layer of the graph neural network combined with the attention weight to obtain the graph convolution feature vector of the MXenes material to be predicted, and then input to the global attention layer of the graph neural network.

[0042] In an example, a multi-head attention graph convolution network containing 20 graph convolution layers is constructed herein, each layer uses 8 attention heads, and the attention head dimension is 37.5 (300 / 8): a) learning the interatomic interaction by the multi-head attention mechanism, and calculating the attention weight by using formula (6); b) adding residual connection and layer normalization after each graph convolution layer to prevent gradient disappearance. Specifically, the attention weight is determined based on the following formula: (6) ; In the formula, is the attention weight of the atom node i to the atom node j in the kth attention head, , are the query and key vectors respectively, is the attention head dimension, is the attention weight matrix, is the neighbor set of the atom node i.

[0043] S402, the global attention layer of the graph neural network receives the graph convolution feature vector and performs correlation analysis to generate a global attention enhanced feature. The global attention enhanced feature is processed through a feature reweighting mechanism to obtain a global correlation feature vector of the to-be-predicted MXenes material, which is input to the fusion prediction layer of the graph neural network.

[0044] Specifically, the covalent bond interaction of atoms in the same layer is reweighted by combining a multi-head attention function; The van der Waals force interaction of atoms across layers is reweighted by combining a cross-layer attention function; The multi-head attention function is: (7) ; The cross-layer attention function is: (8) ; In the formula, is the intra-layer attention feature of the lth layer, is the inter-layer attention feature of the lth layer, is the multi-head attention function, is the cross-layer attention function.

[0045] S404, the fusion prediction layer of the graph neural network acquires the global feature vector and the global correlation feature vector input by the global attention layer, and performs fusion prediction on the global feature vector and the global correlation feature vector to output a predicted target performance index of the MXenes material.

[0046] Specifically, the fusion prediction of the global feature vector and the global correlation feature vector to output the predicted target performance index of the MXenes material includes: The global correlation feature vector is transformed by a multi-layer MLP to obtain a global transformed feature vector; The global transformed feature vector and the global feature vector are fused and predicted to output the target performance index; In one example, the MLP includes three layers: a) The first layer: 512-dimensional input is mapped to 256-dimensional, using ELU activation function and 0.1 dropout; b) Second layer: 256 dimensions are mapped to 128 dimensions, using the ELU activation function and dropout of 0.1; c) Output layer: 128 dimensions are mapped to 1-dimensional predicted values, and the calculation is performed using formulas (9)-(11).

[0047] The fully connected network layer is based on the following formula: (9); (10); (11); In the formula, The first layer output is r, where r is the fused feature vector. Let y be the output of layer l, and y be the final predicted value. and These are the weight matrix and bias term for layer 0, respectively. and These are the weight matrix and bias term of the l-th layer, respectively. and These are the weight matrix and bias terms for the output layer.

[0048] Furthermore, in this example, the graph neural network is pre-trained, and the training parameters are set as follows: • Optimizer: Adam, initial learning rate 1×10⁻³, weight decay 0; • Learning rate scheduling: Reduced to 1×10⁻ after 50 epochs 4 ; • Loss function: Mean Absolute Error (MAE); • Number of training epochs: 100 epochs; • Early stopping strategy: Stop if the validation set loss does not improve for 10 consecutive epochs.

[0049] The prediction results on the test set of 2,386 samples are shown in Table 1: Table 1 Performance metrics MAE RMSE R² Magnetic moment 0.0807 μΒ 0.1123 μΒ 0.9513 Bandgap 0.0097 eV 0.0134 eV 0.9748 Lattice constant 0.0032 Å 0.0045 Å 0.9986 According to the evaluation criteria, the prediction is considered accurate when MAE ≤ 0.1 μB (magnetic moment), 0.02 eV (band gap), and 0.005 Å (lattice constant) and RMSE ≤ 0.15 μB (magnetic moment), 0.03 eV (band gap), and 0.008 Å (lattice constant). All indicators in this embodiment meet the accurate prediction criteria, and the verification graph of the prediction results is shown below. Figure 5 As shown in the figure, (a) is the result of the lattice constant; (b) is the result of the band gap; and (c) is the result of the magnetic moment.

[0050] The above examples are only used to illustrate the technical solutions of the present application, and are not intended to limit the present application; although the present application has been described in detail with reference to the foregoing examples, those skilled in the art should understand that the technical solutions recorded in the foregoing examples can be modified, or some technical features can be replaced by equivalent features; and these modifications or replacements do not cause the essence of the corresponding technical solutions to deviate from the spirit and scope of the technical solutions of the embodiments of the present application.

Claims

1. A method for predicting the properties of MXenes materials based on graph neural networks, characterized in that, The method includes: Obtain the structural feature parameters of the MXenes material to be predicted. The structural feature parameters include atomic coordinates, atomic type, lattice parameters, and edge feature information. The edge feature information includes interatomic distance, angle information, and chemical bond type. Based on the structural feature parameters, node-derived information, edge-derived information, and global-derived information are matched respectively. Based on the node-derived information, edge-derived information, and global-derived information, an initial feature vector set is constructed, wherein the initial feature vector set includes node feature vectors constructed based on the node-derived information, edge feature vectors constructed based on the edge-derived information, and global feature vectors constructed based on the global-derived information. The initial feature vector set is processed layer by layer by the graph neural network to obtain the target performance index of the MXenes material to be predicted.

2. The method for predicting MXenes material properties based on graph neural networks according to claim 1, characterized in that, The process of obtaining the structural feature parameters of the MXenes material to be predicted includes: Retrieve the MXenes materials to be predicted from the database; The structural files of the MXenes material to be predicted are parsed using density functional theory, performance data is calculated, and secondary verification is performed with the performance data measured by the model to obtain the structural characteristic parameters.

3. The method for predicting MXenes material properties based on graph neural networks according to claim 1, characterized in that, The matching of node-derived information, edge-derived information, and global-derived information based on the structural feature parameters includes: Based on the atom type, the atomic attribute information of each atom is matched to obtain the node-derived information. The atomic attribute information includes at least atomic number, chemical valence, and atomic radius. Based on the edge feature information, the interaction features between atoms are encoded to obtain edge-derived information used to describe the interaction information of atom pairs; Determine the target extraction dimension of the global derived information; Based on the atomic coordinates and the atomic type, information of each dimension that matches the target extraction dimension is extracted to obtain the global derived information.

4. The method for predicting MXenes material properties based on graph neural networks according to claim 1, characterized in that, Constructing the initial node feature vector based on the node-derived information includes: Based on the mapping function, the node-derived information is input into a multilayer perceptron to map the node-derived information into a high-dimensional hidden space vector, thereby obtaining the initial node feature vector represented by the high-dimensional hidden space vector. The mapping function is: （1）； In the formula, Let be the node feature vector of atomic node i, and ELU be the exponential linear unit activation function. The original feature vector of atomic node i, To embed the weight matrix, This is a bias term.

5. The method for predicting MXenes material properties based on graph neural networks according to claim 1, characterized in that, Constructing the initial edge feature vector based on the edge-derived information includes: The edge-derived information is combined and encoded using Gaussian radial basis functions and spherical harmonic functions to generate edge feature vectors characterizing the interactions between atomic pairs. These edge feature vectors are characterized as follows: （2）； In the formula, Let be the eigenvector of edge (i,j), RBF be the Gaussian radial basis function, and SBF be the spherical harmonic function. It is distance encoding, referring to the Euclidean distance between atomic node i and atomic node j: （3）； It's angular encoding, which is key for the layered structure of MXenes: （4）。 6. The method for predicting MXenes material properties based on graph neural networks according to claim 1, characterized in that, Constructing a global feature vector based on the aforementioned global derived information includes: The global derived information is aggregated and mapped using a multi-layer sensing mechanism or a non-linear transformation structure to obtain the global feature vector, wherein the global feature vector is characterized as follows: （5）； In the formula, This represents the global feature vector of the MXenes material to be predicted. It is an aggregate function. It is the node feature vector of the i-th atomic node.

7. The method for predicting MXenes material properties based on graph neural networks according to claim 1, characterized in that, The process of performing hierarchical processing on each feature vector in the initial feature vector set through each layer of the graph neural network to obtain the target performance index of the MXenes material to be predicted includes: The graph neural network uses a multi-head attention graph convolutional layer with attention weights to process the input node feature vector and edge feature vector, and then inputs the graph convolutional feature vector of the MXenes material to be predicted into the global attention layer of the graph neural network. The global attention layer of the graph neural network receives the graph convolution feature vector and performs correlation analysis to generate global attention enhancement features. The global attention enhancement features are then processed through a feature reweighting mechanism to obtain the global correlation feature vector of the MXenes material to be predicted, which is then input into the fusion prediction layer of the graph neural network. The fusion prediction layer of the graph neural network obtains the global feature vector and receives the global associated feature vector input from the global attention layer, performs fusion prediction on the global feature vector and the global associated feature vector, and outputs the predicted target performance index of the MXenes material.

8. The method for predicting MXenes material properties based on graph neural networks according to claim 7, characterized in that, Attention weights are determined based on the following formula: （6）； In the formula, Let be the attention weight of atomic node i to atomic node j in the k-th attention head. , and are the query and key vectors, respectively. For the attention head dimension, Here is the attention weight matrix. Let i be the set of neighbors of atomic node i.

9. The method for predicting MXenes material properties based on graph neural networks according to claim 7, characterized in that, The process of processing the global attention-enhanced features through a feature reweighting mechanism to obtain the global association feature vector of the MXenes material to be predicted includes: By incorporating a multi-head attention function, the covalent interactions of atoms in the same layer are reweighted; By combining the cross-layer attention function, the van der Waals interactions of cross-layer atoms are reweighted; The multi-head attention function is as follows: （7）； The cross-layer attention function is: （8）； In the formula, The intra-layer attention feature of the l-th layer, For the inter-layer attention features of layer l, For multi-head attention functions, This is a cross-layer attention function.

10. The method for predicting MXenes material properties based on graph neural networks according to claim 7, characterized in that, The global feature vector and the globally correlated feature vector are fused and predicted to output the predicted target performance index of the MXenes material, including: Attention-weighted pooling and fully connected network layers are used to nonlinearly transform the global correlation feature vector through a multi-layer MLP to obtain the global transformed feature vector. The global transformation feature vector and the global feature vector are fused and predicted to output the target performance index; The fully connected network layer is based on the following formula: （9）； （10）； （11）； In the formula, The first layer output is r, where r is the fused feature vector. Let y be the output of layer l, and y be the final predicted value. and These are the weight matrix and bias term for layer 0, respectively. and These are the weight matrix and bias term of the l-th layer, respectively. and These are the weight matrix and bias terms for the output layer.

Images