Odd–Even Effects in the Dynamics of Liquid Crystalline Thin Films on

Jul 30, 2015 - Study of the Dynamics of 5CB Thin Layer Placed on the Fullerene Wall: Computer Simulations. Przemysław Raczyński , Zygmunt Gburski. 2...
0 downloads 0 Views 2MB Size
Article pubs.acs.org/JPCC

Odd−Even Effects in the Dynamics of Liquid Crystalline Thin Films on the Surface of Single Walled Carbon and Silicon Carbide Nanotubes: Computer Simulation Study Krzysztof Górny,* Przemysław Raczyński, Zbigniew Dendzik, and Zygmunt Gburski Institute of Physics, University of Silesia, Uniwersytecka 4, 40-007 Katowice, Poland ABSTRACT: Translational and reorientational dynamics of monolayer films of ncyanobiphenyl (with n = 5−8) series of mesogens on the surface of single walled carbon nanotubes and heterogeneous silicon carbide nanotubes has been studied using molecular dynamics computer simulation. Activation characteristic of diffusion of the mesogens on the surface of the heterogeneous nanotube is selective with respect to the number of carbon atoms in the aliphatic tail of the mesogen, opposite to the case of the mesogens on the surface of the homogeneous carbon nanotube, which reveals nonselective activation. The orientational order parameter ⟨P2⟩ and thermal activation energy of dipolar relaxation of the studied systems also exhibit a characteristic odd−even pattern with respect to the length of the hydrocarbon chain of the mesogen molecule.



INTRODUCTION Molecular systems encapsulated in carbon nanotubes1 or physisorbed on their surface2 have been studied for both scientific reasons and intersting potential applications.3,4 Besides homogeneous carbon nanotubes (CNTs), also a variety of heterogeneous nanotubes have been studied in this context. Molecular systems based on silicon carbide nanotubes (SiCNTs) have been studied because of interesting potential applications, such as desalination of water5 or removal of cholesterol molecules from a protein surface,6 to name but a few. Liquid crystals (LCs) are a phase of matter that combines the mobility of an isotropic liquid with the long-range orientational order, which is normally associated with crystalline solids, and have been used for applications that rely on their anisotropic electrooptical properties. The first cyanobiphenyl, 5CB, was synthesized in 1972 as an outcome of research aiming at finding a mesogen with the specific intention of using them in liquid crystal displays. The series of n-cyanobiphenyls (nCB; 4-n-alkyl4′-cyanobiphenyl) consists of important liquid crystals, with their liquid crystalline phases at the temperatures near room temperature, which have intensively been studied through experiment and with use of computer simulation methods.7 Because the ordering of molecules is known to be essentially influenced by interfaces and dopants,8 anchoring of the liquid crystals and other functional organic molecules at a variety of interfaces has been studied for fundamental scientific reasons but also because the molecular organization of liquid crystals at their interfaces is essential for designing and optimizing a variety of devices ranging from displays to organic field effect transistors and implementing new promising concepts in nanoand biotechnology.9 © 2015 American Chemical Society

Anchoring of liquid crystals at their interfaces with such substrates as graphene,10 carbon nanotubes,11 or polymer nanofibers12 proved to essentially affect their orientational order and their electrooptic and dielectric properties. In a series of publications, Zannoni and co-workers have investigated orientation and ordering of mesogen films at their interfaces with a variety of substrates, such as crystalline silicon,13 cristobalite and glassy silica,14 and fullerene planar structure15 as well as studying the impact of the cylindrical geometric confinement.16,17 Alignment of carbon nanotubes, which can be viewed as highly anisotropic rigid particles, is crucial for properties of nanotube based materials. Carbon nanotubes dispersed in liquid crystalline medium may form spontaneously ordered aggregates.18−21 Reorientational dynamics of polar molecular systems in confined geometries is often studied in terms of dipolar relaxation.1,2,22 Broadband dielectric relaxation spectroscopy provides the link between the methods which probe the properties of the individual molecules and techniques characterizing the bulk properties of the sample.23 Relaxation spectra are sensitive to intermolecular interactions and provide information on cooperative processes, revealing features mutual to different types of condensed matter systems.24,25



SIMULATION DETAILS The models of 5CB, 6CB, 7CB, and 8CB liquid crystals molecules have been adopted from the CHARMM-type united atom potential developed by Tiberio and co-workers,7 and have Received: June 22, 2015 Revised: July 29, 2015 Published: July 30, 2015 19266

DOI: 10.1021/acs.jpcc.5b05961 J. Phys. Chem. C 2015, 119, 19266−19271

Article

The Journal of Physical Chemistry C been fitted against experimental data and proved to well describe the isotropic−nematic phase transition temperatures and the appropriate values of the order parameter assigned to these phases. The homogeneous carbon nanotube as well as the silicon carbide nanonube have been modeled using all atoms potentials described elswhere.1,5,6 Interactions between the mesogens have been described taking into account electrostatics and van der Waals forces modeled with Lennard-Jones 12−6 potential with Lorentz−Berthelot mixing rules and cutoff of 20 Å. Interactions between SiCNTs or CNTs and mesogens have been described with a model that accouts for van der Waals interactions modeled with Lennard-Jones 12−6 potential. Periodic boundary conditions were applied and long-range interactions were calculated using particle mesh Ewald (PME)26 summation technique with grid spacing of 1.5 Å. Equations of motion were integrated using the Brunger− Brooks−Karplus (BBK) scheme,27 implemented in NAMD 2.9,28 with the time step of integration of 1 fs. The simulation was performed in NVT ensemble, for temperatures 270, 290, 310, 330, and 350 K. Temperature was controlled using the Langevin thermostat with damping coefficient γ = 5.0 ps−1. In order to generate the initial configurations of mesogen monolayers physisorbed on the surface of each nanotube, we performed a series of preliminary NPT simulations of nanotubes immersed in bulk liquid crystals. Subsequently, we imposed cylindrical geometrical constraints to remove mesogens which do not belong to the immediate vicinity of the surface of the nanotube. Then we performed auxiliary NVT simulations to test whether the mesogen molecules form a uniformly distributed and continuous monolayer, repeating the whole procedure if necessary. Depending on the specific mesogen and the nanotube, the simulation cell contained monolayers that consisted of different numbers of molecules. For the system based on the homogeneous CNT, the simulation cell of 100 × 100 × 57 Å3 contained a single (10,10) carbon nanotube (with diameter of approximately 13 Å) covered with the monolayer of 40 mesogens in the case of 5CB, 38 mesogens in the case of 6CB, 36 mesogens in the case of 7CB, and up to 34 mesogens in the case of 8CB. In the case of the heterogeneous SiCNT, the simulation cell of 100 × 100 × 74 Å3 contained a single (10,10) silicon carbide nanotube (with diameter approximately 16 Å) covered with the monolayer of 59 mesogens in the case of 5CB, 56 mesogens in the case of of 6CB, 53 mesogens in the case of 7CB, and up to 50 mesogens in the case of 8CB. A snapshot of the monolayer of 6CB mesogens on the surface of the silicon carbide nanotube is shown in Figure 1. The initial states have been then used to perform the simulations of the mesogen monolayers on the surface of the nanotubes. The simulated system was thermalized for 5 ns before each trajectory production. The trajectories were produced over 20 000 000 time steps (20 ns).

Figure 1. 6CB mesogens on the SiCNT surface: a snapshot of the simulated system.

from the long-time limit of the mean square displacement defined as ⟨Δr 2(t )⟩ =

1 N

N

∑ |ri(t ) − ri(0)|2 i=1

(2)

which represents a canonical average over all intervals of the MD trajectory corresponding to a time t; N is the total number of mesogens, and ri is the position vector of the ith mesogen, determined from two frames separated by an interval t. Figure 2 shows the thermal activation plot of diffusion of liquid crystal molecules on the surface of carbon and silicon carbide nanotubes determined from calculated trajectories. We observed that silicon carbide nanotube affects the diffusion of nCB molecules in a selective way. Activation energy of nCB molecules increases with the number of carbon atoms (n) in the aliphatic tail of the molecule, from approximately EA = 13.5



RESULTS AND DISCUSSION We have studied dynamical properties of 5CB, 6CB, 7CB, and 8CB liquid crystal molecules on the surface of single walled carbon and silicon carbide nanotubes. First, we studied the diffusion characteristics of the studied liquid crystal layers. To do so, we calculated the diffusion coefficients, D, using the Einstein relation ⟨Δr 2(t )⟩ = 6Dt

Figure 2. Thermal activation characteristic of diffusion coefficient of 5CB (□), 6CB (○), 7CB (△), and 8CB (▽) molecules on the surface of the carbon nanotube and on the surface of the silicon carbide nanotube (analogous solid symbols), plotted versus the number of carbon atoms (n) in the aliphatic tail of the nCB molecule. The inset shows activation energies.

(1) 19267

DOI: 10.1021/acs.jpcc.5b05961 J. Phys. Chem. C 2015, 119, 19266−19271

Article

The Journal of Physical Chemistry C kJ/mol for 5CB, up to EA = 16.5 kJ/mol for 8CB mesogen molecule. On the other hand, the nCB series on the carbon nanotube interface does not exhibit the selective diffusion pattern. In the latter case, the activation energy of diffusion is approximately EA = 11 kJ/mol and does not depend on the number of carbon atoms in the aliphatic tail of the nCB molecule. Considerably higher values of activation energy in the case of the silicon carbide nanotube interface, as compared to the homogeneous carbon nanotube, suggest that in this case the diffusion mechanism is based on hopping between the potential minima imposed by the structure of the surface of the silicon carbide nanotube. This feature does not exist in the case of nCB on the surface of the carbon nanotube.

N

Q=

∑i = 1 [3û i(t ) ⊗ û i(t ) − I] 2N

(5)

where û(t) is an arbitrary selected molecular axis, I is the identity matrix; the sum runs over all the molecules forming a single layer on the surface of the nanotube, and the time averaging is calculated over all recorded time steps in the MD trajectories. The principal axis of inertia has been selected for û(t). For the systems with spherical and cylindrical geometries modifications of the order parameter were proposed,16,17 which corresponds to an arbitrary direction which accounts for the specific geometry confinements. In the case of both, CNT and SiCNT based nCB layers, the mesogens anchoring were planar and aligned along the nanotube axis in the whole range of temperatures. Figure 4

Figure 3. Lindemann index of 5CB (□), 6CB (○), 7CB (△), and 8CB (▽) molecules on the surface of the carbon nanotube and on the surface of the silicon carbide nanotube (analogous solid symbols), plotted versus the number of carbon atoms (n) in the aliphatic tail of the nCB molecule.

Figure 3 shows the temperature dependence of the Lindemann index, a measure of thermally driven disorder, defined as δi =

δ=

1 N−1 1 N

N

∑ j = 1, j ≠ i

⟨rij 2⟩ − ⟨rij⟩2 ⟨rij⟩

(3)

N

∑ δi i=1

(4)

calculated for all studied systems, where rij = |rij| is the distance between the ith and jth mesogen and N is the total number of mesogens. It does not exhibit discontinuities, which might suggest changes in the diffusion mechanism of a given system in the whole range of temperatures studied and indicates greater mobility of the liquid crystal molecules on the surface of the cabon nanotube, compared to the systems based on its silicon carbide counterpart. Orientational order of the system is usually characterized in terms of the order parameter ⟨P2⟩, expressing ordering of the liquid crystal molecules with respect to local director n(z). The values of ⟨P2⟩ and n can be calculated from the largest eigenvalue and corresponding eigenvector of the ordering matrix Q defined as14

Figure 4. Director components of 5CB mesogens on the surface of SiCNT at temperature T = 350 K.

shows the director components of the 5CB mesogens on the surface of the silicon carbide nanotube at the temperature T = 350 K, in Cartesian coordinates in which the z axis is oriented along the axis of the nanotube, which shows that the averaged director of the sample is parallel to the axis of the nanotube (only very slight precession of the director around the geometrical axis of the nanotube can be discerned). As a consequence, calculating the order parameter with respect to the geometrical axis of the nanotube would lead to only a very small decrease of the value of the order parameter, compared to the value obtained using the normally used definition of the 19268

DOI: 10.1021/acs.jpcc.5b05961 J. Phys. Chem. C 2015, 119, 19266−19271

Article

The Journal of Physical Chemistry C order parameter (eq 5). Because of the preceding, we decided to plot the order parameter corresponding to the instantaneous director of the mesogen system rather than that corresponding to the fixed axis aligned along the axis of the nanotube. Figure 5 presents the temperature dependence of ⟨P2⟩ for a series of nCB liquid crystal molecules on the surface of a single

Figure 5. Temperature dependence of ⟨P2⟩ for 5CB (□), 6CB (○), 7CB (△), and 8CB (▽) molecules on the surface of the carbon nanotube and on the surface of the silicon carbide nanotube (analogous solid symbols), plotted versus the number of carbon atoms (n) in the aliphatic tail of the nCB molecule.

walled carbon nanotube and its silicon carbide counterpart. For all studied systems, ⟨P2⟩ decreases with increasing temperature from 0.8 to 0.65. The liquid crystal molecules on the surface of the silicon carbide nanotube generally exhibit a higher degree of order than the molecules on the surface of the carbon nanotube. To gain deeper insight into the degree of order of different molecules of the nCB series, we replotted the values of the ⟨P2⟩ order parameter against the number of carbon atoms (n) in the aliphatic tail of a given nCB series molecule. Figure 6 shows ⟨P2⟩ versus n. For all temperatures studied, the odd− even effect occursthe degree of order which can be assigned to the layers composed of nCB molecules with an odd value of n are higher than that which can be assigned to their even counterparts. This feature agrees well with the experimentally and computationally observed odd−even pattern of nematic− isotropic transition temperatures.7 The odd−even effect is more pronounced for the mesogens on the surface of the silicon carbide nanotube than in the case of the systems based on the carbon nanotube (Figure 6). This can be attributed to the difference in the diffusion mechanism in the case of the mesogens on the surfaces of the SiCNT and CNT (Figure 2). In the system based on SiCNT, the adsorption energy landscape is more corrugated, compared to the case of CNT based systems, leading to the mechanism of diffusion which is based on hopping. Dipolar relaxation characteristic was studied in terms of normalized total dipole moment autocorrelation function (TDCF) Φ(t ) =

M(0) M(t ) M(0) M(0)

Figure 6. Values of ⟨P2⟩ order parameter for 5CB (□), 6CB (○), 7CB (△), and 8CB (▽) molecules on the surface of the carbon nanotube and on the surface of the silicon carbide nanotube (analogous solid symbols), plotted versus the number of carbon atoms (n) in the aliphatic tail of the nCB molecule, at the temperatures T = 270 K, T = 310 K, and T = 350 K. N

M(t ) =

∑ μi

(7)

i=1

using the procedure described in our previous works.1,29 The dipole moment of each mesogen was recalculated for each trajectory frame, with respect to the charge distribution and position of atoms in the molecule. TDCF was then plotted, − ln(− ln(Φ(t))) against ln(t), and fitted with Kohlrausch− Williams−Watts (KWW) decay function Φ(t ) = −(t /τKWW )β

(8)

where τKWW is a measure of the characteristic relaxation time and 0 < β < 1 is interpreted as a measure of distribution of the relaxation time or cooperativity of the relaxation process. The mean relaxation time, τ, has been calculated using the following relation 0

τ=

∫∞ Φ(t ) dt =

τKWW ⎛ 1 ⎞ Γ⎜ ⎟ β ⎝β⎠

(9)

where Γ denotes the gamma Euler function. The calculated values of the mean relaxation time were then plotted on the activation plot ln(τ) against 1/T. Figure 7 shows the dipolar relaxation thermal activation energies, EA′ , determined from the

(6)

where M is the total dipole moment of the system, defined as a sum of dipole moments μi of mesogens 19269

DOI: 10.1021/acs.jpcc.5b05961 J. Phys. Chem. C 2015, 119, 19266−19271

The Journal of Physical Chemistry C



REFERENCES

(1) Dendzik, Z.; Górny, K.; Gburski, Z. Cooperative Dipolar Relaxation of a Glycerol Molecular Cluster in Nanoscale Confinement - A Computer Simulation Study. J. Phys.: Condens. Matter 2009, 21, 425101. (2) Marti, J.; Gordillo, M. C. Structure and Dynamics of Liquid Water Adsorbed on the External Walls of Carbon Nanotubes. J. Chem. Phys. 2003, 119, 12540. (3) Dappe, Y. J. Encapsulation of Organic Molecules in Carbon Nanotubes: Role of the van der Waals Interactions. J. Phys. D: Appl. Phys. 2014, 47, 083001. (4) Rahman, M.; Lee, W. Scientific Duo of Carbon Nanotubes and Nematic Liquid Crystals. J. Phys. D: Appl. Phys. 2009, 42, 063001. (5) Garakani, F. M.; Kalantarinejad, R. A Molecular Dynamics Simulation of Water Transport through C and SiC Nanotubes: Application for Desalination. Int. J. Nano Dimens. 2012, 2, 151−157. (6) Raczyński, P.; Górny, K.; Samios, J.; Gburski, Z. Interaction Between Silicon-Carbide Nanotube and Cholesterol Domain. A Molecular Dynamics Simulation Study. J. Phys. Chem. C 2014, 118, 30115−30119. (7) Tiberio, G.; Muccioli, L.; Berardi, R.; Zannoni, C. Towards in Silico Liquid Crystals. Realistic Transition Temperatures and Physical Properties for n-Cyanobiphenyls via Molecular Dynamics Simulations. ChemPhysChem 2009, 10, 125−136. (8) Lin, I.-H.; Miller, D. S.; Bertics, P. J.; Murphy, C. J.; de Pablo, J. J.; Abbott, N. L. Endotoxin-Induced Structural Transformations in Liquid Crystalline Droplets. Science 2011, 332, 1297−1300. (9) Lagerwall, J. P.; Scalia, G. A New Era for Liquid Crystal Research: Applications of Liquid Crystals in Soft Matter Nano-, Bio- and Microtechnology. Curr. Appl. Phys. 2012, 12, 1387−1412. (10) Alam, T. M.; Pearce, C. J. Impact of Graphene Incorporation on the Orientational Order of Graphene/Liquid Crystal Composites. Chem. Phys. Lett. 2014, 592, 7−13. (11) Manjuladevi, V.; Gupta, R. K.; Kumar, S. Effect of Functionalized Carbon Nanotube on Electro-Optic and Dielectric Properties of a Liquid Crystal. J. Mol. Liq. 2012, 171, 60−63. (12) Rasna, M. V.; Zuhail, K. P.; Manda, R.; Paik, P.; Haase, W.; Dhara, S. Discontinuous Anchoring Transition and Photothermal Switching in Composites of Liquid Crystals and Conducting Polymer Nanofibers. Phys. Rev. E: Stat., Nonlinear, Soft Matter Phys. 2014, 89, 052503. (13) Pizzirusso, A.; Berardi, R.; Muccioli, L.; Ricci, M.; Zannoni, C. Predicting Surface Anchoring: Molecular Organization Across a Thin Film of 5CB Liquid Crystal on Silicon. Chem. Sci. 2012, 3, 573−579. (14) Roscioni, O. M.; Muccioli, L.; Della Valle, R. G. D.; Pizzirusso, A.; Ricci, M.; Zannoni, C. Predicting the Anchoring of Liquid Crystals at a Solid Surface: 5-Cyanobiphenyl on Cristobalite and Glassy Silica Surfaces of Increasing Roughness. Langmuir 2013, 29, 8950−8958. (15) D'Avino, G.; Muccioli, L.; Zannoni, C. From Chiral Islands to Smectic Layers: A Computational Journey Across Sexithiophene Morphologies on C60. Adv. Funct. Mater. 2015, 25, 1985−1995. (16) Chiccoli, C.; Pasini, P.; Evangelista, L. R.; Teixeira De Souza, R.; Zannoni, C. Nematics in Hybrid Cylindrical Cells. Mol. Cryst. Liq. Cryst. 2013, 576, 42−52. (17) Chiccoli, C.; Pasini, P.; Evangelista, L. R.; Teixeira-Souza, R. T.; Zannoni, C. Molecular Organization of Nematic Liquid Crystals Between Concentric Cylinders: Role of the Elastic Anisotropy. Phys. Rev. E: Stat., Nonlinear, Soft Matter Phys. 2015, 91, 022501. (18) van der Schoot, P.; Popa-Nita, V.; Kralj, S. Alignment of Carbon Nanotubes in Nematic Liquid Crystals. J. Phys. Chem. B 2008, 112, 4512−4518. (19) Ji, Y.; Huang, Y. Y.; Terentjev, E. M. Dissolving and Aligning Carbon Nanotubes in Thermotropic Liquid Crystals. Langmuir 2011, 27, 13254−13260. (20) Goncharuk, A. I.; Lebovka, N. I.; Lisetski, L. N.; Minenko, S. S. Aggregation, Percolation and Phase Transitions in Nematic Liquid Crystal EBBA Doped with Carbon Nanotubes. J. Phys. D: Appl. Phys. 2009, 42, 165411.

Figure 7. Thermal activation energies of dipolar relaxation of 5CB (□), 6CB (○), 7CB (△), and 8CB (▽) molecules on the surface of the carbon nanotube and on the surface of the silicon carbide nanotube (analogous solid symbols), plotted versus the number of carbon atoms (n) in the aliphatic tail of the nCB molecule.

activation plot for simulated mesogen monolayers, for different values of the length of the aliphatic tail of nCB molecules. The odd−even effect which occurs in the activation energy of the dipolar relaxation of the studied mesogen monolayers exhibits a different characteristic than in the analogous dependence of ⟨P2⟩. The activation energies of dipolar relaxation of mesogen systems composed of nCB molecules with an odd value of n are lower than that with even values of n. This can be attributed to the fact that, in the case of strongly ordered systems, the influence of thermal disorder on the reorientation is lower.



CONCLUSIONS We studied translational and reorientational dynamics of monolayers of n-cyanobiphenyl molecules with n ranging from 5 to 8, on the surface of single walled carbon and silicon carbide nanotubes. Contrary to the homogeneous carbon nanotube, the silicon carbide nanotube affects the diffusion of nCB mesogens in a selective waythe activation energy of diffusion increases with the increase of the number of carbon atoms in the nCB aliphatic tail. In the case of the mesogen monolayers on the surface of the silicon carbide nanotube, as well as in the case of the systems based on the homogeneous carbon nanotubes, ordering of the mesogens considerably depends on n and the characteristic odd−even effect occurs. The monolayers composed of the mesogens with an odd value of n exhibit a higher degree of order than the mesogens with even values of n. Dipolar relaxation exhibits a different characteristic of dependency of thermal activation energy on the value of n. Mesogens with an odd value of n exhibit lower activation energies than their counterparts with even values of n.



Article

AUTHOR INFORMATION

Corresponding Author

*E-mail: [email protected]. Notes

The authors declare no competing financial interest. 19270

DOI: 10.1021/acs.jpcc.5b05961 J. Phys. Chem. C 2015, 119, 19266−19271

Article

The Journal of Physical Chemistry C (21) Lisetski, L. N.; Lebovka, N. I.; Naydenov, S. V.; Soskin, M. S. Dispersions of Multi-Walled Carbon Nanotubes in Liquid Crystals: A Physical Picture of Aggregation. J. Mol. Liq. 2011, 164, 143−147. (22) Całus, S.; Kityk, A. V.; Eich, M.; Huber, P. Inhomogeneous Relaxation Dynamics and Phase Behaviour of a Liquid Crystal Confined in a Nanoporous Solid. Soft Matter 2015, 11, 3176−3187. (23) Paluch, M.; Dendzik, Z.; Rzoska, S. J. Scaling of High-Pressure Viscosity Data in Low-Molecular-Weight Glass-Forming Liquids. Phys. Rev. B: Condens. Matter Mater. Phys. 1999, 60, 2979−2982. (24) Dendzik, Z.; Paluch, M.; Gburski, Z.; Ziolo, J. On the Universal Scaling of the Dielectric Relaxation in Dense Media. J. Phys.: Condens. Matter 1997, 9, L339−L346. (25) Schneider, U.; Brand, R.; Lunkenheimer, P.; Loidl, A. Scaling of Broadband Dielectric Data of Glass-Forming Liquids and Plastic Crystals. Eur. Phys. J. E: Soft Matter Biol. Phys. 2000, 2, 67−73. (26) Darden, T.; York, D.; Pedersen, L. Particle Mesh Ewald: An N log(N) Method for Ewald Sums in Large Systems. J. Chem. Phys. 1993, 98, 10089−10092. (27) Brunger, A.; Brooks, C. L., III; Karplus, M. Stochastic Boundary Conditions for Molecular Dynamics Simulations of ST2 Water. Chem. Phys. Lett. 1984, 105, 495−500. (28) Phillips, J. C.; Braun, R.; Wang, W.; Gumbart, J.; Tajkhorshid, E.; Villa, E.; Chipot, C.; Skeel, R. D.; Kalé, L.; Schulten, K. Scalable Molecular Dynamics with NAMD. J. Comput. Chem. 2005, 26, 1781− 1802. (29) Górny, K.; Dendzik, Z.; Pabiszczak, M.; Gburski, Z. Non-Debye Dipolar Relaxation of Ethylene Glycol Embedded in ZSM-5 Zeolite Host Matrix - Computer Simulation Study. J. Non-Cryst. Solids 2013, 364, 28−33.

19271

DOI: 10.1021/acs.jpcc.5b05961 J. Phys. Chem. C 2015, 119, 19266−19271