Polyethylene (PE) and polyethylene glycol (PEG) are exogenously depolymerizable. In exogenous type depolymerization processes, molecules liberate small units from their terminals. Microbial depolymerization processes of PEG have been documented in detail. Studies of PEG biodegradation include utilization of PEG of average molecular weight 20,000 by Psedonomas aeruginosa , degradation of PEG 20,000 by anaerobic bacteria isolated from sludge of a municipal anaerobic digester , and confirmation of efficient biodegradation of PEG by Pseudomonas stutzeri . Microbial depolymerization processes of PE and PEG have also been studied mathematically. A mathematical model was proposed in numerical studies of biodegradation of PE . The mathematical techniques devised for PE biodegradation processes were applied to biodegradation processes of PEG .
Molecules are broken down randomly in an endogenous type depolymerization process. Polyvinyl alcohol (PVA) and polylactic acid (PLA) are depolymerizable in endogenous type depolymerization processes. A mathematical model was proposed for enzymatic degradation process of PVA . The mathematical techniques devised for PVA were applied to an enzymatic hydrolysis of polylactic acid (PLA) . Techniques originally devised for endogenous type microbial depolymerization processes were applied to exogenous type depolymerization processes of PE and PEG .
In this study, mathematical techniques devised in previous studies were reapplied to a depolymerization process of PEG 4000. Weight distribution of PEG 4000 before and after cultivation of Sphingopyxis macrogoltabidus strain 103 were applied in inverse analyses for a molecular factor and a time factor of a degradation rate.
2. Description of Mathematical Formulation for Exogenous Type Depolymerization Process
Let w(t,M) [mg] be the weight distribution of a polymer with respect to the molecular weight M at time t, and let v(t) [mg] be the total weight of polymer molecules with molecular weight between A and B at time t. The total weight v(t) over the interval [A, B] is expressed in terms of the integral of w(t,M) from A to B,
In particular, the total weight v(t) of the entire residual polymer at time t is expressed in terms of the integral of w(t,M) over the interval [0,∞)
Integral (1) is a reasonable approximation of the integral (2) with suitable values of A and B. In this study, integral (2) over the infinite interval (2) was replaced by an integral over a finite interval (1) with A = 103.5 and B = 103.75. Similarly, an integral with the lower limit 0 was replaced by an integral with the lower limit A = 103.5, and an integral with the upper limit ∞ was replaced by an integral with the upper limit B = 103.75.
Suppose that is the population of viable cells at time t in a medium in which the polymer is the sole carbon source. The system of Equations (3), (4) for the weight distribution w(t,M) and the microbial population was proposed in previous studies -.
The parameter L is the molecular weight of a monomer unit, e.g. PE: L = 28 (CH2CH2), PEG: L = 44 (CH2CH2O), and k and h are positive parameters. The function is the molecular factor of degradation rate, and the function is the time factor of degradation rate. Equations (2) and (3) lead to the expression for the derivative of v(t),
System of Equations (3), (4) for w(t,M) and is associated with the initial conditions
where and are the initial weight distribution and the initial microbial population, respectively.
3. Inverse Problems for Molecular Factor and Time Factor of Degradation Rate and Their Solutions
Once the molecular factor and values of parameters , k, and h are set, initial value problem (3), (4), (6), (7) becomes solvable a problem for unknown functions w(t,M) and . In order to obtain the function and values of parameters , k, and h, the change of variables from to (8) was applied.
Suppose that , , and are functions which correspond to w(t,M), , and v(t), respectively, according to the change of variables (8), that is
, , .
In view of the expression
equations (3) and (4) lead to
respectively. The expression
leads to the equation
Expressions (5) and (10) lead to
Let F1(M) be the weight distribution for so that
and F2(M) be the weight distribution for so that
Equation (9), the initial condition (12), and the terminal condition (13) form an inverse problem for , for which the solution of the initial value problem (9), (12) also satisfies the terminal condition (13). Numerical techniques for the inverse problem (9), (12), (13) were devised in previous studies. Those techniques were applied to weight distributions of PEG 4000 before and after cultivation of Sphingopyxis macrogoltabidus strain 103 for one day, three days, five days, and seven days . Figure 1 shows the molecular factor of degradation rate based on weight distributions before after cultivation of the microorganism for three days.
Once the inverse problem for was solved, Equation (9) was solved for with the initial condition , where . Figure 2 shows the curve . The numerical results show that the value of t, , , , , and correspond to the value of , , , , , and .
Figure 1. This figure shows the molecular factor of degradation rate based on weight distributions before and after cultivation of Sphingopyxis macrogoltabidus strain 103 for three days.
4. Application of Numerical Techniques to Inverse Problem for Time Factor of Degradation Rate
According to a previous study, function is well approximated by an exponential function
so that . In this study, function was replaced with the exponential function (14) with a value of the parameter approximately equal to 0.574436, and equations were solved numerically.
Let be the solution of the equation (10) with initial value . The change of variables (8) leads to the expression , where
Given m pairs of values of t and , , equations for , k, and h,
where , were set. Each of intervals , , and , for , k, and h, respectively, was divided into one hundred equally spaced subintervals, and the error
was evaluated at the 1013 combinations of the nodal points, among which
Figure 2. CURVE .
gave the minimum value to the error. Figure 3 shows the curve for those values of , k, and h. In the previous study, initial value problem (3), (6) was solved numerically, once the time factor was obtained. In this study, initial value problem (3), (4), (6), (7) was solved numerically. Figure 4 shows experimental and numerical results for the transition of the total weight distributions of PEG 4000 during the cultivation of Sphingopyxis macrogoltabidus strain 103 for seven days.
Figure 3. Curve for .
Figure 4. Transition of the total weight of PEG 4000 during cultivation of Sphingopyxis macrogoltabidus strain 103 for seven days.
Mathematical techniques for depolymerization processes were demonstrated. In particular, exogenous type depolymerization process of PEG 4000 was analyzed. The weight distributions of PEG 4000 before and after cultivation of Sphingopyxis macrogoltabidus strain 103 for one day, three days, five days, and seven days are introduced into the analysis. Figure 1 shows the molecular factor based on the weight distributions before and after cultivation for three days. The inverse problem for the time factor of the degradation rate was solved (Figure 3). Once the molecular factor and values of the parameters that determine the time factor of the degradation rate were obtained, the initial value problem (3), (4), (6), (7) was solved numerically, and the exogenous type depolymerization process of PEG 4000 was simulated.
The authors thank Ms. Y. Shimizu and Dr. S. Jimbo for their technical support.
 Obradors, N. and Aguilar, J. (1991) Efficient Biodgradation of High-Molecular-Weight Polyethylene Glycols by Pure Cultures of Pseudomonas Stutzeri. Applied and Environmental Microbiology, 2383-2388.
 Watanabe, M., Kawai, F., Shibata, M., Yokoyama, S., Sudate, Y. and Hayashi, S. (2004) Analytical and Computational Techniques for Exogenous Depolymerization of Xenobiotic Polymers. Mathematical Biosciences, 192, 19-37. https://doi.org/10.1016/j.mbs.2004.06.006
 Watanabe, M. and Kawai, F. (2009) Effects of Microbial Population in Degradation Process of Xenobiotic Polymers. In: Howlett, P., Nelson, M. and Roberts, A.J., Eds., Proceedings of the 9th Biennial Engineering Mathematics and Applications Conference, EMAC-2009, Vol. 51 ANZIAM J, C682-C96. http://journal.austms.org.au/ojs/index.php/ANZIAMJ/article/view/2433 https://doi.org/10.21914/anziamj.v51i0.2433
 Watanabe, M. and Kawai, F. (2006) Mathematical Modelling and Computational Analysis of Enzymatic Degradation of Xenobiotic Polymers. Applied Mathematical Modelling, 30, 1497-1514. https://doi.org/10.1016/j.apm.2005.12.011
 Watanabe, M., Kawai, F., Tsuboi, S., Nakatsu, S. and Ohara, H. (2007) Study on Enzymatic Hydrolysis of Polylactic Acid by Endogenous Depolymerization Model. Macromolecular Theory and Simulations, 16, 619-626. https://doi.org/10.1002/mats.200700015
 Watanabe, M. and Kawai, F. (2014) Simulation for Microbial Depolymerization Processes of Polyethylene Glycol. In: Gaol, F.L., Soewito, B., Bououdina, P., Chen, M.-S., Eds., Proceedings of the 2013 International Conference on Advances in Intelligent Systems in Bioinformatics, intel-13, Advances in Intelligent Systems Research, Atlantis Press, February 2014, 1951-6851. https://www.atlantis-press.com/proceedings/intel-13/11361 https://doi.org/10.2495/intelsys130111
 Watanabe, M. and Kawai, F. (2016) Simulation of Microbial Depolymerization Process with Exponential Consumption of Carbon Source. 2016 International Conference on Computer, Mechatronics and Electronic Engineering (CMEE 2016), DEStech Publications, Inc., Lancaster, Pennsylvania, 426-431. https://doi.org/10.12783/dtcse/cmee2016/5356
 Watanabe, M. and Kawai, F. (2017) Computational Analysis for Biodegradation of Exogenously Depolymerizable Polymer. 2017 The 5th International Conference on Mechanical Engineering, Material Science and Civil Engineering, December 2017, Kuala Lumpur, Malaysia, 012015, IOP Conference Series: Materials Science and Engineering, Published Online: 6 April 2018. IOP Publishing. https://doi.org/10.1088/1757-899x/324/1/012015
 Watanabe, M. and Kawai, F. (2018) Computational Study of Bacterial Depolymerization Process of Xenobiotic Polymers. 2018 2nd International Conference on Material Engineering and Advanced Manufacturing Technology (MEAMT 2018), Article Number 02007, p. 7, Section Material Simulation. https://doi.org/10.1051/matecconf/201818902007
 Watanabe, M. and Kawai, F. (2018) Numerical Simulation of Exogenous Type Microbial Depolymerization Process with Weight Distributions and Microbial Population. Proc. of The Seventh Intl. Conf. on Advances in Civil, Structural and Environmental Engineering-ACSEE 2018, 19-22.
 Watanabe, M. and Kawai, F. (2018) Exponential Approximation of Consumption Rate of Carbon Source in Numerical Simulation of Microbial Depolymerization Process. International Journal of Advances in Science Engineering and Technology, 7.
 Watanabe, M. and Kawai, F. (2019) Mathematical Modeling and Inverse Analysis for Microbial Depolymerization Processes of Xenobiotic Polymers. Journal of Materials Science and Chemical Engineering, 7, 39-46.
 Watanabe, M. and Kawai, F. (2018) Inverse Analysis for Microbial Depolymerization Process of Polyethylen Glycol 4000. 2018 International Conference on Computer, Communications and Mechanical Engineering (CCME 2018). https://doi.org/10.12783/dtcse/ccme2018/28689