KINETICS OF THE CONFORMATIONAL CHANGES OF THE NITROGEN MATRIX ISOLATED DISUBSTITUTED SILACYCLOHEXANES: EXPERIMENT MEETS THEORY AT FINITE TEMPERATURE RANGE

Joanna Stocka1, Pawel Rodziewicz2, Valdas Sablinskas1, Gamil A. Guirgis3

1 Institute of Chemical Physics, Vilnius University, Saulėtekio al. 3, 10257 Vilnius, Lithuania

2 Institute of Chemistry, Jan Kochanowski University, 15G Świętokrzyska St., 25-406, Kielce, Poland

3 Department of Chemistry and Biochemistry, College of Charleston, Charleston, SC 29424, USA

[email protected]

Organosilicon compounds find a usage in the field of surface science due to its good surface adhesion properties which are mostly related to a $\pi$ electron orbital [1]. The substitution of the carbon atom in the ring with a silicon atom enhances adhesion since it acts as hydrolytically sensitive center that can react with inorganic substrates such as glass to form stable covalent bonds [2]. 1-chloromethyl-1-fluorosilacyclohexane and 1-chloro-1-chloromethylsilacyclohexane were newly synthesized molecular compounds with unknown structural parameters and conformational diversity. In our previous works the detailed conformational, structural and spectroscopic analysis of the latter mentioned molecules is presented. [3, 4] The next step is theoretical studies of the dynamical structural changes and lifetime of different conformers at finite temperature with the help of Car-Parrinello molecular dynamics simulations.

Car-Parrinello molecular dynamics simulations have been performed with PBE functionals, DFT (PW), Nosé–Hoover chain thermostats, temperature was set to 20, 50 and 70 K (to mimic experimental conditions), the production run of 200 ps was derived utilizing the CPMD program package.

Figure 1
Fig. 1. Snapshots from the CMPD simulations at 20K for 1-Chloro-1-Chloromethyl-1-silacyclohexane.

In the matrix the molecule remains in the global energy minimum structure (chair axial trans), the structural parameters and the interaction with the matrix were traced by Radial Pair Distribution Function and bond distance analysis.

The biggest advantage of the Molecular Dynamics Simulations is that they explicitly describe the experimental conditions, due to the temperature influence. The static calculations in most cases, describe the isolated molecules without taking into account the temperature impact. This is a huge step in the computational description of the experimental conditions. Despite the fact the molecule does not undergo any conformational changes, the molecule – matrix interactions were traced.


[1] H. Butt, K Graf, M. Kappl, Physics and Chemistry of Interfaces, Wiley, 2003

[2] J. Ceponkus, V. Sablinskas, V.Aleksa, M.Pucetaite, R.Platakyte, C.W.Reed, C.Cotter, G.Guirgis, Vibrational Spectroscopy, 81, 136-143, 2015

[3] T.M.C McFadden, R. Platakyte, J. Stocka, J. Ceponkus, V. Aleksa, T. Carrigan-Broda, V. Sablinskas, P. Rodziewicz, G.A. Guirgis, Experimental (Raman and IR) and computational (DFT, MP2) studies of conformational diversity of 1-chloromethyl-1-fluorosilacyclohexane, Journal of Molecular Structure, 1221 (2020).

[4] J. Stocka, R. Platakyte, T.M.C. McFadden, J. Ceponkus, V. Aleksa, A.G. Hanna, V. Sablinskas, P. Rodziewicz, G.A. Guirgis, Conformational diversity of 1-chloro-1-chloromethylsilacyclohexane with experimental (Raman and IR) and computational (DFT, MP2) methods, Journal of Molecular Structure, 1249 (2022).