ARTICLE   Open Access    

A bio-based hyperbranched flame retardant towards the fire-safety and smoke-suppression epoxy composite

  • These authors contributed equally: Zhiqian Lin, Wangbin Zhang

More Information
  • The design of fully biological flame retardants is vital for environment-friendly and sustainability development. Herein, a fully biological hyperbranched flame retardant named PA-DAD was synthesized successfully through a simple neutralization reaction between 1,10-Diaminodecane (DAD) and phytic acid (PA). PA-DAD as a kind of reactive flame retardant was subsequently employed to composite with epoxy resin (EP). The obtained EP composite, comprising 25 wt% PA-DAD, exhibited excellent flame resistance with an appreciative limiting oxygen index (LOI) of 28.0%, and a desired V-0 rating of UL-94. The favorable attributes stem from the evident flame-retardant characteristics of PA-DAD, manifesting particularly within the gaseous and condensed phases of combustion. Benefitting from the PA-DAD with the synergetic effect of smoke suppression and flame resistance, the EP/25% PA-DAD composite displayed highlighted reductions both in smoke production and heat release rate. In contrast with neat EP, total smoke production (TSP), peak smoke production release (pSPR), the peak heat release rate (pHRR), and the rate of fire growth (FIGRA) of the EP/25% PA-DAD composite were decreased by 49.5%, 57.0%, 72.2% and 77.8%, respectively. Moreover, the EP/25% PA-DAD composite resulted in well-preserved mechanical properties, especially enhanced toughness, compared with the neat EP. The strategies in our work provided a facile, green, and highly efficient way for fabricating high fire-retardant EP composites.
  • Corydalis DC. is the largest and most diverse genus of Papaveraceae, comprising more than 500 species, which are mainly divided into four subgenera (subg. Bipapillatae, subg. Cremnocapnos, subg. Sophorocapnos, and subg. Corydalis) and 39 sections[1]. Corydalis is widely distributed in the Northern Hemisphere with a few species extending into East Africa[2]. Moreover, this genus demonstrates remarkable diversity in China, particularly in the Himalaya–Hengduan Mountains and the adjacent regions[2,3]. As a genus with spectacular radiation, Corydalis exhibits an extremely high level of diversity in its leaves, subterranean organs, fruits, seeds, flower color, and the length of spurs, and can adapt to diverse habitats, such as riversides, forests, shrubs, grasslands, screes, or even cliffs. Previous studies speculated that Corydalis had an ancient origin, and underwent rapid radiation after the middle Miocene, which was most likely promoted by the continuous orogenesis and climate change associated with the uplift of the Qinghai-Tibetan Plateau (QTP)[47].

    The sequencing of plant genomes has greatly advanced our understanding of the underlying diversification mechanisms[810]. Although considerable progress has been made recently in the classification and phylogeny of Corydalis, the underlying diversification mechanisms remain poorly understood. Particularly, polyploidy is common in this genus, suggesting that polyploidization possibly contributed to the diversification of Corydalis. Polyploidy could drive dramatic changes in the genome landscapes, such as genome size and structural variation, providing a genetic basis for the adaptation and diversification of species-rich groups. While more and more studies link polyploidization events with speciation, fewer studies have documented chromosomal variation in plant groups that underwent rapid radiations, thus the extent to which it may have contributed to radiation is still elusive, especially in biodiversity hotspots. In Corydalis, only the genomes of C. tomentella and C. yanhusuo have been released to date[11,12], which has hindered our understanding of its diversification and adaptation.

    Moreover, a large number of Corydalis species have been frequently used in folk medicine due to their antibacterial, antiviral, and anticancer activities. For instance, C. yanhusuo, C. bungeana, and C. decumbens are the most famous medicinal plants recorded in the Pharmacopoeia of China (http://db.ouryao.com/yd2020/). High medicinal value of C. sheareri, C. hendersonii, C. incisa, C. repens, C. edulis, C. racemosa, and C. pallida, also have been reported[2]. Previous phytochemical investigations have isolated various components from Corydalis, including alkaloids, coumarins, flavonoids, anthraquinones, triterpenes, steroids, and organic acids[13]. Of particular interest are alkaloids, especially benzylisoquinoline alkaloids (BIAs), which have been reported to play a crucial role in sedation, releasing pain, promoting blood circulation, and inhibiting cancer cells[13]. Compared with other genera in Papaveraceae, Corydalis can produce multiple types of BIAs, such as cavidines, apocavidine, tetrahydropalmatine, corydalis, protopine, dehydroapocavidine, and dehydrocavidine. Notably, it was also reported that species-specific BIAs were also isolated in different species of Corydalis[1316], which indicates the remarkable diversity of BIAs in these plant groups. However, the limited genomic resources hampered an in-depth understanding of the diversity of BIAs in Corydalis. Fortunately, biosynthesis pathways of BIAs have been proposed in Corydalis[12,17], as well as in other Papaveraceae species, such as Macleaya cordata[18], Eschscholzia californica[19], and Papaver somniferum[2022]. The biosynthesis of BIAs in Corydalis involves a series of enzymes, namely berberine bridge enzyme (BBE), berberine bridge enzyme-like (BBEL), C-methyltransferase (CMT), cytochrome P450 (CYP), norcoclaurine synthase (NCS), N-methylcoclaurine 30-hydroxylase (NMCH), coclaurine N-methyltransferase (CNMT), 4-hydroxyphenylpyruvate decarboxylase (HPPDC), O-methyltransferase (OMT), tyrosine aminotransferase (TAT), tetrahydroprotoberberine N-methyltransferase (TNMT), tyrosine decarboxylase (TYDC), and tyramine 3-hydroxylase (TYR).

    In this study, PacBio long read sequencing, chromosome conformation capture (3C)-based Hi-C sequencing, and Illumina short-read sequencing were used to assemble a high-quality genome of Corydalis sheareri, one species from subg. Corydalis, the largest and most diverse lineages of Corydalis. Genomic comparison was then carried out between C. sheareri and C. tomentella, one previously reported diploid genome from subg. Sophorocapnos. Additionally, we identified the candidate BIAs biosynthesis genes in representative species in Papaveraceae and traced their evolution history by combing phylogenetic reconstruction, chromosomal location, and gene duplication analyses. Our study will provide more insights into the genome evolution as well as the BIAs diversity in Corydalis.

    Fresh leaves, stems, rhizomes, and flowers at different developmental stages of C. sheareri were collected from Zhongshan Botanical Garden (Nanjing, China) (Fig. 1a). After collection, these samples were immediately frozen in liquid nitrogen or dried in silica gel followed by preservation at −80 °C in the laboratory. The silica gel-dried leaves were used for the flow cytometry measurement, and the material stored in liquid nitrogen was used for genome and transcriptome sequencing. Genomic DNA was extracted using a modified CTAB method. Total RNA was extracted from leaves, stems, rhizomes, and flowers at different developmental stages using RNAprep Pure Plant Kit (Tiangen, Beijing).

    Figure 1.  Overview of Corydalis sheareri genome. (a) Photo of C. sheareri. (b) Genome survey of C. sheareri. (c) Heatmap for coverage pattern of heterozygous k-mer pairs, in which X axis indicates the normalized minor k-mer coverage, while Y axis indicates the total k-mer pairs coverage. (d) Genomic features of eight pseudochromosome. The outermost circle (blue) represents each chromosome of the genome. The bar charts of the second to fifth circles suggest gene density, LTR density, Copia, and Gypsy density, respectively. The inner circular shows inter-chromosomal synteny.

    For the genome survey, the 350 bp paired-end library was constructed according to the Illumina protocols and sequenced on the Illumina Novaseq 6000 platform. For PacBio HiFi sequencing, the DNA SMRT bell library with an insert size of approximately 15 kb was prepared using SMRTbell® Express Template Prep Kit 2.0 (Pacific Biosciences, PN 101-853-100), and subsequently sequenced on the PacBio Sequel II platform (Pacific Biosciences, USA). For high-throughput 3C-based Hi-C sequencing, the library was generated using the standard procedures and sequenced on the Illumina Novaseq 6000 platform. For transcriptome sequencing, RNA from different tissues was pooled equally for library construction to obtain more expressed genes. Thereafter, the cDNA library with an insert size of 300−500 bp was prepared using VAHTS mRNA-seq v2 Library Prep Kit for Illumina (Vazyme) and sequenced on the Illumina Novaseq 6000 platform to generate paired-ends reads. All sequencing was carried out in Berry Genomics Company (www.berrygenomics.com), Beijing, China. The raw Illumina data for genome survey, Hi-C, and RNA-seq reads, was trimmed using Fastp v0.23.2[23] to remove the adaptors and low-quality paired reads. The HiFi reads were generated through CCS software (https://github.com/PacificBiosciences/ccs) with the following parameters: --min-passes = 3 --min-rq = 0.99.

    To guide genome sequencing and assembly, flow cytometry[24], and genome survey were used to estimate the genome size. For flow cytometry, the dried material was chopped with a razor blade and then the DNA content was measured following the protocol of CyStain®PI Absolute P (Sysmex-Partec, Germany) with an Elite flow cytometer (BD FACSCalibur, USA) at the Institute of Botany, Chinese Academy of Sciences (China). Finally, the genome size was inferred based on the external reference standards (Glycine max = 1,100 Mb). For genome survey, the obtained clean reads were used to estimate the genome size based on the k-mer method with KMC v3[25], and Genomescope v2.0[26]. Given that multiple peaks were detected in the k-mer spectrum, Smudgeplot[26] was further employed to visualize and evaluate the ploidy and genome structure through the analysis of heterozygous k-mer pairs.

    Hifiasm v0.14.2[27] was used for de novo assembly of the draft contig genome of C. sheareri, which contained unphased contigs from two homolog genomes (csh v1.0). Purge_dups v1.2.3[28] was used to improve the assembly by removing duplications. The filtered Hi-C reads were aligned to the draft genome (csh v1.0) using Juicer v1.6.2[29], and scaffolding was performed on the contigs with 3D-DNA v180922[30]. Juicebox v1.9.8[31] was used to visualize the Hi-C heatmap, and the scaffolds were manually adjusted to get the chromosomal-level assembly of C. sheareri (csh v2.0). BWA v0.7.15[32] and minimap2[33] were used to assess the quality of the genome by aligning the Illumina short reads and the long HiFi reads to the genome, respectively. BUSCO v4.14[34] was performed to evaluate the integrity of the assembled genome by searching against the 1,614 conserved single-copy genes obtained from the embryophyta_odb10 database.

    LTR_Finder v1.07[35], MITE-Hunter v1.0[36], and Repeat-Masker v4.1.0 (www.repeatmasker.org) were used to predict the repeat sequences. Protein-coding genes were predicted by combining the results of ab initio-based, homology-based, and RNAseq-based predictions. For ab initio prediction, Augustus v3.2.2[37], Snap v6.0[38], Glimmer hmm v3.0.4[39], and GeneMark-ET v4.57[40] were utilized to predict the gene structure in the repeat-masked genome. GeMoMa v1.7.1[41] was used to perform homology prediction with Arabidopsis thaliana, Oryza sativa, Macleaya cordata, and Papaver rhoeas as references. For RNAseq-based prediction, transcriptomic data from different tissues were assembled de novo using trinity v2.2.0[42]. PASA r20140417[43] was used to predict the gene structure based on the obtained transcript. All predicted protein-coding genes were annotated by blast against five databases, including KEGG (www.genome.jp/kegg/brite.html), Gene Ontology (GO) terms, NR (https://ftp.ncbi.nlm.nih.gov/blast/db/FASTA/nr.gz), SwissProt (https://ftp.uniprot.org/pub/databases/uniprot/current_release/knowledgebase/complete/uniprot_sprot.fasta.gz), and eggNOG (http://eggnog5.embl.de/#/app/home). The tRNA (transfer RNA) was predicted by tRNAscan-SE v2.0[44]. Other types of noncoding RNAs (ncRNAs), including rRNA (ribosomal RNA), miRNA (microRNA), and snRNA (small nuclear RNA), were annotated by BLAST against the Rfam database (https://ftp.ebi.ac.uk/pub/databases/Rfam/14.1/).

    Eight species from Papaveraceae were selected for phylogenomic analysis, including C. sheareri, C. tomentella, Eschscholzia californica, Macleaya cordata, Papaver somniferum, Capnoides sempervirens, Ceratocapnos vesicaria, and Hypecoum procumbens, with Aquilegia coerulea (Ranunculaceae) as outgroup. The genomes or transcriptomes were directly retrieved from Genome Warehouse in National Genomics Data Center (https://ngdc.cncb.ac.cn/gwh), National Center for Biotechnology Information (www.ncbi.nlm.nih.gov/datasets/genome), and ONEKP database (https://db.cngb.org/onekp). Gene families were clustered using protein sequences by OrthoFinder v2.2.7[45] with the parameter '-S diamond' . Each gene family was aligned with MAFFT[46]. Single or low-copy gene families (one or two copies in the polyploid species Papaver somniferum, while one copy in the other eight species) with more than 50 amino acids were retained to reconstruct the phylogenetic tree by RAxML v8.1.17[47] with 1,000 bootstrap replicates under PROTCATWAG model using protein sequences. Divergence time was estimated under a relaxed molecular clock model by the MCMCTree implement in the PAML v4.8[48]. Calibration points retrieved from TimeTree (http://timetree.org) were used as priors in divergence time estimation, including the split of Papaveraceae and Ranunculaceae (103~117 Ma), the crown age of Papaveraceae (65~111 Ma), the split age of Hypecoideae and Fumarioideae (63~96 Ma) and the crown age of Fumarioideae (38~44 Ma). Gene family expansion and contraction were inferred by CAFE v4.1[49], with an input species tree constructed from the single or low-copy orthologs. Meanwhile, we performed functional enrichment analysis for the expanded gene families by BLAST against the KEGG and GO databases.

    The whole-genome duplication (WGD) and whole-genome triplication (WGT) events were identified by Ks method, syntenic analyses, and phylogenomic methods. For Ks method, the homologous gene pairs were firstly identified by the all-against-all BLASTP search (e-value cutoff < 1e-5). Then, YN00 in PAML v4.8[48] was called by WGDI[50] to calculate the synonymous substitution rate (Ks) of each gene pair between two species or within a single species. For syntenic analyses, collinear blocks for intra- and interspecies comparisons were detected using MCScanX v0.8[51] with '-s 15', meaning that each block contained at least 15 collinear gene pairs. JCVI v0.8.12[52] was used to draw dotplots of C. sheareri, C. tomentella, V. vinifera, and A. trichopoda with the default parameters. TBtools[53] was used to visualize the synteny between C. sheareri and C. tomentella. Additionally, structural variants between C. sheareri and C. tomentella, i.e., inversion, translocation, and duplication, were identified with SyRI v1.6.3[54]. For phylogenomic methods, gene families were firstly identified by OrthoFinder v2.2.7[45] and multiple sequence alignments were performed by MAFFT[46]. Then, Maximum-Likelihood (ML) trees were constructed using RAxML[47], with bootstrap values estimated from 100 replicates using the PROTCATWAG model. The WGD/WGT event was identified by tree2GD (https://sourceforge.net/projects/tree2gd/) with the default parameter and WGD/WGT events were considered to have occurred according to any of the following conditions: (1) gene duplication (GD) > 500, of which the number of (AB)(AB) type is over 250; (2) GD > 1,500, of which (AB)(AB) type is over 100, and at the same time, the sum of (AB)(AB) type and (AB)A or (AB)B type are over 1,000[55].

    To gain more insights into the diversity of BIAs in Corydalis, 12 gene families involved in BIAs biosynthesis, i.e., BBE, BBEL, CMT, CYP719, CYP80B, CYP82N, NCS, NMT, OMT, TAT, TYDC, and TYR, were identified. All BIA biosynthesis genes reported in C. tomentella genome were firstly retrieved by the gene ID reported previously[12]. Then, TBLASTN was conducted to identify the candidate BIA-related genes in C. sheareri and seven other species of Papaveraceae, including Chelidonium majus, Capnoides sempervirens, Ceratocapnos vesicaria, Hypecoum procumbens, Eschscholzia californica, Macleaya cordata, Papaver somniferum, with Aquilegia coerulea as the outgroup. A sequence is regarded as a candidate gene if it encompasses the entire domain region and the pairwise amino acid identity between the queries and the targets exceed 40%. The annotation for each candidate gene was manually checked according to the blast result. All candidate sequences were confirmed by BLAST against the InterPro (www.ebi.ac.uk/interpro/) and the NCBI Conserved Domain Database (CDD, www.ncbi.nlm.nih.gov/Structure/cdd/wrpsb.cgi) with e < 1e-10.

    To infer the expansion mechanism of the BIAs biosynthesis genes, phylogenetic reconstruction, chromosomal location, and gene duplication analyses were performed. ML analyses were conducted based on the protein sequences by IQ-TREE v1.6.8[56] using the JTT model and 100 bootstrap replications. Gff files, gene files, and targeted BIA biosynthesis gene IDs of C. sheareri and C. tomentella were downloaded or extracted. The chromosome location information was obtained from the genome annotation file, and visualized by TBtools[53]. Gene duplication events were analyzed using MCScanX[51] with a BLASTp search (e < 1e-10). The synonymous (Ks) and nonsynonymous (Ka) values of the identified species-specific tandem duplicated gene pairs were calculated by the M0 model in PAML v4.8[48].

    Flow cytometry indicated that C. sheareri had a genome size of approximately 580 Mb (Supplementary Fig. S1a). Based on the obtained 66.29 Gb short paired-end reads (Supplementary Table S1), the 23-mer distribution showed that the estimated genome size of C. sheareri was 261 Mb (Fig. 1b), and the k-mer spectrum exhibited four distinct peaks at ~40, 80, 120, 160, which was highly similar to the results for autotetraploids (M. sativa and S. spontaneum)[26]. Moreover, nucleotide heterozygosity analysis showed 1.73% aaab and 1% aabb, which was consistent with the expectation that the heterozygous rate of autotetraploid AAAB would be greater than that of AABB[26]. The Smudgeplot analysis also revealed that more heterozygous k-mer pairs concentrated at 1/4 for the normalized coverage of minor k-mer and 4n for the total coverage of k-mer pairs, and the prevalence of AAAB (53%) was considerably greater than AABB (21%) (Fig. 1c). All these results suggested that the genome of C. sheareri exhibited complex genome structure and possibly an autotetraploid.

    PacBio Sequel II sequencing yielded approximately 33.02 Gb of high quality Hifi reads, and the average length and N50 of filtered subreads (1,990,989) were 16,584 bp and 16,614 bp, respectively (Supplementary Table S1). De novo assembly generated the draft genome of C. sheareri (csh v1.0, ~550 Mb) with contig N50 of 9.18 Mb, and the longest contig was approximately 25.09 Mb (Table 1). The assembly was further scaffolded with 67.96 Gb Hi-C data (Supplementary Table S1). Finally, the high-quality genome of C. sheareri (csh v2.0, ~282 Mb) comprised eight pseudochromosomes and 36 contigs (Supplementary Fig. S1b, Supplementary Table S2), with contig N50 of 11.39 Mb (Table 1). The minimum length of the chromosome was greater than 30 Mb (Supplementary Table S2). Approximately 97.62% of DNA reads and 88.32% of RNA-seq reads could be mapped to the assembly genome. BUSCO analysis revealed that 97.1% (1,567/1,614) of the core eukaryotic genes were completely present in the C. sheareri genome, of which 91.4% were single copy (1,475) and 5.7% were duplicated (92), while 0.9% (15), and 2.0% (32) were partially present or missing, respectively.

    Table 1.  Genome assembly and annotation of Corydalis sheareri.
    Analytical process Characteristic C. sheareri
    Genome survey Genome size (flow cytometry) (Mb) 580
    Genome size (k-mer spectrum) (Mb) 261
    Assembly_csh v1.0 Number of contigs 2,713
    Assembly size (Mb) 550
    Contig N50 (Mb) 9.18
    Shortest contig (bp) 16,037
    Largest contig (bp) 25,090,490
    Assembly_csh v2.0 Total number of contigs 179
    Assembly size (Mb) 282
    Contig N50 (Mb) 11.39
    Number of pseudochromosomes 8
    GC content 37.11%
    Annotation Number of protein-coding genes 26,287
    Mean gene length (bp) 5,092
    Mean CDS length (bp) 1,440
    Complete BUSCOs (C) 1,567
    Percentage of repeat sequences (%) 44.47
     | Show Table
    DownLoad: CSV

    Approximately 44.47% (125,404,725 bp) of the C. sheareri genome was annotated as transposable elements (TEs), of which 26.47% were retrotransposons and 7.37% were transposons (Supplementary Table S3). For retrotransposons, a total of 84,322 LTR elements were identified, of which 53,097 (7.27%) belonged to the Copia superfamily and 29,747 (5.05%) belonged to the Gypsy superfamily (Supplementary Table S3). Based on a combination of homology search, de novo prediction, and RNA-seq based prediction, a total of 26,287 protein-coding genes were confidently annotated for C. sheareri, and the mean lengths of the predicted gene and coding sequence were 5,092 and 1,440 bp, respectively (Table 1). Approximately, 92.39% (24,286/26,287) of the genes were functionally annotated, of which, 24,257, 19,162, 8,566, 5,531, and 23,330 genes showed high similarity to known proteins in the NR, SwissProt, KEGG, GO, and eggNOG databases, respectively (Supplementary Fig. S2). In addition, a total of 2,746 noncoding RNA genes were identified in the genome of C. sheareri, including 1,014 tRNA genes, 919 rRNA genes, 609 snRNA genes, and 86 miRNA genes (Supplementary Table S4).

    A total of 28,545 orthologous groups (OGs) were identified in nine selected species (Fig. 2b, Supplementary Table S5). Of them, 7,769 gene families were shared by all species, while 1,714 gene families were specific to Corydalis. Additionally, 1,025 single or low-copy nuclear gene families were identified in these species. After removing the OGs less than 50 aa, 1,003 single or low-copy nuclear gene families were retained to infer a high-confidence species tree. As expected, three subfamilies were recovered in the phylogenomic analysis, and Hypecoideae was strongly supported as a sister to Fumarioideae. Within Fumarioideae, C. sheareri and C. tomentella formed a highly supported clade, while the relationship of Ceratocapnos vesicaria, Capnoides sempervirens, and Corydalis was not fully resolved. Within Papaveroideae, Eschscholzia californica diverged firstly, and Papaver somniferum was sister to Macleaya cordata (Fig. 2a). Molecular dating indicated that the divergence of C. sheareri and C. tomentella was dated to 24.92 Ma, with a 95% confidence interval (95% CI) of 16.35~33.22 Ma (Fig. 2a). Gene family expansion and contraction analyses revealed that 674 and 566 gene families expanded and contracted in Corydalis, respectively (Fig. 2a). GO and KEGG enrichment analyses showed that the significantly expanded gene families were mainly related to response to stimulus, membrane, cell periphery, response to chemical, and cellular response to stimulus, which primarily enriched in secondary metabolites pathways such as phenylpropanoid biosynthesis, pentose, and glucuronate interconversions, photosynthesis, flavonoid biosynthesis, and monoterpenoid biosynthesis (Supplementary Fig. S3). Specifically, cytochrome P450 (CYP) and photosynthesis proteins were also enriched in the KEGG analysis (Supplementary Fig. S3).

    Figure 2.  Genome evolution of Corydalis sheareri. (a) Phylogenetic tree and divergence time estimation. Pal., Palaeogene; Neo., Neogene; Hyp., Hypecooideae. (b) Gene families identified in each species. (c) Collinearity within the genome of C. sheareri. (d) Ks distributions of anchor pairs for the paralogous genes of C. sheareri, C. tomentella and Aquilegia coerulea, and for orthologous genes between C. sheareri and C. tomentella, A. coerulea, respectively. (e) Collinearity analysis between C. sheareri and C. tomentella. (f) Structural variation detection between C. sheareri and C. tomentella performed by SyRI.

    WGD events play a crucial role in the duplication and retention of genes. Intra-genomic colinearity analyses uncovered remnants of one WGD event in C. sheareri (Fig. 2c). Obviously, one signature peak of synonymous substitutions per synonymous site (Ks) distribution was detected for the C. sheareri genome at approximately 1.0 (Fig. 2d), indicating an ancient WGD event. Similarly, previously sequenced genomes of Ranunculales species, including C. tomentella, and A. coerulea, also showed a signature peak at 1.0~1.2 in their genomes (Fig. 2d), suggesting this WGD event was probably shared by all Ranunculales species. Comparison of the C. sheareri paralogue Ks distribution between Corydalis and Aquilegia coerulea also indicated a WGD occurred in Ranunculales (Fig. 2d). Previous studies reported that the Vitis vinifera genome had an ancestral hexaploidization[57] and the Amborella trichopoda genome shows no evidence of more recent lineage-specific genome duplications[58]. The inter-genomic comparison and syntenic depths analyses among Corydalis sheareri, Vitis vinifera, and Amborella trichopoda further confirmed that this ancient WGD event, in which two paralogous segments in the Corydalis sheareri genome correspond to one and three orthologous regions in the Amborella trichopoda (with no WGD), and Vitis vinifera (with one WGT) genomes, respectively (Supplementary Fig. S4). The tree2GD analyses revealed 1,281 GDs with 934 (73%) (AB)(AB) retention type, which also strongly suggested the identified WGD event was shared by the ancestor of Ranunculales (node 5) (Supplementary Fig. S5). The genome of Corydalis sheareri and C. tomentella exhibited good collinearity (Fig. 2e, f, Supplementary Fig. S6). A number of large structural variations were detected, including multiple invasions across the genome, and a small fraction of translocations on chromosomes 6 (Fig. 2f). To further verify the authenticity of these structural variations, several specific inversion regions were randomly checked by mapping the HiFi reads to the assembled genome with Integrative Genomics Viewer (IGV)[59]. Our results confirmed that these inversions are genuine and not assembly artifacts (Supplementary Fig. S7).

    A total of 172 genes involved in the BIAs biosynthesis have been identified in Corydalis sheareri, which is notably larger than that found in C. tomentella (125), Macleaya cordata (128), and Aquilegia coerulea (138), while smaller than that observed in Papaver somniferum (364) and Eschscholzia californica (197) (Fig. 3b; Supplementary Table S6), which has undergone two WGD events[12]. For each gene family, the upstream BIA biosynthetic genes, such as TAT, TYDC, TYR, BBE, and CYP80B (NMCH), generally possess small copy numbers in Corydalis. Conversely, the downstream BIA biosynthetic genes, including OMT, CMT, NMT, and BBEL, typically have a higher copy number (Fig. 3b; Supplementary Table S6). The only exception is NCS, which encodes enzymes to catalyze the condensation of dopamine and 4-hydroxyphenylacetaldehyde (4-HPAA) to produce norcoclaurine, in the upstream of BIA biosynthesis, which also has a greater number of copies (Fig. 3a, b).

    Figure 3.  BIAs biosynthesis related gene families identified in this study. (a) The biosynthetic pathway of benzylisoquinoline alkaloids (BIAs) in Corydalis. The species-specific tandem duplicated gene pairs or clusters in Corydalis identified in Fig. 5 and Table 2 were highlighted with red color. (b) The copy number of BIAs biosynthesis related gene families. Four species with transcriptome sequence are represented by asterisks. (c) Gene duplication type identified in C. sheareri (left) and C. tomentella (right). Abbreviations: BBE, berberine bridge enzyme; BBEL, berberine bridge enzyme-like; CMT, C-methyltransferase; CNMT, coclaurine N-methyltransferase; CYP, Cytochrome P450; HPPDC, 4-hydroxyphenylpyruvate decarboxylase; NCS, norcoclaurine synthase; NMT, N-methyltransferase; OMT, O-methyltransferase; TAT, tyrosine aminotransferase; TNMT, tetrahydroprotoberberine N-methyltransferase; TYDC, tyrosine decarboxylase; TYR, tyramine 3-hydroxylase. Aco, Aquilegia coerulea; Cma, Chelidonium majus; Cse, Capnoides sempervirens; Csh, Corydalis sheareri; Cto, C. tomentella; Cve, Ceratocapnos vesicaria; Eca, Eschscholzia californica; Hpr, Hypecoum procumbens; Mco, Macleaya cordata; Pso, Papaver somniferum. SD, Singleton; DD, dispersed duplications; PD, proximal duplications; TD, Tandem duplications; WGD, segmental/WGD duplications.

    Chromosomal location analysis indicated that some BIA biosynthetic genes are unevenly distributed among chromosomes. For instance, all 22 NCS genes of C. sheareri are located on chromosome 6, whereas, in C. tomentella, all members are mapped on the chromosome 5. Additionally, in C. sheareri, 34 out of the 39 BBEL genes are located on chromosome 4, 17 out of the 22 NMT genes are mapped on the chromosome 6, and nine out of 15 CYP82N genes are found on chromosome 3 (Fig. 4a). Similarly, in C. tomentella, 25 out of 31 BBEL genes are found on chromosome 6, and 15 out of the 20 NMT genes are located on chromosome 5 (Fig. 4b). On the contrary, some other genes are widely distributed throughout genomes but are uneven among chromosomes in both two species. For instance, seven out of eight chromosomes (except for chr02) harbor the 22 CMT genes in C. sheareri, and 11 of them are located on chromosome 4. Furthermore, five chromosomes (1, 2, 4, 6, and 8) contain the OMT genes, and both chromosome 2 and 4 harbor nine members (Fig. 4a). In C. tomentella, 20 OMT genes are distributed across five chromosomes (1, 2, 4, 5, and 6). and chromosome 6 harbors eight members, chromosome 1 harbors four (Fig. 4b). Nineteen CMT genes are located on seven chromosomes (except for chr01), and 9 of them are located on chromosome 8 (Fig. 4b). Interestingly, a gene cluster including seven NCS genes and ten NMT genes, was identified in chromosome 6 within a 270-kb region of C. sheareri genome (Fig. 4a).

    Figure 4.  Chromosome location and duplication event analyses of BIA biosynthesis genes in (a) C. sheareri and (b) C. tomentella. The chromosome number is indicated at the top of each chromosome. Tandem duplicated genes are indicated with orange color. The scale bar on the left indicates the length (mb) of chromosomes. The gene and species abbreviations were the same as Fig. 3. Chr, chromosome.

    Gene duplication analysis revealed that the majority of BIAs biosynthesis genes were generated through gene duplication events in Corydalis. Over half of these genes were identified as tandem duplications, with 55.2% in C. tomentella and 64.5% in C. sheareri, respectively (Fig. 3c; Supplementary Table S6). Additionally, 18% to 27.2% of genes were identified as dispersed duplications, and one-tenth of genes, corresponding to 10.5% to 12.0%, were recognized as proximal duplications (Fig. 3c; Supplementary Table S6). Meanwhile, a small fraction, ranging from 4.7% to 5.6%, were identified as segmental/WGD (Fig. 3c; Supplementary Table S6). Specifically, BBEL genes exhibit an extremely high probability of tandem duplication in both C. sheareri and C. tomentella. For instance, 31 BBEL genes (79.4%) of C. sheareri were found to be tandemly duplicated on chromosome 4, while 24 members (80%) of C. tomentella tandemly duplicated on chromosome 6 (Figs 3c, 4; Supplementary Table S6). Moreover, NCS, CYP719, and TYR genes also demonstrated an extremely high rate of tandem duplication in C. sheareri, with 19 out of 22 NCS genes (86%), five out of six CYP719 genes (83%) and nine out of 11 TYR genes (82%) were observed to be tandemly duplicated (Figs 3c, 4a; Supplementary Table S6).

    Among the genes associated with the BIAs biosynthesis, the largest is the BBEL family, with 39 members in Corydalis sheareri and 30 members in C. tomentella, respectively (Fig. 3b; Supplementary Table S6). Phylogenetic analysis of BBE and BBEL revealed three monophyletic clades (BBE, I, and II). Clade II has undergone a significant expansion in Corydalis, encompassing 34 members from C. sheareri and 23 from C. tomentella. Seven C. sheareri or C. tomentella-specific monophyletic groups were identified in clade II, and one of them comprised a notably high number of 14 members from C. sheareri (Fig. 5a).

    Figure 5.  Phylogenetic trees of BIAs biosynthesis genes. (a) BBE and BBEL, (b) CYP719, CYP80B (NMCH), and CYP82N, (c) NCS, (d) CMT, (e) NMT, (f) OMT, (g) TAT, (h) TYDC, and (i) TYR. Corydalis sheareri and C. tomentella – specific monophyletic groups are highlighted with green and pink color, respectively. The gene and species abbreviations were the same as Fig. 3.

    CYPs are important for determining chemical diversity in metabolism, in which, CYP719, CYP80B (NMCH), and CYP82N, have been identified as key components of BIA biosynthesis[12,60,61]. Phylogenetic analysis of these three CYP subfamily suggested that 13 CYP82N genes of C. sheareri (a total of 15 members) were clustered into four monophyletic groups, and two members of C. tomentella (CtoCYP82N5 and CtoCYP82N6) were clustered into one highly supported monophyletic group (Fig. 5b). For CYP719 genes, two C. sheareri specific monophyletic groups were identified, consisting of three (CshCYP719_01, CshCYP719_02, and CshCYP719_03) and two members (CshCYP719_04 and CshCYP719_06), respectively. Similarly, four members of C. tomentella clustered into two monophyletic groups in the phylogenetic tree (Fig. 5b).

    Phylogenetic analysis of the NCS genes revealed four highly supported monophyletic clades (I−IV). In clade III, two C. sheareri specific monophyletic groups were identified, comprising seven (CshNCS02, CshNCS03, CshNCS04, CshNCS05, CshNCS06, CshNCS16 and CshNCS17), and three (CshNCS19, CshNCS20 and CshNCS21) members, respectively (Fig. 5c). In clade IV, one C. sheareri specific monophyletic group, and one C. tomentella specific monophyletic group were identified, containing three (CshNCS07, CshNCS08 and CshNCS09) and two members (CtoNCS4 and CtoNCS9), respectively (Fig. 5c).

    Phylogeny the CMT genes revealed seven monophyletic clades (I−VII), and the majority of Corydalis members are scattered throughout the phylogenetic tree (Fig. 5d). Conversely, phylogenetic tree of the NMT genes revealed five highly supported clades (I−V) and eight species-specific monophyletic groups were recognized in Corydalis (Fig. 5e). In clade II, three species-specific monophyletic groups were identified, comprising four C. sheareri members (CshNMT04, CshNMT05, CshNMT07, and CshNMT08), four C. tomentella members (CtoNMT14, CtoNMT15, CtoNMT16, and CtoNMT17), and two C. sheareri members (CshNMT19 and CshNMT20), respectively. In clade IV, four species-specific monophyletic groups, with three C. sheareri members (CshNMT17, CshNMT21, and CshNMT22), three C. tomentella members (CtoNMT4, CtoNMT5, and CtoNMT6), two C. sheareri members (CshNMT14 and CshNMT16) and two C. tomentella members (CtoNMT2 and CtoNMT3), respectively. In clade V, one monophyletic group, containing five C. sheareri members (CshNMT10, CshNMT11, CshNMT12, CshNMT13, and CshNMT15), was identified. Similarly, phylogenetic analysis of the NMT genes revealed eight well-supported clades (I−VIII). In clade III, three C. tomentella members (Cto9OMT, CtoOMT10, CtoOMT11) formed one well-supported monophyletic group. In clade VIII, another three C. tomentella members (Cto6OMT, CtoOMT7, and CtoOMT8) also formed one well-supported monophyletic group. While in clade V, two C. sheareri - specific monophyletic groups were identified, containing four (CshOMT03, CshOMT04, CshOMT05, and CshOMT09), and three (CshOMT08, CshOMT10, and CshOMT11) members, respectively.

    Notably, the close relationship of the species-specific tandem duplicated gene pairs or clusters with high sequence similarity was confirmed in the phylogenetic analyses, such as CshBBEL03 and its paralogs (CshBBEL04 and CshBBEL05), CshBBEL30 and its paralog CshBBEL31, CshCYP719_01 and its paralogs (CshCYP719_02 and CshCYP719_03), CshNCS07 and its paralogs (CshNCS08 and CshNCS09), CshNCS19 and its paralogs (CshNCS20 and CshNCS21), CshNMT19 and its paralog CshNMT20, CshTYR03 and its paralogs (CshTYR04, CshTYR05, and CshTYR06), CtoCYP82N5 and its paralog CtoCYP82N6, and CtoCYP719A2 and its paralog CtoCYP719A2-2 (Fig. 5). The Ka/Ks ratio of these gene pairs ranged from 0.15543 to 0.43823 with an average of 0.292909 (Table 2), suggesting that purifying selection was the primary evolutionary force on these species-specific tandem duplicated gene pairs or clusters.

    Table 2.  Selective pressure and sequence similarity of the species-specific tandem duplicated gene pairs or clusters in Corydalis.
    Gene pairs or clusters Ka Ks Ka/Ks Similarity
    CshBBEL03, CshBBEL04, CshBBEL05 0.1656 0.4321 0.38315 97.54%
    CshBBEL30, CshBBEL31 0.1090 0.3349 0.32543 98.73%
    CshCYP719_01, CshCYP719_02,
    CshCYP719_03
    0.0841 0.3688 0.22803 98.81%
    CshNCS07, CshNCS08, CshNCS09 0.0546 0.2236 0.24430 99.11%
    CshNCS19, CshNCS20, CshNCS21 0.1325 0.8527 0.15543 94.61%
    CshNMT19, CshNMT20 0.1435 0.4799 0.29900 98.68%
    CshTYR03, CshTYR04, CshTYR05, CshTYR06 0.1635 0.373 0.43823 97.79%
    CtCYP719A2, CtCYP719A2-2 0.0841 0.3808 0.22096 97.73%
    CtCYP82N5, CtCYP82N6 0.0913 0.2672 0.34165 97.79%
     | Show Table
    DownLoad: CSV

    In this study, by combining data from PacBio long-read sequencing, 3C-based Hi-C sequencing, and Illumina short-read sequencing, we assembled the genome of Corydalis sheareri, one species from subg. Corydalis, the largest and most diverse lineages of Corydalis. Genome survey showed that its genome is extremely complex (Fig. 1b). Both the GenomeScope and Smudgeplot analyses implied that the genome structure of C. sheareri might be an autotetraploid (a special tetraploid with three homologous chromosomes and one non-homologous chromosome) with a special AAAB karyotype (Fig. 1b, c). Intriguingly, the karyotype of C. sheareri is remarkably similar to that of C. yanhusuo[11], which also belongs to subg. Corydalis.

    As previously reported, the diploid C. tomentella has a genome size of 258 Mb[12], and the estimated tetraploid C. sheareri genome size was 580 Mb (Supplementary Fig. S1a), which is nearly twice that of C. tomentella. This indicates that polyploidization might have played a significant role in the genome evolution within this genus. It is noteworthy that both C. sheareri and C. yanhusuo are tetraploid, yet their genome sizes differ by more than threefold. In Corydalis, both the genome size (https://cvalues.science.kew.org/search/angiosperm) and the chromosome number (https://ccdb.tau.ac.il/; http://legacy.tropicos.org/Project/IPCN) varied considerably. We deduce that diploidization following polyploidization, chromosomal rearrangements, including inversions, translocations, and changes in chromosome number via fusion and fission, as well as gene loss might be common and could thus trigger the genome size diversity in Corydalis. In our study, the detection of large-scale chromosomal structural variants, especially multiple inversions, between the genomes of C. sheareri and C. tomentella (Fig. 2e, f), seems to provide strong evidence in support of this hypothesis.

    Polyploidization and structural variations might not only lead to changes in genome size but also substantially affect the gene content. In this study, we identified 172 candidate genes involved in the BIAs biosynthesis within the C. sheareri genome and traced their evolutionary history in Papaveraceae. As previously reported, BBELs, CMTs, NMTs, OMTs, and CYPs (CYP719, CYP82N) are the key enzymes downstream of the BIA biosynthetic pathway (Fig. 3a). In the Corydalis genome, the genes that encode these downstream enzymes, such as BBEL, CMT, NMT, and OMT, tend to possess a larger number of copies than the upstream genes of the BIA biosynthetic pathway (Fig. 3a, b), which enables plants to synthesize various BIAs. Interestingly, NCS, one upstream gene family, also has a relatively high copy number (Fig. 3b). NCS catalyzes the condensation of dopamine and 4-HPAA to produce norcoclaurine, which was identified as one of rate - limiting enzymes in BIA biosynthesis in opium poppy[62].

    Furthermore, the chromosomal location and gene duplication analyses demonstrated that these genes frequently appeared as tandem duplications scattered throughout the genome (Figs 3c, 4). All these findings corroborated that tandem duplications might play a key role in the diversity of BIA biosynthetic genes in Corydalis. What's more, our comprehensive phylogenetic analyses, which cover representative species throughout Papaveraceae, revealed that clade I and II of NCS, clades I, II, and VI of CMT, clades I and III of NMT, clades II, VI, and VII of OMT, two members of TAT, four members of TYDC, and clade I of TYR, each retain a single copy in Corydalis (Fig. 5). In addition, the fact that their homologs are shared by the majority of species suggests that they likely originate before the divergence of Papaveraceae or even Ranunculales. By contrast, the identification of C. sheareri or C. tomentella - specific monophyletic groups in BBEL, CYP719, CYP82N, clades III and IV of NCS, clades V and VI of CMT, clades II, IV and V of NMT, clades V and VIII of OMT, clade III of TYR (Fig. 5) indicate that recent gene duplications might occur frequently for these members in Corydalis. Particularly, the close relationship of some tandem duplicated gene pairs or clusters was confirmed in our phylogenetic analyses (Fig. 5), strongly suggesting that these genes may be generated from the recent duplication events that occurred within C. sheareri, or C. tomentella. Additionally, more C. sheareri-specific tandem duplication events are identified compared to those in C. tomentella, which is likely associated with the significant expansion of BIAs biosynthetic genes in C. sheareri.

    Until recently, with the burst of plant genome sequencing projects and the advancement of bioinformatic tools, biosynthetic pathways for many natural products have been elucidated[6365]. Numerous studies have stated that tandem duplication is a major factor contributing to the diversity of secondary metabolites biosynthesis, by recruiting novel genes and potentially introducing new metabolic pathways. For instance, Papaver somniferum has undergone significant tandem duplication events, which result in the emergence of morphinan and noscapine biosynthesis pathways[22,66]. The divergence and expansion of CYP genes strongly contribute to the alkaloid diversity in Coptis[67]. Tandem duplications are also common for the triterpene biosynthetic genes in Aralia elata, especially for CYP72A, CSLM, and UGT73, which may drive the diversity of triterpenoids[68]. In the case of Scutellaria, tandem duplications of the CYP82D subfamily shape the flavonoid diversification[69]. In Corydalis, each species was found to contain a particular set of BIAs, some of which are common to other species but not in the same combinations. BBEL genes have been reported to be expanded and considered to be related to the appearance of cavidines and coptisine in C. tomentella[12]. In this study, tandem duplications are also found to be evident in other genes involved in the biosynthesis of BIAs, particularly NCS, CMT, NMT, OMT, and CYP. Despite not conducting the functional experiment, we still have reason to believe that tandem duplications may drive the expansion of BIA biosynthesis genes, which is conducive to the complexity of biosynthesis pathway and further contributes the BIAs diversity in Corydalis. The incorporation of sequencing data from more Corydalis species, in conjunction with the conduction of multi-omics landscape investigations and the performance of functional experiment studies related to these recently duplicated genes in the future, may lead to an updated canvas for illustrating the genetic mechanisms underlying the diversity of BIAs.

    In this study, we proposed that C. sheareri might be a complex autotetraploid with a special AAAB karyotype. Genomic comparison detected large syntenic blocks between C. sheareri and its relatives C. tomentella, and also uncovered large-scale chromosomal structural variations (particularly inversions) between these two genomes, which might have profound effects on the divergence of Corydalis. Furthermore, we identified 172 candidate genes involved in BIAs biosynthesis in C. sheareri and traced their evolution history in Papaveraceae. We deduce that tandem duplication has played a prominent role in the expansion of the BIAs biosynthesis genes, especially for BBEL, CMT, NMT, OMT, and NCS. Moreover, the identification of the species-specific tandem duplication events implies that the growth in the number of gene members is likely to complicate the metabolic pathway, and consequently, has further contributed to the diversification of BIAs biosynthesis, ultimately leading to the diversity of BIAs in Corydalis. Our study provides more insights into the genome evolution for species-rich taxa with radiation, as well as the mechanisms underlying the diversity of the BIA biosynthetic pathway. It is also of great value for future genetic studies and medicinal applications of Corydalis.

  • The authors confirm contribution to the paper as follows: study conception and design: Liu YY, Peng D, Li JM; data collection: Liu YY, Peng D, Yu CL, Liu YJ, Chen M, Kan SL, Cao YN; analysis and interpretation of results: Liu YY, Peng D, Yu CL, Liu YJ, Kan SL, Cao YN; draft manuscript preparation: Liu YY, Peng D, Wang HW, Li JM, Peng D. All authors reviewed the results and approved the final version of the manuscript.

  • The whole-genome sequence data, including Illumina short reads, PacBio HiFi reads, Hi-C interaction reads, transcriptome data, and genome annotation files, have been deposited in The National Genomics Data Center (NGDC), under the project number: PRJCA035358.

  • We thank Dr Yuan-Yuan Feng (Institute of Botany, Chinese Academy of Sciences) for her help in the flow cytometry measurement. We also thank Mr Hai-Kuan Zhang (Berry Genomics Company) for his help in the genome assembly and annotation. This research was funded by National Natural Science Foundation of China (32000170) and Xinyang Academy of Ecological Research Open Foundation (2023XYMS05).

  • The authors declare that they have no conflict of interest.

  • Supplemental Table S1 The Various Proportions of the Epoxy Resin Systems.
    Supplemental Table S2 TGA data of PA-DAD and EP composites.
    Supplemental Table S3 Kinetic parameters of the curing reaction.
    Supplemental Table S4 The fire performances of PA-DAD.
    Supplemental Table S5 Structural assignments for the main fragments observed in the Py-GC/MS of EP/25 %PA-DAD.
    Supplemental Table S6 Mechanical properties of neat EP and EP composites.
    Supplemental Fig. S1 1H NMR spectra of (a) PA, (b) DAD, and (c) PA-DAD.
    Supplemental Fig. S2 The DSC curves of (a) neat EP, (b) EP/10% PA-DAD, (c) EP/20% PA-DAD, (d) EP/25% PA-DAD at different heating rates and the fitting curves based on Kissinger, and Ozawa methods for neat EP and EP composites.
    Supplemental Fig. S3 The UL-94 test for neat EP and EP composites.
    Supplemental Fig. S4 The THR curves vs time for EP composites at a flux of 35 kW/m2.
    Supplemental Fig. S5 (a) FTIR spectrogram of the composition of residual char; (b) Raman spectrogram of the composition of residual char; (c)XPS spectrogram of the composition of residual char.
    Supplemental Fig. S6 (a)The C1s spectra of the composition of residual char; (b) The O1s spectra of the composition of residual char; (c) The N1s spectra of the composition of residual char; (d) The P2p spectra of the composition of residual char.
    Supplemental Fig. S7 (a) The Gram-Schmidt curves of neat EP and EP/25% PA-DAD; (b) The intensity of hydrocarbons versus the temperature curves of neat EP and EP/25% PA-DAD; (c) The intensity of CO2 versus the temperature curves of neat EP and EP/25% PA-DAD; (d) The intensity of CO versus the temperature curves of neat EP and EP/25% PA-DAD; (e) The intensity of carbonyl compounds versus the temperature curves of neat EP and EP/25% PA-DAD; (f) The intensity of aromatic compounds versus the temperature curves of neat EP and EP/25% PA-DAD.
    Supplemental Fig. S8 Py-GC/MS spectra of EP/25% PA-DAD.
    Supporting Information Details of characterizations.
  • [1]

    Kong Q, Wu T, Zhang J, Wang DY. 2018. Simultaneously improving flame retardancy and dynamic mechanical properties of epoxy resin nanocomposites through layered copper phenylphosphate. Composites Science and Technology 154:136−44

    doi: 10.1016/j.compscitech.2017.10.013

    CrossRef   Google Scholar

    [2]

    Gholipour-Mahmoudalilou M, Roghani-Mamaqani H, Azimi R, Abdollahi A. 2018. Preparation of hyperbranched poly (amidoamine)-grafted graphene nanolayers as a composite and curing agent for epoxy resin. Applied Surface Science 428:1061−69

    doi: 10.1016/j.apsusc.2017.09.237

    CrossRef   Google Scholar

    [3]

    Fang F, Song P, Ran S, Guo Z, Wang H, et al. 2018. A facile way to prepare phosphorus-nitrogen-functionalized graphene oxide for enhancing the flame retardancy of epoxy resin. Composites Communications 10:97−102

    doi: 10.1016/j.coco.2018.08.001

    CrossRef   Google Scholar

    [4]

    Xu B, Zhang Q, Zhou H, Qian L, Zhao S. 2023. Small change, big impact: Simply changing the substitute on Si atom towards significant improvement of flame retardancy and toughness of epoxy resins. Composites Part B: Engineering 263:110832

    doi: 10.1016/j.compositesb.2023.110832

    CrossRef   Google Scholar

    [5]

    Xu B, Wei S, Liu Y, Wu M. 2023. Effects of different metal ions in phosphonitrile-modified organometallic complexes on flame retardancy of epoxy resin. Polymer Degradation and Stability 214:110407

    doi: 10.1016/j.polymdegradstab.2023.110407

    CrossRef   Google Scholar

    [6]

    Li D, Zhang Z, Wang S, Xu M, Li B. 2022. A monomolecular intumescent flame retardant for improvement simultaneously of fire safety, smoke suppression, and mechanical properties of epoxy resin. Journal of Applied Polymer Science 139:52104

    doi: 10.1002/app.52104

    CrossRef   Google Scholar

    [7]

    Luo F, Wu K, Wang S, Lu M. 2017. Melamine resin/graphite nanoflakes hybrids and its vacuum-assisted prepared epoxy composites with anisotropic thermal conductivity and improved flame retardancy. Composites Science and Technology 144:100−6

    doi: 10.1016/j.compscitech.2017.03.025

    CrossRef   Google Scholar

    [8]

    Huo S, Yang S, Wang J, Cheng J, Zhang Q, et al. 2020. A liquid phosphorus-containing imidazole derivative as flame-retardant curing agent for epoxy resin with enhanced thermal latency, mechanical, and flame-retardant performances. Journal of Hazardous Materials 386:121984

    doi: 10.1016/j.jhazmat.2019.121984

    CrossRef   Google Scholar

    [9]

    Huo S, Yang S, Wang J, Cheng J, Zhang Q, et al. 2020. A Liquid Phosphaphenanthrene-Derived Imidazole for Improved Flame Retardancy and Smoke Suppression of Epoxy Resin. ACS Applied Polymer Materials 2:3566−75

    doi: 10.1021/acsapm.0c00577

    CrossRef   Google Scholar

    [10]

    Xu B, Wu M, Liu Y, Wei S. 2023. Study on flame retardancy behavior of epoxy resin with phosphaphenanthrene triazine compound and organic zinc complexes based on phosphonitrile. Molecules 28:3069

    doi: 10.3390/molecules28073069

    CrossRef   Google Scholar

    [11]

    Xu B, Wei S, Liu Y, Zhao S, Qian L. 2022. Preparation of an organometallic complex based on phosphonitrile and its flame retardant application in epoxy resin. Journal of Materials Research and Technology 21:4921−39

    doi: 10.1016/j.jmrt.2022.11.099

    CrossRef   Google Scholar

    [12]

    Lou G, Rao Q, Li Q, Bai Z, He X, et al. 2023. Novel ionic complex with flame retardancy and ultrastrong toughening effect on epoxy resin. Chemical Engineering Journal 455:139334

    doi: 10.1016/j.cej.2022.139334

    CrossRef   Google Scholar

    [13]

    Xu B, Liu Y, Wei S, Zhao S, Qian L, et al. 2022. A Phosphorous-Based Bi-Functional Flame Retardant Based on Phosphaphenanthrene and Aluminum Hypophosphite for an Epoxy Thermoset. International Journal of Molecular Sciences 23:11256

    doi: 10.3390/ijms231911256

    CrossRef   Google Scholar

    [14]

    Ma W, Xu B, Shao L, Liu Y, Chen Y, et al. 2019. Synthesis of (1,4-Methylenephenylphosphinic acid) Piperazine and Its Application as a Flame Retardant in Epoxy Thermosets. Macromolecular Materials and Engineering 304:1900419

    doi: 10.1002/mame.201900419

    CrossRef   Google Scholar

    [15]

    Huo S, Song P, Yu B, Ran S, Chevali VS, et al. 2021. Phosphorus-containing flame retardant epoxy thermosets: Recent advances and future perspectives. Progress in Polymer Science 114:101366

    doi: 10.1016/j.progpolymsci.2021.101366

    CrossRef   Google Scholar

    [16]

    He W, Song P, Yu B, Fang Z, Wang H. 2020. Flame retardant polymeric nanocomposites through the combination of nanomaterials and conventional flame retardants. Progress in Materials Science 114:100687

    doi: 10.1016/j.pmatsci.2020.100687

    CrossRef   Google Scholar

    [17]

    Xue Y, Shen M, Zeng S, Zhang W, Hao L, et al. 2019. A novel strategy for enhancing the flame resistance, dynamic mechanical and the thermal degradation properties of epoxy nanocomposites. Materials Research Express 6(12):125003

    doi: 10.1088/2053-1591/ab537f

    CrossRef   Google Scholar

    [18]

    Huang G, Chen W, Wu T, Guo H, Fu C, et al. 2021. Multifunctional graphene-based nano-additives toward high-performance polymer nanocomposites with enhanced mechanical, thermal, flame retardancy and smoke suppressive properties. Chemical Engineering Journal 410:127590

    doi: 10.1016/j.cej.2020.127590

    CrossRef   Google Scholar

    [19]

    Xue Y, Feng J, Huo S, Song P, Yu B, et al. 2020. Polyphosphoramide-intercalated MXene for simultaneously enhancing thermal stability, flame retardancy and mechanical properties of polylactide. Chemical Engineering Journal 397:125336

    doi: 10.1016/j.cej.2020.125336

    CrossRef   Google Scholar

    [20]

    Zhang Y, Jing J, Liu T, Xi L, Sai T, et al. 2021. A molecularly engineered bioderived polyphosphate for enhanced flame retardant, UV-blocking and mechanical properties of poly(lactic acid). Chemical Engineering Journal 411:128493

    doi: 10.1016/j.cej.2021.128493

    CrossRef   Google Scholar

    [21]

    Xu Y, Liu L, Yan C, Hong Y, Xu M, et al. 2021. Eco-friendly phosphonic acid piperazine salt toward high-efficiency smoke suppression and flame retardancy for epoxy resins. Journal of Materials Science 56:16999−7010

    doi: 10.1007/s10853-021-06384-1

    CrossRef   Google Scholar

    [22]

    Feng J, Sun Y, Song P, Lei W, Wu Q, et al. 2017. Fire-resistant, strong, and green polymer nanocomposites based on poly(lactic acid) and core-shell nanofibrous flame retardants. ACS Sustainable Chemistry & Engineering 5:7894−904

    doi: 10.1021/acssuschemeng.7b01430

    CrossRef   Google Scholar

    [23]

    Yin W, Chen L, Lu F, Song P, Dai J, et al. 2018. Mechanically robust, flame-retardant poly(lactic acid) biocomposites via combining cellulose nanofibers and ammonium polyphosphate. ACS Omega 3:5615−26

    doi: 10.1021/acsomega.8b00540

    CrossRef   Google Scholar

    [24]

    Yang H, Shi B, Xue Y, Ma Z, Liu L, et al. 2021. Molecularly engineered lignin-derived additives enable fire-retardant, UV-shielding, and mechanically strong polylactide biocomposites. Biomacromolecules 22:1432−44

    doi: 10.1021/acs.biomac.0c01656

    CrossRef   Google Scholar

    [25]

    Liu L, Qian M, Song Pa, Huang G, Yu Y, et al. 2016. Fabrication of green lignin-based flame retardants for enhancing the thermal and fire retardancy properties of polypropylene/wood composites. ACS Sustainable Chemistry & Engineering 4:2422−31

    doi: 10.1021/acssuschemeng.6b00112

    CrossRef   Google Scholar

    [26]

    Liu L, Huang G, Song P, Yu Y, Fu S. 2016. Converting industrial alkali lignin to biobased functional additives for improving fire behavior and smoke suppression of polybutylene succinate. ACS Sustainable Chemistry & Engineering 4:4732−42

    doi: 10.1021/acssuschemeng.6b00955

    CrossRef   Google Scholar

    [27]

    Yang H, Yu B, Xu X, Bourbigot S, Wang H, et al. 2020. Lignin-derived bio-based flame retardants toward high-performance sustainable polymeric materials. Green Chemistry 22:2129−61

    doi: 10.1039/D0GC00449A

    CrossRef   Google Scholar

    [28]

    Lou G, Ma Z, Dai J, Bai Z, Fu S, et al. 2021. Fully biobased surface-functionalized microcrystalline cellulose via green self-assembly toward fire-retardant, strong, and tough epoxy biocomposites. ACS Sustainable Chemistry & Engineering 9:13595−605

    doi: 10.1021/acssuschemeng.1c04718

    CrossRef   Google Scholar

    [29]

    Mendis GP, Weiss SG, Korey M, Boardman CR, Dietenberger M, et al. 2016. Phosphorylated lignin as a halogen-free flame retardant additive for epoxy composites. Green Materials 4:150−59

    doi: 10.1680/jgrma.16.00008

    CrossRef   Google Scholar

    [30]

    Guo W, Wang X, Gangireddy CSR, Wang J, Pan Y, et al. 2019. Cardanol derived benzoxazine in combination with boron-doped graphene toward simultaneously improved toughening and flame retardant epoxy composites. Composites Part A: Applied Science and Manufacturing 116:13−23

    doi: 10.1016/j.compositesa.2018.10.010

    CrossRef   Google Scholar

    [31]

    Wang X, Zhou S, Guo WW, Wang PL, Xing W, et al. 2017. Renewable cardanol-based phosphate as a flame retardant toughening agent for epoxy resins. ACS Sustainable Chemistry & Engineering 5:3409−16

    doi: 10.1021/acssuschemeng.7b00062

    CrossRef   Google Scholar

    [32]

    Li C, Wan J, Kalali EN, Fan H, Wang DY. 2015. Synthesis and characterization of functional eugenol derivative based layered double hydroxide and its use as a nanoflame-retardant in epoxy resin. Journal of Materials Chemistry A 3:3471−79

    doi: 10.1039/C4TA05740F

    CrossRef   Google Scholar

    [33]

    Fang F, Huo S, Shen H, Ran S, Wang H, et al. 2020. A bio-based ionic complex with different oxidation states of phosphorus for reducing flammability and smoke release of epoxy resins. Composites Communications 17:104−8

    doi: 10.1016/j.coco.2019.11.011

    CrossRef   Google Scholar

    [34]

    Zhu ZM, Shang K, Wang LX, Wang JS. 2019. Synthesis of an effective bio-based flame-retardant curing agent and its application in epoxy resin: Curing behavior, thermal stability and flame retardancy. Polymer Degradation and Stability 167:179−88

    doi: 10.1016/j.polymdegradstab.2019.07.005

    CrossRef   Google Scholar

    [35]

    Wang P, Liao D, Hu X, Pan N, Li W, et al. 2019. Facile fabrication of biobased P−N−C-containing nano-layered hybrid: Preparation, growth mechanism and its efficient fire retardancy in epoxy. Polymer Degradation and Stability 159:153−62

    doi: 10.1016/j.polymdegradstab.2018.11.024

    CrossRef   Google Scholar

    [36]

    Zhang J, Mi X, Chen S, Xu Z, Zhang D, et al. 2020. A bio-based hyperbranched flame retardant for epoxy resins. Chemical Engineering Journal 381:122719

    doi: 10.1016/j.cej.2019.122719

    CrossRef   Google Scholar

    [37]

    Yang D, Dong L, Hou X, Zheng W, Xiao J, et al. 2020. Synthesis of bio-based poly (cyclotriphosphazene-resveratrol) microspheres acting as both flame retardant and reinforcing agent to epoxy resin. Polymers for Advanced Technologies 31:135−45

    doi: 10.1002/pat.4755

    CrossRef   Google Scholar

    [38]

    Li Z, Liu Z, Zhang J, Fu C, Wagenknecht U, et al. 2019. Bio-based layered double hydroxide nanocarrier toward fire-retardant epoxy resin with efficiently improved smoke suppression. Chemical Engineering Journal 378:122046

    doi: 10.1016/j.cej.2019.122046

    CrossRef   Google Scholar

    [39]

    Zhao X, Xiao D, Alonso JP, Wang DY. 2017. Inclusion complex between beta-cyclodextrin and phenylphosphonicdiamide as novel bio-based flame retardant to epoxy: Inclusion behavior, characterization and flammability. Materials & Design 114:623−32

    doi: 10.1016/j.matdes.2016.11.093

    CrossRef   Google Scholar

    [40]

    Gao YY, Deng C, Du YY, Huang SC, Wang YZ. 2019. A novel bio-based flame retardant for polypropylene from phytic acid. Polymer Degradation and Stability 161:298−308

    doi: 10.1016/j.polymdegradstab.2019.02.005

    CrossRef   Google Scholar

    [41]

    Wang D, Wang Y, Li T, Zhang S, Ma P, et al. 2020. A bio-based flame-retardant starch based on phytic acid. ACS Sustainable Chemistry & Engineering 8:10265−74

    doi: 10.1021/acssuschemeng.0c03277

    CrossRef   Google Scholar

    [42]

    Wang WZ, Zhang YH. 2010. Synthesis of semiaromatic polyamides based on decanediamine. Chinese Journal of Polymer Science 28:467−73

    doi: 10.1007/s10118-010-9049-2

    CrossRef   Google Scholar

    [43]

    Zhang T, Yan H, Shen L, Fang Z, Zhang X, et al. 2014. A phosphorus-, nitrogen- and carbon-containing polyelectrolyte complex: preparation, characterization and its flame retardant performance on polypropylene. RSC Advances 4:48285−92

    doi: 10.1039/C4RA09243K

    CrossRef   Google Scholar

    [44]

    Zhu ZM, Xu YJ, Liao W, Xu S, Wang YZ. 2017. Highly flame retardant expanded polystyrene foams from Phosphorus–Nitrogen–Silicon synergistic adhesives. Industrial & Engineering Chemistry Research 56:4649−58

    doi: 10.1021/acs.iecr.6b05065

    CrossRef   Google Scholar

    [45]

    Tan Y, Shao ZB, Yu LX, Long JW, Qi M, et al. 2016. Piperazine-modified ammonium polyphosphate as monocomponent flame-retardant hardener for epoxy resin: flame retardance, curing behavior and mechanical property. Polymer Chemistry 7:3003−12

    doi: 10.1039/C6PY00434B

    CrossRef   Google Scholar

    [46]

    Ali Lakho D, Yao D, Cho K, Ishaq M, Wang Y. 2017. Study of the curing kinetics toward development of fast-curing epoxy resins. Polymer-Plastics Technology and Engineering 56:161−70

    doi: 10.1080/03602559.2016.1185623

    CrossRef   Google Scholar

    [47]

    Laufer G, Kirkland C, Morgan AB, Grunlan JC. 2012. Intumescent multilayer nanocoating, made with renewable polyelectrolytes, for flame-retardant cotton. Biomacromolecules 13:2843−48

    doi: 10.1021/bm300873b

    CrossRef   Google Scholar

    [48]

    Gao M, Wu W, Yan Y. 2009. Thermal degradation and flame retardancy of epoxy resins containing intumescent flame retardant. Journal of Thermal Analysis and Calorimetry 95:605−4

    doi: 10.1007/s10973-008-9766-8

    CrossRef   Google Scholar

    [49]

    Yan YW, Chen L, Jian RK, Kong S, Wang YZ. 2012. Intumescence: An effect way to flame retardance and smoke suppression for polystryene. Polymer Degradation and Stability 97:1423−31

    doi: 10.1016/j.polymdegradstab.2012.05.013

    CrossRef   Google Scholar

    [50]

    Breulet H, Steenhuizen T. 2005. Fire testing of cables: comparison of SBI with FIPEC/Europacable tests. Polymer Degradation and Stability 88:150−58

    doi: 10.1016/j.polymdegradstab.2004.01.031

    CrossRef   Google Scholar

    [51]

    Xing W, Zhang P, Song L, Wang X, Hu Y. 2014. Effects of alpha-zirconium phosphate on thermal degradation and flame retardancy of transparent intumescent fire protective coating. Materials Research Bulletin 49:1−6

    doi: 10.1016/j.materresbull.2013.08.033

    CrossRef   Google Scholar

    [52]

    Yang G, Wu WH, Wang YH, Jiao YH, Lu LY, et al. 2019. Synthesis of a novel phosphazene-based flame retardant with active amine groups and its application in reducing the fire hazard of Epoxy Resin. Journal of Hazardous Materials 366:78−87

    doi: 10.1016/j.jhazmat.2018.11.093

    CrossRef   Google Scholar

    [53]

    Lou G, Pei G, Wu Y, Lu Y, Wu Y, et al. 2021. Combustion conversion of wood to N, O co-doped 2D carbon nanosheets for zinc-ion hybrid supercapacitors. Chemical Engineering Journal 413:127502

    doi: 10.1016/j.cej.2020.127502

    CrossRef   Google Scholar

    [54]

    Zhu W, Kim D, Han M, Jang J, Choi H, et al. 2023. Fibrous cellulose nanoarchitectonics on N-doped Carbon-based Metal-Free catalytic nanofilter for highly efficient advanced oxidation process. Chemical Engineering Journal 460:141593

    doi: 10.1016/j.cej.2023.141593

    CrossRef   Google Scholar

    [55]

    Luo H, Rao W, Zhao P, Wang L, Liu Y, et al. 2020. An efficient organic/inorganic phosphorus-nitrogen-silicon flame retardant towards low-flammability epoxy resin. Polymer Degradation and Stability 178:109195

    doi: 10.1016/j.polymdegradstab.2020.109195

    CrossRef   Google Scholar

    [56]

    Song J, Chen C, Zhu S, Zhu M, Dai J, et al. 2018. Processing bulk natural wood into a high-performance structural material. Nature 554:224−28

    doi: 10.1038/nature25476

    CrossRef   Google Scholar

    [57]

    Xiao L, Huang J, Wang Y, Chen J, Liu Z, et al. 2019. Tung Oil-based modifier toughening epoxy resin by sacrificial bonds. ACS Sustainable Chemistry & Engineering 7:17344−53

    doi: 10.1021/acssuschemeng.9b04284

    CrossRef   Google Scholar

    [58]

    Wang X, Peng J, Zhang Y, Li M, Saiz E, et al. 2018. Ultratough bioinspired graphene fiber via sequential toughening of hydrogen and ionic bonding. ACS Nano 12:12638−45

    doi: 10.1021/acsnano.8b07392

    CrossRef   Google Scholar

  • Cite this article

    Lin Z, Zhang W, Lou G, Bai Z, Xu J, et al. 2023. A bio-based hyperbranched flame retardant towards the fire-safety and smoke-suppression epoxy composite. Emergency Management Science and Technology 3:21 doi: 10.48130/EMST-2023-0021
    Lin Z, Zhang W, Lou G, Bai Z, Xu J, et al. 2023. A bio-based hyperbranched flame retardant towards the fire-safety and smoke-suppression epoxy composite. Emergency Management Science and Technology 3:21 doi: 10.48130/EMST-2023-0021

Figures(8)  /  Tables(2)

Article Metrics

Article views(4742) PDF downloads(597)

ARTICLE   Open Access    

A bio-based hyperbranched flame retardant towards the fire-safety and smoke-suppression epoxy composite

Emergency Management Science and Technology  3 Article number: 21  (2023)  |  Cite this article

Abstract: The design of fully biological flame retardants is vital for environment-friendly and sustainability development. Herein, a fully biological hyperbranched flame retardant named PA-DAD was synthesized successfully through a simple neutralization reaction between 1,10-Diaminodecane (DAD) and phytic acid (PA). PA-DAD as a kind of reactive flame retardant was subsequently employed to composite with epoxy resin (EP). The obtained EP composite, comprising 25 wt% PA-DAD, exhibited excellent flame resistance with an appreciative limiting oxygen index (LOI) of 28.0%, and a desired V-0 rating of UL-94. The favorable attributes stem from the evident flame-retardant characteristics of PA-DAD, manifesting particularly within the gaseous and condensed phases of combustion. Benefitting from the PA-DAD with the synergetic effect of smoke suppression and flame resistance, the EP/25% PA-DAD composite displayed highlighted reductions both in smoke production and heat release rate. In contrast with neat EP, total smoke production (TSP), peak smoke production release (pSPR), the peak heat release rate (pHRR), and the rate of fire growth (FIGRA) of the EP/25% PA-DAD composite were decreased by 49.5%, 57.0%, 72.2% and 77.8%, respectively. Moreover, the EP/25% PA-DAD composite resulted in well-preserved mechanical properties, especially enhanced toughness, compared with the neat EP. The strategies in our work provided a facile, green, and highly efficient way for fabricating high fire-retardant EP composites.

    • Epoxy resin (EP) plays a vital role in human daily life, which has been employed in different areas of adhesives, coatings, insulation materials, and structural composites for its outstanding chemical resistance, adhesion, high mechanical properties, and electrical insulation[16]. However, the highly intrinsic flammability and dense smoke of EP during combustion poses some threats to human safety and property[711]. The modification of flame retardancy for EP is imperative and indispensable[1214]. In this regard, many strategies have been developed to improve flame retardancy of EP composites. Incorporation of flame retardants (FRs) into polymeric matrix is generally one of the most efficient and economical ways to pursue excellent flame retardancy for EP[15,16]. Thus, numerous FRs have been designed to hinder fire hazards of polymers in the past decades[1721]. Nevertheless, a significant proportion of FRs for the EP matrix are originated from non-renewable petrochemical sources, which goes against the principles of global sustainability and environmental biodegradability. On this account, biological flame retardants have been attracting increasing attention for their inherent renewability, abundance, and biodegradability[2228].

      So far, numerous bio-based FRs have been prepared and investigated for imparting flame retardancy to EP, such as lignin[29], cardanol[30,31], eugenol[32], phytic acid[3335], itaconic acid[36], resveratrol[37] and cyclodextrin[38,39], and so on. Among them, phytic acid (PA) possesses high flame-retardant efficiency because of its high phosphorus content of 28 wt%. PA is a phosphorus-rich organic acid, derived from the seeds and roots of the plant, which could provide an acid source for intumescent FRs[40,41]. Wang et al.[35] prepared a nano-layered hybrid as a kind of intumescent flame retardant for EP by a neutralization reaction of melamine and PA, and the obtained composites showed both good smoke suppression and flame resistance. Fang et al.[33] improved the flame resistance of EP via synthesizing one unique phosphorus compound by the neutralization of phenyl phosphonate and PA, affecting condensed and gas phases. Zhu et al.[34] synthesized a macromolecule ammonium phytic acid via a neutralization reaction of PA and N-aminoethyl piperazine for reducing smoke emissions and enhancing flame retardancy as a curing agent for EP. However, it is noted that the amines chelated with PA reported in the literature were petroleum-based products, which did not achieve all bio-based FRs.

      In this work, a fully biological flame retardant (PA-DAD) was designed via a neutralization reaction between 1,10-diaminodecane (DAD) and PA that from renewable castor oil. The resulting PA-DAD exhibited both favourable smoke suppression and flame retardancy due to the generated bi-phase effect from PA-DAD in EP composites. The existing mechanism of flame retardancy for PA-DAD was explored using SEM-EDS, TG-IR, and Py-GC/MS methods. The transferring for heat and oxygen can be declined by the formed compact char, protecting the residual matrix from additional thermal degradation and decreasing emissions of toxic gases. The mechanical performances of EP/PA-DAD composites with a satisfactory toughness were methodically investigated. PA-DAD was demonstrated as a high-efficient green smoke suppressant as well as flame retardant for EP, expecting to contribute to low-carbon programs.

    • Dicyandiamide (DCD, 98%), Anhydrous ethanol (99.5%), and 2-methylimidazole (2-MI, 98%) all originated from Aladdin Chemistry Co., Ltd (Shanghai, China). Phytic acid aqueous solution (PA, 70wt%) was purchased by Shanghai Yuanye Bio-Technology Co., Ltd (Shanghai, China). Shanghai Macklin Biochemical Co., Ltd (Shanghai, China) supplied the 1,10-Diaminodecane (DAD, 97%). Anhui Shanfu New Material Technology Co., Ltd (Huangshan, China) supplied solid epoxy resin E12. All reagents were applied with no additional dealing.

    • As displayed in Fig. 1a, the PA-DAD was prepared in ethanol-water solution by neutralization reaction. In brief, 18.86 g of PA solution (0.02 mol) and 42.64 g of DAD (0.12 mol) were dissolved into 50 ml of deionized water and 400 ml of anhydrous ethanol, respectively. Then, the PA solution was dropwise added into the stirred DAD solution, yielding in the formation of white precipitate. Subsequently, the final product was acquired by filtration, being washed multiple times by anhydrous ethanol, and dried in a 60 °C vacuum oven for 8 h.

      Figure 1. 

      (a) Schematic diagram for the prepared processing of PA-DAD, (b) FTIR spectra of DAD, PA and PA-DAD, (c) 31P NMR curve of PA-DAD and PA, (d) TG curves of PA-DAD.

    • The different EP/PA-DAD composites were prepared as depicted in the following. The solid epoxy resin E12, along with the prepared PA-DAD, promoter 2-MI and curing agent DCD, were thoroughly mixed using a disintegrator to keep the distribution of mixture uniform. The neutralization reaction occurred at a temperature of 25 °C. The different mixing conditions are displayed in Supplemental Table S1. The above mixture was put in a mold made by stainless steel, and then cured with a press vulcanizer by hot pressing under a pressure of 10 MPa for 25 min at 145 °C. When the temperature was cooled to room temperature, the EP/PA-DAD composites to be used for different tests were obtained.

    • The characterizations employed in this work mainly include: a nuclear magnetic resonance spectrometer (NMR), Fourier Transform infrared spectroscopy (FTIR), Energy-dispersive X-ray spectroscopy (EDS), Scanning Electron Microscope (SEM), Raman spectra, X-ray photoelectron spectroscopy (XPS), Differential scanning calorimetry (DSC), Thermogravimetry Analysis (TG), Limiting Oxygen Index (LOI), Pyrolysis gas chromatography-mass spectrometry (Py-GC/MS), a cone calorimeter, the tensile strength, the impact strength, and the flexural strength and so on. The corresponding details of characterizations are supplied in the Supporting Information.

    • The FTIR curves of DAD, PA and PA-DAD were carried out as depicted in Fig. 1b. Regarding PA, the 2,400 cm−1 peak was identified as arising from the (P=O)-OH chemical bond[34]. In the curve of DAD, the 3,329 cm−1 absorption peak corresponded to the stretching vibration of N-H, while the scope for 2,800−3,000 cm−1 was related to C-H bonds[42]. By contrast, the peak assigned to N-H and (P=O)-OH of PA-DAD disappeared. Instead, the new appeared peak of 1,554 cm−1 and 1,075 cm−1 were ascribed to NH+3 and PO23, respectively[35,43]. The above change suggests that the bio-based ionic complex of PA-DAD was successfully prepared.

      Additionally, the chemical structure of PA-DAD was demonstrated via 31P NMR,1H NMR, and FTIR. As presented in Supplemental Fig. S1aS1c and Fig. 1c, the signal peak at 4.8 ppm in all spectra were ascribed to the deuterated water. In the context of PA, two distinct peaks appeared at approximately 4 ppm in the spectrum. Concerning DAD, peaks at 2.5, 1.4, and 1.2 ppm were related to the -CH2- protons in various chemical environments. When it comes to PA-DAD, the two peaks at 2.5 and 1.4 ppm exhibited a noticeable shift towards the left, indicating that a neutralization reaction occurred between DAD and PA. Moreover, in the spectra of 31P NMR, as depicted in Fig. 1b, the signal of phosphorus shifted from −0.32 ppm of PA to 1.47 ppm of PA-DAD.

      The thermal-degradation behaviors of PA-DAD are also shown in Fig. 1d and Supplemental Table 2. The process of thermal degradation for PA-DAD involved two primary stages, with the temperatures of maximum weight loss rates recorded at 227 and 399 °C. The initial degradation stage can be attributed to and the dehydration reaction from phytic acid, and the destruction of ionic bonds within PA-DAD[44]. The second degradation stage was associated with the releasing of NH3, H2O, and also the generated poly-/pyro-/ultra-phosphoric acids as PA-DAD further decomposed[40,45].

    • The curing behaviour from effects of PA-DAD on EP for was explored via non isothermal DSC testing. Supplemental Fig. S2aS2d displayed the DSC curves of EP/PA-DAD composites with different heating rates of 5, 10, 15, 20, 25 °C min−1 separately. As the heating rate increased, the temperatures of the initial peaks were elevated, and the range of curing temperature became broader because of the intensified thermal effects per unit time[46]. The reaction activation energy (Ea) had been calculated from Kissinger and Ozawa equation, as shown in Fig. 2a &b and as well summarized in Supplemental Table S3. The Ea increased significantly when the PA-DAD was incorporated into the EP. The increased presence of active sites, stemming from the amino group of PA-DAD, needed higher energy to be input for chemical reactions during the initial stages of curing. Consequently, compared to pure EP, the EP/PA-DAD systems required higher curing temperatures.

      Figure 2. 

      The calculated results of (a) Kissinger, and (b) Ozawa methods for neat EP and its different composites, (c) TG and (d) DTG curves for EP composites.

      The behaviors of thermal degradation from EP/PA-DAD composites were inspected via TG with a nitrogen atmosphere in Fig. 2c & d, and related data were summarized in Supplemental Table S2. The residual mass for EP composites was improved after imparting the PA-DAD. The yields of char for neat EP, EP/10% PA-DAD, EP/20% PA-DAD, and EP/25% PA-DAD composites increased from 11.1% (for neat EP) to 13.6%, 17.1%, and 17.2% at 800 °C, respectively. A higher residual mass indicated an optimistic impact on high flame resistance of composites. Compared with neat EP, the Tonset (the temperature at 5 wt% mass loss) and Tmax (the maximum degradation temperature) of EP/PA-DAD decreased continuously as the increased addition amount of PA-DAD. The initial thermal decomposition of EP/PA-DAD primarily resulted from the premature breakdown of PA-DAD. Consequently, PA-DAD can be decomposed as poly-/pyro-/ultra-phosphoric acids at lower temperatures, accelerating the decomposition of EP matrix. The generated poly-/pyro-/ultra-phosphoric acids facilitated crosslinking reactions and catalyzed the carbonization of EP molecules. The process could lead to the formatting for a protective and stable char layer during high-temperature stages. The protective layer significantly enhanced the fire retardancy for EP composites during combustion, providing a safeguard for the underlying matrix[47,48].

    • Both LOI values as well as the UL-94 rating were initially employed to assess the flame retardancy for EP composites, and corresponding data and results are summarized in Table 1. Neat EP, which is inherently flammable, exhibited a quite low LOI value of 20.4% and did not meet V-0 rating in the test of UL-94. Upon adding PA-DAD into EP, there were significant increases in LOI values of EP composites. When the content of PA-DAD reached 25 wt%, a notably high LOI value with 28.0% for the EP/25% PA-DAD composite was achieved and successfully satisfied V-0 rating in the test of UL-94 (as indicated in Supplemental Fig. S3). The above findings demonstrated that PA-DAD served as an efficient biological flame retardant for EP, enhancing fire resistance of composites.

      Table 1.  Detailed results corresponding to UL-94 and LOI for EP composites.

      SampleUL-94LOI (%)
      RatingDrippingt1/t2(s)
      Neat EPNRYES> 3020.4
      EP/10%PA-DADNRNO> 3023.6
      EP/20%PA-DADV-1NO8.0/2.425.3
      EP/25%PA-DADV-0NO1.1/0.928.0

      Cone calorimetry tests were performed to assess the actual fire performance for EP composites, and the associated results are reviewed in Table 2. EP/PA-DAD composites' time to ignition (TTI) was found to be shorter than that for pure EP. Because an earlier decomposition for PA-DAD would form a stable, and protective char layer, keeping the EP matrix from further thermal degradation. As depicted in Fig. 3a, the highly flammable EP would burn out rapidly after ignition, and its peak heat release rate (pHRR) was 1,007 kW/m2. Reversely, after introducing 25 wt% PA-DAD to EP, the pHRR of the composites has been reduced to 280 kW/m2, decreasing by 72.2% than that of neat EP. Moreover, as shown in Supplemental Fig. S4, the total heat release (THR) reduced significantly from 99.2 MJ/m2 of pure EP to 71.2 MJ/m2 of EP/25% PA-DAD, which had decreased by 28.3%. Besides, the rate of fire growth (FIGRA) was the ratio between pHRR and TpHRR (time to pHRR) according to the HRR curves to evaluate the actual fire hazard for material qualitatively[49,50]. The FIGRA value (1.38 kW/m2·s) of EP/25% PA-DAD was 77.8% lower than that of neat EP with 6.21 kW/m2·s, underscoring the superior fire safety enhancement achieved by PA-DAD in the EP composites. Additionally, the smoke released during combustion was generally recognized as another critical factor in evaluating the fire retardancy of the material. In Fig. 3b & c, the total smoke product (TSP) and peak of smoke product rate (pSPR) values of EP/PA-DAD composites were inhibited significantly. The values of THR and pHRR for the EP/25% PA-DAD composite experienced noteworthy reductions of 49.5% and 57.0%, respectively, in contrast with neat EP. In detail, as presented in Fig. 3d & e, the peak of CO production named pCOP, and the peak of CO2 production named pCO2P for EP/25% PA-DAD decreased by 66.5%, and 78.9%, respectively. In the event of a fire, fatalities can occur due to suffocation when the concentrations of carbon dioxide (CO2), and carbon monoxide (CO) in the surrounding air reach critical levels. Thus, the introduced PA-DAD to EP dramatically improved the fire safety and provided more escaped time for people. Moreover, lower production of CO and CO2 disclosed the more carbon fixed char residue during combustion, contributing to exert condensed phase effect for flame-retardant mechanism. To underscore the superior attributes of the newly formulated PA-DAD, the TSP reduction, pCOP reduction of PA-DAD and previously reported flame retardant for EP are displayed in Fig. 3f. The detailed fire performances are also summarized as in Supplemental Table S4. The biological PA-DAD in this work showed a significant reduction of TSP and pCOP than those of previous reported FRs.

      Table 2.  Data from Cone calorimeter for EP and its composites.

      SamplesTTI
      (s)
      TPHRR
      (s)
      pHRR
      (kW/m2)
      THR
      (MJ/m2)
      pSPR
      (m2/s)
      TSP
      (m2)
      pCOP
      (g/s)
      pCO2P
      (g/s)
      FIGRA
      (kW/m2·s)
      Neat EP86 ± 3162 ± 51007 ± 3299.2 ± 1.90.26 ± 0.01027.7 ± 2.10.0258 ± 0.00070.586 ± 0.0136.21
      EP/10% PA-DAD71 ± 2121 ± 4740 ± 2696.1 ± 2.10.28 ± 0.01225.1 ± 1.70.0274 ± 0.00060.343 ± 0.0096.12
      EP/20% PA-DAD63 ± 2184 ± 4403 ± 2471.9 ± 0.70.12 ± 0.00616.6 ± 0.90.0125 ± 0.00050.176 ± 0.0072.19
      EP/25% PA-DAD34 ± 1203 ± 3280 ± 1971.2 ± 0.90.11 ± 0.00714.0 ± 1.10.0086 ± 0.00030.124 ± 0.0081.38

      Figure 3. 

      (a) HRR, (b) SPR, (c) TSP, (d) COP, (e) CO2P curves vs time for EP composites under a flux with 35 kW/m2, and (f) the TSP reduction vs pCOP reduction of PA-DAD compared with other works.

    • To disclose the mechanisms from flame retardants within EP matrix, the char residues of EP composites from cone calorimetry were further studied. As displayed in Fig. 4, the digital photos indicated that EP/PA-DAD composites had a dense and intumescent carbon layer. With an increasing loading of PA-DAD, the volumes of residues increased as well expanded gradually. The dense and intumescent carbon layer provided a desired shield to reduce the heat and mass transferring between condensed and gas phases. SEM-EDS was applied to evaluate the quality of char residues. As depicted in Fig. 4a1d1 & 4e, there were four elements including C, O, N and P in the char residues of EP/PA-DAD composites. And P content of the char layer was increased gradually with the increased loading of PA-DAD. The element P was introduced by PA-DAD, accelerating the catalysis and carbonization for EP to produce a good char layer in dense, continuous, and phosphorous-containing. Thus, the char layers of EP/PA-DAD composites in Fig. 4a2d2 were relatively compact compared with neat EP, enhancing the flame retardancy via condensed phase.

      Figure 4. 

      Digital photos of residual chars from cone calorimetry for (a) neat EP, (b) EP/10% PA-DAD, (c) EP/20% PA-DAD, (d) EP/25% PA-DAD. (a1)−(d1) EDS results and (a2)−(d2) SEM photos for corresponding residue char. (e) The images of elemental mapping for residue chars of EP/25% PA-DAD.

      The compositions of char residues were also analyzed by FTIR methods. In Supplemental Fig. S5a, the absorption peaks at around 1,618, and 2,852−2,958 cm−1 were related to C=C, and C-H' stretching vibration from CH2 group, respectively[33,43]. For EP/25% PA-DAD, a new 1,160 cm−1 absorption peak indicated the presence of P-O bonds in compound. Two typical peaks of G band at 1,580 cm−1 and D band at 1,350 cm−1 were presented in Raman spectra of char residues for EP and EP/25% PA-DAD composite, which corresponded to graphite char and the disordered char, respectively (see Supplemental Fig. S5b)[51]. The value of ID/IG for EP/25% PA-DAD char was 1.11, higher than that of char from neat EP (0.89), suggesting a smaller microstructure size that enhances fire resistance[52]. As displayed in Supplemental Fig. S5c, the composition of elements from the surface of EP/25% PA-DAD char was analyzed by XPS. The elements including C, N, O, and P existed in the chemical structure of EP/25% PA-DAD char, which was consistent with the results of EDS characterization. The high-solution of elements was shown in Supplemental Fig. S6aS6d. Three peaks of C-C in aliphatic and aromatic species at 284.5 eV, C=O at 288.5 eV, and C-O in ether and/or hydroxyl group at 286.0 eV were deconvoluted from C1s spectra. In the high-solution O1s, three kinds of oxygen configurations appear at 535.1, 532.1 and 530.4 eV, corresponding to -COOH, -O- of C-O-C or C-O-P and =O from phosphate or carbonyl groups, respectively. The N1s spectra could be divided into two peaks of 400.2 and 399.2 eV, assigning to nitrogen functionality in the pyrrolic group and oxidized nitrogen compounds[53,54]. The only 133.5 eV peak from P2p spectra was ascribed to pyrophosphate and polyphosphate. The results of XPS demonstrated that poly-/pyro-/ultra-phosphoric acids were generated while the combustion occurred, forming a strong char layer via dehydration and esterification of EP.

      The TG-IR was carried out to explore the pyrolysis products of PA-DAD, pure EP and EP/25%PA-DAD composite under a heating rate with 20 °C/min and N2 atmosphere. Figure 5a showed the FTIR 3D and 2D for decomposed fragments from PA-DAD at 200−500 °C when thermal degradation. The primary gaseous products of PA-DAD were hydrocarbons at 2970 cm−1, CO2 at 2360 cm−1, PO· at 1258 cm−1, P-O-C at 1097 cm−1 and NH3 at both 965 and 930 cm−1. The fragments with P-O-C and PO· were demonstrated by a strong signal at 450−500 °C, being cable of quenching active radicals with OH· and H· during burning, thereby suppressing the combustion. The absorption peak for NH3 was displayed between 400−500 °C, which contributed to the dilution of flue gases and a reduction in the rate of burning. In Fig. 5b, both neat EP and EP/25% PA-DAD displayed the analogous characteristic peaks from gaseous products when degradation, such as aromatic compounds at 1,605, 1,510, 834 and 689 cm−1, hydrocarbons at 2,800−3,000 cm−1, CO2 at 2,360 cm−1, H2O at 3,600−4,000 cm−1.

      Figure 5. 

      (a) FTIR 3D and 2D spectra from pyrolysis products for PA-DAD. (b) FTIR 3D and 2D spectra from the pyrolysis products for pure EP and EP/25% PA-DAD with the maximum rate of decomposition.

      Nevertheless, the EP/25% PA-DAD showed the weaker peak intensities compared with pure EP, indicating the incorporation for PA-DAD could obstacle some productions of volatiles effectively during burning. In addition, the curves of total volatile gas intensity versus temperatures, and the Gram-Schmidt (GS) curves for pure EP and EP/25% PA-DAD composite are depicted in Supplemental Fig. S7. The EP/25% PA-DAD composite' intensities of pyrolytic products were lower than that of neat EP consistently. Some reasons for the dramatical reduction of volatile pyrolysis products from EP/25% PA-DAD are as follows. At the beginning of combustion, the earlier thermal decomposition for PA-DAD could lead to the carbonization layers with nano-scale, acting as a block to inhibit the transferring of heat and mass. When combustion continues, the PA-DAD generated poly-/pyro-/ultra-phosphoric acids for catalyzing the degradation of EP, forming intumescent and compact char for protecting the matrix from further thermal degradation.

    • The pyrolysis products obtained from EP/PA-DAD composites were investigated by the method of Py-GC/MS further. The main pyrolysis products from the pyrogram (in Supplemental Fig. S8) for EP/25% PA-DAD, were also listed in Supplemental Table S5. The 2, 3, 5, 6, and 10 peaks in the pyrogram were attributed to phenol, p-Cresol, 4-isopropyl phenol, 4-isopropenylphenol, and 2,2-bis(4-hydroxyphenyl) propane. The aforementioned substances stemmed from the reactions involving the opening of aromatic rings and the subsequent rearrangement reactions of the bisphenol A chains. The nitrogen-containing compounds of (E)-Alpha-cyanocinnamamide at peak 8, have the capacity to decompose into non-flammable, small-molecule gas NH3. Above decomposition process effectively reduced the concentration of the flammable gases[55]. The detailed pyrolysis routes and speculative flame-retardant mechanism for EP/25% PA-DAD are shown in Fig. 6. During the burning process, the decomposition of PA-DAD can yield products with P-O-C, PO·, and NH3, thereby capturing the active radicals of H·, OH· and decreasing the flammable gases. When combustion continues, the residual phosphorus from condensed generated the poly-/pyro-/ultra-phosphoric acids, catalyzing the degradation of EP for promoting the producing of the intumescent and compact chars. Herein, the formation of intumescent and compact chars reduced the oxygen, and toxic gas emissions to prevent the matrix from continuously degrading. Thus, the above results about incorporated PA-DAD further demonstrate the excellent flame retardancy and enhanced residual char of composites.

      Figure 6. 

      Speculated pyrolysis and flame-retardant mechanisms for composites with PA-DAD.

    • Mechanical performances including the tensile, flexural and impact strength for neat EP as well EP/PA-DAD composites are displayed in Figs. 7 & 8. The above relatively detailed data are reviewed in Supplemental Table S6. In Fig. 7a & d, the typical strain-stress for neat EP and its corresponding composites are presented. As demonstrated in Fig. 7c & f, the toughness values of pure EP and its corresponding composites were calculated by the integration values from stress-strain curves[56,57]. Obviously, the tensile, flexural strength and toughness were firstly increased and then decreased with the increasing contents of PA-DAD. Partial increased modulus can be attributed to the capacity of PA-DAD to enhance crosslink density. The increased toughness was assigned to the ion bonds formed by PA-DAD, which have the effect on dissipating energy[12,58]. The partial agglomeration of PA-DAD will be caused after the added amount reaches a certain level, leading to a noticeable decline in mechanical properties of EP composites[44].

      Figure 7. 

      (a) Tensile stress-strain curves, (b) tensile strength and modulus, (c) tensile toughness, (d) flexural stress-strain curves (e) flexural strength, (f) flexural toughness for EP and its corresponding composites.

      Figure 8. 

      SEM images from the fractured surfaces of tensile samples for (a) neat EP, (b) EP/10% PA-DAD, (c) EP/20% PA-DAD, and (d) EP/25% PA-DAD. (e) EDS mapping photograph of tensile fractured surfaces for EP/25% PA-DAD.

      The PA-DAD were agglomerated in the EP matrix as seen in the tensile fracture when the addition reached 25 wt% (see Fig. 8). The mapping photograph also proved the phosphorus agglomeration in the fracture. When the PA-DAD loading was 10 wt%, the tensile, and flexural strength of composites owned higher values compared with neat EP, which were 17.8% and 26.0% respectively. In particular, the tensile, and flexural toughness were greatly enhanced by 466% and 178% compared with neat EP. When PA-DAD loading reached 20 wt% and 25 wt%, although the tensile, and flexural strength of EP/PA-DAD composites exhibited lower values, the tensile, and flexural toughness still higher than those of neat EP. Some agglomerations caused by the relatively high filler contents led to the tensile and flexural strength decreased, but the ion bonds formed by PA-DAD have the function of dissipating energy for increasing tensile and flexural toughness. In comparison with neat EP, the impact strength (Supplemental Table S6) of EP/10% PA-DAD, EP/20% PA-DAD, and EP/25% PA-DAD were enhanced to 14.61, 13.00 and 8.53 KJ·m−2, and were also enhanced by 99.0%, 77.1%, and 16.2%, respectively.

    • In conclusion, a fully biological flame retardant named PA-DAD has been successfully prepared through a straightforward neutralization reaction of DAD and PA. When the content of PA-DAD achieved 25 wt%, LOI value for EP/25% PA-DAD composite reached up to 28.0%, and also satisfied a V-0 rating of the UL-94 test. In comparison with neat EP, the EP/25% PA-DAD displayed dramatic reductions in pHRR (72.2%), pSPR (57.0%), THR (28.3%), TSP (49.5%) and FIGRA (77.8%). The imparting of PA-DAD can also inhibit the production of volatiles sufficiently. The synergy between condensed and gas phase contributed to the outstanding flame retardancy from EP/ PA-DAD composites, beneficial to proposed flame retardant mechanisms. For the condensed phases, PA-DAD generated poly-/pyro-/ultra-phosphoric acids and facilitated the production of a protective char shield to insulate the transferring of heat and mass. While in the gas phases, the PA-DAD released the PO·, P-O-C, and NH3 and then quenched the active radicals of H· and OH·. The mechanical performances for EP composites have been enhanced to certain degrees after adding the flame retardant of PA-DAD. Our work supplied a friendly and facile approach for preparing biological flame retardant and smoke suppressant materials for high fire-safety required EP composites.

    • The authors confirm contribution to the paper as follows: study conception and design: Lin Z, Zhang W, Dai J; data collection: Lin Z, Zhang W, Lou G, Bai Z; analysis and interpretation of results: Lou G, Bai Z, Xu J, Li H, Zong Y, Chen F, Song P, Liu L, Dai J; draft manuscript preparation: Lin Z, Lou G, Chen F, Dai J. All authors reviewed the results and approved the final version of the manuscript.

    • The datasets generated and analyzed during the current study are available from the corresponding author on reasonable request.

      • This work was financially supported by the National Natural Science Foundation of Zhejiang Province (No. LY21E030001), National Natural Science Foundation of China (No. 51903222), National College Students Innovation and Entrepreneurship Training Program of China (No. 202210341055), the 'Leading geese' projects of Zhejiang Provincial Department of Science and Technology (No. 2022C03128), and the College Student Science and Technology Innovation Activity Plan of Zhejiang Province (No. 2023R412019).

      • The authors declare that they have no conflict of interest.

      • These authors contributed equally: Zhiqian Lin, Wangbin Zhang

      • Supplemental Table S1 The Various Proportions of the Epoxy Resin Systems.
      • Supplemental Table S2 TGA data of PA-DAD and EP composites.
      • Supplemental Table S3 Kinetic parameters of the curing reaction.
      • Supplemental Table S4 The fire performances of PA-DAD.
      • Supplemental Table S5 Structural assignments for the main fragments observed in the Py-GC/MS of EP/25 %PA-DAD.
      • Supplemental Table S6 Mechanical properties of neat EP and EP composites.
      • Supplemental Fig. S1 1H NMR spectra of (a) PA, (b) DAD, and (c) PA-DAD.
      • Supplemental Fig. S2 The DSC curves of (a) neat EP, (b) EP/10% PA-DAD, (c) EP/20% PA-DAD, (d) EP/25% PA-DAD at different heating rates and the fitting curves based on Kissinger, and Ozawa methods for neat EP and EP composites.
      • Supplemental Fig. S3 The UL-94 test for neat EP and EP composites.
      • Supplemental Fig. S4 The THR curves vs time for EP composites at a flux of 35 kW/m2.
      • Supplemental Fig. S5 (a) FTIR spectrogram of the composition of residual char; (b) Raman spectrogram of the composition of residual char; (c)XPS spectrogram of the composition of residual char.
      • Supplemental Fig. S6 (a)The C1s spectra of the composition of residual char; (b) The O1s spectra of the composition of residual char; (c) The N1s spectra of the composition of residual char; (d) The P2p spectra of the composition of residual char.
      • Supplemental Fig. S7 (a) The Gram-Schmidt curves of neat EP and EP/25% PA-DAD; (b) The intensity of hydrocarbons versus the temperature curves of neat EP and EP/25% PA-DAD; (c) The intensity of CO2 versus the temperature curves of neat EP and EP/25% PA-DAD; (d) The intensity of CO versus the temperature curves of neat EP and EP/25% PA-DAD; (e) The intensity of carbonyl compounds versus the temperature curves of neat EP and EP/25% PA-DAD; (f) The intensity of aromatic compounds versus the temperature curves of neat EP and EP/25% PA-DAD.
      • Supplemental Fig. S8 Py-GC/MS spectra of EP/25% PA-DAD.
      • Supporting Information Details of characterizations.
      • Copyright: © 2023 by the author(s). Published by Maximum Academic Press on behalf of Nanjing Tech University. This article is an open access article distributed under Creative Commons Attribution License (CC BY 4.0), visit https://creativecommons.org/licenses/by/4.0/.
    Figure (8)  Table (2) References (58)
  • About this article
    Cite this article
    Lin Z, Zhang W, Lou G, Bai Z, Xu J, et al. 2023. A bio-based hyperbranched flame retardant towards the fire-safety and smoke-suppression epoxy composite. Emergency Management Science and Technology 3:21 doi: 10.48130/EMST-2023-0021
    Lin Z, Zhang W, Lou G, Bai Z, Xu J, et al. 2023. A bio-based hyperbranched flame retardant towards the fire-safety and smoke-suppression epoxy composite. Emergency Management Science and Technology 3:21 doi: 10.48130/EMST-2023-0021

Catalog

  • About this article

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return