[
    {
        "id": 1,
        "chunk": "# SMiPoly: Generation of a Synthesizable Polymer Virtual Library Using Rule-Based Polymerization Reactions  \n\nMitsuru Ohno,\\* Yoshihiro Hayashi, Qi Zhang, Yu Kaneko, and Ryo Yoshida\\*  \n\nCite This: J. Chem. Inf. Model. 2023, 63, 5539−5548",
        "category": " Abstract"
    },
    {
        "id": 2,
        "chunk": "# ACCESS  \n\nMetrics & More  \n\nArticle Recommendations  \n\nSupporting Information  \n\nABSTRACT: Recent advances in machine learning have led to the rapid adoption of various computational methods for de novo molecular design in polymer research, including high-throughput virtual screening and inverse molecular design. In such workflows, molecular generators play an essential role in creation or sequential modification of candidate polymer structures. Machine learningassisted molecular design has made great technical progress over the past few years. However, the difficulty of identifying synthetic routes to such designed polymers remains unresolved. To address this technical limitation, we present Small Molecules into Polymers (SMiPoly), a Python library for virtual polymer generation that implements 22 chemical rules for commonly applied polymerization reactions. For given small organic molecules to form a candidate monomer set, the SMiPoly generator conducts possible polymerization reactions to generate an exhaustive list of potentially synthesizable polymers. In this study, using 1083 readily available monomers, we generated 169,347 unique polymers forming seven different molecular types: polyolefin, polyester, polyether, polyamide, polyimide, polyurethane, and polyoxazolidone. By comparing the distribution of the virtually created polymers with approximately 16,000 real polymers synthesized so far, it was found that the coverage and novelty of the SMiPoly-generated polymers can reach 48 and $53\\%$ , respectively. Incorporating the SMiPoly library into a molecular design workflow will accelerate the process of de novo polymer synthesis by shortening the step to select synthesizable candidate polymers.  \n\n![](images/c6bb04d570f5ef74cd3ac4b847c054e12983cf4fffb330e5e854a9a0bb03b4d2.jpg)",
        "category": " Abstract"
    },
    {
        "id": 3,
        "chunk": "# INTRODUCTION  \n\nIn recent years, data-driven and computer-aided molecular design technologies have spread rapidly1,2 across various areas of materials research, such as drug discovery,3,4 catalyst design,5−8 and polymer design.9−12 The task of molecular design boils down to the forward and backward predictions.13 The objective of the forward prediction is to obtain a forward mapping from any given chemical structure to its properties. In a conventional scenario, supervised learning is performed to learn the structure−property relationships. In the backward prediction, the inverse mapping of the forward model is obtained to identify promising chemical structures that exhibit desired properties. The most traditional approach to solving the inverse problem is to perform high-throughput virtual screening; millions or billions of candidate molecules are computationally created, and the forward model is used to determine whether their properties fall into the target region or not.14−16 Alternatively, adaptive heuristic search algorithms, such as evolutionary computatio n,17−19 Bayesian optimization,20 particle swarm optimization,21 and Monte Carlo computation11,13,22 have also been used to iteratively modify the candidate molecules to achieve the desired properties.  \n\nHere, the performance of molecular generators plays a key role in determining the success or failure of molecular design.  \n\nTraditionally, virtual libraries have been constructed by probabilistically recombining chemical fragments from which substituent groups or ring structures of existing molecules are pre-extracted.17,23 Recently, molecular generators using statistically trained generative models have also become widely used.24 For example, a probabilistic language model is trained to generate simplified molecular input line entry system (SMILES)25 strings by converting existing molecules into the training string set.11,13,15,22 Alternatively, deep generative models based on the graph representation of molecules have been studied extensivel y.26−28 While the traditional methods limit the searchable space to a combination of predefined fragment sets, generative machine-learning models can create novel molecules from a wider chemical space. However, it is not guaranteed that such computationally designed virtual molecules are feasible to synthesize.29 One way to overcome this drawback is to use synthetic reaction models.30,31 From a large set of synthetic reactions of organic compounds in a public database, a deep neural network is trained to predict the synthetic products of a given set of reactants.32 In addition to such generative machine-learning models, rule-based models compiling prior knowledge of chemical reactions,14 hybrid models of machine learning and rule-based methods,16,24,33 and retrosynthetic analysis34,35 have been proposed. By feeding commercially available reactants to such a prediction model, a virtual library can be constructed in which synthesizability and synthetic routes are taken into account.  \n\nThe aforementioned methodologies have been developed primarily for designing small molecules. However, the development of general-purpose polymer generators has lagged far behind that of small molecule generators. Thus far, polymer repeating units with their chemical structures have been generated by statistical generative models trained from existing polymers in the same manner as the generation of small molecules. Wu et al.11 designed and synthesized three amorphous polymers with high thermal conductivity by solving an inverse problem using a probabilistic language model trained on synthetic polymers recorded in the polymer properties database PoLyInfo.36 In addition, several applications of polymer generative models have been reported for inverse design for the dielectric constant and solubility of polymers.16,37,38 Ma and Luo15 constructed a virtual library of more than one million polymers using a deep language generative model trained on the synthetic polymers in PoLyInfo as well. Methodological and practical research on the computational design of polymeric materials is less advanced than for small molecules.39,40 The lack of data resources for polymeric materials is a significant bottleneck. Furthermore, the synthesis of computationally generated polymers presents a greater technical hurdle than for small materials; in order to estimate the synthesizability of the generated polymers, it is necessary to break a designed polymer down into starting monomers and investigate their availability. Computational methodologies for this have not been well studied.  \n\nRecently, Kim et al.16 provided a generative model of synthetically accessible polymer repeating units with a rulebased polymerization reaction algorithm. Using this system, they constructed a database called the Open Macromolecular Genome (OMG) that contains highly synthesizable virtual polymers. The OMG will be an important database for datadriven polymer research, but there is room for improvement in the definition of rule sets. From the viewpoint of synthetic organic chemistry, the reactivity of a substrate is influenced by the steric and electrical effects of the substituents in the reaction center. Furthermore, as pointed out in that paper, the selectivity of the reaction is influenced by coexisting functional groups in the reactant molecule. Therefore, it is necessary to provide reaction rules that take these factors into account.  \n\nIn this study, we built a virtual library generator for polymers that implements a comprehensive rule set for practically applied polymerization reactions with a Python open-source library, Small Molecules into Polymers (SMiPoly). The generators implement 22 reaction rules that consist of six chain polymerization reactions and 16 step-growth polymerization reactions. The types of polymerization reactions used in actual polymer synthesis are quite limited compared to the overall organic synthetic reactions. Our rule set covers a significant portion of common polymerization reactions considered in previous studies.41−44  \n\nUsing 1083 input monomers carefully selected based on availability, cost, safety, and legal compliance, seven classes of polymers, i.e., polyolefins, polyesters, polyethers, polyamides, polyimides, polyurethanes, and polyoxazolidones, were polymerized in silico. The generated virtual polymers have significantly promising synthesizability with the navigation of polymerization routes. The distribution of the 169,347 generated polymers and nearly 16,000 synthesized polymers forming the 21 polymer classes in PoLyInfo were compared to systematically and quantitatively evaluate the coverage and novelty of the virtual library with respect to the existing polymers. Our library generator covered approximately $50\\%$ of the chemical space of the polymers synthesized to date, and approximately $50\\%$ of the polymers produced were outside the distribution of the existing polymers.",
        "category": " Introduction"
    },
    {
        "id": 4,
        "chunk": "# METHODS  \n\nThe SMiPoly library implements a rule set of 22 polymerization reactions, which consists of two submodules, “monc.py” and “polg.py” (Figure 1). The submodule “monc.py” extracts monomers from the given set of starting molecules and classifies them into 19 different monomer classes by checking whether or not to include polymerizable functional groups of the 22 polymerization reactions. The submodule “polg.py” applies the functional group transformations programmed into an applicable polymerization reaction to generate the chemical structure of a polymer repeating unit from the given one or two starting monomers. The current implementation is capable of generating seven classes of polymers. In the following, each step is described in detail.  \n\n![](images/43a5dcb639b689f55ac275d7f0612a1a9baae28837132154768f292478ac639a.jpg)  \nFigure 1. Workflow of the monomer classification and the polymer generation with SMiPoly.  \n\nPolymerization Reactions. The 22 polymerization reactions are classified in accordance with domain knowledge as summarized in Figure 2.41,42,44 The definition of the terms follows the chemical terminology of the International Union of Pure and Applied Chemistry (IUPAC)43,45 except for “stepgrowth polymerization” and “addition-condensation” as described below.  \n\nAt first, we classified common polymerization reactions into two types: chain polymerization and step-growth polymerization. These two types are classified according to the reaction mechanism of the polymer chain growth. The polymerization reactions in which the growth of a polymer chain proceeds exclusively by a chain reaction are categorized as chain polymerization, and others are defined as step-growth polymerization. The step-growth polymerization here indicates a polymerization reaction in which the growth of polymer chains proceeds by reactions between molecules of all degrees of polymerization. The chain polymerization can be further classified into addition chain polymerization and ring-opening chain polymerization on the basis of the type of chemical reactions involved in the growth step. The step-growth polymerization can be classified into three types: polycondensation, polyaddition, and addition-condensation. The additioncondensation here means a polymerization reaction in which addition and condensation reactions are repeated alternately.  \n\n![](images/7c99c5520db5b138060cd7c4435b4fd39702489c8ed349cc6382832ac4386640.jpg)  \nFigure 2. Classification of polymerization reactions implemented in the SMiPoly library. The definition of the terms is in accordance with the chemical terminology of the International Union of Pure and Applied Chemistry (IUPAC)43,45 except for “step-growth polymerization” and “additioncondensation”.  \n\nEach reaction class consists of several polymerization reactions defined by specific functional group transformations. A polymerization reaction is expressed as a set of the starting monomer(s) and the functional group transformation represented by reaction SMiles ARbitrary Target Specification (SMARTS).46 The 22 polymerization reactions currently implemented and their functional group transformations represented by reaction SMARTS are summarized in Table 1, and illustrative examples of the rule-specific functional group transformations are shown in Table ??. The current rule set can generate the seven polymer classes: polyolefin, polyester, polyether, polyamide, polyimide, polyurethane, and polyoxazolidone.  \n\nChain polymerization is applied to a single starting monomer having at least one polymerizable functional group. Step-growth polymerization reactions are applied to two starting monomers with at least two polymerizable functional groups within each molecular entity or a single monomer with two different functional groups for the same polymerization reaction.  \n\nSelection and Mapping Monomers to Polymerization Reactions. First, from the given set of starting molecules, we select those that can form monomers in the polymerization reaction rule set according to the presence or absence of the class-specific functional groups. The monomer or monomer pair is then mapped to an applicable polymerization reaction depending on the presence or absence of polymerizable functional groups in the reaction rule. The polymerizable functional groups in the rule set include vicinal and/or geminal functional groups as summarized in Table 1. Here, starting molecules containing incompatible functional groups are automatically excluded from the calculation.  \n\nIn this study, 1083 starting molecules were manually extracted from the literature and lists of commercially available compounds. Their chemical structures were expressed in SMILES strings.25 SMiPoly’s submodule “monc.py” classified the 1083 compounds into 19 monomer classes according to the presence or absence of the class-specific functional groups. In the step-growth polymerization, the upper and lower limits for the occurrence number of a polymerizable functional group in a given monomer were limited to the default values of 2 and 4, respectively. The classification result is summarized in Tables 2 and 3, respectively.  \n\nLibrary Generation. To generate a virtual polymer library, an applicable set of polymerization reactions is performed according to the classification of a given monomer set. Polymerizable functional groups identified by “polg.py” are denoted by asterisks. If necessary, an option can be specified to limit the class of polymers generated. The generator lists all applicable monomer combinations for a given polymerization reaction according to the pre-computed monomer class and converts one or two starting monomers to products. The polymerization of polyolefins by addition chain polymerization and polyesters and polyamides composed of monomers with two different functional groups for the same polymerization reaction within the same molecular entity produces bipolymers as well as homopolymers. Polymers with regioisomeric structures were treated as different compounds. The relation between monomer classes, polymerization reactions, and polymer classes is illustrated in Figure S1.",
        "category": " Materials and methods"
    },
    {
        "id": 5,
        "chunk": "# RESULTS AND DISCUSSION  \n\nAs summarized in Table 2, the 1083 starting molecules were classified into one or more of the 19 different monomer classes according to the presence or absence of the class-specific reaction sites. From these monomers, a total of 574,837 virtual polymers were generated by applying the 22 reaction rules. Here, polymers with the same repeating units but different starting monomers were distinguished. These virtual polymers constituted the seven polymer classes (polyolefin, polyester, polyether, polyamide, polyimide, polyurethane, and polyoxazolidone) as summarized in Table 3. Each of these polymers is defined by a single repeating unit because the current virtual library consists of either homopolymers or alternation copolymers. The class of polyolefins is given by homopolymers and binary alternating copolymers that are generated from all combinations of the olefin and/or cyclic olefin monomers. The  \n\nTable 1. 22 Polymerization Reactions with Monomer Classes To Be Applied, and Their Functional Group Transformations   \n\n\n<html><body><table><tr><td colspan=\"6\">ReactionSMARTS</td></tr><tr><td colspan=\"6\">polymer class monomer class 1</td></tr><tr><td rowspan=\"2\">1</td><td colspan=\"2\"> polymerization reaction</td><td></td><td rowspan=\"2\"> reaction type</td><td rowspan=\"2\"> addition chain polymerization</td></tr><tr><td colspan=\"3\">polyolefin vinyl</td></tr><tr><td colspan=\"5\"></td></tr><tr><td rowspan=\"2\">2</td><td>[CX3;H2,H1,H0:1]=[C;H2,H1,H0:2] >*-[C:1][C:2]-* polyolefin</td><td>cyclic olefin</td><td></td><td rowspan=\"2\">addition chain polymerization</td></tr><tr><td rowspan=\"2\"></td><td colspan=\"3\">[CX3;R:1]=[CX3;R:2]>*-[CX4;R:1][CX4;R:2]-*</td></tr><tr><td colspan=\"5\">polyolefin</td></tr><tr><td rowspan=\"3\">3b</td><td>vinyl</td><td>vinyl [ CX3;H2,H1,H0;!R:1]=[CX3;H2,H1,H0;!R:2].[CX3;H2,H1,H0;!R:3]=[CX3;H2,H1,H0;!R:4]></td><td> addition chain polymerization</td><td rowspan=\"3\"></td></tr><tr><td>*-[CX4:1][CX4:2][CX4:3][CX4:4]-*</td><td></td><td></td></tr><tr><td colspan=\"3\">polyolefin vinyl</td></tr><tr><td rowspan=\"3\">4b</td><td></td><td>cyclic olefin [CX3;H2,H1,H0;!R:1]=[CX3;H2,H1,HO;!R:2].[CX3;H1,H0;R:3]=[CX3;H1,H0;R:4]></td><td></td><td rowspan=\"3\"> addition chain polymerization</td></tr><tr><td>*-[CX4:1][CX4:2][CX4:3][CX4:4]-*</td><td></td><td></td></tr><tr><td colspan=\"3\">polyolefin cyclic olefin</td></tr><tr><td rowspan=\"2\">sb 6</td><td></td><td>[CX3;H1,H0;R:1]=[CX3;H1,H0;R:2].[CX3;H1,H0;R:3]=[CX3;R:4]>*-[CX4:1][CX4:2][CX4:3][CX4:4]-*</td><td>cyclic olefin</td><td rowspan=\"2\"> addition chain polymerization</td></tr><tr><td rowspan=\"2\"> polyester</td><td colspan=\"3\">lactone ring-opening chain polymerization</td></tr><tr><td rowspan=\"2\">7</td><td>[CX3;R:1](=[OX1])[OX2;R:2]>(*-[CX3:1](=[OX1]).[OX2:2]-*)</td><td></td><td rowspan=\"2\"> polycondensation</td></tr><tr><td rowspan=\"2\"></td><td colspan=\"3\">polyester hydroxy carboxylic acidc</td></tr><tr><td rowspan=\"2\">8b</td><td>polyester hydroxy carboxylic acid</td><td>([OX2H1;$(OC=*):1].[CX3:2](=[O])[OX2H1])>(*-[OX2:1].[CX3:2](=[O])-*）</td><td> polycondensation</td></tr><tr><td colspan=\"3\">hydroxy carboxylic acidc ([OX2H1;!$(OC=*):1].[CX3:2](=[O])[OX2H1]}).([OX2H1;!$(OC=*):3].[CX3:4](=[O])[OX2H1])></td></tr><tr><td rowspan=\"3\">9</td><td>(*-[OX2:1].[CX3:2](=[O])[OX2:3].[CX3:4](=[O])-*)</td><td></td><td></td></tr><tr><td></td><td>di/polycarboxylic acid di/polyol</td><td rowspan=\"2\">polycondensation</td></tr><tr><td colspan=\"3\">polyester</td></tr><tr><td rowspan=\"2\"></td><td></td><td></td><td>([CX3:1](=[O])[XH1,CBr][CX3:2](=[O])[OXH1CBr]).([O,S;X;H1;$([OS]C=*):3].[O,S;X;H1;$([O,S]C=*):4])></td><td rowspan=\"2\"></td></tr><tr><td colspan=\"3\">(*-[CX3:1](=[O]).[CX3:2](=[O])-[O,S;X2;!$([O,S]C=*):3].[O,S;X2;!$([O,S]C=*):4]-*)</td></tr><tr><td rowspan=\"3\">10</td><td>polyesterd ([O,S;X2;H1;!$([O,S]C=*):1].[O,S;X2;H1;!$([O,S]C=*):2]).[C-]#[O+]></td><td>di/polyol carbon monoxidee</td><td rowspan=\"2\"> polycondensationf</td></tr><tr><td>(*-[O,S;X2;!$([O,S]C=*):1].[O,S;X2;!$([O,S]C=*):2][CX3](=[O])-*)</td><td></td></tr><tr><td>polyester cyclic anhydride</td><td>epoxide</td><td> ring-opening chain polymerization</td></tr><tr><td rowspan=\"2\">11</td><td></td><td>[C,c;R:1][CX3,c;R](=[OX1])[OX2,o;R][CX3,c;R](=[OX1])[C,c;R:2].[CX4;R:3]1[OX2;R:4][CX4;R:5]1></td><td rowspan=\"2\"></td></tr><tr><td colspan=\"3\">([C,c:1][CX3](=[OX1])(-*).[C,c:2][CX3](=[OX1])[OX2][CX4:3][CX4:5][OX2:4]-*)</td></tr><tr><td>12</td><td>polyether epoxide</td><td></td><td> ring-opening chain polymerization</td></tr><tr><td colspan=\"2\"></td><td>[CX4;H2,H1,H0;R:1]1[O;R][C;R:2]1>*-[CX4:1][CX4:2][O]-*</td><td></td></tr><tr><td colspan=\"2\">13</td><td>hindered phenol</td><td> polycondensationg</td></tr><tr><td colspan=\"2\"></td><td>[c]1([OH1:1])[c:2][c:3][c;H1:4][c:5][c:6]1>[c]1([OX2:1]-[*])[c:2][c:3][c:4](-*)[c:5][c:6]1</td><td> polycondensation</td></tr><tr><td colspan=\"2\">14 polyetherh</td><td>bis(p-halogenated aryl)sulfone di/polyol (without thiol)</td><td></td></tr><tr><td colspan=\"2\"></td><td>[c:1]1[c:2]p[c:3]([F,CLBr,1]})[c:4][c:][c:6]1[SX4](=[OX1])(=[OX1])[c:7]2[c:8][c:9][c:10]([F,CLBr,1])[c:11][c:12]2. ([OX2;H1;!$([O,S]C=*):13].[OX2;H1;$([O,S]C=*):14])></td><td></td></tr><tr><td colspan=\"2\"></td><td>[c:1]1[c:2][c:3](-[*])[c:4][c:5][c:6]1[SX4](=[OX1])(=[OX1])[c:7]2[c:8][c:9][c:10]([OX2;1$([O,s]C=*):13].</td><td></td></tr><tr><td colspan=\"2\">15 polyetheri</td><td>[OX2;!$([O,S]C=*):14]-[*])[c:11][c:12]2</td><td> polycondensation</td></tr><tr><td colspan=\"2\"></td><td>bis(p-fluoroaryl)ketone di/polyol (without thiol)</td><td>[c:1]1[c:2][c:3]([F])[c:4][c:5][c:6]1[CX3](=[OX1])[c:7]2[c:8][c:9][c:10]([F])[c:11][c:12]2.([OX2;H1;!$([O,S]C=*):13].</td></tr><tr><td colspan=\"2\">[OX2;H1;!$([O,s]C=*):14])></td><td></td><td></td></tr><tr><td colspan=\"2\">[OX2;!$([O,S]C=*):14]-[*])[c:11][c:12]2</td><td>[c:1]1[c:2][c:3](-[*])[c:4][c:5][c:6]1[CX3](=[OX1])[c:7]2[c:8][c:9][c:10]([OX2;!$([O,S]C=*):13].</td><td></td></tr><tr><td colspan=\"2\">16 polyamide</td><td>lactam</td><td> ring-opening chain polymerization</td></tr><tr><td colspan=\"2\"></td><td>[CX3;R:1](=[OX1])[NX3;R:2]>(*-[CX3:1](=[OX1])[NX3:2]-*)</td><td> polycondensation</td></tr><tr><td colspan=\"2\">17</td><td>polyamide amino acidc ([NX3;H2,H1;!$(OC=*):1].[CX3:2](=[O])[OX2H1])>(*-[NX3:1].[CX3:2](=[O])-*)</td><td></td></tr><tr><td colspan=\"2\">18b</td><td>amino acidc amino acidc</td><td> polycondensation</td></tr><tr><td colspan=\"2\"> polyamide</td><td>([N&X3;H2,H1;!$(NC=*):1].[CX3:2](=[O])[OX2H1]).([N&X3;H2,H1;!$(NC=*):3].[CX3:4](=[O])[OX2H1])></td><td></td></tr><tr><td colspan=\"2\"></td><td></td><td></td></tr><tr><td colspan=\"2\"> polyamide</td><td>(*-[NX3;!$(NC=*):1].[CX3:2](=[O])[NX3;!$(NC=*):3][CX3:4](=[O])-*)</td><td></td></tr><tr><td colspan=\"2\">19</td><td>di/polycarboxylic acid di/polyamine</td><td> polycondensation</td></tr><tr><td colspan=\"2\"></td><td>([CX3:1](=[O])[OXH1,CBr][CX3:2](=[O])[OX2H1,CI,Br]).([N&X3;H2,H1;$(NC=*):3].[N&X3;H2,H1;!$(NC=*):4])></td><td></td></tr><tr><td colspan=\"2\"></td><td>(*-[CX3:1](=[O]).[CX3:2](=[O])-[NX3;!$(NC=*):3].[NX3;!$(NC=*):4]-*)</td><td></td></tr><tr><td colspan=\"2\">20</td><td>di/polycyclic anhydride primary di/polyamine</td><td> polycondensation</td></tr><tr><td colspan=\"2\">polyimide</td><td>([CX3,c;R:1](=[OX1])[OX2,o;R][CX3,c;R:2](=[OX1]).[CX3,c;R:3](=[OX1])[OX2,0;R][CX3,c;R:4](=[OX1])).</td><td></td></tr><tr><td colspan=\"2\"></td><td></td><td></td></tr><tr><td colspan=\"2\"></td><td>([C,c:5][NX3;H2;!$(N[C,S]=*)].[C,c:6][NX3;H2;!$(N[C,S]=*)])></td><td></td></tr><tr><td colspan=\"2\"></td><td></td><td></td></tr><tr><td colspan=\"2\">polyurethane</td><td>([CX3,c;R:1](=[OX1])[NX3;R]([C,c:5][C,c:6]-*)[CX3,c;R:2](=[OX1]),[CX3,c;R:3](=[OX1])[NX3;R](-*)[CX3;R:4](=[OX1]))</td><td></td></tr><tr><td colspan=\"2\">21 （[NX2:1]=[CX2]=[OX1,SX1:2].[NX2:3]=[CX2:4]=[OX1,SX1:5]).</td><td>di/polyisocyanate di/polyol</td><td> polyaddition</td></tr></table></body></html>  \n\nTable 1. continued   \n\n\n<html><body><table><tr><td></td><td>polymer class polymerization reaction</td><td> monomer class 1</td><td>monomer class 2</td><td>reaction type</td></tr><tr><td colspan=\"5\">22</td></tr><tr><td rowspan=\"4\"></td><td colspan=\"3\">([OX2,SX2;H1;!$([O,S]C=*):6].[OX2,SX2;H1;!$([O,S]C=*):7])> (*-[CX3](=[OX1,X1:2])[NX3:1][X3:3][CX3:4](=[OX1,SX1:5])[OX2,SX2;$([O,S]C=*):6][X2,SX2;!$([O,S]C=*):]-*)</td></tr><tr><td colspan=\"3\">polyoxazolidone di/polyepoxide di/polyisocyanate polyaddition</td></tr><tr><td colspan=\"3\">([CX4;H2,H1,H0;R:1]1[OX2;R:2][CX4;H1,H0;R:3]1.[CX4;H2,H1,H0;R:4]2[OX2;R:5][CX4;H1,H0;R:6]2). ([OX1,SX1:7]=[CX2:8]=[NX2:9].[OX1,SX1:10]=[CX2:11]=[NX2:12])></td></tr><tr><td colspan=\"3\">(*-[OX2:2][CX4:3]([CX4:1]-*).[CX4;R:6]1[OX2;R:5][CX2;R:8](=[OX1,SX1:7])[NX3;R:9][CX4;R:4]1.</td></tr></table></body></html>  \n\naThe 7 different polymer classes can be generated by the 22 polymerization rules. bBipolymer. $^c\\mathrm{A}$ monomer with two different functional groups for the same polymerization reaction within the same molecular entity. dPolycarbonate. eSynthetic equivalent of phosgene. fOxidative carbonylation in the case of CO. gUsually oxidative polymerization was applied. hPolyethersulfone (PES) was generated. iPolyetherketone (PEK) or polyetheretherketone was generated.  \n\nTable 2. Result of Monomer Classification   \n\n\n<html><body><table><tr><td>monomer class</td><td>no. of compounds</td></tr><tr><td>vinyl</td><td>462</td></tr><tr><td>cyclic olefin</td><td>52</td></tr><tr><td>epoxid</td><td>71</td></tr><tr><td>lactone</td><td>13</td></tr><tr><td>lactam</td><td>1</td></tr><tr><td>hydroxy carboxylic acid</td><td>7</td></tr><tr><td>amino acid</td><td>1</td></tr><tr><td>hindered phenol</td><td>8</td></tr><tr><td>bis(p-halogenated aryl) sulfone</td><td>4</td></tr><tr><td>bis(p-fluoroaryl)ketone</td><td>1</td></tr><tr><td>di/poly epoxide</td><td>14</td></tr><tr><td>di/polycarboxylic acid</td><td>85</td></tr><tr><td>di/polyol (include thiol)</td><td>162</td></tr><tr><td>di/polyol (without thiol)</td><td>132</td></tr><tr><td>di/polyamine</td><td>149</td></tr><tr><td>primary di/polyamine</td><td>143</td></tr><tr><td>di/polyisocyanate</td><td>19</td></tr><tr><td>cyclic carboxylic acid ahnydride</td><td>33</td></tr><tr><td>di/polycyclic catboxylic acid anhydride</td><td>28</td></tr></table></body></html>  \n\nTable 3. Number of SMiPoly-Generated Polymers Classified into Seven Polymer Classes   \n\n\n<html><body><table><tr><td> polymer class</td><td> no. of generated polymers</td></tr><tr><td>polyolefin</td><td>523,080</td></tr><tr><td>polyester</td><td>25,204</td></tr><tr><td>polyether</td><td>868</td></tr><tr><td>polyamide</td><td>15,554</td></tr><tr><td>polyimide</td><td>5,445</td></tr><tr><td>polyurethane</td><td>4,200</td></tr><tr><td>polyoxazolidone</td><td>486</td></tr><tr><td>total</td><td>574,837</td></tr></table></body></html>  \n\nclasses of polyesters and polyamides consist of homopolymers and binary alternating copolymers generated from the combinations of hydroxy or amino carboxylic acid. The other polymer classes are all composed of homopolymers. Figure 3 depicts several examples of generated polymers for each polymer class along with their polymerization reactions. The Supporting Information provides more randomly selected examples that belong to the seven different polymer classes (Figure S2).  \n\nWe investigated the degree of overlap or coverage and novelty between the distributions of polymers synthesized to date and the virtually created polymers. As the set of existing polymers, 16,223 polymers were extracted from PoLyInfo, the world’s largest database of synthetic polymers. The uniform manifold approximation and projection (UMAP)47 of all polymers with their chemical structures in PoLyInfo and the SMiPoly virtual library is displayed in Figure 4, which shows that SMiPoly is generally able to cover the distribution of existing polymers. A more detailed, quantitative evaluation was performed for each of the seven polymer classes as described below.  \n\nPoLyInfo defines 21 different polymer classes based on its own classification criteria. The polymer classes are identified by the first three letters of the polymer ID (referred to as PID in PoLyInfo), such as P01, P02, P03, etc. Here, the definitions of polymer classes in SMiPoly and PoLyInfo do not exactly match. To eliminate ambiguity in the classification criteria, the polymers generated by SMiPoly and those recorded in PoLyInfo were reclassified to 21 classes using the classifier function (poly.polyinfo_classifier) in the Python library RadonPy,48 which mimics the classification criterion of PoLyInfo’s 21 classes. The six classes of RadonPy that were consistent with SMiPoly and PoLyInfo were used for comparison. The correspondence of the six classes with the RadonPy’s 21 classes, SMiPoly’s seven classes, and PoLyInfo’s 21 classes is summarized in Table 4. In addition, only polymers with two asterisks in the SMILES string of the repeating unit, i.e., linear homopolymers and alternating copolymers, were considered for comparison because it is difficult to represent the branching and cross-linking polymers with SMILES. Further elimination of structural redundancy reduced the number of polymers produced by SMiPoly to 169,347 using the “poly.full_match_smiles_listself” function of RadonPy.48  \n\nHere, we denote the sets of actually synthesized and virtual polymers belonging to a certain polymer class by $S_{\\mathrm{{R}}}$ and $s_{\\mathrm{v}},$ respectively. We assessed the coverage and novelty of $S_{\\mathrm{v}}$ with respect to $S_{\\mathrm{{R}}}$ as follows:  \n\n(1) Evaluate the pairwise similarity between all polymers in $S_{\\mathrm{{R}}}$ and $S_{\\mathrm{v}}$ . As the similarity measure, we employed the Tanimoto coefficient, in which the chemical structure of each polymer was transformed into a 2048-bit vector with an extended connectivity fingerprint with a radius of three atoms (ECFP6).49 To consider the repeating structure of polymers, the ECFP6 descriptor was constructed after generating the macrocyclic oligomer with 10-mer of the repeating unit using RadonPy.4   \n(2) Set the threshold values of the Tanimoto coefficient as $\\gamma\\in$ $\\{0,0.1,...,0.9,1\\}$ .   \n(3) Coverage: calculate the percentage of polymers in $S_{\\mathrm{{R}}}$ with a similarity greater than γ to those in $S_{\\mathrm{v}}$ .   \n(4) Novelty: calculate the percentage of polymers in $S_{\\mathrm{v}}$ with a similarity less than $\\gamma$ to those in $S_{\\mathrm{{R}}}$ .   \n(5) Vary the threshold $\\gamma$ from 0 to 1, and draw a curve representing the balance between (1-coverage) and novelty (coverage−novelty (CN) curve) as shown in Figures 5 and 6.  \n\n![](images/2dfe89588a7a23d3b6229ae69c38de9784bd49b84badc1993196f9e7289de4a4.jpg)  \n\n![](images/06d24df4e9d7284ec0a412af653583244aefc2f4864ba3ba145c051b6202a148.jpg)  \ndifferent polymer classes and their polymerization reactions: (a) polyolefin, (b) polyester, (c) polyether, (d) polyamide, (e) polyimide, (f) polyurethane, and $\\mathbf{\\eta}(\\mathbf{g})$ polyoxazolidone.  \n\nThe CN curve shows an upward or downward convex pattern depending on the inclusive relationship of the distributions of $S_{\\mathrm{{R}}}$ and $\\dot{S}_{\\mathrm{v}}$ . If the two distributions are nearly coincident, the CN curve is drawn on the $45^{\\circ}$ line. If the distribution of $S_{\\mathrm{v}}$ encompasses $s_{\\mathrm{{R}}},$ the CN curve deviates slightly from the $45^{\\circ}$ line and shows an upward convex pattern. This is the most desirable virtual library, given that it contains reasonably novel and diverse polymers while maintaining a high coverage of existing ones. Conversely, if the virtual library set fails to reproduce a part of the synthetic polymers, the CN curve shows downward convexity. This is visually illustrated in Figure S3.  \n\nFigure 5 summarizes the CN curves for each of the six polymer classes and the distributions of $S_{\\mathrm{R}}$ and $S_{\\mathrm{v}}$ visualized on a two-dimensional plane by performing UMAP.47 The virtual libraries for polyamide, polyimide, and polyurethane, which exhibited moderately upward convex CN curves, achieved wellbalanced coverage and novelty. When the threshold of Tanimoto similarity was set to $\\gamma\\geq0.5.$ , the coverage and novelty were approximately 67 and $45\\%$ for polyamide, 70 and $59\\%$ for polyimide, and 70 and $52\\%$ for polyurethane, respectively. The CN curves for polyether showed a slightly downward convex pattern. The coverage and novelty were 51 and $38\\%$ when the threshold was set to $\\gamma\\geq0.5$ .  \n\nWe also investigated the coverage with respect to all 16,223 polymers in PoLyInfo (Figure 6). The coverage and novelty of the 169,347 virtual polymers for all these polymers were 48 and $53\\%$ under $\\gamma\\geq0.5.$ , respectively.  \n\nThe performance of the rule-based virtual polymer generator depends on the diversity of the starting monomer set and the comprehensiveness of polymerization reaction rules. As shown by the UMAP projection in Figure 6, the current SMIPoly virtual library clearly failed to cover some regions of the existing polymer distribution in PoLyInfo, which are denoted by regions A−C (see example polymers belonging to each region in Figure S4 in the Supporting Information). The presence of region A was due to the fact that the current SMiPoly does not implement the condensation reaction of hydroxycarboxylic acids and amino acids even though both types of monomers were included when the library was created. Region B arose because the starting monomer set did not include the relevant molecules; as shown in Table 2, the starting monomer set contains very few amino acid molecules. The polymers belonging to region C were mainly polyarylenes. The current SMiPoly could not cover this area because polymerization reactions with aryl coupling reactions and their monomers were undefined. This is the main reason for the low coverage and novelty in the generation of the class of polyethers. The class of polyethers defined by PoLyInfo (P07 and P37) includes polyoxides with chain structures as well as polymers with cyclic structures including $\\scriptstyle{\\mathrm{C-O-C}}$ bonds and acetal structures in the main chain, but the current SMiPoly does not support these polymerization reactions or monomers.",
        "category": " Results and discussion"
    },
    {
        "id": 6,
        "chunk": "# CONCLUSIONS  \n\nTo generate synthesizable virtual polymer candidates, we have developed the Python library SMiPoly that implements 22 chemical rules for polymerization reactions. For any given set of starting molecules, the generator creates synthetic products by conducting automatically selected polymerization reactions applicable to the input monomers. All polymers produced are potentially synthesizable.  \n\nIn this study, using 1083 available monomers, we generated 169,347 unique linear polymers forming seven different types of polymers, such as polyolefin, polyester, polyether, polyamide, polyimide, polyurethane, and polyoxazolidone. The reaction space that can be represented by the current release of SMiPoly covers approximately $50\\%$ $(\\gamma\\ge0.5)$ of the entire chemical space of polymers synthesized to date. When focusing on nitrogen containing polymers defined by SMIPoly, it covered $70\\%$ of the corresponding chemical space with around $50\\%$ novelty $(\\gamma\\geq$ 0.5). The coverage can be monotonically increased by further expanding the list of raw monomers and the set of polymerization reactions.  \n\n![](images/fe46c2736fc416b4402a52d5d858dd34d27026369707a6c716070e50a26bd7d3.jpg)  \n· Polyolefins (hydrocarbons) Polystyrenes Polyvinyls · Polyacrylics · Polyhalo-olefins · Polydienes Polyethers Polysulfides Polyesters Polyamides ? Polyurethanes Polyureas Polyimides Polyanhydrides Polycarbonates Polyimines Polysilanes Polyphosphazenes Polyketones Polysulfones · Polyphenylenes   \nFigure 4. UMAP projection of all polymers with their chemical structures in (a) PoLyInfo (16,223 polymers) and (b) SMiPoly virtual library (169,347 polymers). ECFP6 descriptor, which was constructed after generating the macrocyclic oligomer with 10-mer of the repeating unit, was used as input of the UMAP. The 21 classes of the polymer backbones are color-coded according to the classification of RadonPy.  \n\nTable 4. Correspondence of Polymer Classes between Existing Polymers in PoLyInfo, 21 Classes of RadonPy, and Seven Classes of SMiPoly-Generated Polymers   \n\n\n<html><body><table><tr><td> polymer class</td><td> RadonPy</td><td>PoLyInfo PIDa</td><td> SMiPoly</td><td>no. of polymers in PoLyInfo</td><td> no. of polymers in SMiPoly</td></tr><tr><td>polyolefin</td><td>hydrocarbon of polyolefin</td><td>P01, P31</td><td>polyolefin</td><td>3,139</td><td>139,990</td></tr><tr><td rowspan=\"6\"></td><td>polystyrene</td><td>P02, P32</td><td></td><td></td><td></td></tr><tr><td> polyvinyl</td><td>P03, P33</td><td></td><td></td><td></td></tr><tr><td> polyacrylate</td><td>P04,P34</td><td></td><td></td><td></td></tr><tr><td> halogenated polyolefin</td><td>P05, P35</td><td></td><td></td><td></td></tr><tr><td>polydiene</td><td>P06,P36</td><td></td><td></td><td></td></tr><tr><td> polyester</td><td>P09, P39</td><td>polyester</td><td>2,775</td><td>9,903</td></tr><tr><td>polyester polyether</td><td>polyether</td><td>P07, P37</td><td> polyether</td><td>4,692</td><td>5,547</td></tr><tr><td> polyamide</td><td> polyamide</td><td>P10, P40</td><td> polyamide</td><td>5,059</td><td>14,402</td></tr><tr><td> polyimide</td><td>polyimide</td><td>P13, P43</td><td>polyimide</td><td>2,755</td><td>6,553</td></tr><tr><td> polyurethane</td><td> polyurethane</td><td>P11, P41</td><td>polyurethane polyoxazolidone</td><td>695</td><td>2,772</td></tr><tr><td colspan=\"4\"> no. of unique polymers</td><td>16,223</td><td>169,347</td></tr></table></body></html>  \n\nHere, we remark on a comparison with the OMG virtual polymer database, which was created using a rule-based polymerization algorithm similar to SMiPoly; OMG contains approximately 12 million virtual polymers generated from 17 different polymerization reaction rules using approximately 80,000 small molecules as starting materials. Table S2 in the Supporting Information summarizes whether each of the 21 polymer classes of PoLyInfo can be generated by SMiPoly and OMG. Currently, SMiPoly implements $59\\%$ (10/17) of the reaction rules defined in the OMG, and the OMG implements $36\\%$ (8/22) of the reaction rules defined in SMiPoly.  \n\nFinally, the limitations of the current version of SMiPoly are summarized. For example, it is necessary to extend the rule set of polymerization reactions, such as ring-opening metathesis polymerization (ROMP), addition-condensation polymerization to generate phenol and/or melamine resin, and main-chain growth with aryl coupling reaction. The difficulty in defining these polymerization reactions arises from the complexity of defining the monomer structures to be reacted. In addition, in order to improve the accuracy of monomer selection and polymerization reaction prediction, olefin monomers need to be further classified in detail. More importantly, although branched polymers play an important role in the development of thermosetting polymers, the generation of branched polymers cannot be handled by the current version of SMiPoly. Since this limitation is due to SMILES-based modeling, it is necessary to introduce a structural representation that can handle stochastic polymer structures such as BigSMILES.50  \n\n![](images/4317311d37f81eba74120272d6bb153a1bd82e63095033aabb1809ccc418fd0f.jpg)  \nFigure 5. For each class of polymers, the coverage and novelty of the virtual library relative to existing polymers are evaluated based on CN curves and visual inspection of UMAP projections of the two sets of chemical structures: (a) polyolefin, (b) polyester, (c) polyether, (d) polyamide, (e) polyimide, and (f) polyurethane.",
        "category": " Conclusions"
    },
    {
        "id": 7,
        "chunk": "# ASSOCIATED CONTENT",
        "category": " References"
    },
    {
        "id": 8,
        "chunk": "# Data Availability Statement  \n\nThe source code of SMiPoly is available from the GitHub website(https://github.com/PEJpOhno/SMiPoly). A sample script to generate the polymer library is in the directory named “sample_script”. The list of starting molecules is stored in the directory “sample_data”. The list of the selected monomers and generated polymers is also available.",
        "category": " Results and discussion"
    },
    {
        "id": 9,
        "chunk": "# $\\bullet$ Supporting Information  \n\nThe Supporting Information is available free of charge at https://pubs.acs.org/doi/10.1021/acs.jcim.3c00329.  \n\nAdditional experimental details, materials, and methods, including (1) functional group transformations in the polymerization reaction rules (Table S1), (2) monomer class, polymerization reaction, and polymer class (Figure S1), (3) generated polymers belonging to seven different polymer classes and their polymerization reactions (Figure S2), (4) drawing CN curve (Figure S3), (5) examples of polymers belonging to Regions $_{\\mathrm{A-C}}$ in Figure 6 (Figure S4), and (6) comparison of polymerization reaction rule sets between SMiPoly and OMG (Table S2) (PDF)",
        "category": " Materials and methods"
    },
    {
        "id": 10,
        "chunk": "# AUTHOR INFORMATION  \n\nCorresponding Authors Mitsuru Ohno − Daicel Corporation, 530-0011 Osaka, Japan; $\\circledcirc$ orcid.org/0000-0002-4588-1927; Email: mt_ohno@ jp.daicel.com Ryo Yoshida − The Institute of Statistical Mathematics, Research Organization of Information and Systems, Tachikawa, Tokyo 190-8562, Japan; The Graduate University for Advanced Studies, SOKENDAI, Tachikawa, Tokyo 190-8562, Japan; National Institute for Materials Science, 305-0047 Ibaraki, Japan; orcid.org/0000-0001- 8092-0162; Email: yoshidar@ism.ac.jp   \nAuthors Yoshihiro Hayashi − The Institute of Statistical Mathematics, Research Organization of Information and Systems, Tachikawa, Tokyo 190-8562, Japan; The Graduate University for Advanced Studies, SOKENDAI, Tachikawa, Tokyo 190-8562, Japan; $\\circledcirc$ orcid.org/0000-0002-7650- 4083 Qi Zhang − The Institute of Statistical Mathematics, Research Organization of Information and Systems, Tachikawa, Tokyo 190-8562, Japan Yu Kaneko − Daicel Corporation, 530-0011 Osaka, Japan; orcid.org/0000-0001-8617-0876  \n\nComplete contact information is available at: https://pubs.acs.org/10.1021/acs.jcim.3c00329",
        "category": " Abstract"
    },
    {
        "id": 11,
        "chunk": "# Author Contributions  \n\nM.O. and R.Y. designed and conceived the project, and wrote a preliminary draft of the paper. M.O. designed and developed the SMiPoly library, and performed the experiments with assistance from Y.H., ${\\bf{\\cal Z}}.{\\bf{\\cal Q}}.,$ and Y.K. M.O., Y.H., and R.Y. wrote and revised the manuscript. All authors discussed the result and commented on the manuscript.",
        "category": " References"
    },
    {
        "id": 12,
        "chunk": "# Notes  \n\nThe authors declare no competing financial interest.  \n\n![](images/aff97ee7e75321f40a66f30f25b4e17fb536b1fd678223ee70d40ac333132876.jpg)  \nFigure 6. Coverage and novelty of virtual libraries with respect to existing polymers are evaluated based on the CN curve and visual inspections of th distribution of polymers’ chemical structures constituting the two sets by UMAP projections.",
        "category": " Results and discussion"
    },
    {
        "id": 13,
        "chunk": "# ACKNOWLEDGMENTS  \n\nThis work was supported in part by MEXT under “Program for Promoting Researches on the Supercomputer Fugaku” (Project ID: hp210264 to R.Y.), JST CREST Grant Number JPMJCR19I3, and the Grant-in-Aid for Scientific Research (A) 19H01132 from the Japan Society for the Promotion of Science (JSPS).",
        "category": " Acknowledgments"
    },
    {
        "id": 14,
        "chunk": "# REFERENCES  \n\n(1) Matsubara, S. Digitization of Organic Synthesis \u0001 How Synthetic Organic Chemists Use AI Technology \u0001. Chem. Lett. 2021, 50, 475− 481.   \n(2) Williams, W. L.; Zeng, L.; Gensch, T.; Sigman, M. S.; Doyle, A. G.; Anslyn, E. V. The Evolution of Data-Driven Modeling in Organic Chemistry. ACS Cent. Sci. 2021, 7, 1622−1637.   \n(3) Gómez-Bombarelli, R.; Wei, J. N.; Duvenaud, D.; HernándezLobato, J. M.; Sánchez-Lengeling, B.; Sheberla, D.; AguileraIparraguirre, J.; Hirzel, T. D.; Adams, R. P.; Aspuru-Guzik, A. Automatic Chemical Design Using a Data-Driven Continuous Representation of Molecules. ACS Cent. Sci. 2018, 4, 268−276. PMID: 29532027. (4) Segler, M. H. S.; Kogej, T.; Tyrchan, C.; Waller, M. P. Generating Focused Molecule Libraries for Drug Discovery with Recurrent Neural Networks. ACS Cent. Sci. 2018, 4, 120−131. PMID: 29392184. (5) Santiago, C. B.; Guo, J.-Y.; Sigman, M. S. Predictive and mechanistic multivariate linear regression models for reaction development. Chem. Sci. 2018, 9, 2398−2412.   \n(6) Toyao, T.; Maeno, Z.; Takakusagi, S.; Kamachi, T.; Takigawa, I.; Shimizu, K. I. Machine Learning for Catalysis Informatics: Recent Applications and Prospects. ACS Catal. 2020, 10, 2260−2297.   \n(7) Reid, J. P.; Sigman, M. S. Holistic prediction of enantioselectivity in asymmetric catalysis. Nature 2019, 571, 343−348.   \n(8) Yada, A.; Nagata, K.; Ando, Y.; Matsumura, T.; Ichinoseki, S.; Sato, K. Machine Learning Approach for Prediction of Reaction Yield with Simulated Catalyst Parameters. Chem. Lett. 2018, 47, 284−287. (9) Ramprasad, R.; Batra, R.; Pilania, G.; Mannodi-Kanakkithodi, A.; Kim, C. Machine learning in materials informatics: recent applications and prospects. npj Comput. Mater. 2017, 3, 54.   \n(10) Sanchez-Lengeling, B.; Aspuru-Guzik, A. Inverse molecular design using machine learning: Generative models for matter engineering. Science 2018, 361, 360−365.   \n(11) Wu, S.; Kondo, Y.; Kakimoto, M.-a.; Yang, B.; Yamada, H.; Kuwajima, I.; Lambard, G.; Hongo, K.; Xu, Y.; Shiomi, J.; Schick, C.; Morikawa, J.; Yoshida, R. Machine-learning-assisted discovery of polymers with high thermal conductivity using a molecular design algorithm. npj Comput. Mater. 2019, 5, 66.   \n(12) Sattari, K.; Xie, Y.; Lin, J. Data-driven algorithms for inverse design of polymers. Soft Matter 2021, 17, 7607−7622.   \n(13) Ikebata, H.; Hongo, K.; Isomura, T.; Maezono, R.; Yoshida, R. Bayesian molecular design with a chemical language model. J. Comput.- Aided Mol. Des. 2017, 31, 379−391.   \n(14) Lyu, J.; Wang, S.; Balius, T. E.; Singh, I.; Levit, A.; Moroz, Y. S.; O’Meara, M. J.; Che, T.; Algaa, E.; Tolmachova, K.; Tolmachev, A. A.; Shoichet, B. K.; Roth, B. L.; Irwin, J. J. Ultra-large library docking for discovering new chemotypes. Nature 2019, 566, 224−229.   \n(15) Ma, R.; Luo, T. PI1M: A Benchmark Database for Polymer Informatics. J. Chem. Inf. Model. 2020, 60, 4684−4690. PMID: 32986418.   \n(16) Kim, S.; Schroeder, C. M.; Jackson, N. E. Open Macromolecular Genome: Generative Design of Synthetically Accessible Polymers. ACS Polym. Au 2023, DOI: 10.1021/acspolymersau.3c00003. in press. (17) Huan, T. D.; Mannodi-Kanakkithodi, A.; Kim, C.; Sharma, V.; Pilania, G.; Ramprasad, R. A polymer dataset for accelerated property prediction and design. Sci. Data 2016, 3, 160012.   \n(18) Pilania, G.; Iverson, C. N.; Lookman, T.; Marrone, B. L. Machine-Learning-Based Predictive Modeling of Glass Transition Temperatures: A Case of Polyhydroxyalkanoate Homopolymers and Copolymers. J. Chem. Inf. Model. 2019, 59, 5013−5025. PMID: 31697891.   \n(19) Schadler, L.; Chen, W.; Brinson, L.; Sundararaman, R.; Gupta, P.; Prabhune, P.; Iyer, A.; Wang, Y.; Shandilya, A. A perspective on the data-driven design of polymer nanodielectrics. J. Phys. D: Appl. Phys. 2020, 53, 333001.   \n(20) Yang, X.; Zhang, J.; Yoshizoe, K.; Terayama, K.; Tsuda, K. ChemTS: an efficient python library for de novo molecular generation. Sci. Technol. Adv. Mater. 2017, 18, 972−976. PMID: 29435094. (21) Kumar, J. N.; Li, $\\mathrm{Q.;}$ Tang, K. Y. T.; Buonassisi, T.; GonzalezOyarce, A. L.; Ye, J. Machine learning enables polymer cloud-point engineering via inverse design. npj Comput. Mater. 2019, 5, 73.   \n(22) Wu, S.; Lambard, G.; Liu, C.; Yamada, H.; Yoshida, R. iQSPR in XenonPy: A Bayesian Molecular Design Algorithm. Mol. Inf. 2020, 39, No. 1900107.   \n(23) St. John, P. C.; Phillips, C.; Kemper, T. W.; Wilson, A. N.; Guan, Y.; Crowley, M. F.; Nimlos, M. R.; Larsen, R. E. Message-passing neural networks for high-throughput polymer screening. J. Chem. Phys. 2019, 150, 234111.   \n(24) Yang, J.; Tao, L.; He, J.; McCutcheon, J. R.; Li, Y. Machine learning enables interpretable discovery of innovative polymers for gas separation membranes. Sci. Adv. 2022, 8, No. eabn9545.   \n(25) Weininger, D. SMILES, a chemical language and information system. 1. Introduction to methodology and encoding rules. J. Chem. Inf. Comput. Sci. 1988, 28, 31−36.   \n(26) De Cao, N.; Kipf, T. MolGAN: An implicit generative model for small molecular graphs. arXiv preprint arXiv:1805.11973 2018,   \n(27) You, J.; Liu, B.; Ying, Z.; Pande, V.; Leskovec, J. Graph Convolutional Policy Network for Goal-Directed Molecular Graph Generation. In Advances in Neural Information Processing Systems. 2018. (28) Popova, M.; Isayev, O.; Tropsha, A. Deep reinforcement learning for de novo drug design. Sci. Adv. 2018, 4, No. eaap7885.   \n(29) Gao, W.; Coley, C. W. The Synthesizability of Molecules Proposed by Generative Models. J. Chem. Inf. Model. 2020, 60, 5714− 5723. PMID: 32250616.   \n(30) Guo, Z.; Wu, S.; Ohno, M.; Yoshida, R. Bayesian Algorithm for Retrosynthesis. J. Chem. Inf. Model. 2020, 60, 4474−4486. PMID: 32975943   \n(31) Zhang, $\\mathrm{Q.;}$ Liu, C.; Wu, S.; Yoshida, R. Bayesian Sequential Stacking Algorithm for Concurrently Designing Molecules and Synthetic Reaction Networks. arXiv preprint arXiv:2204.01847 2022, (32) Mikulak-Klucznik, B.; et al. Computational planning of the synthesis of complex natural products. Nature 2020, 588, 83−88. (33) Badowski, T.; Gajewska, E. P.; Molga, K.; Grzybowski, B. A. Synergy Between Expert and Machine-Learning Approaches Allows for Improved Retrosynthetic Planning. Angew. Chem., Int. Ed. 2020, 59, 725−730.   \n(34) Chen, L.; Kern, J.; Lightstone, J. P.; Ramprasad, R. Data-assisted polymer retrosynthesis planning. Appl. Phys. Rev. 2021, 8, No. 031405. (35) Taniwaki, H.; Kaneko, H. Molecular design of monomers by considering the dielectric constant and stability of the polymer. Polym. Eng. Sci. 2022, 62, 2750.   \n(36) Otsuka, S.; Kuwajima, I.; Hosoya, J.; Xu, Y.; Yamazaki, M. PoLyInfo: Polymer Database for Polymeric Materials Design. In 2011 International Conference on Emerging Intelligent Data and Web Technologies, 2011; pp 22 −29, DOI: 10.1109/EIDWT.2011.13. (37) Liu, D.-F.; Feng, Q.-K.; Zhang, Y.-X.; Zhong, S.-L.; Dang, Z.-M. Prediction of high-temperature polymer dielectrics using a Bayesian molecular design model. J. Appl. Phys. 2022, 132, No. 014901.   \n(38) Gurnani, R.; Kamal, D.; Tran, H.; Sahu, H.; Scharm, K.; Ashraf, U.; Ramprasad, R. polyG2G: A Novel Machine Learning Algorithm Applied to the Generative Design of Polymer Dielectrics. Chem. Mater. 2021, 33, 7008−7016.   \n(39) Audus, D. J.; de Pablo, J. J. Polymer informatics: opportunities and challenges. ACS Macro Lett. 2017, 6, 1078−1082.   \n(40) Chen, L.; Pilania, G.; Batra, R.; Huan, T. D.; Kim, C.; Kuenneth, C.; Ramprasad, R. Polymer informatics: current status and critical next steps. Mater. Sci. Eng., R 2021, 144, No. 100595.   \n(41) Chan, C. H.; Chen, J.-T.; Farrell, W. S.; Fellows, C. M.; Keddie, D. J.; Luscombe, C. K.; Matson, J. B.; Merna, J.; Moad, G.; Russell, G. T.; Théato, P.; Topham, P. D.; Sosa Vargas, L. Reconsidering terms for mechanisms of polymer growth: the “step-growth” and “chain-growth” dilemma. Polym. Chem. 2022, 13, 2262−2270.   \n(42) Odian, G. Principles of polymerization; John Wiley & Sons, 2004. (43) Jenkins, A. D.; Kratochvíl, P.; Stepto, R. F. T.; Suter, U. W. Glossary of basic terms in polymer science (IUPAC Recommendations 1996). Pure Appl. Chem. 1996, 68, 2287−2311.   \n(44) Kagiya, T. A Classification of Polymerization Reaction. Kobunshi 1964, 13, 855−858. 854,   \n(45) Compendium of Chemical Terminology [Internet]. https:// goldbook.iupac.org, cited: 2023-04-08.   \n(46) Daylight. SMARTS: a language for describing molecular patterns [Internet]. https://www.daylight.com/dayhtml/doc/theory/theory. smarts.html, cited: 2022-10-14.   \n(47) McInnes, L.; Healy, J.; Melville, J. Umap: Uniform manifold approximation and projection for dimension reduction. arXiv preprint arXiv:1802.03426 2018,   \n(48) Hayashi, Y.; Shiomi, J.; Morikawa, J.; Yoshida, R. RadonPy: automated physical property calculation using all-atom classical molecular dynamics simulations for polymer informatics. npj Comput. Mater. 2022, 8, 222.   \n(49) Rogers, D.; Hahn, M. Extended-Connectivity Fingerprints. J. Chem. Inf. Model. 2010, 50, 742−754. PMID: 20426451.   \n(50) Lin, T.-S.; Coley, C. W.; Mochigase, H.; Beech, H. K.; Wang, W.; Wang, Z.; Woods, E.; Craig, S. L.; Johnson, J. A.; Kalow, J. A.; Jensen, K. F.; Olsen, B. D. BigSMILES: A Structurally-Based Line Notation for Describing Macromolecules. ACS Cent. Sci. 2019, 5, 1523−1531. PMID: 31572779.",
        "category": " References"
    }
]