1
|
Li K, Kikugawa G, Kawagoe Y, Zhao Y, Okabe T. Determination of interaction parameters in a bottom-up approach employed in reactive dissipative particle dynamics simulations for thermosetting polymers. SOFT MATTER 2024; 20:4591-4607. [PMID: 38805009 DOI: 10.1039/d3sm01743e] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/29/2024]
Abstract
The limitations in previous dissipative particle dynamics (DPD) studies confined simulations to a narrow resin range. This study refines DPD parameter calculation methodology, extending its application to diverse polymer materials. Using a bottom-up approach with molecular dynamics (MD) simulations, we evaluated solubility parameters and bead number density governing nonbonded interactions via the Flory-Huggins parameter and covalent-bonded interactions. Two solubility parameter methods, Hildebrand and Krevelen-Hoftyzer, were compared for DPD simulations. The Hildebrand method, utilizing MD simulations, demonstrates higher consistency and broader applicability in determining solubility parameters for all DPD particles. The DPD/MD curing reaction process was examined in three epoxy systems: DGEBA/4,4'-DDS, DGEBA/MPDA and DGEBA/DETA. Calculations for the curing profile, gelation point, radial distribution function and branch ratio were performed. Compared to MD data for DGEBA/4,4'-DDS, the maximum deviation in secondary reactions between epoxy and amine groups according to DPD simulations with Krevelen-Hoftyzer was 14.8%, while with the Hildebrand method, it was 1.7%. The accuracy of the DPD curing reaction in reproducing the structural properties verifies its expanded application to general polymeric material simulations. The proposed curing DPD simulations, with a short run time and minimal computational resources, contributes to high-throughput screening for optimal resins and investigates mesoscopic inhomogeneous structures in large resin systems.
Collapse
Affiliation(s)
- Kaiwen Li
- Department of Finemechanics, Graduate School of Engineering, Tohoku University, 2-1-1 Katahira, Aoba-ku, Sendai 980-8577, Japan
- Institute of Fluid Science, Tohoku University, 2-1-1 Katahira, Aoba-ku, Sendai 980-8577, Japan.
| | - Gota Kikugawa
- Institute of Fluid Science, Tohoku University, 2-1-1 Katahira, Aoba-ku, Sendai 980-8577, Japan.
| | - Yoshiaki Kawagoe
- Department of Aerospace Engineering, Graduate School of Engineering, Tohoku University, 6-6-01, Aramaki Aza Aoba, Aoba-ku, Sendai, Miyagi, 980-8597, Japan
| | - Yinbo Zhao
- School of Aerospace Engineering and Applied Mechanics, Tongji University, Shanghai, 200092, PR China
| | - Tomonaga Okabe
- Department of Aerospace Engineering, Graduate School of Engineering, Tohoku University, 6-6-01, Aramaki Aza Aoba, Aoba-ku, Sendai, Miyagi, 980-8597, Japan
- Department of Materials Science and Engineering, University of Washington, BOX 352120, Seattle, WA 98195-1750, USA
- Research Center for Structural Materials, Polymer Matrix Hybrid Composite Materials Group, National Institute for Materials Science, 1-2-1 Sengen, Tsukuba, Ibaraki 305-0047, Japan
| |
Collapse
|
2
|
Kawagoe Y, Kikugawa G, Shirasu K, Kinugawa Y, Okabe T. Dissipative Particle Dynamics Simulation for Reaction-Induced Phase Separation of Thermoset/Thermoplastic Blends. J Phys Chem B 2024; 128:2018-2027. [PMID: 38373192 PMCID: PMC10911110 DOI: 10.1021/acs.jpcb.3c07756] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/25/2023] [Revised: 02/05/2024] [Accepted: 02/08/2024] [Indexed: 02/21/2024]
Abstract
Reaction-induced phase separation occurs during the curing reaction when a thermoplastic resin is dissolved in a thermoset resin, which enables toughening of the thermoset resin. As resin properties vary significantly depending on the morphology of the phase-separated structure, controlling the morphology formation is of critical importance. Reaction-induced phase separation is a phenomenon that ranges from the chemical reaction scale to the mesoscale dynamics of polymer molecules. In this study, we performed curing simulations using dissipative particle dynamics (DPD) coupled with a reaction model to reproduce reaction-induced phase separation. The curing reaction properties of the thermoset resin were determined by ab initio quantum chemical calculations, and the DPD parameters were determined by all-atom molecular dynamics simulations. This enabled mesoscopic simulations, including reactions that reflect the intrinsic material properties. The effects of the thermoplastic resin concentration, molecular weight, and curing conditions on the phase-separation morphology were evaluated, and the cure shrinkage and stiffness of each cured resin were confirmed to be consistent with the experimental trends. Furthermore, the local strain field under tensile deformation was visualized, and the inhomogeneous strain field caused by the phase-separated structures of two resins with different stiffnesses was revealed. These results can aid in understanding the toughening properties of thermoplastic additives at the molecular level.
Collapse
Affiliation(s)
- Yoshiaki Kawagoe
- Department
of Aerospace Engineering, Tohoku University, 6-6-01, Aoba, Aramaki, Aoba-ku, Sendai, Miyagi 980-8579, Japan
| | - Gota Kikugawa
- Institute
of Fluid Science, Tohoku University, 2-1-1, Katahira, Aoba-ku, Sendai, Miyagi 980-8577, Japan
| | - Keiichi Shirasu
- Department
of Finemechanics, Tohoku University, 6-6-01, Aoba, Aramaki, Aoba-ku, Sendai, Miyagi 980-8579, Japan
| | - Yuuki Kinugawa
- Department
of Aerospace Engineering, Tohoku University, 6-6-01, Aoba, Aramaki, Aoba-ku, Sendai, Miyagi 980-8579, Japan
| | - Tomonaga Okabe
- Department
of Aerospace Engineering, Tohoku University, 6-6-01, Aoba, Aramaki, Aoba-ku, Sendai, Miyagi 980-8579, Japan
- Department
of Materials Science and Engineering, University
of Washington, P.O. Box 352120, Seattle, Washington 98195-1750, United States
- Research
Center for Structural Materials, Polymer Matrix Hybrid Composite Materials
Group, National Institute for Materials
Science, 1-2-1 Sengen, Tsukuba, Ibaraki 305-0047, Japan
| |
Collapse
|
3
|
Tomiyoshi Y, Oya Y, Kawakatsu T, Okabe T. Reaction-induced morphological transitions in a blend of diblock copolymers and reactive monomers: dissipative particle dynamics simulation. SOFT MATTER 2023; 20:124-132. [PMID: 38054239 DOI: 10.1039/d3sm00959a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/07/2023]
Abstract
The dissipative particle dynamics (DPD) method is applied to the morphological transitions of microphase-separated domains in a mixture of symmetric AB-diblock copolymers and reactive C-monomers, where polymerization and cross-linking reactions take place among C-monomers. The initial structure for the DPD simulation is an equilibrated cylindrical domain structure prepared by the density-biased Monte Carlo method with density profiles obtained from the self-consistent field theory. By introducing a cross-linking reaction among reactive C-monomers, we confirmed that the DPD simulation reproduces the morphological transitions observed in experiments, where the domain morphology changes due to segregation between A-blocks of diblock copolymers and cross-linking networks of C-monomers. When the cross-linking reaction of C-monomers is sufficiently fast compared to the deformation of the domains, the initial cylindrical domains are preserved, while the distance between the domains increases. On the other hand, when the formation of the cross-linking network is slow, the domains can deform and reconnect with each other in the developing cross-linking network. In this case, we observe morphological transitions from the initial domain morphology with a large-curvature interface to another domain morphology with a smaller-curvature interface, such as the transition from the cylindrical phase to the lamellar phase. We calculated the spatial correlations in the microphase-separated domains and found that such correlations are affected by the speed of the formation of the cross-linking network depending on whether the bridging between microphase-separated domains occurs in a nucleation and growth process or in a spinodal decomposition process.
Collapse
Affiliation(s)
- Yoshinori Tomiyoshi
- Center for Soft Matter Physics, Ochanomizu University, Bunkyo-ku, Tokyo 112-8610, Japan.
| | - Yutaka Oya
- Department of Materials Science and Technology, Tokyo University of Science, Katsushika-Ku, 125-8585, Tokyo, Japan
| | - Toshihiro Kawakatsu
- Department of Physics, Graduate School of Science, Tohoku University, Sendai 980-8578, Japan
| | - Tomonaga Okabe
- Department of Aerospace Engineering, Graduate School of Engineering, Tohoku University, Sendai 980-8578, Japan
| |
Collapse
|
4
|
Shu T, Mitra G, Alberts J, Viana MP, Levy ED, Hocky GM, Holt LJ. Mesoscale molecular assembly is favored by the active, crowded cytoplasm. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.09.19.558334. [PMID: 37781612 PMCID: PMC10541124 DOI: 10.1101/2023.09.19.558334] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/03/2023]
Abstract
The mesoscale organization of molecules into membraneless biomolecular condensates is emerging as a key mechanism of rapid spatiotemporal control in cells1. Principles of biomolecular condensation have been revealed through in vitro reconstitution2. However, intracellular environments are much more complex than test-tube environments: They are viscoelastic, highly crowded at the mesoscale, and are far from thermodynamic equilibrium due to the constant action of energy-consuming processes3. We developed synDrops, a synthetic phase separation system, to study how the cellular environment affects condensate formation. Three key features enable physical analysis: synDrops are inducible, bioorthogonal, and have well-defined geometry. This design allows kinetic analysis of synDrop assembly and facilitates computational simulation of the process. We compared experiments and simulations to determine that macromolecular crowding promotes condensate nucleation but inhibits droplet growth through coalescence. ATP-dependent cellular activities help overcome the frustration of growth. In particular, actomyosin dynamics potentiate droplet growth by reducing confinement and elasticity in the mammalian cytoplasm, thereby enabling synDrop coarsening. Our results demonstrate that mesoscale molecular assembly is favored by the combined effects of crowding and active matter in the cytoplasm. These results move toward a better predictive understanding of condensate formation in vivo.
Collapse
Affiliation(s)
- Tong Shu
- Institute for Systems Genetics, NYU Langone Medical Center, 435 E 30th Street, New York, NY 10016, USA
| | - Gaurav Mitra
- Department of Chemistry, New York University, New York, New York, USA
| | | | | | - Emmanuel D. Levy
- Department of Chemical and Structural Biology, Weizmann Institute of Science, Rehovot, Israel
| | - Glen M. Hocky
- Department of Chemistry, New York University, New York, New York, USA
- Simons Center for Computational Physical Chemistry, New York University, New York, New York, USA
| | - Liam J. Holt
- Institute for Systems Genetics, NYU Langone Medical Center, 435 E 30th Street, New York, NY 10016, USA
| |
Collapse
|
5
|
Livraghi M, Pahi S, Nowakowski P, Smith DM, Wick CR, Smith AS. Block Chemistry for Accurate Modeling of Epoxy Resins. J Phys Chem B 2023; 127:7648-7662. [PMID: 37616478 PMCID: PMC10493980 DOI: 10.1021/acs.jpcb.3c04724] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/14/2023] [Revised: 07/26/2023] [Indexed: 08/26/2023]
Abstract
Accurate molecular modeling of the physical and chemical behavior of highly cross-linked epoxy resins at the atomistic scale is important for the design of new property-optimized materials. However, a systematic approach to parametrizing and characterizing these systems in molecular dynamics is missing. We therefore present a unified scheme to derive atomic charges for amine-based epoxy resins, in agreement with the AMBER force field, based on defining reactive fragments─blocks─building the network. The approach is applicable to all stages of curing from pure liquid to gelation to fully cured glass. We utilize this approach to study DGEBA/DDS epoxy systems, incorporating dynamic topology changes into atomistic molecular dynamics simulations of the curing reaction with 127,000 atoms. We study size effects in our simulations and predict the gel point utilizing a rigorous percolation theory to recover accurately the experimental data. Furthermore, we observe excellent agreement between the estimated and the experimentally determined glass transition temperatures as a function of curing rate. Finally, we demonstrate the quality of our model by the prediction of the elastic modulus based on uniaxial tensile tests. The presented scheme paves the way for a broadly consistent approach for modeling and characterizing all amine-based epoxy resins.
Collapse
Affiliation(s)
- Mattia Livraghi
- Friedrich-Alexander-Universität
Erlangen-Nürnberg (FAU), Institute for Theoretical Physics,
PULS Group, Interdisciplinary Center for Nanostructured Films (IZNF), Cauerstrasse 3, Erlangen 91058, Germany
| | - Sampanna Pahi
- Friedrich-Alexander-Universität
Erlangen-Nürnberg (FAU), Institute for Theoretical Physics,
PULS Group, Interdisciplinary Center for Nanostructured Films (IZNF), Cauerstrasse 3, Erlangen 91058, Germany
| | - Piotr Nowakowski
- Group
for Computational Life Sciences, Division of Physical Chemistry, Ruđer Bošković Institute, Bijenička cesta 54, Zagreb 10000, Croatia
| | - David M. Smith
- Group
for Computational Life Sciences, Division of Physical Chemistry, Ruđer Bošković Institute, Bijenička cesta 54, Zagreb 10000, Croatia
| | - Christian R. Wick
- Friedrich-Alexander-Universität
Erlangen-Nürnberg (FAU), Institute for Theoretical Physics,
PULS Group, Interdisciplinary Center for Nanostructured Films (IZNF), Cauerstrasse 3, Erlangen 91058, Germany
| | - Ana-Sunčana Smith
- Friedrich-Alexander-Universität
Erlangen-Nürnberg (FAU), Institute for Theoretical Physics,
PULS Group, Interdisciplinary Center for Nanostructured Films (IZNF), Cauerstrasse 3, Erlangen 91058, Germany
- Group
for Computational Life Sciences, Division of Physical Chemistry, Ruđer Bošković Institute, Bijenička cesta 54, Zagreb 10000, Croatia
| |
Collapse
|
6
|
Addressing diffusion behavior and impact in an epoxy-amine cure system using molecular dynamics simulations. Sci Rep 2023; 13:138. [PMID: 36599868 PMCID: PMC9813372 DOI: 10.1038/s41598-022-26835-2] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/22/2022] [Accepted: 12/21/2022] [Indexed: 01/06/2023] Open
Abstract
To deepen understanding of diffusion-controlled crosslinking, molecular dynamics (MD) simulations are carried out by taking the diffusion image of 3,3'-diamino diphenyl sulfone (3,3'-DDS) and polyethersulfone (PES) with epoxy resin varying temperatures from 393.15 to 473.15 K over crosslinking conversion of 0-85%. The diffusion of PES and 3,3'-DDS into the bulk increased with increasing the temperature as a result of enhanced mobility of the molecules when the difference between the glass-transition temperature (Tg) and the curing temperature. Beyond the onset points of the converged crosslinking conversion ratio of 3,3'-DDS and PES, their diffusion properties are obviously restricted with crosslinking conversion ratio. At low crosslinking conversion ratios (> 10%), the diffusion coefficients of triglycidyl p-aminophenol (TGAP) were 1.1 times higher than those of diglycidyl ether of bisphenol F (DGEBF) because of the lower molecular weight of TGAP. On the other hand, the diffusion coefficients of TGAP decreased when the crosslinking ratio was up to ~ 60% because, compared with DGEBF, it had more functional groups available to react with the curing agent. At higher crosslinking ratios, the diffusion coefficients of both resins converged to zero as a result of their highly crosslinked structures.
Collapse
|
7
|
Kawagoe Y, Kikugawa G, Shirasu K, Okabe T. Thermoset resin curing simulation using quantum-chemical reaction path calculation and dissipative particle dynamics. SOFT MATTER 2021; 17:6707-6717. [PMID: 34169305 DOI: 10.1039/d1sm00600b] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
Thermoset resin, which is commonly used as a matrix in carbon-fiber-reinforced plastic, requires curing procedures. We propose a curing simulation technique involving a dissipative particle dynamics (DPD) simulation, which can simulate a larger system and longer time period than those of conventional all-atom molecular dynamics (AA-MD) simulations. The proposed curing DPD simulation can represent the thermoset resin exothermic reaction process precisely by considering each reactivity according to the reaction types calculated via quantum-chemical reaction path calculations. The cure reaction process given by the curing DPD simulation agrees well with that given by a conventional curing AA-MD simulation, but with run-time and computational-resource reductions of 1/480 and 1/10 times, respectively. We also conduct reverse mapping, through which the AA-MD system can be reconstructed from the DPD system, to evaluate the structural and thermomechanical properties. The X-ray diffraction pattern and thermomechanical properties of the reconstructed system agree well with those of the systems derived from the curing AA-MD simulation and experimental setup. Therefore, a cured-resin AA-MD system can be obtained from a curing DPD simulation at an extremely low computational cost, and the thermomechanical properties can be evaluated precisely using this system. The proposed curing simulation technique can be applied in high-throughput screening for better materials properties and in large system calculations.
Collapse
Affiliation(s)
- Yoshiaki Kawagoe
- Department of Aerospace Engineering, Tohoku University, Sendai 980-8579, Japan.
| | - Gota Kikugawa
- Institute of Fluid Science, Tohoku University, Sendai 980-8577, Japan.
| | - Keiichi Shirasu
- Department of Aerospace Engineering, Tohoku University, Sendai 980-8579, Japan.
| | - Tomonaga Okabe
- Department of Aerospace Engineering, Tohoku University, Sendai 980-8579, Japan. and Department of Materials Science and Engineering, University of Washington, Seattle, Washington, USA
| |
Collapse
|
8
|
van der Haven DLH, Köhler S, Schreiner E, In 't Veld PJ. Closed-Form Coexistence Equation for Phase Separation of Polymeric Mixtures in Dissipative Particle Dynamics. J Phys Chem B 2021; 125:7485-7498. [PMID: 34196184 DOI: 10.1021/acs.jpcb.0c11274] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
To date, no extensive study of the phase diagram for binary fluid mixtures in dissipative particle dynamics (DPD) has been published. This is especially pertinent for newer parameterization schemes where the self-self interaction, or the effective volume, of different particle types is varied. This work presents an exhaustive study of the parameter space concerning DPD particles with soft interaction potentials. Moreover, we propose a closed-form coexistence equation or binodal curve that is inspired by the Flory-Huggins model. This equation describes the phase diagram of all binary mixtures made up out of monomers, homopolymers, and the mixtures thereof when self-self interactions are varied. The mean absolute percentage error (MAPE) of the equation on simulated data, including validation simulations, is 1.02%. The equation can a priori predict the phase separation of mixtures using only DPD interaction parameters. The proposed coexistence equation can therefore be used to directly validate interaction parameters resulting from novel parameterization schemes, including coarse graining and equations of state, without the need for additional simulations. Finally, it is shown that the choice of bond potential markedly influences phase behavior.
Collapse
Affiliation(s)
- Dingeman L H van der Haven
- Department of Chemical Engineering and Chemistry, Eindhoven University of Technology, Eindhoven 5600 MB, The Netherlands
| | - Stephan Köhler
- Polymer Physics, BASF SE, Ludwigshafen am Rhein 67056, Germany
| | | | | |
Collapse
|
9
|
Henry MM, Thomas S, Alberts M, Estridge CE, Farmer B, McNair O, Jankowski E. General-Purpose Coarse-Grained Toughened Thermoset Model for 44DDS/DGEBA/PES. Polymers (Basel) 2020; 12:polym12112547. [PMID: 33143261 PMCID: PMC7693565 DOI: 10.3390/polym12112547] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/30/2020] [Revised: 10/23/2020] [Accepted: 10/23/2020] [Indexed: 12/26/2022] Open
Abstract
The objective of this work is to predict the morphology and material properties of crosslinking polymers used in aerospace applications. We extend the open-source dybond plugin for HOOMD-Blue to implement a new coarse-grained model of reacting epoxy thermosets and use the 44DDS/DGEBA/PES system as a case study for calibration and validation. We parameterize the coarse-grained model from atomistic solubility data, calibrate reaction dynamics against experiments, and check for size-dependent artifacts. We validate model predictions by comparing glass transition temperatures measurements at arbitrary degree of cure, gel-points, and morphology predictions against experiments. We demonstrate for the first time in molecular simulations the cure-path dependence of toughened thermoset morphologies.
Collapse
Affiliation(s)
- Michael M. Henry
- Micron School of Materials Science and Engineering, Boise State University, Boise, ID 83725, USA; (M.M.H.); (S.T.); (M.A.)
| | - Stephen Thomas
- Micron School of Materials Science and Engineering, Boise State University, Boise, ID 83725, USA; (M.M.H.); (S.T.); (M.A.)
| | - Mone’t Alberts
- Micron School of Materials Science and Engineering, Boise State University, Boise, ID 83725, USA; (M.M.H.); (S.T.); (M.A.)
| | | | | | - Olivia McNair
- The Boeing Company, St. Louis, MO 63134, USA; (C.E.E.); (O.M.)
| | - Eric Jankowski
- Micron School of Materials Science and Engineering, Boise State University, Boise, ID 83725, USA; (M.M.H.); (S.T.); (M.A.)
- Correspondence:
| |
Collapse
|
10
|
Pervaje AK, Tilly JC, Detwiler AT, Spontak RJ, Khan SA, Santiso EE. Molecular Simulations of Thermoset Polymers Implementing Theoretical Kinetics with Top-Down Coarse-Grained Models. Macromolecules 2020. [DOI: 10.1021/acs.macromol.9b02255] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Amulya K. Pervaje
- Department of Chemical and Biomolecular Engineering, North Carolina State University, Raleigh, North Carolina 27606, United States
| | - Joseph C. Tilly
- Department of Chemical and Biomolecular Engineering, North Carolina State University, Raleigh, North Carolina 27606, United States
| | | | - Richard J. Spontak
- Department of Chemical and Biomolecular Engineering, North Carolina State University, Raleigh, North Carolina 27606, United States
- Department of Materials Science and Engineering, North Carolina State University, Raleigh, North Carolina 27606, United States
| | - Saad A. Khan
- Department of Chemical and Biomolecular Engineering, North Carolina State University, Raleigh, North Carolina 27606, United States
| | - Erik E. Santiso
- Department of Chemical and Biomolecular Engineering, North Carolina State University, Raleigh, North Carolina 27606, United States
| |
Collapse
|
11
|
Sarupria S. Editorial: “Introduction to the special issue on advanced molecular simulations: Methods and applications”. JOURNAL OF THEORETICAL & COMPUTATIONAL CHEMISTRY 2018. [DOI: 10.1142/s0219633618020017] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
Affiliation(s)
- Sapna Sarupria
- Department of Chemical & Biomolecular Engineering, Clemson University, Clemson, SC 29634, USA
| |
Collapse
|