Invariant Tokenization of Crystalline Materials for Language Model Enabled Generation

We consider the problem of crystal materials generation using language models (LMs). A key step is to convert 3D crystal structures into 1D sequences to be processed by LMs. Prior studies used the crystallographic information framework (CIF) file stream, which fails to ensure SE(3) and periodic invariance and may not lead to unique sequence representations for a given crystal structure. Here, we propose a novel method, known as Mat2Seq, to tackle this challenge. Mat2Seq converts 3D crystal structures into 1D sequences and ensures that different mathematical descriptions of the same crystal are represented in a single unique sequence, thereby provably achieving SE(3) and periodic invariance. Experimental results show that, with language models, Mat2Seq achieves promising performance in crystal structure generation as compared with prior methods.

Paper

References (42)

Scroll for more · 30 remaining

Similar papers

Peer review

Reviewer vCYe6/10 · confidence 3/52024-07-09

Summary

This paper presents a novel approach for generating crystal materials using language models. The key innovation lies in the Mat2Seq method, which converts 3D crystal structures into 1D sequences while ensuring SE(3) and periodic invariance. This approach addresses the challenge of representing crystal structures in a way that is unique and invariant under different mathematical descriptions.

Strengths

1. The paper is well-written, providing clear descriptions of the problem background and proposed methods. 2. It tackles an interesting problem by using language models to generate crystal structures and introduces the consideration of SE(3) and periodic invariance for the first time. 3. The experiments are comprehensive, covering standard benchmarks and attempting to discover crystals with specific properties, demonstrating the robustness and versatility of the method. 4. Although the paper ensures the invariance of the generated structures, enforcing this invariance at the tokenizer level might limit the diversity of the generated structures to a small subspace. To obtain structures with different rotations, additional post-processing may be required, unlike models based directly in the spatial domain which can generate various rotational variants naturally.

Weaknesses

1. The paper should discuss the computational complexity differences between language models and specialized crystal models, as language models generally have higher parameter counts and computational demands. 2. The comparison of language models is limited. The experiments only compare with CrystalLLM, despite mentioning other methods. 3. The rationale for comparing only with CrystalLLM in Section 4.2 and not with other methods from Section 4.1 should be clarified.

Questions

1. For crystalline materials, isn't the original data already the most primitive unit? If we first apply Niggli cell reduction before using CIF files to describe unit cells, can we ensure a consistent representation? 2. While generating continuous crystal structures in 3D space seems natural, I am curious whether this token generation method can ensure that the reconstructed 3D crystal structure is continuous. Specifically, can the generated sequences always be converted back into a meaningful crystal representation?

Rating

6

Confidence

3

Soundness

3

Presentation

4

Contribution

3

Limitations

1. When performing conditional generation, it appears that generating structures for each specific property requires fine-tuning under specific conditions. This could result in high complexity in practical applications.

Reviewer Sr8M6/10 · confidence 5/52024-07-11

Summary

In this paper, the authors focus on the application of language models in the field of material generation. Starting from CIF files that represent crystal unit cell structures, they primarily utilize the Niggli reduction method to organize structures under different translations, rotations, and unit cell expansions into a unique representation. This representation is then used for crystal generation and conditional generation tasks. Compared to CrystaLLM, the proposed method demonstrates better generalization capabilities across multiple tasks.

Strengths

1. The authors introduce a new process, Mat2Seq, which uniquely represents crystal structures using Niggli reduction. This approach ensures SE(3) invariance, periodic invariance, and completeness, avoiding data augmentation problem and shortening token length in language model. 2. The paper showcases Mat2Seq’s ability to generalize to novel crystal structures that were not seen during training. This is a crucial aspect for practical applications and demonstrates the robustness of the method.

Weaknesses

1. The authors utilized a unique representation method for crystal structures, which theoretically promises to become an essential step in data preprocessing. This representation method primarily aims to improve the model's ability to recognize equivariant transformations in materials. However, the evaluation does not demonstrate this method's generalization capability to equivariant structures. It remains unclear from the authors' validation whether this method help language model learns from SE(3) equivariance or if the improved performance is due to the new token representation. 2. In the experimental validation in Section 4.2 and the conditional generation in Section 4.3, the evaluation of capabilities is very limited. The authors only provided the proportion of generated structures without detailing whether new structures can be generated, how many of the generated structures are valid, and accuracy of these generated structures. The generated structures need detailed explanation and evaluation.

Questions

1. In Section 3.1, the authors list two conditions for uniqueness and three types of transformations: translation, rotation, and unit cell expansion. It is suggested that Figure 1 should clearly illustrate all these transformations. Additionally, the "Change Lattice" operation in Figure 1 is confusing. Which type of transformation does this correspond to? Simply changing the length in the a-direction should not alter the atomic positions. Clarification is needed on whether the authors intended to illustrate unit cell expansion or some other operation. 2. In Table 1, the authors demonstrate Mat2Seq's ability to recognize crystal uniqueness, but the specific evaluation details are not provided. How were the values obtained, especially the 30% proportion? The methodology and criteria for this evaluation should be clearly explained. 3. In Section 4.3, while band gap is indeed an important property for semiconductors, it should be noted that whether the band gap is zero is a classification task distinguishing metals from non-metals. The authors' use of <0.5 as a threshold is inappropriate, especially when they mention that "values from 0 to 0.5 are marked as 0," which exacerbates the issue. This oversimplification could lead to misleading results. 4. In Section 4.2, in evaluating the ability to generate structures of new materials, the authors only indicate whether the generated structures match the given chemical formulas. It is recommended to include an evaluation of the atomic distance differences between the real and generated structures using RMSE to better assess the generation accuracy. 5. In Section 4.3, compared to commonly used machine learning models for material screening, generative models have the potential to produce out-of-distribution (OOD) structures, aiding researchers in discovering new structures with desired properties. However, the authors do not inform readers about the chemical formula repetitions, structure repetitions, or generation errors in the generated structures. A standard should be defined to evaluate the generated structures based on validity, uniqueness, and novelty.

Rating

6

Confidence

5

Soundness

3

Presentation

3

Contribution

2

Limitations

The authors utilize an existing unique unit cell representation method in the data preprocessing stage for language models. This method holds promise as a necessary step for future work. However, the evaluation does not provide sufficient evidence to demonstrate its application value. Additionally, the conditional generation applications mentioned by the authors are also quite limited, based on the evaluations provided in the paper.

Reviewer 5Kst4/10 · confidence 3/52024-07-12

Summary

There are several challenges when developing LMs for materials: 1) each crystal structure consists of an infinite number of atoms and a unique and invariant unit cell must therefore be selected for each crystal 2) the unit cell can be represented in a one-dimensional (1D) sequence that maintains invariance under arbitrary rotations and ensures completeness. To tackle these challenges, the paper proposed a method called Mat2Seq that systematically transforms 3D crystal structures into 1D sequences. This is achieved by first identifying SO(3) equivariant unit cells and subsequently converting these into SE(3) invariant sequences. The experimental results in crystal structure prediction and crystal discovery with desired properties validate the efficacy of Mat2Seq.

Strengths

- I understand why it is important for the inherent symmetry of materials to be reflected in the text when designing a language model for material generation. - The importance of unique representation and completeness from the structure of the material to the sequence is emphasized.

Weaknesses

- Overall, the novelty and contribution of the paper are low. Instead of proposing new algorithms, they simply apply technologies that are already widely used in materials science to the process of representing materials as 1d sequence text. More specifically, the method for determining the equivariant unit cell is derived from Niggli cell reduction rather than being an original contribution of this paper. This concept is already widely used in the field of materials science, and the authors’ application of it to the text representation of materials does not seem like a significant contribution. - It is unclear what Section 3.3 is trying to convey. The authors claim to achieve text representation reflecting SE(3) symmetry following SO(3) symmetry, but it is doubtful whether this is well conveyed in this section. - The evaluation metrics used are limited. To demonstrate efficiency more effectively, the paper should present values for various metrics such as Validity, Stability, and S.U.N., which are commonly addressed in other studies.

Questions

- Could you provide a more detailed explanation with examples regarding the input for the language model (LM)? - Why is the target property value inserted at the beginning of the text input when generating materials for a specific target property? - Does “irreducible atom sets” refer to the atom sets contained within the minimum unit cell?

Rating

4

Confidence

3

Soundness

2

Presentation

2

Contribution

2

Limitations

Please refer to the Weakness section.

Reviewer JJwt8/10 · confidence 3/52024-07-15

Summary

This article propose a novel method, known as Mat2Seq, to tackle this challenge. Mat2Seq to converts 3D crystal structures into 1D sequences and ensures that different mathematical descriptions of the same crystal are represented in a single unique sequence, thereby probably achieving SE(3) and periodic invariance. Experimental results show that, with language models, Mat2Seq achieves promising performance in crystal structure generation as compared with prior methods. Overall this work gives new insights and should be accepted for the conference after minor revisions.

Strengths

1. Mat2Seq to converts 3D crystal structures into 1D sequences 2. it reports the framework for creating unique and complete crystal sequence representations, followed by the construction of a material LLM capable of generating novel crystal structure with desired properties of interest.

Weaknesses

1. Lack of experimental data validation for the generated data.

Questions

1. Author needs to verify the generated data with more experimental results as currently it have multiple limitations. 2. Author needs to compare the accuracy with already existing CIF data of crystals.

Rating

8

Confidence

3

Soundness

3

Presentation

3

Contribution

3

Limitations

The limitations of the current Mat2Seq include: (1) it cannot be directly used for other atomic systems, like molecules and proteins; (2) the extension to model disordered materials remains a challenging frontier; and (3) large-scale training with more stable crystal structures can potentially enhance the robustness and performance when more computational resources are available.

Reviewer vCYe2024-08-10

The authors' rebuttal has addressed my concerns, and I have no further questions.

Authorsrebuttal2024-08-10

Thank you for your responses

Dear Reviewer vCYe, Thank you for your recognition and responses! Your suggestions undoubtedly help us enhance the clarity of this work, and we are glad that your concerns have been addressed. Yours sincerely, Authors

Reviewer Sr8M2024-08-12

Thank you to the authors for the detailed response and the additional experimental evaluations. These supplements have enriched the paper's experimental assessment, particularly the added evaluations on the validity, uniqueness, and novelty of the generated structures, which have improved the overall quality of the paper. As a result, I have decided to adjust my score to 6. However, there is still a point worth discussing. Regarding the authors' evaluation of conditional generation in the last reply, although the authors have explained their approach to merge materials with zero and near-zero band gaps, I still have some concerns. In the MP database, the ratio of materials with a zero band gap to those with a non-zero band gap is approximately 1:1, and in the MP20 dataset, the proportion of materials with a zero band gap is even higher. When evaluating the conditional generation capability, the authors chose to merge materials with Eg = 0 and near-zero band gaps into the range 0 ≤ Eg < 0.5 eV. This approach, which is based on a dataset where the proportion of Eg = 0 is very high and has very low proportion of 0 < Eg < 0.5 eV, seems confusing to me and does not adequately validate the conditional generation capability. I recommend a separate discussion on this aspect to provide a clearer validation of the conditional generation capability. Once again, thank you to the authors for their hard work and for the thorough response to our feedback.

Authorsrebuttal2024-08-13

Further responses from authors

Dear Reviewer Sr8M, Thank you for your responses and recognition of our work. We appreciate your insightful questions regarding the evaluation of the conditional generation ability of the proposed method. Below, we provide further clarifications and discussions which we will also include in the manuscript. From our understanding, your question mainly concerns the high proportion of materials with zero band gap and the low proportion of materials with a band gap value between 0 and 0.5 eV. You suggest that this may not adequately validate the model's conditional generative ability, as the proportion of materials with 0 ≤ Eg < 0.5 eV is substantial. In other words, if the ratio of a group of materials in a dataset is large enough, an unconditional model might also generate a significant proportion of materials that satisfy the given group condition. To clarify this, let us begin with the dataset distribution. The training set we used for generating materials with high or low band gap values contains 61,541 crystals in total, with 46,933 crystals (77.9%) having 0 ≤ Eg < 0.5 eV, 9,814 crystals (15.9%) with 0.5 ≤ Eg ≤ 3 eV, and 4,794 crystals (7.8%) with Eg > 3 eV. As you mentioned in your comments, the ratio of crystals with 0 ≤ Eg < 0.5 eV is indeed large in the training set. Therefore, to better demonstrate the model's conditional generation capacity, we not only show the success rate of generating crystals with low band gap values (<0.5 eV) but also the success rate of generating crystals with high band gap values (> 3 eV). This approach highlights the ability of the proposed method to significantly alter the distribution from the training set. As shown in Table 4 of the paper and Table 9 in the general response, Mat2Seq can significantly change the generated crystal distribution from 7.8% with Eg > 3 eV in the training data to more than 90% in the conditional generation distribution, with a remarkable 92.2% uniqueness ratio and 89.8% validity ratio. Therefore, we believe that the conditional generation capacity of Mat2Seq is well demonstrated. Furthermore, another reason we do not group crystals with a specific band gap value (e.g., 0 eV) into a single group but rather use a range of band gap values for this task is that current state-of-the-art machine learning-based band gap predictors (such as ComFormer, ALIGNN, and others) treat the prediction of band gap values as a regression task, rather than first classifying materials as metal or non-metal and then performing regression for non-metal materials. As a result, it is challenging to determine whether a generated material satisfies a 0 eV band gap because the mean absolute error (MAE) can be as large as 0.122 eV. Thank you again for your recognition and insightful suggestions and comments. If you have any additional questions, we would be more than happy to answer them. Yours sincerely, The Authors

Authorsrebuttal2024-08-13

Email notification was down on OpenReview

Dear Reviewer Sr8M, The email notification system on OpenReview was down for a period of time, and we posted our further responses to your insightful additional question during that period. Thank you for your recognition, as well as your insightful suggestions and comments. If you have any additional questions, we would be more than happy to answer them. Yours sincerely, The Authors

Authorsrebuttal2024-08-13

We appreciate the chance to answer any additional questions

Dear Reviewer 5Kst, Thank you again for the valuable time you invested in reviewing this work. For your previous questions and concerns, we have provided detailed clarifications along with extensive additional experiments. We tend to believe these clarifications and additional experiments could thoroughly address your concerns and questions. As the author-reviewer discussion period is closing soon, we would greatly appreciate the chance to answer any additional questions you may have. Yours sincerely, The Authors

Reviewer 5Kst2024-08-14

Thanks for the rebuttal

Thanks for the detailed rebuttal and response. 1. I am not arguing that the contribution is weak solely because the proposed method only uses Niggli reduction. However, I still question whether the methods applied after the initial Niggli reduction algorithm—such as determining the origin, assigning unique ordering of atoms, and representing lattice vectors in terms of lengths and angles (which are also widely used in the field of materials science)—are strong enough as the main contribution. 2. This issue has been resolved. 3. Thanks for conducting the experiment. Could you explain why the structural validity is low compared to the other baselines in this case? Additionally, it would be beneficial to include the results of CrystalLLM in the final version. 4, 5, 6: The unclear parts have been clarified. Some of my concerns have been resolved, but some still remain. I will raise my score from 3 to 4. Thanks.

Authorsrebuttal2024-08-14

Thank you for joining the discussion

Dear Reviewer 5Kst, Thank you very much for joining the discussion. Following your responses, we would like to provide further clarification on your remaining questions, specifically questions 1 and 3. For your first question, we respectfully disagree with your point and would like to emphasize that **Mat2Seq, to the best of our knowledge, is the first work that can achieve the unique and SE(3)-invariant conversion from 3D crystal structures to 1D sequences.** If this is not the case, could you please point us to related works that have solved this problem before? Additionally, **if this problem has been widely studied and solved, why have previous LLM-based methods, including your cited CrystalLLM from Meta and the other two works, failed to maintain such uniqueness and SE(3) invariance for the 1D crystal sequence representation?** We also respectfully disagree with your point of viewing the novelty of a method solely by whether the components of the method are easy or not. Rather, **a method can be novel and efficient even if each step is straight forward, as long as it solves challenges that have not been addressed by previous studies, and these challenges are important to be solved rather than ignored.** Furthermore, let's consider the novelty of this work by evaluating what the field will be like with or without it. **Currently, none of the previous studies on LLM-based crystal generation have achieved unique and SE(3)-invariant 1D sequence representations.** Without this work, researchers in this field might remain confused and continue using similar approaches, like CIF files used by CrystalLLM, which are far from unique. With this work, we not only demonstrate how to achieve uniqueness and SE(3) invariance—eliminating the need for extensive data augmentation—but also show that this approach enhances the model's performance in generating crystals. For your third question: 1. As we can see from the results, the structural validity is actually similar to FlowMM when the generation temperature is set to 1.35, and it can be further improved to achieve higher structural validity by lowering the temperature. Additionally, the combined validity, where the generated crystal is both structurally valid and compositionally valid, is 83.4%, which is higher than FlowMM’s 80.6% and DiffCSP’s 83.3%, as measured by the structural validity ratio multiplied by the composition validity ratio. 2. The reason we do not directly compare with CrystalLLM (Meta) for generation tasks is simply that **CrystalLLM from Meta uses pre-trained LLMs trained on a vast amount of text data**, while **Mat2Seq and CrystaLLM from the UK team are trained from scratch solely for crystal sequences.** Thank you again for joining the discussion. We are glad that some of your concerns have been addressed, and we appreciate the opportunity to further discuss the remaining issues. If you have any additional questions, we are more than willing to answer them. Yours sincerely, The Authors

Program Chairsdecision2024-09-25

Decision

Accept (poster)

© 2026 NYSGPT2525 LLC