Authors:Min Ji, Xinlu Cheng, Weidong Wu Abstract: Journal of Theoretical and Computational Chemistry, Ahead of Print. The density functional theory (DFT) was used to investigate some toxic derivatives of aromatic hydrocarbons adsorption on perfect graphene (pG) and graphene-doped with B/Al/Ga (BG/AlG/GaG). And the parallel and vertical adsorptions were considered for the position relation between the adsorbent and adsorbate. The adsorption energy, adsorption distance, charge transfer and density of states (DOS) were discussed in optimized structures. The greater adsorption energy, shorter adsorption distance and more charge transfer were found in AlG by studying the four kinds of molecules (phenol/m-cresol/PCP/p-NP) adsorption on pG/BG/AlG/GaG. Then, 10 derivatives adsorption on AlG were reported, and the adsorption energy increased in the order of pentachlorophenol [math] 2,4,6-trichlorophenol [math] 2,4-dichlorophenol [math] p-cresol [math] m-cresol [math] phenol [math] o-chlorophenol [math] o-cresol [math] 2,4,6-trintrotoluene [math] para-nitrophenol. The interaction between these derivatives and the substrate was chemisorption for AlG and physisorption for pG. The oxygen atom in nitro group was more closer to the substrate than in hydroxyl group about optimized structures. Citation: Journal of Theoretical and Computational Chemistry PubDate: 2017-01-25T08:46:31Z DOI: 10.1142/S0219633617500043

Authors:Amirali Abbasi, Jaber Jahanbin Sardroodi Abstract: Journal of Theoretical and Computational Chemistry, Ahead of Print. Over the past years, an interest has arisen in resolving the problems of the increased carbon monoxide and carbon dioxide emissions, leading to the serious air pollution and many detrimental effects. A convenient solution would be a process that could utilize metal oxide nanoparticles such as TiO2 to control the concentration of atmospheric pollutants. The chemisorption of CO and CO2 molecules over the semiconductor titanium dioxide (TiO[math] is such a process. In this way, density functional theory (DFT) calculations were performed to investigate CO and CO2 adsorptions on undoped and N-doped TiO2 anatase nanoparticles. The supercell approach is conducted to construct the considered nanoparticles and the adsorption of COx molecule was simulated by use of these chosen nanoparticles. By including van der Waals (vdW) interactions between COx molecule and TiO2 nanoparticle, we found that both CO and CO2 molecules can bind strongly to the N-doped nanoparticles. The adsorption on the five-fold coordinated titanium site of TiO2 nanoparticles including the bond lengths, bond angles, adsorption energies, density of states (DOSs), Mulliken population analysis and molecular orbitals has been broadly studied in this work. Based on the obtained results, it can be concluded that the adsorption on the N-doped nanoparticle is more energetically favorable than the adsorption on the pristine one, representing the higher tendency of N-doped nanoparticles for COx detention, compared to the undoped ones. Therefore, the results indicate that the N-doped TiO2 would be an ideal COx gas sensor in the environment. Citation: Journal of Theoretical and Computational Chemistry PubDate: 2017-01-25T08:46:30Z DOI: 10.1142/S0219633617500055

Authors:Lidan Zhang, Jiguang Du, Gang Jiang Abstract: Journal of Theoretical and Computational Chemistry, Ahead of Print. Density functional theory (DFT) calculations were used to investigate the gas phase reaction of U[math] with COS to produce US[math]CO and UO[math]CS. It is shown that the two reactions are exothermic and the formation of UO[math]CS has the lower energy barrier which agrees with the experimental result that UO[math] is the main product. The reaction mechanisms and the potential energy profiles (CPEPs) considering different spin states were presented in detail. Diverse analyses including atoms in molecules, natural bond orbital were used to study the bonding properties of all the involved species. Citation: Journal of Theoretical and Computational Chemistry PubDate: 2017-01-17T11:31:38Z DOI: 10.1142/S0219633617500109

Authors:Marilena Carbone Abstract: Journal of Theoretical and Computational Chemistry, Ahead of Print. In the current study, the NO adsorption on the Si(100) surface was investigated by DFT including van der Waals forces (vdW). Stable molecular and dissociative configurations were found and compared to previous studies. Furthermore, additional states were investigated where NO adsorbs across dimers. The transformation of NO from molecular adsorbate into a dissociated adsorbate migrated into the subsurface was investigated by nudged elastic band. Several pathways were explored, either multi-staged, or direct from molecular into subsurface migrated configurations, both including and excluding vdW forces. The energy barriers of the single steps of multi-staged pathways never exceed 0.15[math]eV and are, in general, smaller when NO is adsorbed across dimers rather than bridged on a single dimer and when including vdW. Furthermore, the oxygen-bridged configurations are kinetically more accessible than the nitrogen-bridged ones. Citation: Journal of Theoretical and Computational Chemistry PubDate: 2017-01-10T06:29:17Z DOI: 10.1142/S0219633617500092

Authors:S. Ali Beyramabadi, Tina Khadivjam, Atoosa Gonabadi, Ali Morsali, Azar Gharib, Maryam Khashi, Mahdi Khorsandi-Chenarboo Abstract: Journal of Theoretical and Computational Chemistry, Ahead of Print. The Mepivacaine drug is used as a local anesthetic in dentistry, which could exist as three different tautomers. Herein, geometry, energy behavior as well as tautomerization of these tautomers have been investigated by employing density functional theory (DFT) and considering the solvent effects with the polarizable continuum model (PCM) model. The most stable tautomer of the Mepivacaine has a carbonyl and an –NH amine groups in its structure. The frontier orbitals and the energy gap of the molecule have been computed using the natural bond orbital analysis (NBO). Also, the armchair (5,5) single wall carbon nanotube (SWCNT) was used for investigation of the noncovalent interactions of the Mepivacaine molecule with the pristine SWCNT and the –COOH functionalized SWCNT in several forms. Geometries of the possible forms have been optimized. The most stable form for noncovalent interactions of the drug with each of the pristine and functionalized SWCNTs have been determined. The intermolecular H-bonds have essential role in energy behavior of the noncovalent interactions between the Mepivacaine drug and the investigated SWCNTs. Citation: Journal of Theoretical and Computational Chemistry PubDate: 2017-01-10T06:29:17Z DOI: 10.1142/S0219633617500080

Authors:Reza Ghiasi, Nasrin Sadeghi Abstract: Journal of Theoretical and Computational Chemistry, Ahead of Print. This study used mpw1pw91 quantum chemical calculations in gas and solution phases to clarify the interaction between C20 and Cr(CO)5 fragment. It also sought to clarify the effects of solvent polarity on dipole moment, structural parameters, and frontier orbital energies of the complex. Energy decomposition analysis (EDA) was applied to analyze the bonding interaction between the C20 and Cr(CO)5 fragment. Percentage composition in terms of the defined groups of frontier orbitals for the complex was evaluated to characterize the metal–ligand bonds. The Cr–C bonds within the complex were examined using quantum theory of atoms in molecules (QTAIM) analysis. In order to determine the back-bonding effects in these bonds, QTAIM analysis was applied to calculate of the quadrupole polarization of the carbon atom. Citation: Journal of Theoretical and Computational Chemistry PubDate: 2017-01-10T06:24:57Z DOI: 10.1142/S0219633617500079

Authors:Xiao-Yong Yang Abstract: Journal of Theoretical and Computational Chemistry, Volume 15, Issue 08, December 2016. By combining group representation theory with density functional theory (DFT) investigations, we reveal that (BCO)[math] systems have three nonoverlapping level sequences and the energy level diagrams in each level sequence have the same general arrangement of level. The stability of monocyclic carbon boronyl–boron carbonyl mixed clusters [(BCO)[math](CBO)n] ([math]–6) systems are theoretically investigated based on DFT. MO analysis reveals that the delocalized nature of delocalized [math] valence orbitals within the –CBO groups is a key factor for the increase in the stability of (BCO)[math](CBO)n([math]–6) with increasing [math]. The order of stability of (BCO)[math](CBO)n([math], 3 and 4) isomers is explained by recursion analysis at charge distribution. The determined order of stability is consistent with the energy order. The stability of (BCO)[math](CBO)n ([math] to 6) increase with increasing n. The charge distribution of a single –CBO substituted (BCO)5(CBO) ([math]) is considered a starting point of our recursion analysis for the order of stability of (BCO)4(CBO)2 ([math]) isomers. The charge distribution in the single –CBO-substituted (BCO)5(CBO) ([math]) shows that meta and para positions have the highest electron density, making them the most favored positions for substitution by highly electronegative –CBO groups. Citation: Journal of Theoretical and Computational Chemistry PubDate: 2016-12-27T08:53:39Z DOI: 10.1142/S021963361650067X

Authors:Mostafa Lashkarbolooki, Ali Zeinolabedini Hezave, Mahdi Bayat, Mohammad Hasan Khademi, Behzad Vaferi Abstract: Journal of Theoretical and Computational Chemistry, Volume 15, Issue 08, December 2016. Thermal conductivity is one of the most important properties of materials especially liquids. Thermal conductivity is highly interesting since it is the main parameter required for most heat-transfer calculations or design of industrial equipments such as heat exchangers. Unfortunately, thermal conductivity is extremely hard to be experimentally measured due to some operational issues which shifted the researchers’ interest toward proposing a correlation to predict this property. In the light of this limitation, in the current study, a four-parameter correlation based on critical temperature ([math], critical pressure ([math] and acentric factor ([math] of substances is proposed to predict thermal conductivity of liquids. In this way, 956 experimental data points collected from previously published literatures were divided into two different subsets, namely training and testing subsets, in the first stage and then used to find the optimum values of fitting parameters of the proposed correlation. Based on the obtained results of error analysis, it can be concluded that the proposed correlation has capability of both extrapolating and correlating the thermal conductivity of liquids. • Thermal conductivities of liquid hydrocarbons and aromatics are predicted. • A four-parameter correlation is proposed to predict thermal conductivity. • The correlation has capability of extrapolating and correlating the thermal conductivity. Citation: Journal of Theoretical and Computational Chemistry PubDate: 2016-12-27T08:53:23Z DOI: 10.1142/S0219633616500656

Authors:Bikram Nath, Chandan Kumar Mondal Abstract: Journal of Theoretical and Computational Chemistry, Volume 15, Issue 08, December 2016. Zeno and anti-Zeno effects in the evolution of the multi-photonic dissociation dynamics of the diatomic molecule HBr[math] owing to repeated measurements demand if the system in the initial state have been studied. The effects have been calculated numerically for the case of vibrational population transfer and dissociation dynamics of HBr[math] taking it as a model. We use time-dependent Fourier grid Hamiltonian (TDFGH) method as a mathematical tool in presence of intense radiation field as perturbation. The effects have been explored through a probable mechanism of population transfer from the ground vibrational state to the different upper vibrational states which ultimately go to the dissociation continuum. The results show significant differences in the mechanism of population transfer and the significant role of time interval of measurement ([math] in Zeno and anti-Zeno effects. In case of survival probability of ground vibrational states, there is Zeno effect when the frequency of the laser to which the molecule is submitted is near the vibrational [math] to [math] resonance, while there is anti-Zeno effect if it is far from this resonance. Zeno and anti-Zeno effects on the dissociation dynamics of HBr[math] owing to repeated measurements asking if the system is in the initial state have been studied. The study shows the significant differences in the mechanism of population transfer and important role of time interval of measurement in Zeno and anti-Zeno effects. In case of survival probability of ground vibrational states, there is Zeno effect when the frequency of the laser to which the molecule is submitted is near the vibrational [math] to [math] resonance while there is anti-Zeno effect if it is far from this resonance. Citation: Journal of Theoretical and Computational Chemistry PubDate: 2016-12-27T08:53:21Z DOI: 10.1142/S021963361650070X

Authors:Gabriele Milani, Federico Milani Abstract: Journal of Theoretical and Computational Chemistry, Volume 15, Issue 08, December 2016. A robust extrapolation model requiring few input parameters to predict the kinetic constants characterizing the curing behavior of NR vulcanized in presence of sulfur and two accelerants (TBSS and DPG) at different concentrations is discussed. The numerical model based on the reproduction of rheometer curves by means of the well-known Han’s kinetic model, which describes with kinetic base the most important NR vulcanization phases, namely curing initiation, formation of matured crosslinked polymer and reversion. The derived mathematical model is a closed form exponential function depending on only three kinetic constants. The procedure proposed is a two-step one. In the first step, kinetic constants of NR in the presence of single activators (i.e. either only with S and TBBS or S and DPG) are estimated by means of an interactive trial and error optimization software (GURU) that proceeds in approximating more and more strictly normalized experimental rheometer curves with Han’s function. Four different concentrations of S and TBSS (or S and DPG) are assumed as calibration points. In the second step, from the results obtained previously, kinetic constants for NR with any S-TBBS-DPG concentration of technical relevance are deduced by means of standard mathematical extrapolation. The procedure is benchmarked on 16 different S-TBBS-DPG concentrations at two temperatures (150[math]C and 180[math]C), for which both experimental data are available and kinetic constants are previously derived with GURU. Quite good agreement is found, meaning that the approach may be useful for practical purposes, because expensive and cumbersome experimental investigations can be avoided. • A robust simple mathematical extrapolation method to predict the kinetic constants characterizing NR curing is proposed. • A validation on an extensive experimental campaign at different sulphur and two accelerants (TBSS and DPG) concentrations is provided. • The extrapolated kinetic model, based on three kinetic constants and able to characterize reversion, allows performing simulations at temperatures and concentrations different from those experimentally tested. Citation: Journal of Theoretical and Computational Chemistry PubDate: 2016-12-27T08:53:19Z DOI: 10.1142/S0219633616500681

Authors:Juan Frau, Francisco Muñoz, Daniel Glossman-Mitnik Abstract: Journal of Theoretical and Computational Chemistry, Ahead of Print. The validity of the “Koopmans in DFT” (KID) procedure have been assessed by means of the calculation of several Conceptual DFT reactivity descriptors calculated through a [math]SCF procedure compared with the results of the HOMO and LUMO energies of the neutral system. Three resveratrol derivatives were considered: cis- and trans-piceid and resveratrone-6-O-[math] glucoside. The Minnesota latest family of density functionals have been considered for the calculations in connection with water as a solvent simulated with the SMD parametrization. It is shown that the range-separated hybrids MN12SX and N12SX fulfill the KID procedure with great accuracy. Citation: Journal of Theoretical and Computational Chemistry PubDate: 2016-12-23T03:59:10Z DOI: 10.1142/S0219633617500067

Authors:Jitrayut Jitonnom, Wijitra Meelua Abstract: Journal of Theoretical and Computational Chemistry, Ahead of Print. Group 4 metallocene-mediated cationic ring-opening polymerizations of a series of lactones and cyclic carbonates, with different ring sizes ([math]–8) have been theoretically studied. Using the “naked cation” approach in combination with density functional theory, the activated chain-end mechanism and the influence of transition metals, solvent and monomer ring size on the polymerizability were explored in detail. The results showed that the cationic metallocene–monomer complex, [catalyst][monomer][math], is formed, generating cationic (carbocation ion) species responsible for polymer chain growth. We found that poor polymerizability of five-membered lactone and six-membered ring carbonate depends not only on the nature of the monomer ring size but also the relative stability of the complex, which was found to correlate well with the ring strain. Subsequently, several propagation steps take place through an SN2 reaction which involves ring opening of an active monomer, via alkyl–oxygen bond cleavage. Based on the computed activation energies of all metallocene systems, the first propagation was found to be the rate-determining step of the overall propagation and the hafnocene was found to be most active with the energy barrier of 17.6[math]kcal/mol, followed by zirconocene (18.6[math]kcal/mol) and titanocene (19.5[math]kcal/mol), respectively. The mechanistic study may be applicable to the cationic ROP of lactides and other related monomers. Citation: Journal of Theoretical and Computational Chemistry PubDate: 2016-12-22T07:04:57Z DOI: 10.1142/S0219633617500031

Authors:Yuan Zhang, Huanjie Wang, Weiren Xu, Fancui Meng Abstract: Journal of Theoretical and Computational Chemistry, Ahead of Print. Martini coarse-grained force field simulations have been carried out to estimate the free energy profiles of the spontaneous membrane-translocating peptide TP2 and one negative control peptide ONEG with POPC as the model bilayer. The results show that the free energy minimum of TP2 is [math]20[math]kJ/mol lower than that of ONEG. In addition, the minimum of TP2 shifts slightly to the bilayer center compared with ONEG. The translocation barrier height for TP2 and ONEG are 119.0[math]kJ/mol and 155.7[math]kJ/mol, respectively. The lower central energy barrier of TP2 facilitates the transition between two leaflets of POPC. Both translocating peptides induce the formation of funnel-shaped structures at the bilayer center, but TP2 has a more compact structure and brings less perturbation compared with ONEG. Subsequently all atom molecular simulations testify the findings. It is indicated that compared with its negative control ONEG, TP2 binds better with lipid and penetrates deeper into bilayer with less perturbation to the bilayer structure. Our findings may shed light on the design and virtual screening of spontaneous membrane-translocating peptides. Citation: Journal of Theoretical and Computational Chemistry PubDate: 2016-12-05T09:33:23Z DOI: 10.1142/S021963361750002X

Authors:Ahmed Taki Eddine Ardjani, Sidi Mohamed Mekelleche Abstract: Journal of Theoretical and Computational Chemistry, Ahead of Print. A theoretical study of the antioxidant behavior of N[math]-(2-hydroxy-3-methoxy-benzylidene)-4-tert-buty-lenzohydrazide (1), N[math]-(5-bromo-2-hydroxy-benzylidene)-4-tert-butyl benzohydrazide (2) and N[math]-(2-hydroxy-3-methoxybenzylidene)-4-methyl-benzene-sulfonohydrazide (3) and their tautomers 1 [math] –3 [math] have been carried out at B3LYP/6-31[math]G(2d,2p). The numerical values of descriptors, namely, bond dissociation enthalpy, proton affinity (PA), electron transfer enthalpy (ETE), ionization potential, and proton dissociation enthalpy (PDE) have been calculated in gas phase and media solution (EtOH, DMSO and water). The obtained results show that the hydrogen atom transfer (HAT) mechanism is more favored thermodynamically in gas phase, whereas the sequential proton loss electron transfer (SPLET) mechanism is more preferred in solvents. Moreover, the couple (3,3[math]) is found to be the most potent antioxidant as expected experimentally. Furthermore, the BDE values of compound 3 [math] is much lower than that of ascorbic acid (AA), indicating that the tautomerization of compounds 1–3 has great influence on the antioxidant activity of these compounds. The antioxidant power of compounds (3.3[math]) was also rationalized by the calculation of the atomic spin density. In addition, the molecular docking study of compounds 1–3 and 1[math] –3[math] on xanthine oxidase (XO) as the protein target revealed important interactions between active compounds and amino acids. Moreover, compound 3 is predicted to be a potential inhibitor with higher activity. Citation: Journal of Theoretical and Computational Chemistry PubDate: 2016-11-17T02:10:23Z DOI: 10.1142/S0219633617500018

Authors:Marilena Carbone Abstract: Journal of Theoretical and Computational Chemistry, Ahead of Print. The adsorption of [math]-AminoTiophene on Si(100)2[math][math][math]1 was investigated by van der Waals corrected DFT and climbing image nudged-elastic band, in view of potential applications in silicon-based technologies. The overall scenario indicates that dissociative states are more favorable than the molecular ones, the one occurring through N–C bond breakage and Si–N and Si–C bond formation, having the largest adsorption energy (2.71[math]eV). Furthermore, this configuration is also kinetically easily accessible, being connecting to one of the physisorbed states (Phys1) by a nearly barrierless transition. Also the molecular states are relatively easily kinetically accessible, with transition barriers from the corresponding physisorbed states in the 0.05–0.30[math]eV range. At variance with this, the transitions to the dissociative state characterized by N–H bond breakage and Si–N and Si–H bond formation (N–H Diss) either from physisorbed or from molecular states are all significantly higher, i.e. in the 0.63–2.70[math]eV range. Finally, the effects of the coverage on the adsorption energy were evaluated for the N–H Diss configuration and indicating a gain, whose extent depends both on the coverage and on the surface arrangement, i.e. whether cis or trans. The trend is different if the vdW forces are excluded. Citation: Journal of Theoretical and Computational Chemistry PubDate: 2016-11-17T02:10:21Z DOI: 10.1142/S0219633617400016

Authors:Anbang Li, Kaifu Gao Abstract: Journal of Theoretical and Computational Chemistry, Ahead of Print. Poisson–Boltzmann (PB) model is a widely used implicit solvent approximation in biophysical modeling because of its ability to provide accurate and reliable PB electrostatic salvation free energies ([math] as well as electrostatic binding free energy ([math] estimations. However, a recent study has warned that the 0.5[math]Å grid spacing which is normally adopted can produce unacceptable errors in [math] estimation with the solvent excluded surface (SES) (Harris RC, Boschitsch AH and Fenley MO, Influence of grid spacing in Poisson–Boltzmann equation binding energy estimation, J Chem Theory Comput 19: 3677–3685, 2013). In this work, we investigate the grid dependence of the widely used PB solver DelPhi v6.2 with molecular surface (MS) for estimating both electrostatic solvation free energies and electrostatic binding free energies. Our results indicate that, for the molecular complex and components the absolute errors of [math] are smaller than that of [math], and grid spacing of 0.8[math]Å with DelPhi program ensures the accuracy and reliability of [math]; however, the accuracy of [math] largely relies on the order of magnitude of [math] itself rather than that of [math] or [math]. Our findings suggest that grid spacing of 0.5[math]Å is enough to produce accurate [math] for molecules whose [math] are large, but finer grids are needed when [math] is very small. Citation: Journal of Theoretical and Computational Chemistry PubDate: 2016-11-17T02:10:16Z DOI: 10.1142/S0219633616500711

Authors:Ana María Mendoza-Wilson, Francisco Javier Carmelo-Luna, Humberto Astiazarán-García, Bertha I. Pacheco-Moreno, Iván Anduro-Corona, María Lucila Rascón-Durán Abstract: Journal of Theoretical and Computational Chemistry, Ahead of Print. The objective of this study was to determine the physicochemical properties of the oligomers of procyanidins (PCs) including PA1, PA2, PB1, PB2, PC1 and a B-type tetramer, taking into account of their conformations related to the interflavan links using the density functional theory (DFT). This information may provide useful insight into the potential effect of physicochemical properties on the absorption of PCs. The results indicate that A-type and B-type PCs in all of their conformations tend to be more stable in water than in octanol, showing a hydrophilic character due to their negative log [math] values, which increase with the degree of polymerization (DP). The studied PCs, including the B-type tetramer, can achieve an appropriate molecular size (i.e. width and length) that can allow them to pass through the pores in the paracellular route in the human small intestine. The factor that could limit the absorption of the PC oligomers with increases in size is the higher number of hydroxyl groups exposed to the outside of the molecule due to their potential to interact with other molecules, which is based on electrostatic potential maps. Citation: Journal of Theoretical and Computational Chemistry PubDate: 2016-11-03T06:16:29Z DOI: 10.1142/S0219633616500693

Authors:Xin Zhang, Yuan Zhao, Xinli Duan, Hui N. Zhang, Zexing Cao, Yirong Mo Abstract: Journal of Theoretical and Computational Chemistry, Ahead of Print. The deamination reaction of 8-oxoguanine (8-oxoG) catalyzed by 8-oxoguanine deaminase (8-oxoGD) plays a critically important role in the DNA repair activity for oxidative damage. In order to elucidate the complete enzymatic catalysis mechanism at the stages of 8-oxoguanine binding, departure of 2-hydroxy-1H-purine-6,8(7H,9H)-dione from the active site, and formation of 8-oxoxanthine, extensive combined QM(PM3)/MM molecular dynamics simulations have been performed. Computations show that the rate-limiting step corresponds to the nucleophilic attack from zinc-coordinate hydroxide group to free 8-oxoguanine. Through conformational analyses, we demonstrate that Trp115, Trp123 and Leu119 connect to O8@8-oxoguanine with hydrogen bonds, and we suggest that mutations of tryptophan (115 and 123) to histidine or phenylalanine and mutation of leucine (119) to alanine could potentially lead to a mutant with enhanced activity. On this ground, a proton transfer mechanism for the formation of 8-oxoxanthine was further discussed. Both Glu218 and water molecule could be used as proton shuttles, and water molecule plays a major role in proton transfer in substrate. On the other hand, comparative simulations on the deamination of guanine and isocytosine reveal that, for the helping of hydrogen bonds between O8@8-oxoguanine and enzyme, O8@8-oxoguanine is the fastest to be deaminated among the three substrates which are also supported by the experimental kinetic constants. Citation: Journal of Theoretical and Computational Chemistry PubDate: 2016-10-25T03:37:09Z DOI: 10.1142/S0219633616500668

Authors:Dating Wu, Hui Zhang Abstract: Journal of Theoretical and Computational Chemistry, Ahead of Print. We present results of a detailed numerical investigation of the phase separation kinetic process of the macromolecular microsphere composite (MMC) hydrogel. Based on the Flory-Huggins-de Gennes-like reticular free energy, we use the time-dependent Ginzburg–Landau (TDGL) mesoscopic model (called MMC-TDGL model) to simulate the phase separation process. Domain growth is investigated through the pair correlation function. Then we obtain the time-dependent characteristic domain size, which reflects the growth kinetics of the MMC hydrogel. The results indicate that the growth law based on the MMC-TDGL equation is consistent with the modified Lifshitz–Slyozov theory. Citation: Journal of Theoretical and Computational Chemistry PubDate: 2016-10-14T09:52:05Z DOI: 10.1142/S0219633616500644