Special Issue on “Computational Methods for Polymers”

Polymers play a key role in our daily lives. Natural polymers include proteins, cellulose, nucleic acids, lignin, natural rubber, and wood resins. Examples of manmade polymers are synthetic rubbers and plastics. To describe polymer properties, processes and reactions mathematically, phenomena at di ﬀ erent length and time scales should be accounted for. The length scales are from angstrom to meter, and the time scales are from femtosecond to hour. The wide ranges of the scales require the use of an appropriate modeling method or a combination of appropriate methods. Examples of the methods are: quantum chemical / mechanical methods such as density functional theory (electronic scale); molecular dynamics and Monte Carlo (molecular scale); dissipative particle dynamics, Brownian dynamics, and lattice Boltzmann method (microscopic scale); dissipative particle dynamics and ﬁeld theoretic polymer simulation (meso scale); and control volume methods and ﬁnite element (macroscopic scale) [1–3]. Special Issue used

electron-donating methyl with isosteric but electron-withdrawing chlorine substituents results in a significant increase in comonomer incorporation. Using DFT, they explored the steric and electronic properties of the ancillary ligand and found that the relative reactivity ratios are mainly determined by the electrophilicity of the metal center. They report that electronic effects observed in these catalysts affect the relative barriers for insertion and a capture-like transition state.
Lathrop et al. [9] report an experimental and theoretical/computational study on modeling and monitoring reversible addition-fragmentation chain-transfer (RAFT) polymerization of methyl methacrylate. They carry out parameter estimation to determine RAFT reaction rate coefficients and an initiator efficiency from in-situ 1 H NMR and SEC experimental data. Based on their model, they implement a multi-rate multi-delay observer that include an inter-sample predictor and a deadtime compensator. Their results show that the multi-rate multi-delay observer has satisfactory convergence after a few sampling periods.
Salas et al. [10] present a data-driven strategy for the online estimation of important kinetic parameters in copolymerization of ethylene. The kinetic parameters are chosen based on a global sensitivity analysis. The retrospective cost model refinement algorithm is adapted and implemented to estimate the kinetic parameters in real time. The results demonstrate that the estimates of the kinetic parameters converge to theoretical values of the parameters without requiring prior knowledge of the polymerization model or the theoretical parameter values.
Scott et al. [11] revisit five terpolymerization studies and estimate the 'ternary' reactivity ratios. They compare previously determined binary reactivity ratios and newly estimated 'ternary' reactivity ratios for several systems. They discuss the advantages and challenges associated with 'ternary' reactivity ratio estimation.
Garg et al. [12] present a method of uniquely specifying and robustly achieving user-specified product quality in a complex industrial batch process. They demonstrate this method using a lab-scale uni-axial rotational molding process. The method involves data-driven modeling (subspace state-space model) and model predictive control, and is able to reject raw material variation and achieve product quality which is specified through constraints on quality variables. Results from these experimental studies demonstrate the capability of the proposed method in meeting process specifications and rejecting raw material variability.
Bagalkot et al. [13] present case studies on performing polymer rapid tooling inserts and observing different failures over the life of the tool. They identify critical parameters affecting tool life, and the effect of those parameters on different areas of the tool. They categorize modes of the different failures and the underlying mechanisms including the root causes.
Intensified process units such as polymer membrane reactors pose unique challenges pertaining to design and operation that have not been fully addressed. Bishop and Lima [14] present an approach for modeling membrane reactors. Their model allows for the simulation of polymer membrane reactors under nonisothermal and countercurrent operation. This article demonstrates how operability analyses can be used to identify areas of improvement in membrane reactor design.
Topological indices have been computed for various molecular structures. They are numerical invariants associated with molecular structures and are helpful in that they feature many properties. Among these molecular descriptors, the eccentricity connectivity index is of great importance due to its relationship with pharmaceutical properties. Zheng et al. [15] calculate eccentric connectivity, total eccentricity connectivity, augmented eccentric connectivity, first Zagreb eccentricity, modified eccentric connectivity, second Zagreb eccentricity, and the edge version of eccentric connectivity indices for the molecular graph of a polyethyleneamidoamine dendrimer. They also calculate explicit representations of the polynomials associated with some of these indices.
Dendrimers are branched organic macromolecules with successive layers of branch units surrounding a central core. The molecular topology and the irregularity of their structure plays a central role in determining structural properties like enthalpy and entropy. Hussain et al. [16] determine irregularity indices that are based on the imbalance of edges, for the molecular graphs associated with some general classes of dendrimers. They also provide graphical analyses of these indices for the dendrimers.
In many applications, knowledge on the irregularity of a molecular structure is of importance. The irregularity provides quantitative information on structure-property relationships and structure-activity relationships as well as various physical and chemical properties. Zhao et al. [17] present a study on the computation and comparison of the irregularity measures of different classes of dendrimers. Their investigation focuses on four irregularity indices; the σ irregularity index, the irregularity index by Albertson, the variance of vertex degrees, and the total irregularity index.
Tembely et al. [18] present a mathematical model based on physical principles to simulate droplet impact, spreading, and eventually rebound of a viscoelastic droplet. Their simulations are based on the volume of fluid method in conjunction with a dynamic contact model accounting for the hysteresis between the droplet and substrate. They report that while the kinematic phase of droplet spreading seems to be independent of both the substrate and fluid rheology, the recoiling phase are highly influenced by those operating parameters. Their model can be used to optimize 2D/3D printing of complex fluids.
In PLP, the behavior of the degree of branching and backbiting reactions has not been well-understood yet at high laser frequencies and at relatively high reaction temperatures due to the inherent difficulties in the determination of the degree of branching of polymers. Hamzehlou et al. [19] evaluate the validity of different explanations of the recovery of PLP-molar mass distribution at high laser frequencies using a simulation study. They show that the reduction in the backbiting reaction rate at a high laser frequency, and the consequent decrease in the degree of branching, are not a necessary condition for recovering the PLP-molar mass distribution.
I would like to thank all the contributors for their contributions to this Special Issue, and Processes for supporting this issue.

Conflicts of Interest:
The author declares no conflicts of interest.