Next Article in Journal
Formation of a Conducting Polymer by Different Electrochemical Techniques and Their Effect on Obtaining an Immunosensor for Immunoglobulin G
Next Article in Special Issue
Phosphazene Functionalized Silsesquioxane-Based Porous Polymer as Thermally Stable and Reusable Catalyst for Bulk Ring-Opening Polymerization of ε-Caprolactone
Previous Article in Journal
Experimental and Numerical Investigation of the Mechanical Properties of 3D-Printed Hybrid and Non-Hybrid Composites
Previous Article in Special Issue
Polymerization of Hexene-1 and Propylene over Supported Titanium–Magnesium Catalyst: Comparative Data on the Polymerization Kinetics and Molecular Weight Characteristics of Polymers
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Polymers
Volume 15
Issue 5
10.3390/polym15051166
polymers-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

DFT and ONIOM Simulation of 1,3-Butadiene Polymerization Catalyzed by Neodymium-Based Ziegler–Natta System

by
Alexey N. Masliy
1,
Ildar G. Akhmetov
2,
Andrey M. Kuznetsov
1 and
Ilsiya M. Davletbaeva
3,*
1
Department of Inorganic Chemistry, Kazan National Research Technological University, K. Marx Street 68, 420015 Kazan, Russia
2
Nizhnekamsk Chemical and Technological Institute (Branch), Kazan National Research, Technological University, K. Marx Street 68, 420015 Kazan, Russia
3
Technology of Synthetic Rubber Department, Kazan National Research Technological University, K. Marx Street 68, 420015 Kazan, Russia
*
Author to whom correspondence should be addressed.
Polymers 2023, 15(5), 1166; https://doi.org/10.3390/polym15051166
Submission received: 25 January 2023 / Revised: 17 February 2023 / Accepted: 22 February 2023 / Published: 25 February 2023
(This article belongs to the Collection Polymerization and Kinetic Studies)
Browse Figures
Versions Notes

Abstract

:
Using modern methods of quantum chemistry, a theoretical substantiation of the high cis-stereospecificity of 1,3-butadiene polymerization catalyzed by the neodymium-based Ziegler–Natta system was carried out. For DFT and ONIOM simulation, the most cis-stereospecific active site of the catalytic system was used. By analyzing the total energy, as well as the enthalpy and Gibbs free energy of the simulated catalytically active centers, it was found that the coordination of 1,3-butadiene in the trans-form was more favorable than in the cis-form by 11 kJ/mol. However, as a result of π-allylic insertion mechanism modeling, it was found that the activation energy of cis-1,3-butadiene insertion into the π-allylic neodymium–carbon bond of the terminal group on the reactive growing chain was 10–15 kJ/mol lower than the activation energy of trans-1,3-butadiene insertion. The activation energies did not change when both trans-1,4-butadiene and cis-1,4-butadiene were used for modeling. That is, 1,4-cis-regulation was due not to the primary coordination of 1,3-butadiene in its cis-configuration, but to its lower energy of attachment to the active site. The obtained results allowed us to clarify the mechanism of the high cis-stereospecificity of 1,3-butadiene polymerization by the neodymium-based Ziegler–Natta system.

Graphical Abstract

1. Introduction

For more than 75 years, Ziegler–Natta catalysts have continued to be the main type of catalytic system used for the industrial production of polymers with a regular structure [1,2,3,4]. Here, it is necessary to note that the design of Ziegler–Natta catalysts was developed for the polymerization of conjugated dienes and the production of synthetic rubbers [5,6,7,8,9]. The large-scale industrial production of stereoregular rubbers for mass use, including cis-1,4-polyisoprene (an analogue of natural rubber) and cis-1,4-polybutadiene, is one of the outstanding achievements of global science and technology [10]. Over the past three decades, technologies for the production of synthetic rubbers using lanthanide, primarily neodymium-containing catalytic systems, have been actively developed [11,12,13,14,15,16]. Despite the large number of works, research in this area continues to be relevant. This is due to the unique capabilities of neodymium-containing catalytic systems: the high activity, stereoregularity, and controlled molecular parameters of the resulting polymers, and the ability to use various monomers in their synthesis [17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33,34]. In industrial practice, the so-called “binary” [35,36,37,38,39] and “ternary” [40,41,42,43] neodymium-containing catalytic systems have received the greatest development. Despite the differences in the nature of the initial components, we agree that the structure of the active centers, under the conditions of a comparable set of atoms, is similar. This is confirmed by the similar characteristics of polymers obtained using various neodymium-containing catalysts. In our opinion, Yu. B. Monakov’s group advanced most extensively our understanding of the structure of active centers. They showed that Ziegler–Natta catalysts based on neodymium are polycenter and can contain up to six types of active centers in their structure. The set of active centers contained in the catalytic system depends on the synthesis conditions. In this case, the maximum cis-stereospecificity of the centers corresponds to the structure with the highest chlorine content at the neodymium atom and, on the contrary, decreases up to inversion to trans-stereospecificity with a decrease in the chlorine content in the active center [44,45,46,47,48,49,50,51,52,53,54,55]. At the same time, even in these works, the answer to the following question was not obtained: why is the “outstanding” cis-stereospecificity of neodymium-containing catalytic systems observed with more energetically favorable forms of trans-1,3-butadiene, such as monomers and syn-forms of the terminal π-allyl link, as a “matrix” of the polymer chain.
The development of theoretical chemistry methods, primarily quantum chemical methods, has made it possible to clarify various aspects of the ion-coordination polymerization of dienes [53,55,56,57,58,59,60,61]. At the same time, the significant increase in the technical level of theoretical calculations, which occurred due to the increased technical capabilities of modern computers and software optimized for their use, has allowed researchers to not only explain the results of physical experiments, but also offer new solutions in catalysis [62,63]. In this work, using modern methods of quantum chemistry based on structural and thermochemical data, an attempt was made to theoretically substantiate the high stereospecificity of the polymerization of 1,3-butadiene under the action of a neodymium-containing Ziegler–Natta catalyst. The presented materials form the first part of a research project on the revision of the available experimental data [64,65] and the theoretical creation of new efficient industrial neodymium-containing Ziegler–Natta catalysts.

2. Computational Details

Quantum-chemical calculations were carried out using the Orca 5.0 program package [66,67] in the framework of the density functional theory. The hybrid density functional B3LYP [68,69] was used in the calculations in combination with the second-generation split-valence triple-ζ atomic basis set def2-TZVP presented by Ahlrichs et al. [70,71]. The core electrons of the neodymium atom were described using the def2-ECP pseudopotential [72] recommended for use with def2 basis sets. Since the Nd(III) ion has three unpaired electrons, and the total multiplicity of all the complexes under study was four (confirmed by a series of additional calculations), the spin-unrestricted method was used for all calculations.
A feature of the systems under study was that weak Van der Waals interactions played an important role in them, which were taken into account in the framework of the D3 semiempirical model [73] presented by Grimme et al.
All calculations were performed with the full optimization of molecular geometry without any symmetry constrains. Geometry optimization was carried out taking into account the influence of the solvent (hexane) in the self-consistent reaction field (SCRF) conductor-like polarizable continuum model (C-PCM) [74,75].
The calculations of vibrational spectra performed after geometry optimization did not contain imaginary modes. This meant that the structures found corresponded to minima on the total potential energy surface. These calculations were also used to estimate the thermal corrections needed to calculate the total Gibbs free energy of particles (at 298.15 K and 1 atm).
Since it was found in the first part of the work that the results obtained did not allow us to draw unambiguous conclusions (see conclusions in Section 3.1 below), it became clear that to continue the work, larger catalyst systems would need to be considered. Therefore, it was decided to carry out further studies (all results provided in Section 3.2 below) within the framework of the ONIOM two-layer (QM1/QM2) model [76], in which a part was allocated to the system under study, calculated using the high-level quantum chemical method QM1; the system as a whole was calculated using the low-level QM2 method. Then, the results of these two calculations were combined using a special technique.
As a high-level (QM1) method, we used the same approximation as in the first part of the work, B3LYP/def2-TZVP. Based on the results of preliminary calculations, the tight-binding DFT (semi-empirical DFT) method XTB1 presented by Grimme et al. [77,78] was chosen as a low-level QM2 method.
Van der Waals interactions were also taken into account within the semi-empirical D3 model designed by Grimme et al. [73]. All calculations were performed with the full optimization of molecular geometry without any symmetry constrains. Geometry optimization was carried out taking into account the influence of the solvent (hexane) in the continuum model ALPB—a method for the analytical linearization of the Poisson–Boltzmann equation [79].
To calculate the activation energies of the process, we used the standard procedure for searching for the transition state. The calculations of vibrational spectra performed after geometry optimization did not contain imaginary modes for minima and contained one imaginary mode for transition states. Each transition state identified was checked for compliance with the considered chemical process using the internal reaction coordinate (IRC) procedure.
The ONIOM technique involved the selection of a fragment in the system under study, which was calculated within the framework of the high-level QM1 method. This fragment was the same for all systems that were considered in the framework of this study. The active side (AS) fragment chosen for QM1 calculations consisted of 43 atoms (see Figure S1 in Supplementary Materials). Dangling chemical bonds within the framework of the ONIOM model were modeled using hydrogen atoms, which were located along the broken bonds. These hydrogen atoms were completed by the program itself. All other atoms, including the growing polymer chain, were calculated using the low-level QM2 method.

3. Results and Discussion

3.1. Structural and Thermochemical Analysis of AS

It is known that Nd(III) is characterized by coordination numbers from six to nine. According to numerous studies, Nd(III) is capable of coordinating diene monomers. As a result, the generally accepted mechanism of interaction becomes possible, which includes the stage of the coordination of the monomer at the active site of polymerization. It was shown in [55] that the polymerization of butadiene initiated by the neodymium catalytic system can lead to the formation of six types of active sites (ASs). In this work, the active site AS+C1H2C2HC3HC4H2 (Figure 1), which exhibited the highest cis-stereospecificity of action, was chosen as a model structure for calculations.
It is known that the trans-configuration is more favorable for 1,3-butadiene [80]. Our calculations confirmed that the trans-configuration was the most favorable for the existence of 1,3-butadiene both in the gas phase and in hexane. Thus, the difference in Gibbs free energy between the cis- and trans-configurations of 1,3-butadiene was 12.2 kJ/mol (ΔH0298 is 14.5 kJ/mol). Due to the fact that the AS+C1H2C2HC3HC4H2 catalytic system considered in this work exhibited high cis-stereospecificity, when calculating the active center, its complexes with both cis- and trans-forms of 1,3-butadiene were considered (Figure 1).
As a result of minimizing the total energy, the optimal structures of the AS+C1H2C2HC3HC4H2 catalytic complex were found for both configurations of 1,3-butadiene (hereinafter AS+η-cis-C4H6 and AS+η-trans-C4H6) (Figure 2). The most important geometric characteristics of these structures are presented in Table 1 and Table 2.
According to Figure 2, the observed free space near Nd(III) made it possible to judge the absence of any steric hindrances for the coordination interaction of Nd(III) with 1,3-butadiene in both cis- and trans-configurations. In the presented structures, the Nd(III) ion was connected to the isobutyl carbanion and aluminum ions by chloride bridges. In addition, according to Figure 2, the Nd(III) ion in a given ligand environment created coordination structures of distorted octahedral symmetry.
According to Table 1, neodymium, chlorine, and aluminum ions in the common molecular structure lay almost in the same plane. At the same time, the planes of two Al-Cl-Nd-Cl quadrilaterals were almost perpendicular. Scheme 1 shows the possible modes of the coordination interaction of 1,3-butadiene with the active site of polymerization.
It follows from Table 1 that the distance between Nd(III) and the carbon atoms in the 1,3-butadiene in both the cis- and trans-forms was ~3 Å. This indicated that the main role in the coordination bonding of Nd(III) with 1,3-butadiene was played by π-bonds, as well as electrostatic (dispersion and Van der Waals) interactions. The η4-cis/trans-coordination of 1,3-butadiene by Nd(III) ions was more favorable than the η2-cis/trans-coordination.
The orientation of cis- and trans-1,3-butadiene relative to i-C4H9 in the active center had to be considered separately, since the further incorporation of the monomer occurred as a result of the interaction of these two particles. Table 1 shows the shortest distances between the terminal atoms of these particles, that is, those atoms that should interact with each other in the future. As can be seen from Table 1, in the case of cis-1,3-butadiene, this distance was ~0.4 Å shorter, which indicated a greater possibility for the interaction of this form of 1,3-butadiene. In addition, Figure 2 shows that the orientation of cis-1,3-butadiene relative to i-C4H9 contributed to a more efficient reaction.
An analysis of the total energy, as well as the enthalpy and Gibbs free energy, of the AS+η-cis-C4H6 and AS+η-trans-C4H6 complexes showed that the coordination of trans-1,3-butadiene was more favorable than the coordination binding of cis-1,3-butadiene by 11.0 kJ/mol. The changes in enthalpy and Gibbs free energy under standard conditions for trans–cis transformation were 13.7 and 7.8 kJ/mol, respectively. The fact that these values were somewhat lower (in particular, the change in the Gibbs free energy was lower by 5 kJ/mol) in comparison with 1,3-butadiene, which was outside the zone of the polymerization process, was not a sufficient justification for the formation of cis-1,4-polybutadiene.
The result of the next stage, the polymerization initiation reaction involving AS and 1,3-butadiene, was the formation of the i-C4H9-H2CHCγHCβH2Cα-AS structure. The possible variations in the terminal group structures of the reactive growing chain are presented in Scheme 2:
Since the π-allylic neodymium–carbon bond was an integral part of the terminal group on the reactive growing chain, the i-C4H9-π-anti-H2CHCγHCβH2Cα-AS and i-C4H9-π-syn-H2CHCγHCβH2Cα-AS complexes were used for further calculations.
As can be seen from Figure 3, there were noticeable differences for the terminal allylic end unit, which existed as anti- and syn-π-allylic complexes. First of all, in the case of the anti-π-allylic complexes, the distance from three atoms of trans-1,3-butadiene to Nd(III) practically coincided. This indicated the predominance of π- and dispersion interactions when this form was bound to Nd(III). In the case of the syn-π-allylic complexes, monodentate metal–ligand bonds and π-interactions with the double bond of the terminal unit of the polymer chain were formed.
As a result, the ligand environment of the Nd(III) ion in the terminal group on the reactive growing chain underwent changes when the π-allylic structure was repositioned from the anti-π-allylic complex to the syn-π-allylic complex. The ligand environment of Nd(III) in the active syn-π-allylic complex did not change much compared to the AS+trans-C4H9 complex described above. In the case of the anti-π-allylic complex, the chlorine atoms were located in the same plane and formed a rectangle 3.5 × 4 Å in size. The Nd(III) ion was located above the center of this rectangle, and above it was a semiring of cis-1,3-butadiene.
The Gibbs free energy of the formation of i-C4H9-π-anti-H2CHCγHCβH2Cα-AS and i-C4H9-π-syn-H2CHCγHCβH2Cα-AS complexes from the more favorable and, accordingly, more probable complex AS+trans-C4H6 was −84.5 and −84.9 kJ/mol, respectively. That is, we could conclude that the formation of both forms of the complex was almost equally likely, but this fact also could not explain the predominant formation of cis-1,4-polybutadiene.
Next, we studied the complexes simulating the coordination of 1,3-butadiene in the terminal group on the reactive growing chain. Based on the assumption of the cis-stereospecificity of the AS from the whole variety of structures, the complexes that lead to the formation of cis-1,4-polybutadiene were considered: i-C4H9-π-anti-H2CHCγHCβH2Cα-AS+η-cis-C4H6 (A) and i-C4H9-π-anti-H2CHCγHCβH2Cα-AS+η-trans-C4H6 (B) (Figure 3). The coordination of 1,3-butadiene from the volume of the reaction mixture was accompanied by the displacement of the growing polymer chain relative to the active center and the freeing of space for the monomer. According to classical concepts (Scheme 2), at this moment, a σ-cis-allylic (σ-C1-cis) structure in the terminal group should form, in which a covalent neodymium–carbon bond of the reactive growing chain is formed. However, our calculations showed that at this stage, too, the π-allylic structure of the anti-configuration was energetically more favorable than the σ-C1-cis structure. At the same time, the coordination of cis-1,3-butadiene by the Nd(III) ion mainly occurred due to π-interaction with the terminal CH-CH2 group of 1,3-butadiene (the smallest distance to Nd(III) was presented by the C1 and C2 atoms of 1,3-butadiene), while the coordination of trans-1,3-butadiene occurred due to the π-interaction with the central CH-CH group (the smallest distance to Nd(III) was presented by the C2 and C3 atoms of 1,3-butadiene). It should also be noted that the trans-form of 1,3-butadiene remained practically planar, while in the cis-form some distortion was observed, and the terminal CH2 group deviated from the plane formed by the other three carbon atoms by about 15°. As for the coordination of the polymer fragment relative to AS, it was almost identical in these two cases (Table 3).
It should be noted that the coordination of 1,3-butadiene in the cis- and trans-forms significantly changed the geometry of the active site. If in i-C4H9-π-anti-H2CHCγHCβH2Cα-AS all four chlorine atoms lay almost in the same plane (Figure 3a), then when 1,3-butadiene was coordinated, the Cl1-Nd-Cl2 plane was almost perpendicular to the Cl3-Nd-Cl4 plane (the Cl1 and Cl2 atoms were coordinated with one aluminum atom, and Cl3 and Cl4 with the second). Due to the fact that both the terminal π-allylic structure and 1,3-butadiene in the cis- and trans-forms were linked to Nd(III) via π-bonds, it was not possible to reliably determine the coordination number of Nd(III).
The thermochemical parameters for the formation of i-C4H9-π-anti-H2CHCγHCβH2Cα-AS+η-cis-C4H6 and i-C4H9-π-anti-H2CHCγHCβH2Cα-AS+η-trans-C4H6 complexes (Figure 4) were evaluated by the reaction:
i-C4H9-π-anti-H2CHCγHCβH2Cα-AS+trans-C4H6 ⇄
i-C4H9-π-anti-H2CHCγHCβH2Cα-AS+η-trans/cis-C4H6
The change in the Gibbs free energy under standard conditions for the formation of the i-C4H9-π-anti-H2CHCγHCβH2Cα-AS+η-cis-C4H6 complex was 11.1 kJ/mol, and for the i-C4H9-π-anti-H2CHCγHCβH2Cα-AS+η-trans-C4H6 formation it was 18.8 kJ/mol. This meant that under standard conditions, this reaction did not proceed. However, both reactions were exothermic: the ΔH0298 value in the first case was -46.5 kJ/mol, and in the second -43.2 kJ/mol. Accordingly, the ΔS0298 value for i-C4H9-π-anti-H2CHCγHCβH2Cα-AS+η-cis-C4H6 was -193.3 J/K·mol, and for i-C4H9-π-anti-H2CHCγHCβH2Cα-AS+η-trans-C4H6 it was -208.1 J/K·mol. Due to the fact that the ΔH and ΔS had the same signs, these reactions were reversible, and an increase in the concentration of reagents shifted the equilibrium towards the formation of products. Taking into account that the ΔG0298 in both cases had a rather small positive value, under the conditions of a real reaction mixture, where the concentration of 1,3-butadiene is significantly higher than the concentration of active sites, both reactions should proceed in the forward direction. It should also be noted that the Gibbs free energy of the trans–cis transition of 1.3-butadiene upon interaction with the active center remained equal to 7.7 kJ/mol, as in Figure 2. The values of ΔG0298 for these reactions were quite close. This meant that approximately equal concentrations of both the anti- and syn-forms of terminal π-allyl structures should be expected in the reaction system.
Thus, the analysis of the structural characteristics and thermochemical features of the formation of AS complexes with 1,3-butadiene, corresponding to stages such as monomer coordination and the initiation and growth of the polymer chain, was insufficient to justify the predominant formation of cis-1,4-polybutadiene by the neodymium-based Ziegler–Natta system. This required a more detailed study of the process mechanism.

3.2. Study of Polymerization Mechanism

At the initial stage, it was decided to calculate the transition state for the trans–cis transition of 1,3-butadiene in the reaction system, taking into account its interaction with the AS. The resulting transition state structures are presented in the Supplementary Materials (Tables S7–S57). Further, within the framework of this section, we consider only the change of the Gibbs free energy; however, the Supplementary Materials also provide tables similar to Table 4 but containing the values of the change in the total energy E0 (Table S58) and enthalpy (Table S59). The activation energy (ΔG≠) determined for the trans–cis transition for 1,3-butadiene in the free state, associated with the rotation of the CH2 group around the C=C bond, was 27.0 kJ/mol. A similar transition was calculated for 1,3-butadiene interacting with the active site, and the activation energy in this case increased to 43 kJ/mol. That is, the assumption that the active center facilitated the trans–cis transition for 1,3-butadiene was not confirmed. Thus, the initial coordination of 1,3-butadiene at the active center in the trans-form was legitimate for subsequent calculations.
Next, the transition states associated with the insertion of 1,3-butadiene into the polymer chain were simulated. At the stage of initiation, this was ensured by the rotation of 1,3-butadiene in such a way that the terminal CH2 group was oriented toward the i-C4H9 group bound to Nd(III), as well as with some distance between i-C4H9 and Nd(III) and approaching the terminal CH2 group of 1,3-butadiene. At the stage of polymer chain growth, everything should have proceeded in exactly the same way, only the position of i-C4H9 was occupied by a reactive growing chain. Transition states were successfully identified for trans- and cis-1,3-butadiene. The correspondence of the transition states of the studied reaction was confirmed by the calculation of “descents” to the reactants and products in the IRC procedure. The calculated activation energy (ΔG≠) of the polymerization process for trans-1,3-butadiene was 67.0 kJ/mol, and for cis-1,3-butadiene it was 54.3 kJ/mol. It should be noted that the activation energy of depolymerization was twice the activation energy of the direct reaction and was 122.2 kJ/mol for trans-1,3-butadiene and 112.8 kJ/mol for cis-1,3-butadiene. This indicated an extremely low probability of this process occurring.
If we assume that the initial form in this reaction was the active center with coordinated trans-1,3-butadiene (i-C4H9-AS-trans-C4H6), then the insertion of the first trans-unit into the reactive growing chain proceeded in one stage, the activation energy of which was the activation energy of the insertion of 1,3-butadiene in the trans-configuration into the terminal group on the reactive growing chain, that is, ΔG≠trans = 67.0 kJ/mol. Since the insertion of 1,3-butadiene in the cis-configuration occurred in two stages, the activation energy had to be measured from the structure with the lowest energy (Figure 5), and in this case it was ΔG≠cis = 54.3 + 6.6 = 60.9 kJ/mol.
Thus, it was shown that at the stage of initiation, the formation of the terminal π-allylic structure of the anti-configuration proceeded with a slightly lower activation energy compared to the formation of the syn-π-allylic complexes. At the same time, the energy effect of the insertion of trans-1,3-butadiene into the active center from the volume of the reaction mixture for both forms was almost the same, amounting to ΔGcis = 51.7 kJ/mol and ΔGtrans = 55.1 kJ/mol.
At the next stage of the reactive polymer chain growth, there were two almost thermodynamically equiprobable structures, i-C4H9-π-anti-H2CHCγHCβH2Cα-AS+η-trans-C4H6 and i-C4H9-π-syn-H2CHCγHCβH2Cα-AS+η-trans-C4H6, with somewhat lower energy activation during the formation of the first structure. Further, the formation of four different structures was possible, which for brevity we denote as i-C4H9-CC-AS, i-C4H9-CT-AS, i-C4H9-TC-AS, and i-C4H9-TT-AS (where “T” indicates trans-C4H6 or π-syn-H2CHCγHCβH2Cα fragments and “C” indicates cis-C4H6 or π-anti-H2CHCγHCβH2Cα fragments). To calculate the activation energy of i-C4H9-TC-AS and i-C4H9-TT-AS formation, we considered i-C4H9-C-AS-T and i-C4H9-T-AS-T, respectively. However, to form the other two structures, it was first necessary to calculate the activation energy of the trans–cis transition of 1,3-butadiene. The activation energy (ΔG≠) of the transition
i-C4H9-C-AS-T ⇆ i-C4H9-C-AS-C
practically did not change in comparison with the initiating stage and was equal to 41.0 kJ/mol. The activation energy (ΔG≠) of the transition
i-C4H9-T-AS-T ⇆ i-C4H9-T-AS-C
was also almost identical and equal to 41.9 kJ/mol. That is, this activation energy was constant and was the same in all cases.
Next, the transition states of the insertion of the second 1,3-butadiene molecule into the reactive polymer chain were calculated in all four reactions. Generalized detailed information regarding the change in activation energy and Gibbs free energy during this process is presented in Table 4.
From the analysis of the data in Table 4, it can be seen that the activation energy increased in all cases when the second unit was inserted, and when a cis-unit was inserted into the reactive growing polymer chain, it remained lower than when a trans-unit was inserted. The lowest activation energy was observed in the case of the formation of pure cis-1,4-polybutadiene. In this case, the difference in activation energies compared to the insertion of a trans-unit already exceeded 10 kJ/mol. Also interesting is the fact that the insertion of a cis-unit into a trans-unit still presented a lower activation energy, although the difference was not large. The least likely was the insertion of a trans-unit into a cis-unit.
Next, we consider the reaction of the insertion of the third and subsequent links into the reactive growing polymer chain. Here, we assumed that trans-1,3-butadiene was also inserted into the structures obtained at the previous stage from the bulk of the reaction mixture. However, a serious problem then arose, related to which combinations of cis- and trans-segments of the polymer should be taken into account in the calculations. Even if at the second stage it was still not difficult to determine this, since there were only four options, further on the number of possible combinations increased exponentially. Therefore, it was decided to consider only two options—the insertion of a cis- and trans-1,4-butadiene unit into the most probable polymer structure obtained at the previous stage, and the formation of trans-1,4-polybutadiene for comparison, as practice shows that this is possible with a slow process. It was also decided to carry out calculations before the insertion of five 1,3-butadiene molecules into the reactive growing polymer chain. Since, in this case, for cis-1,4-polybutadiene the number of monomer units exceeded the size of the Kuhn segment, which was 2–3 units, it became possible to extrapolate the results obtained to real objects, the polymer chain length of which reached several thousand monomer units.
At the third stage, taking into account the above, we considered (trans-C4H6)2-AS-trans-C4H6 and (cis-C4H6)2-AS-trans-C4H6 as initial structures, and as reaction products (trans-C4H6)3-AS, (cis-C4H6)3-AS and trans-C4H6(cis-C4H6)2-AS were obtained. The formation of (cis-C4H6)3-AS again required a preliminary trans–cis transition step in butadiene according to the reaction:
(cis-C4H6)2-AS-trans-C4H6 ⇆ (cis-C4H6)2-AS-cis-C4H6
The calculated activation energy of this reaction was significantly lower than the reactions of polymerization initiation and the beginning stage of the reactive polymer chain growth, at 24.6 kJ/mol. This value turned out to be even lower than the activation energy of the trans–cis transition of 1,3-butadiene outside the zone of the polymerization process. The change in the Gibbs free energy in the polymerization reaction zone according to Equation (4) was 14.6 kJ/mol, which was somewhat higher than the value obtained outside the polymerization process zone. Changes in the activation energy and Gibbs free energy were calculated in a similar way for the next stages of the 1,3-butadiene polymerization process.
Thus, the analysis of the data presented in Table 4 allowed us to conclude that the activation energy of the trans–cis transition of 1,3-butadiene interacting with the active center practically stabilized at a value approximately 3 kJ/mol lower compared to the activation energy for 1,3-butadiene, which was outside the zone of the polymerization process. It is expected that this value will be maintained in the future. The activation energy of the insertion of cis-1,3-butadiene into the reactive growing chain was 10–15 kJ/mol lower than the activation energy of the insertion of trans-1,3- butadiene. The activation energy values did not change when both trans- and cis-1,4-polybutadiene terminal units were used for calculations. The almost coinciding activation energy values for the insertion of cis-1,3-butadiene into the reactive growing chain at the last two stages indicated that the value of 61 kJ/mol could be taken as the activation energy of the formation of cis-1,4-polybutadiene at the AS. The obtained value agreed with the experimental data on the activation energy of the polymerization of dienes on lanthanide catalysts equal 30-60 kJ/mol [10,31]. If there are some macroscopic factors that cannot be taken into account in the framework of a quantum chemical study, then the activation energy of the insertion of cis-1,3-butadiene into the reactive growing chain will decrease by another ~10 kJ/mol and will be 51 kJ/mol, while the activation energy of the insertion of trans-1,3-butadiene will not change. Thus, in this case, the difference in activation energies will already be about 20–25 kJ/mol.

4. Conclusions

According to the structural and thermochemical study of the polymerization of 1,3-butadiene initiated by the most cis-stereospecific active site of the NdCl3·3TBF–Al(i-C4H9)3 catalytic system, the coordination of 1,3-butadiene by Nd(III) was equally probable in its cis- and trans-configurations. An analysis of the total energy, as well as the enthalpy and Gibbs free energy, of the AS+η-cis-C4H6 and AS+η-trans-C4H6 complexes showed that the coordination of 1,3-butadiene in the trans-form was more favorable than in the cis-form by 11 kJ/mol.
The transition states associated with the insertion of 1,3-butadiene into the reactive growing polymer chain were simulated. It was found that the energy of the insertion of cis-1,3-butadiene into the reactive growing polymer chain was 10-15 kJ/mol lower than the activation energy of the insertion of trans-1,3-butadiene. The activation energy values did not change when both trans- and cis-1,4-polybutadiene terminal units were used for calculations. The almost coinciding values of the activation energy for the insertion of cis-1,3-butadiene into the reactive growing polymer chain at the last two stages indicated that the value of 61 kJ/mol could be taken as the activation energy of the formation of cis-1,4-polybutadiene at the AS.
Based on the results of the study, the following mechanism of 1,3-butadiene polymerization by the neodymium-based Ziegler–Natta system could be proposed: the coordination of trans-1,3-butadiene at the active site; the trans-cis isomerization of 1,3-butadiene; and the insertion of cis-1,3-butadiene into the reactive growing polymer chain.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/polym15051166/s1, Figure S1: Part of AS (43 atoms), which was calculate on a higher level of accuracy QM1 (a), broken bonds are supplemented by hydrogen atoms (b). This part is the same for all systems calculated in the framework of ONIOM, Table S1: Atomic Cartesian coordinates (in Å) for AS+η-cis-C4H6 optimized at the B3LYP/def2-TZVP level (Figure 1A), Table S2: Atomic Cartesian coordinates (in Å) for AS+η-trans-C4H6 optimized at the B3LYP/def2-TZVP level (Figure 1B), Table S3: Atomic Cartesian coordinates (in Å) for i-C4H9-π-anti-H2CHCγHCβH2Cα-AS optimized at the B3LYP/def2-TZVP level (Figure 2A), Table S4: Atomic Cartesian coordinates (in Å) for i-C4H9-π-sin-H2CHCγHCβH2Cα-AS optimized at the B3LYP/def2-TZVP level (Figure 2B), Table S5: Atomic Cartesian coordinates (in Å) for i-C4H9-π-anti-H2CHCγHCβH2Cα-AS+cis-C4H6 optimized at the B3LYP/def2-TZVP level (Figure 3A), Table S6: Atomic Cartesian coordinates (in Å) for i-C4H9-π-anti-H2CHCγHCβH2Cα-AS+trans-C4H6 optimized at the B3LYP/def2-TZVP level (Figure 3B). The following tables show atomic Cartesian coordinates (in Å) for system optimized within framework of ONIOM. Frist 43 atoms calculated at the B3LYP/def2-TZVP level (Q1). For the Q2 level, we used XTB1. All calculations were carried out taking into account the influence of the solvent in the ALPB continuum model and dispersion corrections in the Grimme D3 semi-empirical model, Table S7: i-C4H9-AS-T, Table S8: i-C4H9-AS-C, Table S9: i-C4H9-AS-trans-cis transition state, Table S10: i-C4H9-AS-T-TS inclusion trans-C4H8 in polymer. Transition state, Table S11: i-C4H9-AS-C-TS inclusion cis-C4H8 in polymer. Transition state, Table S12: i-C4H9-T-AS addition trans-C4H8 chain in polymer, Table S13: i-C4H9-T-AS addition cis-C4H8 chain in polymer, Table S14: i-C4H9-C-AS-T, Table S15: i-C4H9-C-AS-C. Table S16: i-C4H9-T-AS-C. Table S16: i-C4H9-T-AS-T, Table S17: i-C4H9-C-AS-trans-cis transition state, Table S18: i-C4H9-T-AS-trans-cis transition state, Table S19: i-C4H9-C-AS-C-TS inclusion cis-C4H8 in polymer. Transition state, Table S20: i-C4H9-C-AS-T-TS inclusion trans-C4H8 in polymer. Transition state, Table S21: i-C4H9-T-AS-C-TS inclusion cis-C4H8 in polymer. Transition state, Table S23: i-C4H9-T-AS-T-TS inclusion trans-C4H8 in polymer. Transition state, Table S24: i-C4H9-CC-AS addition cis-C4H8 chain in polymer, Table S25: i-C4H9-TC-AS addition cis-C4H8 chain in polymer, Table S26: i-C4H9-CT-AS addition trans-C4H8 chain in polymer, Table S27: i-C4H9-TT-AS addition trans-C4H8 chain in polymer, Table S28: i-C4H9-CC-AS-T, Table S29: i-C4H9-CC-AS-C. Table S30: i-C4H9-TT-AS-T, Table S31: i-C4H9-CC-AS-trans-cis transition state, Table S32: i-C4H9-CC-AS-C-TS inclusion cis-C4H8 in polymer. Transition state, Table S33: i-C4H9-CC-AS-T-TS inclusion trans-C4H8 in polymer. Transition state, Table S34: i-C4H9-TT-AS-T-TS inclusion trans-C4H8 in polymer. Transition state, Table S35: i-C4H9-CCC-AS addition cis-C4H8 chain in polymer, Table S36: i-C4H9-CCT-AS addition trans-C4H8 chain in polymer, Table S37: i-C4H9-TTT-AS addition trans-C4H8 chain in polymer, Table S38: i-C4H9-CCC-AS-T, Table S39: i-C4H9-CCC-AS-C, Table S40: i-C4H9-TTT-AS-T, Table S41: i-C4H9-CCC-AS-trans-cis transition state, Table S42: i-C4H9-CCC-AS-C-TS inclusion cis-C4H8 in polymer. Transition state, Table S43: i-C4H9-CCC-AS-T-TS inclusion trans-C4H8 in polymer. Transition state, Table S44: i-C4H9-TTT-AS-T-TS inclusion trans-C4H8 in polymer. Transition state, Table S45: i-C4H9-CCCC-AS addition cis-C4H8 chain in polymer, Table S46: i-C4H9-CCCT-AS addition trans-C4H8 chain in polymer, Table S47: i-C4H9-TTTT-AS addition trans-C4H8 chain in polymer, Table S48: i-C4H9-CCCC-AS-T, Table S49: i-C4H9-CCCC-AS-C, Table S50: i-C4H9-TTTT-AS-T, Table S51: i-C4H9-CCCC-AS-trans-cis transition state, Table S52: i-C4H9-CCCC-AS-C-TS inclusion cis-C4H8 in polymer. Transition state, Table S53: i-C4H9-CCCC-AS-T-TS inclusion trans-C4H8 in polymer. Transition state, Table S54: i-C4H9-TTTT-AS-T-TS inclusion trans-C4H8 in polymer. Transition state, Table S55: i-C4H9-CCCCC-AS addition cis-C4H8 chain in polymer, Table S56: i-C4H9-CCCCT-AS addition trans-C4H8 chain in polymer, Table S57: i-C4H9-TTTTT-AS addition trans-C4H8 chain in polymer, Table S58: Changes in the activation energy and full energy under standard conditions (kJ/mol) as a result of the introduction of 1,3-butadiene into the growing polymer chain, Table S59: Changes in the activation energy and enthalpy under standard conditions (kJ/mol) as a result of the introduction of 1,3-butadiene into the growing polymer chain.

Author Contributions

Conceptualization, I.G.A. and I.M.D.; methodology, A.N.M. and A.M.K.; software, A.N.M.; validation, A.N.M., I.G.A., A.M.K. and I.M.D.; formal analysis, A.N.M. and I.G.A.; investigation, A.N.M. and I.G.A.; resources, A.M.K.; data curation, A.N.M. and I.G.A.; writing—original draft preparation, A.N.M. and I.G.A.; writing—review and editing, A.M.K. and I.M.D.; visualization, A.N.M.; supervision, I.M.D.; project administration, I.G.A. and I.M.D. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by [Russian Science Foundation] grant number [19-19-00136].

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ziegler, K.; Holzkamp, E.; Breil, H.; Martin, H. The Mulheim normal pressure polyethylene process. Angew. Chem.-Int. Ed. 1955, 67, 541–547. [Google Scholar] [CrossRef]
  2. Von Giulio, N.; Corradini, P. “Kristallstruktur des isotaktischen polystyrols”. Die Makromol. Chem. Macromol. Chem. Phys. 1955, 16, 77–80. [Google Scholar] [CrossRef]
  3. Ziegler, K.; Holzkamp, E.; Breil, H.; Martin, H. Das Mülheimer Normaldruck-Polyäthylen-Verfahren. Angew. Chem. 1956, 68, 721. [Google Scholar] [CrossRef]
  4. Ziegler, K.; Holzkamp, E.; Breil, H.; Martin, H. Polymerisation von Äthylen und anderen Olefinen. Angew. Chem. 1955, 67, 541–547. [Google Scholar] [CrossRef]
  5. Natta, C. Stereospecific polymerizations. J. Polym. Sci. 1960, 48, 219–239. [Google Scholar] [CrossRef]
  6. Natta, G.; Porri, L.; Carbonaro, A.; Stoppa, G. “Polymerization of conjugated diolefins by homogeneous aluminum alkyl-titanium alkoxide catalyst systems. I. Cis-1,4 isotactic poly (1,3-pentadiene)”. Die Makromol. Chem. Macromol. Chem. Phys. 1964, 77, 114–125. [Google Scholar] [CrossRef]
  7. Natta, G.; Porri, L.; Carbonaro, A. “Polymerization of conjugated diolefins by homogeneous aluminum alkyl-titanium alkoxide catalyst systems. II. 1,2-polybutadiene and 3, 4-polyisoprene”. Die Makromol. Chem. Macromol. Chem. Phys. 1964, 77, 126–138. [Google Scholar] [CrossRef]
  8. Tinyakova, E.N.; Dolgoplosk, B.A.; Zhuravleva, T.G.; Kovalevskaya, R.N.; Kuren’Gina, T.N. Synthesis of cis and trans polymers of dienes on oxide catalysts and investigation of their structure and properties. J. Polym. Sci. 1961, 52, 159–167. [Google Scholar] [CrossRef]
  9. Cooper, W.; Vaughan, G. Recent developments in the polymerization of conjugated dienes. Prog. Polym. Sci. 1967, 1, 91–160. [Google Scholar] [CrossRef]
  10. Tornqvist, E.G.M. The stereo rubbers, William M. Saltman, Wiley-Interscience, New York, 1977, 897 pp., $l49.50. J. Polym. Sci. Polym. Lett. Ed. 1978, 16, 53–55. [Google Scholar] [CrossRef]
  11. Hsieh, H.L.; Yeh, H.C. Polymerization of butadiene and isoprene with lanthanide catalysts; characterization and properties of homopolymers and copolymers. Rubber Chem. Technol. 1985, 58, 117–145. [Google Scholar] [CrossRef]
  12. Marina, N.G.; Monakov, Y.B.; Sabirov, Z.M.; Tolstikov, G.A. Lanthanide compounds—Catalysts of stereospecific polymerization of diene monomers. Review. Polym. Sci. USSR 1991, 33, 387–417. [Google Scholar] [CrossRef]
  13. Friebe, L.; Nuyken, O.; Obrecht, W. Neodymium-based Ziegler/Natta catalysts and their application in diene polymerization. In Neodymium Based Ziegler Catalysts–Fundamental Chemistry; Springer: Berlin/Heidelberg, Germany, 2006; pp. 1–154. [Google Scholar] [CrossRef]
  14. Fischbach, A.; Anwander, R. Neodymium Based Ziegler Catalysts Fundamental Chemistry; Nuyken, O., Ed.; Springer: Berlin/Heidelberg, Germany, 2006; Volume 204, pp. 155–281. [Google Scholar] [CrossRef]
  15. Zhang, Z.; Cui, D.; Wang, B.; Liu, B.; Yang, Y. Polymerization of 1,3-conjugated dienes with rare-earth metal precursors. In Molecular Catalysis of Rare-Earth Elements; Springer: Berlin/Heidelberg, Germany, 2010; pp. 49–108. [Google Scholar] [CrossRef]
  16. Wang, F.; Liu, H.; Hu, Y.; Zhang, X. Lanthanide complexes mediated coordinative chain transfer polymerization of conjugated dienes. Sci. China Technol. Sci. 2018, 61, 1286–1294. [Google Scholar] [CrossRef]
  17. Fan, C.; Bai, C.; Cai, H.; Dai, Q.; Zhang, X.; Wang, F. Preparation of high cis-1,4 polyisoprene with narrow molecular weight distribution via coordinative chain transfer polymerization. J. Polym. Sci. Part A Polym. Chem. 2010, 48, 4768–4774. [Google Scholar] [CrossRef]
  18. Manuiko, G.V.; Salakhov, I.I.; Aminova, G.A.; Akhmetov, I.G.; Dyakonov, G.S.; Bronskaya, V.V.; Demidova, E.V. Mathematical modeling of 1,3-butadiene polymerization over a neodymium-based catalyst in a batch reactor with account taken of the multisite nature of the catalyst and chain transfer to the polymer. Theor. Found. Chem. Eng. 2010, 44, 139–149. [Google Scholar] [CrossRef]
  19. Wang, F.; Zhang, C.Y.; Hu, Y.M.; Jia, X.Y.; Bai, C.X.; Zhang, X.Q. Reversible coordinative chain transfer polymerization of isoprene and copolymerization with ε-caprolactone by neodymium-based catalyst. Polymer 2012, 53, 6027–6032. [Google Scholar] [CrossRef]
  20. Wang, F.; Liu, H.; Zheng, W.; Guo, J.; Zhang, C.; Zhao, L.; Zhang, X. Fully-reversible and semi-reversible coordinative chain transfer polymerizations of 1,3-butadiene with neodymium-based catalytic systems. Polymer 2013, 54, 6716–6724. [Google Scholar] [CrossRef]
  21. Georges, S.; Touré, A.O.; Visseaux, M.; Zinck, P. Coordinative chain transfer copolymerization and terpolymerization of conjugated dienes. Macromolecules 2014, 47, 4538–4547. [Google Scholar] [CrossRef]
  22. Zheng, W.; Yan, N.; Zhu, Y.; Zhao, W.; Zhang, C.; Zhang, H.; Zhang, X. Highly trans-1,4-stereoselective coordination chain transfer polymerization of 1,3-butadiene and copolymerization with cyclic esters by a neodymium-based catalyst system. Polym. Chem. 2015, 6, 6088–6095. [Google Scholar] [CrossRef]
  23. Tanaka, R.; Yuuya, K.; Sato, H.; Eberhardt, P.; Nakayama, Y.; Shiono, T. Synthesis of stereodiblock polyisoprene consisting of cis-1,4 and trans-1,4 sequences by using a neodymium catalyst: Change of the stereospecificity triggered by an aluminum compound. Polym. Chem. 2016, 7, 1239–1243. [Google Scholar] [CrossRef] [Green Version]
  24. Díaz de León Gómez, R.E.; Enríquez-Medrano, F.J.; Maldonado Textle, H.; Mendoza Carrizales, R.; Reyes Acosta, K.; Lopez Gonzalez, H.R.; Romero, J.L.O.; Lugo Uribe, L.E. Synthesis and characterization of high cis-polymyrcene using neodymium-based catalysts. Can. J. Chem. Eng. 2016, 94, 823–832. [Google Scholar] [CrossRef]
  25. Dai, Q.; Zhang, X.; Hu, Y.; He, J.; Shi, C.; Li, Y.; Bai, C. Regulation of the cis-1,4-and trans-1,4-polybutadiene multiblock copolymers via chain shuttling polymerization using a ternary neodymium organic sulfonate catalyst. Macromolecules 2017, 50, 7887–7894. [Google Scholar] [CrossRef]
  26. Luo, Y.; Gao, Z.; Chen, J. Rare-earth metal bis (silylamide) complexes supported by ferrocene-substituted amidinate and their performance in cis-1,4 selective polymerization of isoprene. J. Organomet. Chem. 2017, 846, 18–23. [Google Scholar] [CrossRef]
  27. Hollfelder, C.O.; Jende, L.N.; Diether, D.; Zelger, T.; Stauder, R.; Maichle-Mössmer, C.; Anwander, R. 1,3-Diene Polymerization Mediated by Homoleptic Tetramethylaluminates of the Rare-Earth Metals. Catalysts 2018, 8, 61. [Google Scholar] [CrossRef] [Green Version]
  28. Göttker-Schnetmann, I.; Kenyon, P.; Mecking, S. Coordinative Chain Transfer Polymerization of Butadiene with Functionalized Aluminum Reagents. Angew. Chem. 2019, 131, 17941–17945. [Google Scholar] [CrossRef]
  29. González-Zapata, J.L.; Enríquez-Medrano, F.J.; González, H.R.L.; Revilla-Vázquez, J.; Carrizales, R.M.; Georgouvelas, D.; Valencia, L.; de León Gómez, R.E.D. Introducing random bio-terpene segments to high cis-polybutadiene: Making elastomeric materials more sustainable. RSC Adv. 2020, 10, 44096–44102. [Google Scholar] [CrossRef]
  30. Zheng, W.; Yang, Q.; Dong, J.; Wang, F.; Luo, F.; Liu, H.; Zhang, X. Neodymium-based one-precatalyst/dual-cocatalyst system for chain shuttling polymerization to access cis-1,4/trans-1,4 multiblock polybutadienes. Mater. Today Commun. 2021, 27, 102453. [Google Scholar] [CrossRef]
  31. Wang, H.; Cue, J.M.O.; Calubaquib, E.L.; Kularatne, R.N.; Taslimy, S.; Miller, J.T.; Stefan, M.C. Neodymium catalysts for polymerization of dienes, vinyl monomers, and ε-caprolactone. Polym. Chem. 2021, 12, 6790–6823. [Google Scholar] [CrossRef]
  32. Cavalcante de Sá, M.C.; Córdova, A.M.T.; Díaz de León Gómez, R.E.; Pinto, J.C. Modeling of Isoprene Solution Coordinative Chain Transfer Polymerization. Macromol. React. Eng. 2021, 15, 2100005. [Google Scholar] [CrossRef]
  33. Córdova, T.; Enríquez-Medrano, F.J.; Cartagena, E.M.; Villanueva, A.B.; Valencia, L.; Álvarez, E.N.C.; González, R.L.; Díaz-de-León, R. Coordinative Chain Transfer Polymerization of Sustainable Terpene Monomers Using a Neodymium-Based Catalyst System. Polymers 2022, 14, 2907. [Google Scholar] [CrossRef]
  34. Tereshchenko, K.A.; Ulitin, N.V.; Bedrina, P.S.; Shiyan, D.A.; Lifanov, A.D.; Madzhidov, T.I.; Rusanova, S.N.; Volfson, S.I. Analysis of the Mechanism of Polybutadiene Synthesis in the Presence of the Neodymium Versatate+ Diisobutylaluminum Hydride+ Ethylaluminum Sesquichloride Catalytic System within the Solution of the Inverse Kinetic Problem. Ind. Eng. Chem. Res. 2022, 61, 15961–15969. [Google Scholar] [CrossRef]
  35. Iovu, H.; Hubca, G.; Simionescu, E.; Badea, E.G.; Dimonie, M. Polymerization of butadiene and isoprene with the NdCl3· 3TBP-TIBA catalyst system. Die Angew. Makromol. Chem. Appl. Macromol. Chem. Phys. 1997, 249, 59–77. [Google Scholar] [CrossRef]
  36. Srinivasa Rao, G.S.; Upadhyay, V.K.; Jain, R.C. Polymerization of 1,3-butadiene using neodymium chloride tripentanolate–triethyl aluminum catalyst systems. J. Appl. Polym. Sci. 1999, 71, 595–602. [Google Scholar] [CrossRef]
  37. Ren, C.; Li, G.; Dong, W.; Jiang, L.; Zhang, X.; Wang, F. Soluble neodymium chloride 2-ethylhexanol complex as a highly active catalyst for controlled isoprene polymerization. Polymer 2007, 48, 2470–2474. [Google Scholar] [CrossRef]
  38. Hu, Y.; Zhang, C.; Liu, X.; Gao, K.; Cao, Y.; Zhang, C.; Zhang, X. Methylaluminoxane-activated neodymium chloride tributylphosphate catalyst for isoprene polymerization. J. Appl. Polym. Sci. 2014, 131. [Google Scholar] [CrossRef]
  39. Kularatne, R.N.; Yang, A.; Nguyen, H.Q.; McCandless, G.T.; Stefan, M.C. Neodymium catalyst for the polymerization of dienes and polar vinyl monomers. Macromol. Rapid Commun. 2017, 38, 1700427. [Google Scholar] [CrossRef]
  40. Oehme, A.; Gebauer, U.; Gehrke, K.; Beyer, P.; Hartmann, B.; Lechner, M.D. The influence of the catalyst preparation on the homo-and copolymerization of butadiene and isoprene. Macromol. Chem. Phys. 1994, 195, 3773–3781. [Google Scholar] [CrossRef]
  41. Boisson, C.; Barbotin, F.; Spitz, R. Polymerization of butadiene with a new catalyst based on a neodymium amide precursor. Macromol. Chem. Phys. 1999, 200, 1163–1166. [Google Scholar] [CrossRef]
  42. Friebe, L.; Nuyken, O.; Windisch, H.; Obrecht, W. Polymerization of 1,3-Butadiene Initiated by Neodymium Versatate/Diisobutylaluminium Hydride/Ethylaluminium Sesquichloride: Kinetics and Conclusions About the Reaction Mechanism. Macromol. Chem. Phys. 2002, 203, 1055–1064. [Google Scholar] [CrossRef]
  43. Salakhov, I.I.; Akhmetov, I.G.; Kozlov, V.G. Polymerization of butadiene during the action of the catalytic system neodymium versatate-diisobutylaluminum hydride-hexachloro-p-xylene. Polym. Sci. Ser. B 2011, 53, 385–390. [Google Scholar] [CrossRef]
  44. Monakov, Y.B.; Sabirov, Z.M.; Marina, N.G. The nature of active centers and key stages in diene polymerization with lanthanide catalytic systems. Vysokomolekulyarnye Soedineniya. Seriya A Seriya B 1995, 38, 407–417. [Google Scholar]
  45. Sigaeva, N.N.; Usmanov, T.S.; Budtov, V.P.; Spivak, S.I.; Monakov, Y.B. Catalytic activity distribution of active centers in lanthanide-based diene polymerization. Polym. Sci. Ser. B 2000, 42, 1–5. [Google Scholar]
  46. Sabirov, Z.M.; Urazbaev, V.N.; Efimov, V.P.; Mullagaliev, I.R.; Monakov, Y.B. On steroespecificity of the effect of Bi-and monometallic active centers in polymerization of butadiene initiated by lanthanide catalysts. In Doklady Physical Chemistry; Pleiades Publishing, Ltd.: New York, NY, USA, 2000; Volume 372, pp. 89–91. [Google Scholar]
  47. Monakov, Y.B.; Sabirov, Z.M. Basic stages of microstructures formation in polydienes upon lanthanide catalytic systems. Int. J. Polym. Mater. 2001, 50, 57–84. [Google Scholar] [CrossRef]
  48. Sigaeva, N.N.; Usmanov, T.S.; Budtov, V.P.; Spivak, S.I.; Monakov, Y.B. Determination of a cumulative distribution function of kinetic activity of lanthanide ion-coordination catalytic systems used for diene polymerization. Int. J. Polym. Mater. 2001, 49, 475–489. [Google Scholar] [CrossRef]
  49. Monakov, B.Y.; Sigaeva, N.N. Active site multiplicity of catalytic systems in polymerization processes. Polym. Sci. Ser. C 2001, 43, 61–80. [Google Scholar]
  50. Monakov, Y.B.; Sabirov, Z.M.; Duvakina, N.V.; Glukhov, E.A.; Ponomareva, O.A. Nature of active centers in diene copolymerization initiated by Ziegler-Natta catalysts: The role of a non-transition-metal organic compound. Polym. Sci. Ser. A 2001, 43, 91–95. [Google Scholar]
  51. Sigaeva, N.N.; Usmanov, T.S.; Budtov, V.P.; Monakov, Y.B. Effect of organoaluminum compound on kinetic nonuniformity and structure of active centers of neodymium catalytic systems in butadiene polymerization. Russ. J. Appl. Chem. 2001, 74, 1141–1146. [Google Scholar] [CrossRef]
  52. Monakov, Y.B.; Sabirov, Z.M.; Urazbaev, V.N.; Efimov, V.P. Relationship between the stereospecificity of lanthanide catalysts and the structures of active sites and dienes, the nature of a cocatalyst, and preparation conditions. Kinet. Catal. 2001, 42, 310–316. [Google Scholar] [CrossRef]
  53. Monakov, Y.B.; Sabirov, Z.M.; Urazbaev, V.N.; Efimov, V.P. Diene polymerization initiated by NdCl3. 3TBP-based catalytic systems. Multiplicity of active centers and their structure and stereospecificity distributions. Polym. Science. Ser. A 2002, 44, 228–231. [Google Scholar]
  54. Sigaeva, N.N.; Usmanov, T.S.; Budtov, V.P.; Spivak, S.I.; Zaikov, G.E.; Monakov, Y.B. The influence of the nature of organoaluminum compound on kinetic heterogeneity of active sites in lanthanide-based diene polymerization. J. Appl. Polym. Sci. 2003, 87, 358–368. [Google Scholar] [CrossRef]
  55. Urazbaev, V.N.; Efimov, V.P.; Sabirov, Z.M.; Monakov, Y.B. Structure of active centers, their stereospecificity distribution, and multiplicity in diene polymerization initiated by NdCl3-based catalytic systems. J. Appl. Polym. Sci. 2003, 89, 601–603. [Google Scholar] [CrossRef]
  56. Tobisch, S. Theoretical investigation of the mechanism of cis−trans regulation for the allylnickel(II)-catalyzed 1,4 polymerization of butadiene. Acc. Chem. Res. 2002, 35, 96–104. [Google Scholar] [CrossRef] [PubMed]
  57. Tobisch, S. Mechanism of the chain termination of the allylnickel(II)-catalyzed polymerization of 1,3-butadiene. A density functional investigation for the cationic [NiII(RC3H4)(cis-C4H6)L]+ active catalyst. Macromolecules 2003, 36, 6235–6244. [Google Scholar] [CrossRef]
  58. Tobisch, S. Mechanistic insight into the selective trans-1,4-polymerization of butadiene by terpyridine–iron(II) complexes—A computational study. Can. J. Chem. 2009, 87, 1392–1405. [Google Scholar] [CrossRef]
  59. Lin, F.; Wang, M.; Pan, Y.; Tang, T.; Cui, D.; Liu, B. Sequence and regularity controlled coordination copolymerization of butadiene and styrene: Strategy and mechanism. Macromolecules 2017, 50, 849–856. [Google Scholar] [CrossRef]
  60. Bahri-Laleh, N.; Hanifpour, A.; Mirmohammadi, S.A.; Poater, A.; Nekoomanesh-Haghighi, M.; Talarico, G.; Cavallo, L. Computational modeling of heterogeneous Ziegler-Natta catalysts for olefins polymerization. Prog. Polym. Sci. 2018, 84, 89–114. [Google Scholar] [CrossRef]
  61. Guo, Z.A.; Xian, J.Y.; Rong, L.R.; Qin, H.; Jie, Z. Theoretical study of metal ion impact on geometric and electronic properties of terbutaline compounds. Mon. Für Chem.-Chem. Mon. 2019, 150, 1355–1364. [Google Scholar] [CrossRef]
  62. Nsiri, H.; Belaid, I.; Larini, P.; Thuilliez, J.; Boisson, C.; Perrin, L. Ethylene–Butadiene Copolymerization by Neodymocene Complexes: A Ligand Structure/Activity/Polymer Microstructure Relationship Based on DFT Calculations. ACS Catal. 2016, 6, 1028–1036. [Google Scholar] [CrossRef]
  63. Romano, E.; Budzelaar, P.H.; De Rosa, C.; Talarico, G. Unconventional Stereoerror Formation Mechanisms in Nonmetallocene Propene Polymerization Systems Revealed by DFT Calculations. J. Phys. Chem. A 2022, 126, 6203–6209. [Google Scholar] [CrossRef]
  64. Manuyko, G.V.; Aminova, G.A.; Dyakonov, G.S.; Ahmetov, I.G.; Bronskaya, V.V. Kinetic heterogeneity of neodymium catalyst systems modified with methylaluminoxane. Theor. Found. Chem. Eng. 2015, 49, 246–251. [Google Scholar] [CrossRef]
  65. Busygin, V.M.; Gil’manov, K.K.; Gil’mutdinov, N.R.; Akhmetov, I.G.; Salakhov, I.I.; Akhmetova, D.R.; Vagizov, A.M. Method of Producing Polymers and Copolymers of Conjugated Dienes (Versions). RU Patent 2422468, 27 June 2011. [Google Scholar]
  66. Neese, F. The ORCA program system. Wiley Interdiscip. Rev. Comput. Mol. Sci. 2012, 2, 73–78. [Google Scholar] [CrossRef]
  67. Neese, F. Software update: The ORCA program system, version 4.0. Wiley Interdiscip. Rev. Comput. Mol. Sci. 2017, 8, e1327. [Google Scholar] [CrossRef]
  68. Becke, A. Density-functional thermochemistry. III. The role of exact exchange. J. Chem. Phys. 1993, 98, 5648. [Google Scholar] [CrossRef] [Green Version]
  69. Lee, C.; Yang, W.; Parr, R.G. Development of the Colle-Salvetti correlation-energy formula into a functional of the electron density. Phys. Rev. B 1988, 37, 785. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  70. Weigend, F.; Ahlrichs, R. Balanced basis sets of split valence, triple zeta valence and quadruple zeta valence quality for H to Rn: Design and assessment of accuracy. Phys. Chem. Chem. Phys. 2005, 7, 3297–3305. [Google Scholar] [CrossRef]
  71. Stoychev, G.L.; Auer, A.A.; Neese, F. Automatic generation of auxiliary basis sets. J. Chem. Theory Comput. 2017, 13, 554–562. [Google Scholar] [CrossRef]
  72. Dolg, M.; Stoll, H.; Preuss, H. Energy-adjusted abinitio pseudopotentials for the rare earth elements. J. Chem. Phys. 1989, 90, 1730–1734. [Google Scholar] [CrossRef]
  73. Grimme, S.; Antony, J.; Ehrlich, S.; Krieg, H. A consistent and accurate ab initio parametrization of density functional dispersion correction (DFT-D) for the 94 elements H-Pu. J. Chem. Phys. 2010, 132, 154104. [Google Scholar] [CrossRef] [Green Version]
  74. Barone, V.; Cossi, M. Quantum calculation of molecular energies and energy gradients in solution by a conductor solvent model. J. Phys. Chem. A 1998, 102, 1995–2001. [Google Scholar] [CrossRef]
  75. Cammi, R.; Mennucci, B.; Tomasi, J. Fast evaluation of geometries and properties of excited molecules in solution: A Tamm-Dancoff model with application to 4-dimethylaminobenzonitrile. J. Phys. Chem. A 2000, 104, 5631–5637. [Google Scholar] [CrossRef]
  76. Dapprich, S.; Komáromi, I.; Byun, K.S.; Morokuma, K.; Frisch, M.J. A new ONIOM implementation in Gaussian98. Part I. The calculation of energies, gradients, vibrational frequencies and electric field derivatives. J. Mol. Struct. THEOCHEM 1999, 461, 1–21. [Google Scholar] [CrossRef]
  77. Grimme, S.; Bannwarth, C.; Shushkov, P. A robust and accurate tight-binding quantum chemical method for structures, vibrational frequencies, and noncovalent interactions of large molecular systems parametrized for all spd-block elements (Z = 1–86). J. Chem. Theory Comput. 2017, 13, 1989–2009. [Google Scholar] [CrossRef] [PubMed]
  78. Bannwarth, C.; Ehlert, S.; Grimme, S. GFN2-xTB—An accurate and broadly parametrized self-consistent tight-binding quantum chemical method with multipole electrostatics and density-dependent dispersion contributions. J. Chem. Theory Comput. 2019, 15, 1652–1671. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  79. Ehlert, S.; Stahn, M.; Spicher, S.; Grimme, S. Robust and efficient implicit solvation model for fast semiempirical methods. J. Chem. Theory Comput. 2021, 17, 4250–4261. [Google Scholar] [CrossRef]
  80. Anno, T. Out-of-Plane Vibrations of a Conjugated Hydrocarbon: Trans-Butadiene. J. Chem. Phys. 1958, 28, 944–949. [Google Scholar] [CrossRef]
Polymers 15 01166 g001 550
Figure 1. Active site (AS+C1H2C2HC3HC4H2) coordinated to the cis- (a) and trans- (b) forms of 1,3-butadiene (P denotes the i-C4H9 or growing polybutadienyl chain).
Figure 1. Active site (AS+C1H2C2HC3HC4H2) coordinated to the cis- (a) and trans- (b) forms of 1,3-butadiene (P denotes the i-C4H9 or growing polybutadienyl chain).
Polymers 15 01166 g001
Polymers 15 01166 g002 550
Figure 2. Structures of AS+η-cis-C4H6 and AS+η-trans-C4H6 optimized at the B3LYP/def2-TZVP level.
Figure 2. Structures of AS+η-cis-C4H6 and AS+η-trans-C4H6 optimized at the B3LYP/def2-TZVP level.
Polymers 15 01166 g002
Polymers 15 01166 sch001 550
Scheme 1. Structurally different modes of butadiene coordination at the metal M in butadiene π-complexes [56].
Scheme 1. Structurally different modes of butadiene coordination at the metal M in butadiene π-complexes [56].
Polymers 15 01166 sch001
Polymers 15 01166 sch002 550
Scheme 2. Different configurations and coordination modes of the butenyl anion together with their interconversion (anti-syn isomerization) via a σ-C3 intermediate (P denotes the growing polybutadienyl chain) [56].
Scheme 2. Different configurations and coordination modes of the butenyl anion together with their interconversion (anti-syn isomerization) via a σ-C3 intermediate (P denotes the growing polybutadienyl chain) [56].
Polymers 15 01166 sch002
Polymers 15 01166 g003 550
Figure 3. Structures of i-C4H9-π-anti-H2CHCγHCβH2Cα-AS and i-C4H9-π-syn-H2CHCγHCβH2Cα-AS optimized at the B3LYP/def2-TZVP level.
Figure 3. Structures of i-C4H9-π-anti-H2CHCγHCβH2Cα-AS and i-C4H9-π-syn-H2CHCγHCβH2Cα-AS optimized at the B3LYP/def2-TZVP level.
Polymers 15 01166 g003
Polymers 15 01166 g004 550
Figure 4. Structures of i-C4H9-π-anti-H2CHCγHCβH2Cα-AS+η-cis-C4H6 and i-C4H9-π-anti-H2CHCγHCβH2Cα-AS+η-trans-C4H6 optimized at the B3LYP/def2-TZVP level.
Figure 4. Structures of i-C4H9-π-anti-H2CHCγHCβH2Cα-AS+η-cis-C4H6 and i-C4H9-π-anti-H2CHCγHCβH2Cα-AS+η-trans-C4H6 optimized at the B3LYP/def2-TZVP level.
Polymers 15 01166 g004
Polymers 15 01166 g005 550
Figure 5. Energy diagram for the polymerization of 1,3-butadiene on AS. The dotted line corresponds to the generation cycle of the cis-1,4-polybutadiene units, the solid line corresponds to the trans-1,4-polybutadiene units. Ini—initial state of P-AS-trans-C4H6, TS1—trans–cis 1,3-butadiene transition, Int—intermediate P-AS-cis-C4H6, TS2—transition state of 1,3-butadene insertion into polymer chain, Fin—final state.
Figure 5. Energy diagram for the polymerization of 1,3-butadiene on AS. The dotted line corresponds to the generation cycle of the cis-1,4-polybutadiene units, the solid line corresponds to the trans-1,4-polybutadiene units. Ini—initial state of P-AS-trans-C4H6, TS1—trans–cis 1,3-butadiene transition, Int—intermediate P-AS-cis-C4H6, TS2—transition state of 1,3-butadene insertion into polymer chain, Fin—final state.
Polymers 15 01166 g005
Table
Table 1. Some geometric characteristics of AS+η-cis-C4H6 and AS+η-trans-C4H6 structures optimized at the B3LYP/def2-TZVP level.
Table 1. Some geometric characteristics of AS+η-cis-C4H6 and AS+η-trans-C4H6 structures optimized at the B3LYP/def2-TZVP level.
AS+η-cis-C4H6AS+η-trans-C4H6
R(Nd-Cl), Å2.852.87
R(Nd-C(i-C4H9)), Å2.422.39
R(Nd-C1(C4H6)), Å3.283.36
R(Nd-C2(C4H6)), Å3.033.08
R(Nd-C3(C4H6)), Å3.093.14
R(Nd-C4(C4H6)), Å3.303.51
R(Al-Cl), Å2.312.31
R(C(C4H6)-C(i-C4H9), Å3.223.64
∠Cl-Nd-Cl, degrees74.0, 84.073.0, 85.0
∠Cl-Al-Cl, degrees95.794.9
∠Nd-Cl-Al, degrees94.094.0
∠Nd-Cl-Cl-Al, degrees171.0164.0
∠Cl-Cl-Nd-Cl, degrees87.087.0
Table
Table 2. Some geometric characteristics of i-C4H9-π-anti-H2CHCγHCβH2Cα-AS and i-C4H9-π-syn-H2CHCγHCβH2Cα-AS structures optimized at the B3LYP/def2-TZVP level.
Table 2. Some geometric characteristics of i-C4H9-π-anti-H2CHCγHCβH2Cα-AS and i-C4H9-π-syn-H2CHCγHCβH2Cα-AS structures optimized at the B3LYP/def2-TZVP level.
i-C4H9-π-anti-H2CHCγHCβH2Cα-ASi-C4H9-π-syn-H2CHCγHCβH2Cα-AS
R(Nd-Cl), Å2.822.82
R(Nd-Cα(C4H6)), Å2.622.50
R(Nd-Cβ(C4H6)), Å2.672.71
R(Nd-Cγ(C4H6)), Å2.642.84
R(Nd-Cδ(C4H6)), Å2.963.95
R(Al-Cl), Å2.322.32
R(C(C4H9)-C(i-C4H9), Å1.541.54
∠Cl-Nd-Cl, degrees74.0, 91.074.0, 83.0
∠Cl-Al-Cl, degrees95.093.0
∠Nd-Cl-Al, degrees93.094.0, 91.0
∠Nd-Cl-Cl-Al, degrees172.0176.0, 146.0
∠Cl-Cl-Nd-Cl, degrees76.077.0
Table
Table 3. Some geometric characteristics of i-C4H9-π-anti-H2CHCγHCβH2Cα-AS+η-cis-C4H6 and i-C4H9-π-anti-H2CHCγHCβH2Cα-AS+η-trans-C4H6 structures optimized at the B3LYP/def2-TZVP level.
Table 3. Some geometric characteristics of i-C4H9-π-anti-H2CHCγHCβH2Cα-AS+η-cis-C4H6 and i-C4H9-π-anti-H2CHCγHCβH2Cα-AS+η-trans-C4H6 structures optimized at the B3LYP/def2-TZVP level.
i-C4H9-π-anti-H2CHCγHCβH2Cα-AS+η-cis-C4H6i-C4H9-π-anti-H2CHCγHCβH2Cα-AS+η-trans-C4H6
R(Nd-Cl), Å2.872.90
R(Nd-C1(C4H6)), Å3.073.42
R(Nd-C2(C4H6)), Å3.253.10
R(Nd-C3(C4H6)), Å3.963.15
R(Nd-C4(C4H6)), Å4.833.51
R(Nd-Cα(C4H6)), Å2.682.67
R(Nd-Cβ(C4H6)), Å2.672.68
R(Nd-Cγ(C4H6)), Å2.602.64
R(Nd-Cδ(C4H6)), Å2.973.06
R(Al-Cl), Å2.302.30
R(C1(C4H6)-Cδ(C4H6), Å3.933.58
∠Cl-Nd-Cl, degrees74.0, 105.073.0, 87.0
∠Cl-Al-Cl, degrees95.095.0
∠Nd-Cl-Al, degrees93.695.0
∠Nd-Cl-Cl-Al, degrees168.0, 150.0175.0, 165.0
∠Cl-Cl-Nd-Cl, degrees77.079.0
Table
Table 4. Changes in the activation energy and Gibbs free energy under standard conditions (kJ/mol) as a result of the introduction of 1,3-butadiene into the growing polymer chain.
Table 4. Changes in the activation energy and Gibbs free energy under standard conditions (kJ/mol) as a result of the introduction of 1,3-butadiene into the growing polymer chain.
Initial StructureΔG≠
(trans–cis)		ΔG
(trans–cis)		ΔG≠
(Polymer)		ΔG≠
(Total)		ΔG
(Polymer)		Product
I stage
i-C4H9-AS-T43.06.654.360.9−51.7i-C4H9-C-AS
i-C4H9-AS-T 67.067.0−55.1i-C4H9-T-AS
II stage
i-C4H9-C-AS-T4119.442.361.7−38.5i-C4H9-CC-AS
i-C4H9-C-AS-T--75.075.0−32.9i-C4H9-CT-AS
i-C4H9-T-AS-T41.929.238.868.0−46.1i-C4H9-TC-AS
i-C4H9-T-AS-T--72.572.5−39.2i-C4H9-TT-AS
III stage
i-C4H9-CC-AS-T24.614.651.766.3−35.0i-C4H9-CCC-AS
i-C4H9-CC-AS-T--102.0102.0−52.8i-C4H9-CCT-AS
i-C4H9-TT-AS-T--77.577.5−40.4i-C4H9-TTT-AS
IV stage
i-C4H9-CCC-AS-T24.76.354.260.5−60.4i-C4H9-CCCC-AS
i-C4H9-CCC-AS-T--80.180.1−51.6i-C4H9-CCCT-AS
i-C4H9-TTT-AS-T--72.472.4−35.7i-C4H9-TTTT-AS
V stage
i-C4H9-CCCC-AS-T24.811.350.261.5−45.0i-C4H9-CCCCC-AS
i-C4H9-CCCC-AS-T--77.377.3−57.2i-C4H9-CCCCT-AS
i-C4H9-TTTT-AS-T--72.472.4−35.7i-C4H9-TTTTT-AS
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Masliy, A.N.; Akhmetov, I.G.; Kuznetsov, A.M.; Davletbaeva, I.M. DFT and ONIOM Simulation of 1,3-Butadiene Polymerization Catalyzed by Neodymium-Based Ziegler–Natta System. Polymers 2023, 15, 1166. https://doi.org/10.3390/polym15051166

AMA Style

Masliy AN, Akhmetov IG, Kuznetsov AM, Davletbaeva IM. DFT and ONIOM Simulation of 1,3-Butadiene Polymerization Catalyzed by Neodymium-Based Ziegler–Natta System. Polymers. 2023; 15(5):1166. https://doi.org/10.3390/polym15051166

Chicago/Turabian Style

Masliy, Alexey N., Ildar G. Akhmetov, Andrey M. Kuznetsov, and Ilsiya M. Davletbaeva. 2023. "DFT and ONIOM Simulation of 1,3-Butadiene Polymerization Catalyzed by Neodymium-Based Ziegler–Natta System" Polymers 15, no. 5: 1166. https://doi.org/10.3390/polym15051166

APA Style

Masliy, A. N., Akhmetov, I. G., Kuznetsov, A. M., & Davletbaeva, I. M. (2023). DFT and ONIOM Simulation of 1,3-Butadiene Polymerization Catalyzed by Neodymium-Based Ziegler–Natta System. Polymers, 15(5), 1166. https://doi.org/10.3390/polym15051166

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop