NMR OF SOLUTES IN NEMATIC LIQUID CRYSTALS: AN INVESTIGATION OF T H E MECHANISMS OF ORIENTATIONAL ORDERING by ARTHUR J. V A N DER EST B.Sc. (U.B.C.) A THESIS SUBMITTED IN PARTIAL FULFILMENT OF THE REQUIREMENTS FOR T H E DEGREE OF DOCTOR OF PHILOSOPHY in THE FACULTY OF GRADUATE STUDIES Chemistry We accept this thesis as conforming to the required standard THE UNIVERSITY OF BRITISH COLUMBIA June 1987 ® Arthur J . van der Est, 1987 In presenting this thesis in partial fulfilment of the requirements for an advanced degree at the University of British Columbia, I agree that the Library shall make it freely available for reference and study. I further agree that permission for extensive copying of this thesis for scholarly purposes may be granted by the head of my department or by his or her representatives. It is understood that copying or publication of this thesis for financial gain shall not be allowed without my written permission. Department of (2Ht=:ns<.l s T The University of British Columbia 1956 Main Mall Vancouver, Canada V6T 1Y3 Date (O iT-lsl- / ABSTRACT Dipolar and quadrupolar couplings measured from NMR experiments, and order parameters obtained from these couplings, are reported for a number of small solutes dissolved in several nematic liquid crystals. These results are discussed in terms of the solute-solvent between the interactions. It has been shown that the interaction solute molecular quadrupole moment and a mean external electric field gradient due to the liquid crystal accounts ordering of molecular hydrogen. The for most but not all of the remaining contribution to ordering is discussed in terms of possible mechanisms. The anisotropic couplings observed for methane and acetylene account the are interaction Excellent agreement obtained. Evidence is discussed between between the given that in terms the of a model vibrations observed and rotations and calculated these two which takes into of the solute. dipolar couplings solutes experience the same is field gradient as molecular hydrogen. In a mixture of 55wt per cent 1132 (Merck ZLI 1132) and 45wt per cent EBBA (N-(4-ethoxybenzylidene)-4'-n-butylaniline 301.4K the deuterons in T>2 experience at no external electric field gradient. The order parameters of a series of solutes in this mixture are calculated in terms of a simple model for the short range interactions which depend on the size and shape of the agreement gradient with is not solute. the zero These calculated experimental it is shown results. that order In the parameters liquid are crystals combination of in very where the good the field short range interaction model and the electric field gradient - molecular quadrupole moment mechanism predicts the order parameters very well. ii T A B L E OF CONTENTS Abstract ii Table of Contents iii List of Figures v List of Tables vii List of Abbreviations ix Acknowledgements x I. Introduction A. Liquid Crystals B. NMR of Partially Oriented Molecules 1. The Order Matrix 2. The Nuclear Spin Hamiltonian a. The Zeeman Interaction b. Indirect or J Coupling c. Direct Dipole-Dipole Coupling d. Quadrupolar Coupling 3. Spectral Analysis 4. NMR Using Liquid Crystalline Solvents C. Outline of Thesis 1 3 6 6 8 8 9 10 11 12 14 14 II. Experimental A. Preparation of liquid crystals 1. General 2. Preparation of EBBA-d2 a. Syntheses 3. Preparation of 1132/EBBA Mixtures B. Preparation of solutes 1. Hydrogens 2. Acetylene / Acetylene-di / Acetylene-d2 C. Preparation of NMR samples 1. Hydrogens 2. Methanes 3. Acetylenes 4. Other Solutes D. NMR spectra 19 20 20 20 20 22 23 23 23 24 24 24 25 25 26 ITI. Results and Discussion A. Hydrogen 1. Spectra of Partially Ordered H 2 , HD and D2 2. Dipolar and Quadrupolar Coupling Constants 3. External Electric Field Gradients 28 28 28 30 33 iii 4. Electric Field Gradient Molecular Quadrupole Moment Interaction 36 5. D2 in 1132/EBBA Mixtures 37 6. Isotope and Temperature Effects 41 a. 60 wt per cent 1132/EBBA 41 b. 50 and 55wt per cent 1132/EBBA 46 7. Summary 49 B. Methane 52 1. Introduction 52 2. Theory 53 a. Solvent-Solute Interaction Potential 53 b. Calculation of Observables 54 3. Methanes in Liquid Crystal Mixtures 59 4. Summary 63 C. Acetylene 64 1. Introduction 64 2. Experimental Results 66 3. Theoretical Considerations 69 4. Calculation of Order Parameters 70 5. Summary 79 D. Size and Shape Effects 80 1. Introduction 80 2. Short Range Interaction Model 83 a. The Interaction Potential 83 b. The Short Range Potential 84 c. Calculation of Order Parameters 88 3. Ordering of solutes in 55wt per cent 1132/EBBA 89 a. Solutes with C3 or Higher Symmetry 89 b. Solutes of Lower Symmetry 96 c. Liquid Crystal Spectrum as an Internal Standard .... I l l d. Summary 120 4. Ordering of solutes in other liquid crystals 120 a. 1132 and EBBA 121 b. 5CB-a,0-d4 1 IV. Conclusions 2 6 137 V. Appendix I ... 139 VI. Appendix U A. Calculation of Maximum and Minimum Circumferences 1. Projection of solute 2. Maximum circumference 3. Minimum circumference 145 145 145 146 146 VTI. Appendix HI A. Source Code for the Calculation of Expression (51) Vffl. Bibliography 148 148 159 iv List of Figures FIGURE 1: Representation FIGURE 2: l H Spectrum 1132/EBBA of a Nematic Phase of 2,4-Hexadiyne 4 Dissolved in 55wt per cent 15 FIGURE 3: Synthesis of EBBA-d2 21 FIGURE 4: l H Spectrum of H2 and HD in a Nematic Phase 29 FIGURE 5A: 2R Spectrum of D2 and HD in a Nematic Phase 30 FIGURE 5B: 2H Spectrum of D2 in a Nematic Phase 31 FIGURE 7: Spectra of D2 in 1132/EBBA Mixtures 38 FIGURE 8: Ordering of Molecular Hydrogen in Nematic Liquid Crystals. I 40 FIGURE 9: Ordering of Molecular Hydrogen in Nematic Liquid Crystals. H .... 42 FIGURE 10: Dipolar Couplings 1132/EBBA FIGURE 11: Quadrupolar 1132/EBBA FIGURE 12: Dipolar Couplings 1132/EBBA for H2, FIGURE 13: Dipolar Couplings 1132/EBBA for H2, FIGURE for Couplings H2, HD and D2 in 60wt per cent 44 for HD and D2 in 60wt per cent 45 HD and D2 in 50wt per cent 47 HD and D2 in 55wt per cent 48 14: Elastic Tube Model for the Short Range Interactions 85 FIGURE 15: Projection of Solute in XY Plane 86 FIGURE 16A: Short Range interaction Model (Minimum Circumference) 92 FIGURE 16B: Short Range interaction Model (Maximum Circumference) 93 FIGURE 17: Polarizability - Electric Field Mechanism 98 FIGURE 18: FIGURE 19: Short Range Interaction Model; Experimental vs Order Parameters in 55wt per cent 1132/EBBA-d2 Polarizability 1132/EBBA-d2 - Electric Field v Interaction in Calculated 103 55wt per cent 105 FIGURE FIGURE 20A: Asymmetry Parameters for Solutes of Symmetry : Short Range Interaction Model C2v and D2h 109 20B: Asymmetry Parameters for Solutes of C2v Symmetry Polarizability - Electric Field Interaction and D2h 110 FIGURE 21: Structure and Molecule Fixed Axes; EBBA-d2 H2 FIGURE 22: Force Constant in 55wt per cent 1132/EBBA-d2 117 FIGURE 23: Electric Field Gradients in 55wt per cent 1132/EBBA-d2 118 FIGURE 24: Solutes in EBBA-d2 Short Range Interaction Model 124 FIGURE 25: Solutes in 1132 Short Range Interaction Model 125 FIGURE 26: Solutes in 5CB-a,/3-d4: Short Range Interaction Model 129 FIGURE 27: Field Gadients in 5CB-a,/3-d4 131 FIGURE 28: Force Constant in 5CB-a,/3-d4 132 FIGURE 29: FIGURE Solutes in 5CB-a,/3-d4: Short 1,3,5-trichlorobenzene as Reference 30: 2,4-hexadiyne in 5CB-a,j3-d4 vi Range Interaction Model: 134 135 List of Tables TABLE 1: Electric Field Gradients in 1132 and EBBA-d2 37 TABLE 2: D2 in Mixtures of 1132 and EBBA 39 TABLE 3: H2, HD and D2 Dipolar Couplings in Various Liquid Crystals 50 TABLE 4: Methane in 1132, EBBA and 61wt per cent 1132/EBBA at 31 OK: Experimental and Calculated Dipolar Couplings and G/3 Parameters 60 TABLE 5A: Methane in 1132, EBBA and 61wt per cent 1132/EBBA at 31 OK: Experimental and Calculated Quadrupolar Couplings (Hz) ...62 TABLE 5B: Parameters q and 9q/9SF Gradient Corrections TABLE TABLE TABLE eq 2 With and Without External Field 63 6: Acetylene in 1132, EBBA and 55wt per cent 3OIK: Experimental Dipolar Couplings (Hz) 1132/EBBA at 67 7: Acetylene in 1132, EBBA and 55wt per cent 1132/EBBA at 30 IK: Calculated and Experimental Dipolar Couplings (Hz) 71 8: Acetylene in 1132, EBBA and 55wt per cent 1132/EBBA at 30IK: Calculated and Experimental Dipolar Couplings and the Anisotropy in the Indirect C-C Coupling (Hz) 76 TABLE 9: Spectral Parameters for Solutes with C3 or Higher Symmetry TABLE 90 10: Experimental and Calculated Order Parameters for Solutes with C3 or Higher Symmetry Dissolved in 55wt per cent 1132/EBBA at 301.4K 94 TABLE 11: Solute Polarizabilities and References 97 TABLE 12: Experimental an Calculated Solute Order Parameters per cent 1132 at 301.4K in 55wt 13: Deuteron TABLE TABLE TABLE 14: 15A: Force Constants, Electric Fields, Liquid Crystal Splittings and Order Matrix Elements(55wt per cent 1132/EBBA-d2 at 301.4K) 99 107 EBBA-d2 AT 301.4K: Liquid Crystal Order Parameters, Deuteron Quadrupolar Splittings and Dipolar Couplings 114 Furan and Thiophene in 55wt per Experimental Solute Order Parameters 115 vii cent 1132/EBBA-d2; TABLE 15B: Furan and Thiophene in 55wt per cent 1132/EBBA-d2:Liquid Crystal Order Parameters, Deuteron Quadrupolar Splittings and Dipolar Coupling Constants 116 TABLE 16: D2 in 55wt per cent 1132/EBBA-d2 119 TABLE 17: Molecular Quadrupole Moments 122 TABLE 18: Solute Order Parameters: 1132 and EBBA-d2 at 301.4K 123 TABLE 19: D2 in 5CB 127 TABLE 20: Order Parameters of Solutes in 5CB at 294K 127 viii LIST OF ABBREVIATIONS NMR Nuclear Magnetic Resonance S order matrix S^j kl element of order matrix J jj indirect spin-spin coupling between spins a and b AJ anisotropy in J a a b & b direct dipole-dipole coupling between nuclei a and b B quadrupolar coupling of nucleus a fl QJJ nuclear quadrupole moment of deuteron q internal electric field gradient tensor at the site of the deuteron EBBA N-(4-ethoxybenzylidene)-4'-n-butylaniline 1132 Merck ZLI 1132: A mixture of three phenylcyclohexanes and one biphenylcyclohexane (see chapter III) 5CB-a,/3-d4 4-(a,/3-d4-rc-pentyl)-4'-cyanobiphenyl. TTF tetrathiofulvalene U solute-solvent interaction potential F^2 external electric field gradient Q zz component of the molecular quadrupole moment tensor 2 Z 7 , H 7 D gyromagnetic ratio of proton (deuteron) G&kl parameter describing so vent-solute interaction strength Q mth normal mode of vibration of solute k M elastic force constant anisotropy in mean square electric field ix ACKNOWLEDGEMENTS I wish to thank Dr. E . E . Burnell for his guidance attitude and the many hours of discussions we and support. His positive have had, have been a tremendous help to me. It has been an especially great pleasure to have worked with J.C.T. Rendell and E.J. Delikatny. I am grateful for their help and friendship. During the course of this work I was fortunate enough to spend a year at the Vrije Universiteit in Amsterdam and I am indebted to Dr. C.A. de Lange and Dr. J.G. Snijders for their hospitality and support during my stay. I especially wish to thank Drs. J.B.S. Barnhoorn, with whom I worked quite closely, for making it a most rewarding experience. My thanks are also extended to Dr. P.B. Barker, Dr. G.S. Bates, Hoatson, M.Y. Kok, Dr. G.N. Patey, Dr. G.L. and A. Weaver for their collaboration in various aspects of this work. I also wish to thank Dr. F.G Herring and Dr. D.C. Walker for their helpful suggestions I must reserve my deepest gratitude for my wife Karin, whose patience and understanding I am very thankful for. I am also very appreciative of her help in typing this thesis and preparing some of the figures. x I. INTRODUCTION When a molecule is placed in an anisotropic environment its reorientation may also become anisotropic. That is, not all orientations of the molecule relative to an external frame of reference are equally probable. This is a result of the fact that most molecular properties are anisotropic so that their interaction with an anisotropic environment anisotropic fluids and is orientation as solutes. The measurement solvents dependent. they Nematic provide liquid anisotropic crystals are environments for of any property which depends on the orientation of the solute relative to a space fixed axes system can provide information about the solvent-solute interactions in such systems [1-3]. NMR is ideally suited for this purpose because quantities such as the dipolar and quadrupolar couplings and the anisotropy in the chemical shifts, solute and which average orientation dependent part which depend on the orientation of the to zero in isotropic solvents, become observable. The of the dipolar and quadrupolar couplings is often written as a 3X3 matrix called the order matrix [1] and the individual elements of the matrix are called the order parameters of the solute. The order parameters of a solute are a measure of its average orientation and its interaction with the liquid crystalline range of nature liquid of the environment. crystalline Despite solvents intermolecular forces the large number of which have been studied which lead to solutes and [2-5], the solute ordering is wide exact still not understood. Past attempts at understanding the ordering of solutes have generally centred around correlating the solute order parameters with molecular properties such shape as [6, the polarizability [6-10], 13]. In a number of moments other of studies inertia the [11, solute described in terms of an unspecified interaction potential [10, 1 12] and molecular ordering has been 14, 15, 16-20]. The Introduction / 2 success of such studies has, however, been limited and special mechanisms are invoked to explain the ordering of solutes which do not correlate well with the chosen molecular property [21]. In this thesis we will reexamine this problem using a slightly different approach [22]. We will begin by studying in detail the NMR spectra of isotopically substituted modifications of some very simple solutes such as molecular hydrogen and methane. Because of the simplicity of these solutes and because a great deal is known about their physical properties, it should be possible to obtain specific information about the mechanisms the liquid crystals used. which lead to their orientational ordering in We will then use the knowledge gained from these studies in order to vary the liquid crystal environment in a controlled fashion. This should allow a more systematic study of the ordering of larger solutes whose molecular properties are not as well known. For these studies we will also use partially deuterated liquid crystals which will allow us to examine the relationship between solute and solvent ordering of the crystal molecules solvent is ordering. A complete complicated by the description of the internal motions of the liquid [23] and we shall restrict ourselves to using the liquid crystal spectrum as a measure of the anisotropic environment experienced by the solutes. As will be seen this approach is quite successful and we will be able to make some useful conclusions concerning the relative importance of various mechanisms to the orientational ordering of a number of solutes in several liquid crystalline solvents. We begin by presenting a brief description of liquid crystalline solvents and the NMR spectra of solutes in such solvents. Introduction / 3 A. LIQUID CRYSTALS The first reported observation of liquid crytalline behaviour was Reinitzer [24] who found that cholesteryl benzoate melted at 145.5 an opaque fluid which cleared suddenly at 178.5 made by ° C to form ° C . He brought this compound, as well as cholesteryl acetate, to the attention of Lehman [25] who studied the opaque fluid under the polarizing microscope and discovered that it had the optical properties of an uniaxial crystal. He proposed the term "fliissige Kristalle" (liquid crystal) to describe this fluid since it displayed properties of both liquids and crystals. Many compounds have since been discovered which show similar behaviour and it is now known that a wide range of intermediate phases exist between the solid systematically crystal has and liquid. by Friedel [26] come These various phases were first described and despite his strong objection, the term liquid to mean those states of matter whose degree of molecular ordering lies between that of a crystalline solid and an isotropic liquid [27]. In this thesis we will be concerned with the least ordered of these phases. Friedel proposed the name nematic from the Greek word vqua (thread) because of the thread-like discontinuities observed under the polarizing microscope for this phase. The molecules which form this phase are, in general, long and rod-like and are arranged such that,as in an isotropic liquid, there is no long range positional ordering of their centres of mass. Unlike an isotropic liquid, however, all molecules in the phase have the same average orientation with respect to a unique direction approximately referred parallel to to the as long the director. axes of the This direction molecules. A usually "snap lies shot" representation of the nematic phase is shown in figure 1. In the absence of any Introduction / 4 FIGURE 1 Representation of a Nematic Phase A "snap shot" representation of a nematic phase. The long rod-like molecules, represented by the rods in the figure, are aligned so that their long axes are on average parallel to the director, Z. There is no positional ordering of the centres of mass of the molecules. Introduction / 5 external fields the director varies randomly throughout the sample, whereas, in the presence of a strong magnetic field, such as is used in NMR experiments, it normally aligns parallel to the field direction. Alignment of the director may also occur in the presence of an electric field and at the surface of prepared glass plates. A second more ordered liquid crystalline phase, was named smectic by Friedel from the Greek word oprjyqa (soap) because it was first discovered in potassium and ammonium oleate. In this phase the molecules are aligned parallel to the director as in the nematic phase but they also have some degree of positional ordering of their centres of mass. Several modifications of this phase exist and they are labeled smectic A, smectic B, etc. depending on the particular arrangement of the molecules. Both the nematic and smectic phases are examples of thermotropic liquid primarily crystals, so named temperature induced. There formed by amphiphilic induced by the presence compounds exists and of water. because a whose Such the liquid second class liquid crystallinity is of liquid crystals crystalline behaviour systems are usually referred to is as lyotropic liquid crystals and they have attracted a great deal of attention because of their similarity to lipid membranes. Because of its simplicity, the nematic phase is much more attractive as a system for studying intermolecular interactions than either the smectic phase or lytropic liquid crystals. The ordering of the molecules in the nematic phase is completely described by the time or ensemble average of the orientation of a molecule fixed axis system with respect to a single space fixed axis. This is also Introduction / 6 true of solutes dissolved in this phase. We will now describe the effect of such ordering on the NMR spectrum of a molecule and discuss how the parameters obtained from the NMR spectrum can be related to the ordering of the molecule. B. N M R OF PARTIALLY ORIENTED M O L E C U L E S 1. The Order Matrix Before discussing the NMR spectrum of partially oriented molecules we will first introduce the order matrix [1] because it will prove useful in much of what follows. Many of the properties of the molecules that we will be concerned with are of second rank tensorial form. Because we can only measure the ensemble or time average of these properties in the laboratory fixed frame of reference, the transformation of the tensor from a molecule fixed axis sj'stem to a space fixed axis system, averaged over the motions of the molecule, is of interest. For a second rank tensor, T, this is given by: <T > = E <cos0 , cos0 > l U where < T J J > k.,1 fc u k l , (1) are the components of T averaged over the molecular motions in the space fixed axis system and are the components in the molecule fixed frame. (We will use upper case labels for space fixed axes and lower case labels for and the molecule fixed axes throughout). is the angle between axes k and I we have assumed that the components of the tensor in the molecule fixed axes are independent of the internal motions of the molecule. For a molecule in a nematic liquid crystal, the orientations 0^1 the individual a n d ~&kl a r values e e ( i u a u y probable but of Introduction / 7 dfri are not. Because of this < c o s 0 £ j > is zero but <cos 0^j> is not. In most cases, when such a system is placed in a strong magnetic field the director is oriented parallel to the field direction. Because of this, and the fact that there is no positional ordering of the molecules, the lab fixed frame axially symmetric. This leads to the result = < TYY : > < T I J > of reference is = 0 for I / J and < T X X > if Z is the symmetry axis of the lab frame. The trace of the tensor must be invariant to rotations so that <T'xx> = < YY> T Thus we need only specify <r ^2,z <T > = -3 E T ^ + ZZ where = ( f T « - <Tzz>)/2 • > a n c ^ ^ e I £<f c o s 0 6^/ is the Kronecker t r a c e < a ( 2 ) Expression (1) then becomes cos0 lz -J 6 >T kl kl , (3) 6. (1/3)Z T^- is often referred to as the isotropic i part of T because it is the value observed if the molecular motion is isotropic. It is convenient to collect the angular dependent part of expressions 3X3 (3) as a matrix S such that S kl ~ 1 < 3 COS0A2COS0H ~ 6 > kl . (4) This matrix is called the order matrix. It is real, symmetric and traceless and has the transformation properties of a second rank tensor under rotations of the molecule fixed axis system. In general it has five independent elements. However, through an appropriate transformation of axes it can be diagonalized leaving only two independent non-zero elements. If the molecule possesses a three fold or higher rotation axis then in an axis system with z parallel to the rotation axis Introduction / 8 S^j = 0, k / element, S , 1 and S = xx S = ~2$ - Thus S has only one independent zz which is referred to as the order parameter of such a molecule. In zz this case S is the average zz of P2(cos# 2)> the second Legendre polynomial of 2 cosfl . 2Z 2. The Nuclear Spin Hamiltonian The Hamiltonian for a spin system in a large external magnetic field consists of contributions from (i) the Zeeman interaction of the nuclear spins with the magnetic field, (ii) the indirect or J coupling between the nuclear spins via the electron spins, moments and (iii) (iv) the direct for spins dipole-dipole greater than coupling 1/2, the of the interaction nuclear of the magnetic nuclear quadrupole moment of each spin with the electric field gradient at the site of the nucleus. The total Hamiltonian may be written as a sum of contributions from each of these terms. K =X Z + #J + #D + #Q ( 5 ) c. The Zeeman Interaction The Zeeman Hamiltonian has the form *Z = where v & "f Ij. . ( 6 ) is the resonance frequency of nucleus a and l^a is the operator for the component of the spin of nucleus a along the field direction. The resonance frequency depends of the local field at the site of the nucleus and is given by Introduction / 9 7a "a = ~ where y n (7) (1 " ^ZZ,a) » 27T is the gyromagnetic ratio of nucleus a, h is Planck's constant, & the applied magnetic field and O"^ of . Z H nucleus a expressions and is is is the ZZ component of the shielding tensor a averaged over the motions of the molecules. From (3) and (4) we have _ °ZZ,s 1 2 ~ 3f ii* a 3 ? , kl°kl,a kl + In an isotropic solution the second term vanishes and (8) • S t <7 = Z Z a \ ^ a iia i b. Indirect or J Coupling The contribution to the Hamiltonian from the indirect spin-spin coupling between spins a and b is given by: #j = S E I • J a a b>a where J a b is the second rank • I , a b b w tensor describing the strength coupling between spins a and b. Again using expressions = E2 5 « . a b I. • I a b>a t #j I 2 a b>a The second (1/3)Z JJ-J- a b term is is the w > a b kl again S zero kl (3 l for Z a an scalar coupling observed l Z b indirect (2-4) this becomes + E J ? J 2 of the b - I. • I ) . (10) b isotropically tumbling molecule and in a normal high resolution N M R i experiment. expression For most ( 1 0 ) is couplings small and involving may be protons ignored the anisotropic [ 5 ] . This has part of important Introduction / 10 consequences because, as we shall see coupling has exactly the same in a moment, the direct dipole-dipole form as this term and thus the two coupling constants cannot be distinguished in the spectrum. c. Direct Dipole-Dipole Coupling The contribution to the Hamiltonian from the dipolar coupling between the spins in a molecule is given by = EE M D D a b l (3 Z a I - I . Z b a ] ), b (ii) a b>a where _ h7 7b, 3 cos 0 2 a Dab = ~ ( — 2 - ) < a b z - 1 3 2 r 4 7 r > • ( 1 2 > a b 6?abZ is the angle between the laboratory Z axis and the vector joining the two nuclei a and b. If the averaging over the vibrations and rotations can be done separately then expression (12) becomes n ^ab where S a b _ _ _ ( V h . 2 47T 7 a 7 b \ < r " 3 > (13) / ^ ab^ ab » r 0S w is the order parameter of the vector between nuclei a and b and may be related to the elements of the order matrix, S of the molecule fixed axis system. This separation of the vibrational and rotational averages is by no means always valid and we shall see evidence of this in the spectrum of methane. For most solutes the effect of separating the two averages is quite small and the dipolar couplings can be used to determine the order matrix. Moreover, if one distance in the molecule is fixed, then the average structure Introduction / 11 may be determined. It is also apparent that for an isotropically tumbling molecule the dipolar couplings average to zero. d. Quadrupolar Coupling In anisotropic environments there is also a contribution to the Hamiltonian from the coupling between the quadrupole moment of nuclei with spin > electric field gradient at the site of the nucleus. 1/2 and the The Hamiltonian for this interaction has the form I. ' B • 1. #Q = £ a . (14 ) The tensor B , which describes the interaction strength of the nuclear quadrupole a moment of nucleus a with the electric field gradient, is symmetric and traceless so that from expressions (3) and (4) we have *Q = 5 J ZZ.a (3I a " I.*) . B Z where a is the ZZ component of tensor B & (15) and is given by 4 h I [ 2 I - l ] kl a Q a a is the quadrupole moment of nucleus a, e is the charge of the electron, I a is the spin of nucleus a and q^j are the elements of the electric field gradient tensor at the site of nucleus a, averaged over the internal motions of the molecule. We have again made the assumtion that the averages over vibrations and rotations are separable. B ^ a is, of course, zero for isotropic motion and its value may be used to determine the elements of the order matrix for anisotropic Introduction / 12 motion if Q and & are known. Unfortunately, the electric field gradient tensor is very difficult to determine accurately. In this thesis with I > 1/2 which will be studied. is the only nucleus For many compounds the electric field gradient tensor at the site of the deuteron may be taken as axially symmetric which allows a considerable simplification of expression (16). Using values of the field gradient tensor matrix at the from single site of the crystal or molecular beam studies, deuteron may the order then be determined fairly accurately from the spectrum. These, then, are the various contributions to the total Hamiltonian which is given by a + I 8 z b>a g,. (».b + -31 J « , b s „ ) (3 + i ? B . (31* - iz. i Z b - i. • i ) b I. ) . 2 3. Spectral Analysis The analysis of the spectrum of a given spin system inolves calculation of the matrix elements of H with the various spin states of the system. The transition frequencies may then be calculated from the eigenvalues of the Hamiltonian matrix and are functions of the various coupling constants and chemical shifts. These values experimental are spectra adjusted is until obtained. agreement This procedure between is the referred calculated to as and fitting a Introduction / 13 spectrum and is usually done with the aid of a computer. For many simple highly symmetric molecules the analysis of the spectrum is trivial and the coupling constants may be obtained directly from the line positions. In this thesis all spectra, whose analysis is not trivial, have been fit using the computer program LEQUOR [28]. The elements of the order matrix in an appropriate molecule fixed axis system may be calculated from the values of the dipolar and quadrupolar couplings. If the dipolar couplings are used, the average geometry of the molecule must be known. However, in general, the average geometry of a molecule, obtained by electron diffraction or microwave spectroscopy, is not consistent with the observed dipolar couplings. This conformational is averaging a of result of either the solute in changes the liquid in the vibrational and crystalline phase or contributions to the dipolar couplings due to the interaction between the vibrations and rotations of the molecule. A proper treatment of the vibrational problem requires some knowledge of the solvent-solute interaction mechanism and, except for the simplest of solutes, is rather complicated [14-16, 29]. Because of this, it is usually assumed that the averages over the vibrations and rotations may be done separately and a least squares fit of the average geometry to the dipolar couplings is done. The order matrix for this fitted geometry is then obtained. The computer program SHAPE [30] has been developed for this purpose and will be used where necessary to obtain the order matrix of the solutes studied here. Introduction / 14 4. NMR Using Liquid Crystalline Solvents In this thesis we will use the procedure just described to obtain the dipolar and quadrupolar coupling constants of solutes dissolved in a variety of nematic liquid crystals and we will attempt to explain these results obtained from them) in terms of the solute-solvent contributions indirect to the measured coupling and we chemical shifts. will be concerned order parameters interactions. We will ignore dipolar coupling due not (or to the anisotropy in the with the anisotropy in Figure 2 shows the spectrum of 2,4-hexadiyne the in the nematic liquid crystal mixture 55wt per cent 1132/EBBA-d2- The sharp lines are due to the solute and the broad featureless peaks in the baseline are due to the nematic solvent. This difference in linewidths arises because the enormous number of lines in the spectrum of the liquid crystal cannot be resolved and because of the slower tumbling of the solvent molecules due to their large size. This is convenient because it eliminates the need to remove the solvent spectrum as is done when using isotropic solvents. C. OUTLINE OF THESIS Our goal is to interpret the dipolar and quadrupolar couplings and order parameters obtained from spectra such as shown in figure 2, in terms of the intermolecular interactions. To do this, it is important to choose the system carefully because of the problems associated with the interpretation of both the dipolar and quadrupolar couplings. The dipolar couplings depend on the average over all of the motions of the molecule so that we need to study systems in which the internal molecular motions may be either safely ignored or described Introduction / 15 FIGURE 2 l H Spectrum of 2,4-Hexadiyne Dissolved in 55wt per cent 1132/EBBA i 1 -6.0 1 1 -2.5 1 1 1 OJO 1 ZJ6 1 1 1 SJO Frequency (kHz) The narrow peaks are due to the solute and the broad peaks in the baseline are due to the liquid crystalline solvent. Solute concentration = ~1 mole per cent, T = 301.4K. very simply. Quadrupolar couplings can only be observed in molecules which contain deuterons or some other quadrupolar nucleus. The interpretation of the couplings requires a knowledge of the nuclear quadrupole moment and the electric field gradient tensor at the site of the nucleus. These requirements can only be met for the simplest of molecules. Thus we are restricted, initially at least, to studying very simple solutes. We will begin by presenting results obtained for molecular hydrogen and its deuterated analogues [10, 16, 31, 32], For this solute both the dipolar and quadrupolar couplings can be related quite accurately to the Introduction / 16 molecular ordering. Moreover, the electronic properties of this solute are very well known so that it should mechanisms involving experiences an be these external possible properties. electric field to It evaluate has gradient the relative importance of been shown that due to liquid its this solute crystalline environment and that the interaction between the molecular quadrupole moment and the anisotropy in this field gradient accounts for most of the solute ordering [31, 32]. Quantum mechanical effects must be taken into account in order to explain the dipolar couplings of H^, HD and D2 [10, 16]. By mixing two liquid crystals where the measured electric field gradients have opposite sign it will be shown that a system intermediate between may be produced where those in the two the field gradient has component liquid crystals values (see also reference [32]). In these liquid crystal mixtures some very unusual isotope effects are observed in the ordering of molecular hydrogen. These results will be discussed in terms of some posible contributions to the mechanisms which lead to orientational ordering. Because of the small number of couplings which can be measured in hydrogen and because very good of its extremely system for the simple structure, this study of the influence solute does not provide a of internal motion on the observed spectrum. Methane and its deuterated analogues this purpose [33]. Because these molecules are ideally suited for are tetrahedral in their equilibrium geometry, any observed spectrum can only arise because of the coupling between the molecular vibrations and rotations and, in the case of the partially deuterated species, deviations of their average geometries from tetrahedral symmetry. These two effects may be calculated [14-16, 34] in terms of the well known vibrational Introduction / 17 force field [35] and the strength of the solvent solute interaction. If the latter interaction is taken to be of second rank tensorial couplings [33] are predicted very accurately [14, form the observed dipolar 15, 36]. It will be shown that this solute experiences the same external electric field gradient as hydrogen, but that the interaction with the molecular quadrupole moment does not account for all of the observed dipolar couplings (see also reference [36]). Methane and hydrogen are unusual solutes because of their small size and, in the case of methane, very high symmetry. Thus we will consider a solute where vibrational effects are expected to play a smaller but still important role. For this purpose we will choose acetylene. We will see that for this solute the electric field gradient - molecular quadrupole moment mechanism plays an important role and that there is a second mechanism present which is roughly proportional to the size and shape of the solute. The results for hydrogen, methane of other solutes [37]. and acetylene provide a basis for the study For all three of these solutes the electric field gradient - molecular quadropole moment mechanism plays an important role, although it is not the only mechanism. From the results on hydrogen it has been shown [32] that a liquid crystal mixture where the field gradient is zero may be produced. This mixture, then, provides a system for studying other mechanisms [38, 39]. For most solutes vibrational effects are small and will be ignored. Because the short range interactions are expected to play an important role for larger solutes, we will develop a simple model for this mechanism and it will be shown that this model can be used to successfully predict the order parameters of a large number of solutes in the zero field gradient mixture (see also reference [38]). Introduction / 18 Finally, we will see that the short range interaction model combined with the electric field gradient - molecular quadrupole moment mechanism predicts the ordering of solutes quite well in liquid crystals where the electric field gradient is non-zero [38-40]. The results discussed above and the relevant theory will be presented following chapter II which deals with the experimental details. Chapter III is divided up into sections dealing with hydrogen, methane and acetylene respectively, as a discussion of the short range interaction model. as well II. E X P E R I M E N T A L The liquid crystals used were; EBBA: N-(4-ethoxybenzylidene)-4'-M-butylaniline, EBBA-d2: N-(4-ethoxybenzylidene)-2,6-dideutero-4'-n-butylaniline. D. 1132: An eutectic mixture of frans-4-n-alkyl-(4-cyanophenyl)-cyclohexane (alkyl propyl, pentyl and heptyl) and frans-4-n-pentyl-(4'-cyanobiphenyl-4)-cyclohexane [5]. 3y C H 3 ( C H ) 2 4 R = C C H 2 C - i C H 2 C D 2 C D (24 per cent) g o n C_H, / 15 , 3 H C-H.,., 5CB-a,/3-d4: 4-(a,p-d4-n-pentyl)-4'-cyanobiphenyl. C H 4 2 19 (36 per cent) (25 per cent) = Experimental / 20 A. P R E P A R A T I O N OF LIQUID CRYSTALS 1. General All transition temperatures reported, were measured using a Gallenkamp melting point uncorrected. apparatus and are Flash evaporation of all solutions was performed using a Buchi rotary evaporator under reduced pressure. 1132 and non deuterated EBBA were obtained compounds were used without from Merck further temperatures: EBBA (solid-nematic 308 ± and K & K respectively. purification. Measured IK, nematic-isotropic 351 Both transition ± IK) 1132 (nematic-isotropic 347 ± IK). 5CB-a,/3-d^ synthesized in high purity was provided by Dr. G. S. Bates [40, 41]. 2. Preparation of E B B A - d Partially deuterated 2 EBBA was prepared according (reference [42]): to the following reaction sequence shown in figure 3. c. Syntheses N,N-2,6-tetradeutero-n-butylaniline An aqueous solution of 4-n-butylaniline hydrochloride (124 g, 0.66 moles) was prepared by adding freshly distilled 4-n-butylaniline was then removed by (100 g, flash 0.66 moles) to 2M HCI (350 evaporation. The hydrochloride ml). The water salt was then dissolved in D 2 O (100 ml) and refluxed for several hours. The resulting solution was neutralized by the dropwise partially deuterated addition of a saturated 4-n-butylaniline was then extracted solution of KOH. The with ether which was subsequently removed by flash evaporation leaving a brown oily liquid. The crude Experimental / 21 FIGURE 3 Synthesis of EBBA-02 Experimental / 22 product was then purified by vacuum distillation. EBBA-d (reference 2 [43]): A moles), 4-ethoxybenzaldehyde 4.4 mmoles) in toluene during the toluene reaction was was solution of N,N-2,6-d -n-butylaniline 4 (60 g, 0.4 moles) and 4-toluenesulfonic (60 g, 0.4 acid (0.75 g, (600 mis) was refluxed for 4 hours. The water formed removed by means then removed by flash of a Dean-Stark apparatus. The evaporation and the crude EBBA-d2 was recrystallized from methanol and dried under vacuum (yield 61 g, 55 per cent, degree of deuteration: NMR 63 per cent, mass spectroscopy 61 per cent). The measured solid-nematic and nematic-isotropic transition temperatures 351 ± IK) are in good agreement with literature values (308 [44] ± IK; for the non-deuterated compound. 3. Preparation Several of 1132/EBBA mixtures of EBBA Mixtures or EBBA-d2 and 1132 were prepared by adding appropriate weights of each of the liquid crystals to a vial or Erlenmeyer flask, heating the mixture to the isotropic phase, Vortex stirrer. and stirring it vigorously using a Experimental / 23 B. P R E P A R A T I O N OF SOLUTES 1. Hydrogens [10] High purity H electrolyzing was 2 99 purchased from Matheson Ltd. and D per cent pure D 0 of H fitted and D 2 purchased from Merck, 2 Statistical mixtures of H 2 / HD / D 2 was produced by Sharp & Dohme. were prepared by mixing varying amounts in a 100 ml round bottom flask attached to a vacuum line and 2 with a platinum wire which could be heated current through it. The wire was produce the H 2 2 / HD / D 2 then heated to ~1300K for approximately by passing a 1 minute to mixture. 2. Acetylene / Acetylene-d^ / Acetylene-d2 Solid calcium carbide (5 g, 0.1 moles) was placed in a 100 ml round bottom flask equipped with a dropping funnel and a magnetic stirrer. 4M H S 0 2 mis) in a 50 calcium per cent mixture of D 0 carbide. 2 The resulting mixture and H 0 2 was stirred was at 4 (50 added slowly to the room temperature for several hours. The gas evolved during the reaction was driven off with a stream of nitrogen, passed through a H S 0 2 4 / CuS0 4 solution, and collected in a cold trap cooled to liquid nitrogen temperature. The cold trap was then attached to a vacuum rack and the temperature was raised to approx. —10 ° C with a salt/ice bath. The acetylene evolved was then condensed in a clean, dry flask cooled to liquid nitrogen temperature. The proton NMR spectrum of the product dissolved in EBBA showed no impurities and an approximately 2:1 ratio of C H D C H . No further analysis of the purity of the product was done. 2 2 2 to Experimental / 24 C. P R E P A R A T I O N OF NMR SAMPLES For all solutes, a sufficient volume of one of the liquid crystals to fill the NMR coils was placed in either a standard or medium walled 5 mm NMR tube or a 9 mm o.d. standard pyrex glass tube. The sample tubes were then attached to a vacuum line via a piece of rubber tubing and the liquid crystal was thoroughly degassed by performing repeated freeze-pump-thaw cycles until no more gas was evolved. Solutes were then dissolved in the degassed liquid crystals as follows: 1. Hydrogens [10, 32] 9 mm o.d. sample tubes containing degassed liquid crystals were placed in a dewar containing liquid helium and a sufficient volume of hydrogen (H2, D2 or a mixture of H 2 / HD / D2) to produce a final pressure of 15 to 20 atm in the head volume of the tube at room temperature, was condensed into each tube. The tubes were then flame sealed and pressure tested in an oven at approximately 42 5K for 10-15 minutes. 2. Methanes [14, 15, 33, 36] 1 3 CH 4 (90 per cent 1 3 C), (98 per cent D), CHDg Dohme 1 3 CH D, 3 (Canada) CH D 2 2 and 13 CHgD (90 per cent 1 3 C , 98 per cent D), C H D 2 2 (98 per cent D) were obtained from Merck, Sharp and used without purification. Samples of * CH 3 4 / CD4, and CH3D were prepared in 1132, EBBA, and 61.3wt per cent 1132/EBBA, in a similar manner to the hydrogen samples except that (i) for the samples in nitrogen was 1132, 5 mm medium walled used to condense the NMR tubes gasses and (iii) the were final used, (ii) pressures liquid were Experimental / 25 between 5-15 atm. 3. Acetylenes 1 3 C H 2 2 (90 per cent 1 3 C ) and 1 3 C D 2 (98 per cent 2 1 3 C , 99 per cent D) were obtained from Merck, Sharp and Dohme (Canada) and used without purification. A mixture of C H , C H D and C D 2 carbide (see substituted 2 2 2 above). Samples acetylenes were was prepared by the hydrolysis of calcium 2 containing a prepared in mixture 1132, of all of the EBBA and isotopically 55wt per cent 1132/EBBA. All samples were prepared in 9 mm o.d. tubes in the same manner as the methane samples. The volume of gas condensed was sufficient to produce a final pressure of 5 atm in the head volume of the sample tube at room temperature. 4. Other Solutes [38-40] The 28 solutes listed in table 12 (chapter III, section D) were obtained from a variety of chemical suppliers and prepared in one of EBBA, EBBA-d 2 used 1132, without purification. Samples were 55wt per cent 1132/EBBA, 55wt per cent 1132/EBBA-d , or 5CB-a,/3-d . For the solutes which are solids or liquids 2 at room temperature 4 a sufficient amount of the solute to produce a final concentration of 1-2 mole per cent was added to the degassed liquid crystal in a standard paraffin 5 mm NMR tube. film, heated until The tubes the liquid were crystal then capped became and sealed isotropic, and stirred vigorously using a Vortex stirrer. All samples containing EBBA or EBBA-d CHgl or CHgBr with 2 and became orange over a period of several hours if allowed to stand at room temperature. In one sample (CH^Br / 55wt per cent 1132/EBBA) Experimental / 26 a considerable amount of another compound was apparent in the NMR spectrum after the sample was allowed to stand at room temperature overnight. No attempt was made to identify this compound. To avoid these problems all iodoand bromomethane samples were stored at —15 °C until their spectra were measured. For the solutes which are gasses at room temperature samples were prepared in the solutes, allene same manner as and propyne, a for the pressure methanes and acetylenes. of approximately 5 atm For two caused a noticeable drop in the nematic to isotropic phase transition temperatures in 1132, EBB A and 55wt per cent 1132 indicating a solute concentration of considerably more than 1-2 mole per cent. New samples of these two solutes were prepared using a pressure of approximately 0.25 atm. At this pressure no depression of the temperatures nematic to isotropic transition was observed and the orientational order parameter of the molecular symmetry axis of each solute was larger than observed at higher solute pressure. D. NMR S P E C T R A H, H and C free induction decays for the samples Bruker WH400 spectrometer operating at 400.1 deuterons and 100.6 were measured on a MHz for protons 61.4 MHz for MHz for carbon. For the samples containing ^ F nuclei 9 either, (i) ^ F and ^H free induction decays were measured on a Bruker CXP 9 200 spectrometer operating at 188.2 iq (ii) MHz for 1 9 F and 200.0 MHz for *H or o F and spectrometer H free induction decays operating at 282.3 MHz for were 1 9 measured on a Varian F and 46.07 MHz for 2 XL300 H . Where more than one nucleus was measured, one of the nuclei was measured through the decoupling channel ( H on the WH400 and CXP200, J 1 9 F on the XL300) or Experimental / 27 through the lock channel ( H) without removing the sample from the magnet. Thus all spectra were recorded at the same temperature. The temperature was controlled by means of a variable temperature gas flow unit and was calibrated to ±0.3K using the proton chemical shift differences in a sample of ethylene glycol. All samples were heated until the liquid crystal became isotropic and were shaken vigorously before being placed in the probe. Approximately thirty minutes was allowed for temperature equilibration before the spectra were measured. The free induction decays obtained on the Bruker WH400 and CXP200 spectrometers were transferred for processing and plotting to a Nicolet 1280 computer equipped with a Zeta 8 digital plotter. DI. RESULTS A N D DISCUSSION A. HYDROGEN In this section we present results obtained for H2, HD and D 2 in a number of nematic liquid crystals [10, 31, 32] and will begin with a discussion the spectra obtained. 1. Spectra of Partially Ordered As HD and D 2 a result of the Pauli exclusion principle [46] both H 2 and D 2 may exist in one of two forms. In H2 the symmetric combinations of nuclear spin states (total spin 1 = 1) may combine only with the antisymmetric, odd numbered rotational states. The antisymmetric combinations of nuclear spin states (1=0) may combine only with the even numbered rotational states. The first form (1=1) is referred to as ortho-H2 and the second (1 = 0) is para-O-2- The para form with 1 = 0 has no observable NMR spectrum. For D2 the orfAo form has total spin 1=2 and J even whereas para-D2 has 1=1 ortho and para forms of D 2 and J odd. In contrast to H 2 or 0 both the have an observable NMR spectrum. For HD, only one form exists because the two nuclei are different and there are no restrictions on which rotational states may combine with which nuclear spin states. Thus in a mixture of H2, HD and D2 the NMR spectra of orf/io-H , HD, ortho-T)2 2 para-D2 may be observed. The * H spectrum of such a mixture dissolved in a nematic phase is shown in figure 4. For ortho-H2 the spectrum is a doublet of splitting | 3 D | and for HD it is a 1:1:1 HH triplet of splitting |2DJJJJ + J^\- The ^H spectrum of the mixture is shown in figure 5. For HD the spectrum is a doublet of doublets with splittings |2B | and |2DgQ+Jgp|. For para-D2 it is a D 28 Results and Discussion / 29 FIGURE 4 l H Spectrum of H£ and HD in a Nematic Phase 3DHH i i -000 i i i i i i 0.0 -400 i 400 i i BOO Frequency (Hz) Proton spectrum of HD and H2 in 50wt per cent 1132/EBBA: The doublet due to H£ has a splitting |3Djj[j| and the 1:1:1 triplet due to the HD has a splitting |2DHD JHDI (see text). T=298.9K + doublet of splitting |2Bp—6Dp | and for ortho-T>2 D frequencies and intensities of the transitions six lines are observed. The for both ortho- and para-^2 are given in reference [10] and the relative positions of some of the lines are shown in figure 6. From these spectra the values of the scalar, dipolar and quadrupolar couplings may be obtained. The sign of the J coupling is known to be positive [47] and the signs of B , D i m and D D D D may also be determined [10]. Results and Discussion / 30 FIGURE 5 2H Spectrum of D2 and HD in a Nematic Phase —1— -400 —I 0.0 1— -200 I 200 400 Frequency (Hz) Deuteron spectrum of D2 and HD in 55wt per cent 1132/EBBA the doublet of doublets due to HD is labeled with *'s in the figure and has splittings |2DJJD + J|jrj| and 2Br> The remaining lines are due to D2- T=309.4K 2. D i p o l a r a n d Q u a d r u p o l a r C o u p l i n g Constants Using expressions (13) and (16) the dipolar and quadrupolar couplings may be related to the order parameter of the internuclear axis z: Dab = - ( 2~) < r b > 8 4n s zz (18) Results and Discussion / 31 FIGURE 6 2H Spectrum of D2 in a Nematic Phase Frequency (Hz) Deuteron spectrum of D2 in EBBA-d.2- B is the quadrupolar coupling, D is the dipolar coupling between the deuterons of the solute, and J is the indirect coupling between the deuterons of the solute. The centre of the spectrum, which contains no peaks, has been removed. T=301.4K ^qQc 3 B D = 4 < —7 > S " (19) n In expression (19) q=q , zz Q n is the deuteron nuclear quadrupole moment and the angled brackets imply the average over the internal molecular motions. Using these expressions the spectra of H liquid crystals parameters are [10] yield 2 interesting solvent dependent. , HD and D results. Second, the 2 dissolved in several nematic First, the longer the signs of the order average internuclear Results and Discussion / 32 distance the smaller the magnitude of S . Third, ortho- and para-D zz same order parameters to within experimental error. Finally the HD and D 2 have the 2 ratio B/D for is solvent dependent and differs from the gas phase value [49]. In order to understand these results it is necessary to describe the solute-solvent interactions which lead to solute ordering. Without specifying the exact nature of the interactions involved, Burnell, de Lange and Snijders [10, 14-16, 50] have proposed a second rank tensorial form for the interaction potential. For hydrogen it may be written 3 U = -A ( g 2 COS 6 - Zz where A is the strength of the solvent-solute 1 (20) 2 ) f interaction. This potential may then be treated as a perturbation on the freely rotating, and vibrating molecule. For the unperturbed wave functions, rotor functions products of the harmonic oscillator and rigid are used. Then by standard perturbation theory expressions for the observables such as the dipolar and quadrupolar couplings may be obtained. For hydrogen the effects of vibrations will be ignored. The order parameter is a function of the matrix elements of the direction cosines for the rigid rotor wave functions summed over all J states and a parameter describing the solvent-solute interaction and is given by [10, 14-16, 50]: S 2A E E = | r S P j <J|P (cos0 )|J> 2 2 2 P 2 Z2 - <JlP (cos0 )lJ'XJ'|P (cos6 )lJ> 2 Zz Ej - 2 Z2 (21) Ey where (2J+1) exp(-BJ(J+l)/kT) E(2J+1) exp(-BJ(J+l)/kT) ' J (22) Results and Discussion / 33 and P (cos0 ) = 2 - i . lcos\ Z2 z (23) A is the interaction parameter and J is the rotational quantum number. Using these expressions along with <r > averaged over the rotational states in the ground vibrational state, the dipolar couplings may be calculated as a function of A. The experimental couplings are predicted quite well by adjusting the value of A in a least squares fit. From this fit, it is immediately isotope dependence of the order parameters is largely apparent that the a quantum mechanical effect since in the classical limit the same order parameter is predicted for all three molecules. The variation in the sign of the order parameter between liquid crystals is explained in terms of a change in the sign of the parameter describing the solvent-solute interaction. 3. External Electric Field Gradients The liquid crystal dependence of the ratio B/D also provides some insight into the ordering of hydrogen [31]. This ratio should be a molecular property to a good liquid approximation. Thus the changes in the electronic properties crystal dependence of the molecule must be a result of or due to an external contribution to either B or D. It is unlikely that the electronic structure will change much with environment. If the separation of the vibrational and rotational averages is valid then the dipolar coupling depends only on the order parameter and <r"°> (where r is the internuclear distance). Thus the variation in B/D most likely arises because of variations in B. This can be explained very simply by proposing an external electric field gradient due to the liquid crystal Results and Discussion / 34 environment. The total electric field gradient at the site of the deuteron, is then the sum of an internal contribution due to the electrons and the other nucleus within the molecule, and a term due to the liquid crystal environment. The observed quadrupolar coupling is due to the interaction between the total electric field gradient and the deuteron quadrupole moment B where the obs = " 4 < ~7~ > n < ZZ - <* zz) F e S • (24) is the mean electric field gradient due to the liquid crystal and q is internal electric field gradient averaged molecule. If the internal contribution to over the internal motions of the can be calculated then an estimate of the external electric field gradient may be obtained. This requires a knowledge of the internal field gradient, q, at the parameter, S . 22 The value of S coupling using expression (18) site of the deuteron and the order can be calculated from the observed dipolar 2 2 and a quantum mechanically averaged value of q <r"°>. q Here, we expand r"° in a Taylor series in terms of the q parameter £ about the equilibrium position, r'g <r" > = r" (1 3 3 e -3<£> + 6<£ > 2 + . . . ) . dimensionless ( 2 5 ) where £ = ( - r.)/r,. (26) r <£> and < £ > are given by [51]: <f> and = -3(1/ + ^)a[B w ] + 4<J(J + l)>[B w ] ! e c 1 e 6 2 (27) Results and Discussion / 35 (28) <i > = 2(1/ + {)[B J ] Z where a, u and Q and B e e are the Dunham coefficient, equilibrium vibrational frequency g equilibrium rotational constant. taken [10] , l v The values from the literature [52]. of these constants have been For the internal electric field gradient we have [53] q(r) =q(r ) r V e The coefficients c and a Reid and Vaida [54], (1 + a ^ > 3 + a < * > + ... ) (29) 2 have been obtained by fitting the calculated values of 2 The value of <B>/<D> in the absence of the external field gradient is then (1 + (oc,-3)<£> + ( a 2+6-30. ,)<{ >) <B> ^ 2 = ( B / D ) " (. -3< > ? + < e » e • ( 3 0 ) The value of (B/D)eq has been obtained by inserting the experimental value [49] of <B>/<D> calculate for the <B>/<D> 01 state of D . 2 This for the 02 state of D value may then be used to and the 01 state of HD. These 2 two ratios have been measured in molecular beam experiments and our calculated values agree essentially exactly with the experimental values [49]. The internal contribution to B is then experimental coupling. dipolar <B>/<D> The from expression difference (30) between multiplied by this value and the the experimental value of B is the contribution due to the external field gradient. The value value of F r , of the experimental z may then be obtained from expression (24) using the known deuteron nuclear couplings quadrupole moment and Q D [54]. values In table 1 of Results and Discussion / 36 F for D z z value of F in 1132 and EBBA-d2 2 z z a r e s n o w n - It interesting to note that the 1 S is of opposite sign in these two liquid crystals. 4. Electric Field Gradient Molecular Quadrupole Moment Interaction The presence orienting moment of this mechanism [31]. electric field gradient immediately through This the interaction interaction is of second with the rank suggests a possible molecular quadrupole tensorial form and the interaction potential is given by [56]: U Q = " J 2 F u Qmol,IJ • (31) Because both the nematic liquid crystal, and hence F, and the hydrogen molecule have cylindrical symmetry, expression (31) becomes U The Q = " 2 F"ZZ QmoUz (f C O S 0 Z - { ) . (32) 2 2 contribution to the dipolar couplings and order parameters due to this mechanism maj' be calculated using expression (21) with A = value of Q has zz calculation for D 2 been obtained in 1132 from reference and EBBA-d 2 [55]. — $ F Q . The ZZ The results £2 of this are also shown in table 1. It is obvious from the good agreement between the experimental and calculated order parameters that this mechanism dominates the ordering of D 2 in these two liquid crystals. The negative order parameter in EBBA is also explained very simply. The discrepancies between the calculated and experimental order parameters are, in part, due to errors in the calculation of F Z 2 various physical constants required in the . However, these errors do not account for all of discrepancy Results and Discussion / 37 TABLE 1 Electric Field Gradients in 1132 and EBBA-d2 Solvent T (K) (esu S(D2) X10-3 t (exp) (calc) FZZ* XIOH) 1132 298 6.33±0.13 EBBA-d2 301 -7.23±0.20 9.28 8.57±0.65 -9.69±0.27 -11.07 * Fzz external electric field gradient calculated from the ratio of BID for D2 (see text and [10, 15, 32, 36]) I S t n e t S(D£) is the order parameter for D£. The experimental value is obtained from the dipolar coupling between the two deuterons. (The dipolar coupling for D2 in 1132 has been taken from [10].) The calculated values are obtained using expression (20) with A = ~±FzzQzzQzz has been calculated from the values given in reference [55]. and it is likely that other mechanisms play a role. 5. D 2 in 1132/EBBA Mixtures The change in sign of F z z in EBBA and 1132 raises the interesting possibility of producing a mixture of the two liquid crystals in which F intermediate between those in the would be a mixture in which F studying other solutes where z z two components = [32]. z z may have values Of particular interest 0 since this could provide a system for mechanisms F involving z Experimental dipolar and quadrupolar couplings [32] for D z 2 could be ignored. in several mixtures of 1132 and EBBA are shown in table 2 and the spectra obtained in three of these mixtures are shown in figure 7. As can be seen the spectral width, which is largely determined by the solute order parameters, depends very strongly on the liquid crystal composition. In 61.3wt per cent 1132/EBBA spectrum collapses to nearly a single line of width 2.4 at 310K the Hz. In fact, there are Results and Discussion / 38 FIGURE 7 Spectra of D2 in 1132/EBBA Mixtures 76j0wt per rant 61.3w I per cent 2&9wt per ent -1000 -«00 Frequency (Hz) 000 WOO 61.3 MHz deuterium NMR spectra of D2 dissolved in nematic mixtures of EBBA and 1132 at 31 OK. two shoulders with splitting 6|D| + 2|B| = and D versus temperature the values D = 4.7 Hz. By interpolating plots of B —0.45 ± 0.1 Hz and B = 1.0 ± 0.3 Hz are obtained. From table 2 it can be seen that F j , varies continuously z with the liquid crystal composition and that the order parameter follows the field gradient very closely. This is an important result since it provides us with a liquid crystal system in which one orienting mechanism can be controlled. From table 2 it can be seen that F z z and S do not pass through zero at exactly the same composition. Some of this discrepancy can be attributed to the experimental error in F ^ . However, it is likely that more than one mechanism is involved in the ordering of hydrogen. To investigate this, experimental and calculated order Results and Discussion / 39 TABLE 2 D 2 in Mixtures of 1132 and EBBA* 1132 D** (Hz) 0.0 23.9 49.3 55.0 61.3 66.1 76.0 100.0 65.95 + 0.50 33.67±0.14 11.15±0.18 5.59±0.07 -0.45±0.10 -5.38 + 0.14 -16.46±0.06 -49.76±0.60 wt per cent S(exp)t (Hz) (X10-3) ZZ* ( X 1 0 esu) 9.65 4.94 1.64 0.82 0.06 0.79 2.42 7.31 5.42 + 0.8 2.6910.2 0.4710.3 0.0110.1 0.6410.2 1.0610.2 2.1010.1 4.8210.9 -1535.22±1.50 -788.26±0.43 -267.6310.55 -138.02±0.29 1.0±0.30 116.06 + 0.42 373.07±0.20 1151.86+1.70 F 11 * All results are from reference [32] and are at 31 OK except for 55.0 wt per cent 1132 which is at 301.4K ** D is the dipolar coupling between the deuterons in D2 *** B is the quadrupolar coupling of the deuterons in D2 f Calculated from the experimental dipolar coupling $ Fzz the external electric field gradient due to the liquid crystal environment and is calculated from the solute couplings (see text) I S parameters for D [10, 2 and HD in 11 different liquid crystals have been obtained 32, 40]. The experimental values couplings using expressions (18) and (25). parameters for H , HD and D 2 are have 2 been calculated from (Experimental couplings the dipolar and order in. these liquid crystals at various temperatures tabulated in in appendix I.) The calculated values have been obtained using expression (21) with A = ~^zz^zz m ^ e v a ^ u e s °f Qzz have been obtained from reference [55]. In figure 8 the calculated values have been plotted against the experimental values. There are two striking features of this plot. First, the agreement between the experimental and calculated order parameters is very good especially since there are no adjustable parameters involved. This strongly Results and Discussion / 40 FIGURE 8 Ordering of Molecular Hydrogen in Nematic Liquid Crystals. -12 -6 0 6 S (exp) xlO 12 3 Experimental values vs calculated values of S , the order parameter of the molecular symmetry axis, of HD and D2 in 11 liquid crystals. The experimental values are obtained from the dipolar coupling and the quantum mechanically averaged value of <r~3>. The calculated values are obtained using expression (28) with A = — i^zzQzz- FzZ * obtained from the deuteron spectrum (see text). Triangles = HD, Circles = D2. a = EBBA, b = phase V, c = 1167, d = 23.9wt per cent 1132/EBBA, e = 49.3wt per cent 1132/EBBA, f = 55wt per cent 1132/EBBA, g = 61.3wt per cent 1132/EBBA, h = 66.1wt per cent 1132/EBBA, i = 76.4wt per cent j = 5CB-a,$-d4, k = 1132. For composition of phase V and 1167 see ref [10]. zz s Results and Discussion / 41 supports the conclusion that the electric Field gradient moment mechanism dominates the ordering of HD and D crystals and suggests that this may be true for most, - molecular quadrupole 2 in all of these liquid if not all, nematics. Second, the points lie on a straight line whose slope is slightly different from one and which does not pass through the origin. There are several reasons why such a plot may be obtained. If we arbitrarily make the value of (BfD)„ n in expression (30), 0.8 per cent larger the plot shown in figure 9 is obtained. As can be seen this small adjustment can account for most of the deviation in both the slope and intercept., Since we are quite confident of the calculated gas phase value of (B/D)eq, this implies a change in the electronic structure of the molecule between the gas and liquid crystal phases. Although it is quite small, it is unlikely that a change of this magnitude in the electronic structure would occur. A more reasonable explanation of these results is that the field gradient at the site of the deuterons is slightly different from the one which causes the orientation of the molecule. It is also interesting to note in figure 9 that all of the points lie very slightly to the positive side of the line of slope 1. 6. Isotope and Temperature Effects a. 60wt per cent 1132/EBBA In order to investigate the possibility of a second mechanism it is useful to consider the ordering of H2, HD and D 2 in a mixture where S solutes passes through zero. Let us assume z z for all three that the ordering is the result of the interaction between a mean field due to the liquid crystal and some property of the solutes. If all three solutes experience the same mean field, then the Results and Discussion / 42 FIGURE 9 Ordering of Molecular Hydrogen in Nematic Liquid Crystals -6 0 S (exp) x l O 6 12 3 The data are the same as figure 8 except that the equilibrium value of BID is 0.8 per cent larger (see text). Results and Discussion / 43 order parameters of the solutes should all pass though zero at the point where the mean field is zero. However, if more than one mechanism contributes to the ordering, then the order parameter may be zero as a result of the sum of the various contributions being zero. In this case it is likely that the order parameters of the three solutes will pass through zero at different temperatures and liquid crystal compositions. Dipolar couplings obtained for H2, HD and D 2 as a function of temperature in 60wt per cent 1132 are shown in figure 10. The couplings for HD and D respectively. D 2 2 have been scaled by the factors (y^y^) a n ^ ^H^D^ Unfortunately over the temperature range studied, the spectrum of is very poorly resolved and it was only possible to obtain dipolar couplings for this solute at three temperatures. It is very surprising to observe that all of the dipolar couplings obtained for D quadrupolar couplings for D 2 lie between the H 2 2 and HD values. The are more easily determined and they are plotted in figure 11 along with the values obtained for HD as a function of temperature. From the figure we see that Bp(D ) temperatures. This implies that S (D ) 22 consistent with D p D is 2 more negative 2 more positive B (HD) D for all is more positive than S (HD) which is 22 than D ^ . From the two figures it is also apparent that the dipolar and quadrupolar couplings do not pass through zero at exactly the same temperature. This is a direct result of the external electric field gradient. When the order parameter of D 2 or HD is zero then the dipolar coupling is zero and the quadrupolar coupling is entirely due to the external field gradient. If the electric field gradient - molecular quadrupole moment mechanism accounts for all of the ordering then the quadrupolar and dipolar couplings will pass through zero at the same point. The fact that they do not suggests that there is a second mechanism. It is also clear that the order parameters of the Results and Discussion / 44 FIGURE 10 Dipolar Couplings for H2, HD and D2 in 60wt per cent 1132/EBBA 290 294 298 T e m p e r a t u r e 302 (K) The dipolar couplings for HD have been scaled by the factor TTH^TD^ d those for D2 are scaled by (yft/yD)2. The temperature has been calibrated using the chemical shift differences in a sample of ethylene glycol. Circles = H2, triangles = HD and crosses = D2an three solutes do not pass through zero at the temperatures where S is negative |SJJJJ| temperatures where S is postive I S ^ I < |S > n same |S | > n | < |S temperature. At higher H |S^ | but at lower j. This also indicates the presence of a second mechanism. Because of the quantum mechanical averaging, second rank tensorial mechanisms such as the electric field gradient - molecular quadrupole moment interaction have | S J > I S ^ I > | S J . Since the contribution D H to the order parameters from a second mechanism must cancel the contribution Results and Discussion / 45 FIGURE 11 Quadrupolar Couplings for HD and D2 in 60wt per cent 1132/EBBA 294 290 302 298 Temperature (K) The temperature has been calibrated as indicated in the caption to figure 10. Triangles = HD and crosses = D2. due to the electric field gradient - molecular quadrupole moment mechanism when S = 0, the than S Sr, H second mechanism would have to have and S . Moreover, for this mechanism SJJJJ would have than to S U . This is an unusual isotope dependence investigating it somewhat further. to D more negative lie closer to and it is worthwhile Results and Discussion / 46 b. 50 and 55wt per cent 1132/EBBA In 50 and 55wt per cent 1132/EBBA the order parameters of the hydrogens are considerably larger than in 60wt per cent 1132/EBBA but are still much smaller than in the component liquid crystals where the electric field gradient - molecular quadrupole moment mechanism dominates the ordering. Thus results in these two mixtures should provide some information about the dependence of the contributions to the solute ordering on the composition of the liquid crystal. In particular we are interested in the isotope dependence of the couplings. In figures 12 and 13 the observed couplings are plotted as a function of temperature. The most striking feature of these two plots is the fact that the dipolar couplings go through a maximum. This is not too surprising, given the positive slope at lower temperatures and the fact that the dipolar couplings must go to zero at the isotropic phase transition temperature. It should be noted that the plots shown in figures 10-13 are quite different from those usually obtained for the liquid crystal solvent and which are predicted quite well by the mean field theory of Maier and Saupe [57]. responsible This shows quite clearly for the ordering of the liquid that the intermolecular crystal molecules forces are different from those responsible for the ordering of the solutes. It is also apparent from figures 12 and 13 that although the values of the scaled dipolar couplings for HD lie between those for H 2 and D2, they lie much closer experimental dipolar couplings for H , HD and D 2 2 the D 2 values. The in 50, 55 and 60wt per cent 1132/EBBA (as well as several other liquid crystals [10]) are shown in table 3. The difference between the observed HD dipolar coupling and the average of the H 2 and D 2 couplings is also given in table 3. Both the HD and D scaled to the proton value in calculating this difference. 2 have been As can be seen the Results and Discussion / 47 FIGURE 12 Dipolar Couplings for H2, HD and D2 in 50wt per cent 295 300 T e m p e r a t u r e 1132/EBBA 310 305 ( K ) The couplings have been scaled and labeled in the same manner as for figure 10. values obtained are remarkably constant in the three mixtures. In the other liquid crystals the values are again always positive and although there is some dependence on temperature and the solvent used, they are also approximately the same as obtained in the mixtures. Thus it appears that there is a roughly constant extra positive contribution to the dipolar couplings of HD, which implies an extra negative contribution to its order parameter. The mass dependence of electronic properties such as the quadrupole moment and polarizability is very Results and Discussion / 48 FIGURE 13 Dipolar Couplings for H2, HD and D2 in 55wt per cent 1132/EBBA N 8 260 — + bXI o o + + A A A O O + + a a, + + + 240 + A A Vi cd o A A A A 1-4 220 — T3 G) co o co O O O (D O 200 1 300 1 303 1 1 306 1 1 309 1 1 312 315 Temperature (K) The couplings have been scaled and labeled in the same manner as for figure 10. nearly linear and cannot account for this effect. HD does have a small dipole moment which could account for this effect through its interaction with a mean electric field. However, the magnitude of the electric field required is quite large (~10^ V cm"l) and a field of this magnitude would produce a very large order parameter through the interaction with the polarizability. Moreover, it is well known [17] that there is no correlation between the dipole moment of a solute and its ordering, thus this explanation is rather unlikely. The main difference Results and Discussion / 49 between HD and the other two isotopes is that its centre of mass does not lie half way between the two nuclei. This means that as the molecule tumbles about its centre of mass it sweeps out a larger area making it effectively larger that H positive 2 and D - If steric interactions are important this would lead to a more 2 order parameter for HD. Exactly the opposite is observed. This is an intriguing problem and its solution would lead to a better understanding of the intermolecular forces involved in the ordering of hydrogen. At present we are unable to explain the results adequately. 7. Summary For molecular hydrogen we have gained a great deal of insight into its ordering, although the picture is not yet complete. quadrupolar to dipolar couplings is We have seen that the ratio of the solvent dependent and that this solvent dependence can be explained in terms of an external electric field gradient due to the liquid crystalline environment. The ordering of the solute can then be explained largely as a result of the interaction between this electric field gradient and the molecular quadrupole moment. We have also seen that by mixing 1132 and EBBA the value of the field gradient may be varied continuously between the values found in the 1132/EBBA mixture at two 301.4K component there is liquid crystals. no measureable In a 55wt per cent external electric field gradient. In principle, this allows us to determine what role, if any, the electric field gradient plays in the orientation of other solutes and also allows us to investigate other mechanisms. In this mixture the order parameter of hydrogen is not zero which indicates the presence of a second mechanism. In mixtures where the electric field gradient is small a very unusual isotope dependence is observed Results and Discussion / 50 TABLE 3 H2, HD and D2 Dipolar Couplings in Various Liquid Crystals DHD- Temp. (K) L>HH DHD L>DD (Hz) (Hz) (Hz) 290.14 293.67 295.42 297.17 298.92 301.65 303.34 305.03 306.72 50wt per cent 1132/EBBA 376.8110.50 10.5010.02 65.3310.14 373.5610.50 65.8710.14 10.6110.02 383.8110.50 66.4210.14 10.6610.02 s.n 66.96 + 0.14 10.7810.02 s.n 67.3910.14 10.8110.02 394.2010.50 67.7110.14 10.9510.07 394.9210.50 68.0410.14 10.8610.07 399.8310.50 67.8110.14 10.8910.07 396.8010.50 67.7110.14 10.8910.07 300.68 302.43 304.18 305.93 307.69 309.44 311.19 312.94 314.69 55wt per cent 1132/EBBA 200.3710.09 35.7810.14 5.5910.14 203.6910.09 36.2810.14 5.7710.13 207.0610.09 36.7310.14 5.8810.10 209.6410.09 37.1510.14 5.9910.12 212.1710.09 37.5010.14 6.1210.05 214.5010.09 37.7910.14 6.1510.05 216.4310.09 38.0210.14 6.1610.05 217.8310.09 38.2010.14 6.1310.05 218.3510.09 38.1510.14 6.0810.05 291.34 293.29 294.26 295.23 296.20 297.17 298.14 299.11 300.08 301.05 302.03 303.00 304.94 60wt per cent 1132/EBBA -42.6110.23 n.m. -3.8210.28 -34.2510.12 -2.7910.14 -0.7210.04 -29.7510.12 -1.4410.11 n.r. -27.4210.24 -1.3010.11 n.r. -24.0010.19 -0.9310.11 n.r. -19.3910.47 -0.0410.11 n.r. -16.0110.21 0.4610.14 n.r. -13.5910.14 0.91 + 0.14 n.r. -7.3410.19 1.6210.14 n.r. -2.7610.38 2.1010.18 n.r. 0.6010.23 2.8010.11 0.2110.04 3.9610.14 0.4510.04 3.4910.11 10.2510.16 n.m. n.m. *(DHH + D D D ) * (Hz) 14.3910.51 17.1910.51 14.5910.51 11.6510.51 15.3510.51 10.7610.51 11.6210.51 14.2910.74 12.0710.70 10.9810.59 10.0910.66 8.3510.43 8.4410.43 8.7610.43 9.8710.43 10.3410.43 14.2310.41 13.4810.37 11.2110.35 Results and Discussion / 51 Table 3 cont. Temp. (K) (Hz) °HD (Hz) ODD (Hz) DHDi(D + D )* (Hz) H H DD 310 1141.2±0.30 1167T 185.85±0.15 30.6310.02 -9.8210.45 (19.6410.90)* 298 320 3059.3±0.50 2632.0 + 0.30 EBB At 506.6±0.60 434.35 + 0.30 83.110.20 71.2110.10 7.2811.82 2.5610.92 298 2670.5 + 0.30 Phase V t 442.710.60 72.510.30 10.3312.0 298 320 -2340.9±0.30 -1796.1±0.30 1132t -382.310.35 -292.0210.80 -63.27 -48.010.80 22.4910.94 14.2014.2 * D{JD and DrjD scaled to proton value. t Couplings in 1167, EBBA, Phase V, and 1132 are from reference [10] t Scaled by the order parameter of the liquid crystal director with rvspect to the magnetic field direction. n.m. not measured. n.r. D2 spectrum not resolved. s.n coupling not measured because of poor signal to noise. and it appears that there is an extra negative contribution to the ordering of HD. In other liquid crystals this effect is also observed, however, most of the isotope dependence is a result of differences in the quantum mechanical averaging for the three isotopes. At present, there is no adequate explanation of the extra contribution to the ordering of HD, although it is possible that it is related to the shift in the centre of mass. Results B. and Discussion / 52 METHANE 1. Introduction In the preceding section an average that electric field the interaction moment apparently have also seen have discovered gradient between that this field and may liquid strictly valid this then quadrupole orientational ordering. We the field gradient of composition between the values observed for crystals. Throughout be done separately. experiment, and the molecular of the observed results we have made the assumption rotations gradient crystalline environment and in nematic mixtures of 1132 and E B B A as a function two component that molecular hydrogen experiences due to its liquid accounts for most varies continuously the we of the hydrogen that the averages over the vibrations and Given the excellent agreement between assumption for molecules the discussion seems with well full justified. tedrahedral If this symmetry theory assumption is no anisotropic couplings should be observed since the order matrix is traceless and invariant to rotations of anisotropic methane [63], the molecule splittings are observed [33, 58], silane tetramethyl separation of vibrations description vibrations and rotations fact that for such in nematic symmetry. solvents and rotations of the effects solutes of the coupling anisotropic as [62, 63], adamantane [59, 65] and it is clear is not valid on the observed Experimentally, for such [59-61], tetramethylenetetramine for these between couplings, that the molecules. the An molecular is complicated by such effects occur because of the dependence of the solute-solvent interactions on the vibrations problem axis tin [64] and neopentane adequate the fixed and rotations is the nature of the solute-solvent of the solute. Thus the crux of the interactions. In this section we will Results and Discussion / 53 outline a general theory developed by Snijders, de Lange and Burnell [14-16, 50] for the orientation of solutes with internal motions in nematic liquid crystals. A summary of the application analogues dissolved in 1132, of the theory to methane and its deuterated EBBA and a 61.3wt per cent 1132/EBBA mixture will be given in order to determine what role, if any, the electric field gradient plays in determining the observed couplings. 2. Theory The details of the theory presented here are given in references [14-16] and we will only outline the main points. a. Solvent-Solute Interaction Potential In order to describe the anisotropic couplings observed for a solute in a liquid crystalline environment, we require the solvent-solute interaction potential. This potential is determined by the intermolecular interactions between the solvent and solute molecules. The observed couplings result from an average over all molecular motions in the presence of the potential. An exact treatment of this problem is prohibitively complex and instead we consider a simple model for the interaction. We will assume that the solute experiences a mean field due to the liquid crystal and that this mean field may be described by a second rank tensor F. The solute then interacts with this field via some second rank tensorial property p\ Thus we have specified the form of the interaction but not its exact nature. If we assume that the liquid crystal environment has axial symmetry about the Z direction, defined by the magnetic field, then the tensor F should also have this symmetry. The tensor may be written as the sum of an isotropic Results and Discussion / 54 part (the trace) and an anisotropic, symmetric, traceless part. The isotropic part cannot lead to anisotropic couplings and, because of the axial symmetry, the anisotropic part has only one independent non zero element which we will call G. The tensor j3 is a function of the positions may be written in terms potential we have of the normal coordinates ~ 3 Q . For the m interaction U = — - E Fjj /SJJ j = of the nuclei of the solute which 3 G ?, 1 M Q m ) (acosflfczcosfljz - - 6 ) 2 kl , Kl,m (oo) where I, J and Z are space fixed axes and k, I are molecule fixed axis. The potential is a function of both the vibrations and the orientation of the solute and will couple these two types of motion. b. Calculation of Observables This potential may then be treated as a perturbation on the freely rotating and vibrating molecule. wavefunctions Using first order perturbation theory, the perturbed may be obtained from the zeroth order functions for which simple products of the harmonic oscillator and rigid rotor functions are used. Expressions for the expectation values of observables such as the dipolar and quadrupolar couplings in terms of the various Gj3^;(Q ) may then be obtained. Because we m have only specified the form of the potential the values of these parameters and their dependences on the normal coordinates are not known. If we assume that only small displacements occur, the 0 tensor may be expanded in a Taylor series about the equilibrium in geometry in terms of the normal coordinates and truncated after the linear term. Thus the unknown parameters are reduced to Results and Discussion / 55 G(t u and G 9 p V 3 Q | m = Dipolar Couplings: In terms of these parameters the following expressions are obtained for the dipolar couplings Vii*" Dab = a kl b 88 G E E E 3 9d v,- a bk (-^) V m aQ m ab n a '* b ° , n _ r i g i d ( 3 4 ) (35) a f non rigid_ n . + D h , <AAi> • (— Qm=0 9 Q Qm=0 m 2 7TQ (36) V m where j a _ ab*Z ~ , 7 a 7 b , cos6> ~ I 2 _ *i - abJfc / 3 ab a b > i . (37) r, 4n i cos0 h 1 3 ^osfl^cosfl^ - -6 2 kl (38) and <& > - S kl kl (39) Results and Discussion / 56 U) m is the vibrational frequency of normal mode m. The rigid term arises from the ordering of the molecule in the absence of any coupling of the vibrations and rotations. The non rigid term arises from the vibration rotation interaction. The average over the orientations in (35) and (36) must mechanically for hydrogen but for other solutes this is be done quantum not necessary. In the absence of quantum effects this is given by / \j exp(-U(n)/k T) dO B / where respect fl to represents the B P a particular orientation space fixed (40) ex (-u(n)/k T) dn ' axes. For of the methane molecule the fixed interaction axes with energy, U, (expression (33)) is sufficiently small that the exponential in expression (40) can be expanded in a Taylor series and truncated after the linear term. The average over the internal motions in expression (35) is obtained by expanding expression (37) in terms of the normal coordinates and retaining only those terms which make a significant contribution. In a purely harmonic potential only the even terms will contribute. For methane the series may be truncated after the linear term and symmetry greatly reduces the number of non-zero parameters. Because the equilibrium geometry is tetrahedral the parameter Gt$ki q does not contribute e to the coupling. Similarly, the derivative of 0 with respect to the symmetric breathing mode, which retains the tetrahedral symmetry, does not contribute. This leaves only the derivatives with respect to the E (bending), F stretching) and F 2 2 (predominantly (predominantly bending) normal modes. Because the derivatives are evaluated in the equilibrium geometry which has tetrahedral symmetry, they are isotope independent. Thus the dipolar couplings of the series of isotopically substituted methanes may be calculated in terms of these three parameters. For Results and Discussion / 57 the symmetrically substituted methanes (CD^ and CH^) the average geometry is also tetrahedral so that only the non rigid term in expression (34) contributes. For the asymmetrically substituted species (CHgD, CH2D2 and CHDg) there is also a rigid contribution which arises because the average geometry, and hence the /3 tensor averaged over the molecular vibrations, does not have tetrahedral symmetry. partially A fit of deuterated derivatives of the (bend) modes 17 experimental dipolar species dissolved in couplings 1132 and [33] EBBA for in 0 tensor with the respect to the E , F2 has been done [14] and excellent the terms (stretch) agreement various of the and F2 between the experimental and calculated values is obtained. It is clear from these results that the form of the interaction potential is essentially correct. From the parameters obtained, it is found that the effects due to the F2 stretching mode dominate the observed couplings. Quadrupolar Couplings: The quadrupolar couplings have essentially the same form as the dipolar couplings and again two terms are obtained: p B D _ R " g i ~ D B . d + B R D non-rigid (41) eQ 2 B ° g d = ~Zh E < q k l > S k l ( 4 2 ) 2 B non rigid !<, = QD Results and Discussion / 58 2h ^ ( — ) x E E E Thus to calculate which determine the (43) (^~) <<&»AM> • quadrupolar couplings we the dipolar couplings, the require the average gradient tensor at the site of the deuteron <q>, value same of the parameters electric field and the derivatives of this tensor with respect to the various normal modes. The G/3 parameters may be obtained from the fit to the dipolar couplings and some simplifying assumptions may be made deviations values of due about the to the electric field electric field gradient vibrations will be gradient tensor ignored. tensor. elements Then, For the from because rigid part the equilibrium of the local C g y symmetry about the C-D bonds in the equilibrium geometry, there is only one independent element q . on Electronic structure calculations [66] for small distortions the local C g of Thus parameters q EBBA agreement the to the E and F quadrupolar couplings and 9 l y'3QF (str)" ^ ^ ( g y e( shown that symmetry is retained and that the derivatives the field gradient with respect neglected. and y have to 2 are ^ e 2 bending modes determined by two may be unknown experimental couplings in 1132 [33] in terms of these two parameters [15] also produces excellent between experiment and theory. Both parameters here should be electronic properties and independent of the liquid crystalline environment. While this is true of q ^ there is a large difference in the values of 3<l ^9QF (str) e obtained [15] which suggests that this description is not complete. 2 Results and Discussion / 59 3. Methanes in Liquid Crystal Mixtures In section molecular A we saw that the ordering of hydrogen is dominated quadrupole moment electric field gradient mechanism by the and that in a 55wt per cent 1132/EBBA mixture the external electric field gradient is zero. It is important to point out that at this composition the order parameter is not exactly zero which suggests the presence of a second orienting mechanism. In a 60wt per cent mixture the order parameter of hydrogen goes through zero as a function of temperature but there present. It seems likely that is the coupling in methane. methanes measured [36] small but appreciable field gradient electric field gradient also determining the were a in To investigate 1132, EBBA this, and plays the a a role in spectra 61.3wt of the per cent 1132/EBBA mixture at 310K. At this temperature the order parameter of as measured from the dipolar coupling, is zero in the mixture. An analysis of the dipolar couplings of the isotopically substituted methanes was done and these results are given in table 4. as described above Initial fits to results gave relatively large deviations. The measurements the EBBA were repeated several times and it was discovered the ^^CHgD sample was slowly changing over time probably because of a slow degradation of the liquid crystal. The cause of this was not investigated and with the 1 3 CHgD couplings removed a satisfactory fit was obtained. It is immediately apparent that the dipolar couplings of methane do not go through zero in the mixture, although they do he between the values in the two component liquid crystals. This suggests that the external electric field gradient which is responsible for the ordering of D for the observed dipolar couplings in methane. external field gradient there 2 is at least partly responsible If methane does experience will be an extra contribution to an the quadrupolar Results and Discussion / 60 TABLE 4 Methane in 1132, EBBA and 61wt per cent 1132/EBBA at 310K: Experimental and Calculated Dipolar Couplings (Hz) and G/3 Parameters* ' (10-5 erg cm-1) 1 Molecule Dipolar coupling Experimental Calculated EBBA 1 3 CH, 1 3 CH D HH CH HH HD CH CD HH HD DD HD DD DD 3 CH D 2 CHD 2 3 CD, G/3 = G/3 = G/3„ = 2 3 0.68 ± 0.05 10.87 ± 0.1 t t t t -9.15 ± 0.05 0.104 ± 0.10 0.26 ± 0.02 -0.59 ± 0.07 0.15 ± 0.10 0.015 ± 0.010 0.69 10.87 -9.15 0.096 0.25 -0.66 0.13 0.016 -0.197 1.498 -0.078 61.3wt per cent 1132/EBBA 1 3 1 3 CH, HH CH HH HD CH CD HH HD DD HD DD DD CH D 3 CH D 2 CHD 2 3 CD, G0 = 2 G/3 = GpY = 3 -0.197 0.919 -0.109 . 0.0 1.17 -3.476 0.60 -1.67 1.425 -7.174 0.022 0.17 -0.51 0.09 0.010 ± 0.18 ± 0.10 ± 0.18 ± 0.03 + 0.08 ± 0.02 ± 0.02 ± 0.003 ± 0.01 ± 0.02 ± 0.01 ± 0.003 0.13 1.17 -3.78 0.62 -1.68 1.49 -6.99 0.015 0.017 -0.53 0.09 0.003 Results and Discussion / 61 Table 4 cont. Molecule Dipolar Experimental Calculated coupling 1132 1 3 1 3 CH HH CH HH HD CH CD HH HD DD HD DD DD 4 CH D 3 CH D 2 CHD 2 3 CD„ G/3 G/3 G/3, 2 3 -5.17 -3.83 0.57 -8.27 0.31 -7.60 -0.04 <, 0.18 ± 0.14 ± 0.05 ± 0.07 ± 0.7 ± 0.4 ± 0.05 ± 0.07 0.15 -5.23 -3.99 0.66 -8.26 0.59 -7.46 0.02 -0.59 ± 0.07 ~ 0.09 ~ 0.0 -0.55 0.09 0.004 = -0.243 = 0.824 = -0.182 * G/3 ,G/3 and G/3q are the derivatives of the /3 tensor with respect to the E, F (stretch) and F (bend) symmetry modes (see text) t Unreliable results, see text 2 3 2 2 couplings given by: B In table ext 5 the experimental = _ 3 ( ^QDF _Z Z ) ( quadrupolar couplings are given 4 4 ) along with the values calculated as described in section 2 above and the uncorrected values of q^ and ^ l q^^^F (str) couplings ( e given 3qg^/3Qp ( ) 2 str a r e 2 in expression is affected obtained. If the contribution to the quadrupolar (44) is taken into account, only the value of [15, 16]. The corrected parameters shown in table 5 have been obtained using the value of from D 2 experiments. As can be Results and Discussion / 62 T A B L E 5A Methane in 1132, EBBA and 61wt per cent 1132/EBBA at 310K: Experimental and Calculated Quadrupolar Couplings (Hz) Molecule Experimental Calculated EBBA CH D CH D CHD CD„ 1 3 3 2 2 3 t -26.15 ± 0 . 1 0 -3.20 ± 0 . 1 0 18.40 ± 0 . 0 1 0 -25.93 -3.64 18.40 1.3wt per cent 1132/EBBA CH D CH D CHD CD, 1 3 3 2 2 3 -65.26 -49.00 -30.5 -15.17 ± 0.02 ± 0.03 ±0.03 ±0.03 -66.45 -46.79 -31.06 -15.63 1132 CH D CH D CHD CD„ 1 3 3 2 2 3 •90.08 •74.95 -55.20 •37.76 ± ± ± ± 0.35 0.10 0.10 0.10 -92.16 -71.62 -55.07 -39.14 t Unreliable results, see text seen the value of 3Qe^9QF (str) * 2 s n o w independent of the liquid crystalline environment as one would expect for a molecular property. This result provides very strong evidence that the deuterons in both D 2 and methane experience the same external electric field gradient. However, it is clear that for methane, the interaction with the molecular quadrupole moment does not account for all of the observed dipolar couplings. Results and Discussion / 63 T A B L E 5B Parameters' q (kHz) and 9q/9SF, (kHz A " ) With and Without External Field Gradient Corrections 1 eq Parameter EBBA 61.3wt per cent 1132/EBBA 1132 193 -61.4 187 -293 184 582 193 -286 86 187 -252 -10.1 184 -243 -72 Uncorrected 9q/9SF 2 Corrected 9q/9SF B (ext)t 2 D *q q is the principle component of the intramolecular field gradient tensor in the equilibrium geometry and 9q/9SF is the derivative with respect to the F e 2 2 stretching symmetry mode, evaluated in the equilibrium geometry, t Obtained from D experiments [10]. 2 4. Summary In summary, then, we have seen that the observed dipolar couplings for the series of deuterated methanes dissolved in nematic liquid crystals is a result of two contributions. non-rigid averaging contribution is of on substituted the of which dependent asymmetrically deviations For all average the is the a isotopically result of vibrational species there geometry is from substituted the state also fact of that the a rigid tetrahedral species the there is a orientational molecule. For the contribution from the symmetry. The dipolar couplings are predicted very well by assuming that the solvent-solute intraction is of second rank tensorial form and by fitting the calculated couplings to the Results and Discussion / 64 experimental values describe the in terms of three parameters, interaction between the liquid G/3 , crystal a n 2 mean d field, which G, and the derivatives of the second rank tensorial molecular property, p\ with respect to the E (bend), F 2 (stretch) and F (bend) 2 normal modes. The quadrupolar couplings may be calculated in terms of these same three parameters along with qgg. the principle component equilibrium geometry of the intramolecular field gradient tensor and 9q ^/8Qp ( ) e 2 gtr in the the derivative with respect to the F 2 stretching mode. Both of these parameters should be liquid crystal independent. However, this is true only when one takes into account the external electric field gradient as measured from correlation between the D 2 experiments. G/J parameters Finally, and the it is apparent from the external electric field gradient that it also plays some role in determining the dipolar couplings. However, unlike hydrogen where the electric field gradient quadrupole moment mechanism dominates, for methane it makes a relatively small contribution. C. ACETYLENE 1. Introduction In for the previous two sections we have seen that the observed dipolar couplings methane and hydrogen can be explained very well in terms of a simple model in which the solvent-solute interaction is of second rank tensorial form. In hydrogen we moment have mechanism also seen the electric field gradient - molecular quadrupole plays an important role. For methane, the interaction mechanism has not been specified, but there is strong evidence suggesting that this solute also experiences the same electric field gradient as hydrogen. These Results and Discussion / 65 two solutes have thus provided some insight into the solvent-solute responsible for the observed spectra. They are, however, interactions unusual solutes: hydrogen, because of its extremely small size and methane, because of its very high symmetry. In order to see whether the conclusions we have drawn here are generally applicable it is worthwhile considering a somewhat more typical solute. In choosing such a solute there are some important considerations. First, in order to perform the vibrational analysis in a manner similar to that used to analyse the methane small as spectra the possible number of normal modes and the of vibration should be number of dipolar couplings as large as as possible. Second, both the harmonic and anharmonic force fields must be known. A solute which fits these criteria is acetylene. It has 5 normal modes of vibration but by symmetry the derivatives molecular property, will of be the zero /3 tensor, with which describes respect to the some ungerade unspecified (asymmetric) modes. There are two of these which leaves 3 non zero derivatives. A total of 29 dipolar couplings isotopically substituted may be species measured are if all of the available. Thus possible the ratio C and of H unknown parameters to measured couplings is quite favourable for this solute. Diehl et acetylene-1,2- al. [21] have studied the spectra of acetylehe-l-^C and C£ in EBBA and Merck phase IV. The observed dipolar couplings where analysed in terms of the vibrationally averaged structure, but no account of the vibration-rotation interaction was taken. The ratios of the C-H and the C-C bond lengths, as determined from the dipolar couplings were found to be strongly solvent and temperature dependent. This result was explained in terms Results and Discussion / 66 of a two site model. In the model, the acetylene molecules were envisaged as occupying two different environments. In these two environments, the geometry and orientation of the molecules were allowed to vary. The ratio of the order parameters, bond lengths and enthalpies for the two sites were then taken as adjustable parameters and the experimental results were predicted. Although this model correctly predicts the observed behaviour, no evidence was given that the two sites actually exist. In light of the results for methane it seems likely that the solvent and temperature dependence of the dipolar coupling ratios is a result of the coupling between vibrations and rotations. In this section results for the analysis of various isotopically substituted species of acetylene in terms of the model used for methane presented. As will be seen very good agreement couplings is obtained. However, there are and hydrogen will be with the experimental dipolar some small remaining discrepancies. These discrepancies will be discussed in terms of a possible contribution to the C-C dipolar couplings from the anisotropy in the indirect coupling, as well as anharmonic corrections to the calculated rigid dipolar couplings. 2. Experimental Results The experimental dipolar couplings for the various isotopically substituted species in 1132, EBBA and 55wt per cent 1132/EBBA at 301.4K are shown in table 6. These couplings have been LEQUOR where necessary. obtained by fitting the spectra using the program Unfortunately, the quadrupolar coupling constants of the deuterons could not be determined accurately. Because of the large magnitude of the observed quadrupolar splittings (eg. ~120 kHz in 1132 at 30IK), the Results and Discussion / 67 TABLE 6 Acetylene in 1132, EBBA and 55wt per cent 1132/EBBA at 301.4K: Experimental Dipolar Couplings (Hz) Nucleus Measured Molecule Coupling 'H 13 C EBBA H- 1 3 C = 1 3 C-H D H C L>CH' L>HH DCC H- 1 3 C = 1 2 C-H D D D- C 1 3 = 1 3 C-D H C C H ' 1386.471±0.008 147.652±0.008 189.415±0.008 236.929±0.018 1386.3910.04 147.6510.04 189.3510.09 236.6110.05 1386.4±0.2 147.6±0.2 189.410.2 1386.3610.04 147.5010.04 218.7110.05 22.9710.05 4.4010.11 239.0210.12 DCD D CD' DDD HH- 1 2 C 1 2 = C = 1 2 1 2 C-H 189.410.2 C-D 29.210.2 DHD 55wt per cent 1132/EBBA H _ 1 3 C _ 1 3 C _ H D c D D D H - 1 3 C e 1 2 C - H -2412.9710.07 H . C H H H C C D D H' D H C C H -265.6810.07 -337.7810.07 -442.3810.16 H -2412.9110.16 -265.7110.16 -337.8210.11 -2413.5410.11 -265.6810.11 -337.8210.26 -440.5510.12 Results and Discussion / 68 Table 6 cont. Nucleus Measured Molecule Coupling 1 H 13 C 55wt per cent 1132/EBBA D- 1 C 3 _1 3 -380.5710.06 -41.1910.06 -7.9110.09 -444.0810.10 C-D D < CD DD D DCC H - 1 C _1 C - H DH H - 1 C _1 C - D DHD 2 2 2 2 -337.8 + 0.7 H -52.1±0.7 1132 H - 1 C 3 _1 C - H 3 DH C DcH 1 DH H Dec H - i C 3 _1 2 C - H • D C H DCH' DH H D- 1 C 3 _1 C - D 3 -5770.96±0.18 -632.5810.18 -804.9810.19 -1043.3310.42 -5770.7210.40 -632.3110.40 -805.1310.27 -902.2810.10 -97.3710.10 -20.1410.16 -1042.9910.16 DD C CD' DD D D D C C H - 1 C _1 2 C - H DHH -804.610.3 H - i C _1 2 C-D DHD -124.210.3 2 2 -5770.7710.04 -632.3010.04 -804.2110.10 -1043.4710.05 Results and Discussion / 69 contributions to the linewidths due to, e.g. temperature inhomogeneity, will also be large. In the best deuteron spectrum obtained the linewidths in the solute spectrum were ~100 Hz and very little fine structure could be resolved. As can be seen from the results in table 6, the signs of the dipolar couplings change between 1132 and EBBA in the same manner as those for hydrogen. However, the composition where the order parameter of acetylene is zero is quite different from that for hydrogen. In a mixture of approximately 16wt per cent 1132/EBBA the proton spectrum of 1 2 C H 2 2 is a single line at 340K (uncalibrated). Unfortunately this mixture has a fairly narrow nematic range and at lower temperatures a more ordered phase (probably smectic) Because of problems with temperature stability and homogeneity is formed. in the nematic range, a good set of couplings has not yet been obtained in this mixture. 3. Theoretical Considerations Before calculating the dipolar couplings for acetylene in terms of our second rank tensorial interaction model, there are some points which should be emphasized. First, unlike methane, the B tensor will not be isotropic and the largest contribution to the couplings will be the rigid effect which will be dominated by the molecule in its equilibrium geometry. Second, we need to consider the isotope dependence of the derivatives of the B tensor with respect to the various normal coordinates [14, structure of the 16, 67]. The B tensor is entirely determined by the electronic molecule. cartesian displacements Thus, the derivatives of B with are isotope independent. However, because respect to the of the Eckart conditions of zero angular momentum [68], the relationship between the cartesian Results and Discussion / 70 coordinates and internal coordinates parameters 9/3/3Q m | _ Q -0 is mass dependent. Thus, in general, the will be isotope dependent. Crawford [69] has shown m that the derivatives with respect to the internal coordinates of a second rank tensorial property, such as the polarizability, for the various isotopic species can be related to one another through and appropriate rotation of the molecule fixed axes. For methane the |3 tensor is invariant to such a rotation and thus the derivatives are isotope independent. This is not the case for acetylene. It is convenient, here, to relate the derivatives for the isotopically substituted species to a hypothetical molecule with atoms of unit mass [70]. The dipolar couplings may then be calculated in terms of parameters for this hypothetical molecule. Finally, we also note the order parameters for acetylene will be considerably larger than for the asymmetrically substituted methanes and therefore more terms in the expansion of expression (37) in terms of the normal coordinates will be required. 4. Calculation of Order Parameters Initial attempts at fitting the experimental couplings from table 6 did not converge and the values of the derivatives of the |3 tensor with respect to the two stretching modes were unreasonably large. This is a result of the fact that the dependence of the dipolar coupling on these two modes is effectively identical and of opposite sign. Removal of the derivative of the /J tensor with respect to either stretching mode yielded much better results which were independent of which parameter was removed. In table 7, the results where the derivative of the P tensor with respect to the C-C stretch has been set to zero, are presented. The value obtained for the derivative with respect to the symmetric Results and Discussion / 71 TABLE 7 Acetylene in 1132, EBBA and 55wt per cent 1132/EBBA at 30IK: Calculated and Experimental Dipolar Couplings (Hz) Nonrigid Total 1289.88 144.12 181.09 244.20 96.81 2.90 6.47 -5.54 1386.69 147.03 187.56 238.66 1386.43 + 0.04 147.65±0.04 189.38±0.09 236.7710.05 1291.04 144.25 181.25 95.62 2.79 6.33 1386.66 147.04 187.58 1386.4010.22 147.5710.16 189.4010.22 Dee 207.37 22.59 4.37 246.48 11.33 0.19 0.10 -7.31 218.70 22.78 4.46 239.17 218.7110.05 22.9710.05 4.4010.11 239.0210.12 Coupling Molecule Rigid Experiment * EBBA H- 13 C _1 C - H 3 D H C Dor D H H Dec H- 13 C _1 2 C-H D H C D C H ' D H H D- 1 3 C _1 C-D 3 D D C DCD' D D D H- 1 H- 1 2 2 C _1 C - H D H 181.41 6.19 187.59 189.410.2 C _1 C-D DHD 28.10 0.75 28.85 29.210.2 2 H 2 G/3 parameters 1" GB = -0.2546X10" erg G/3, = -1.318X10" erg c m ' G/3 = 0.0 G/3 = - 0 . 2 2 5 4 X 1 0 ' erg c m " rms error = 1.05 Hz* 1 3 eq 5 1 2 5 3 1 Results and Discussion / 72 Table 7 cont. Molecule Coupling Rigid Nonrigid Total Experiment * 55wt per cent 1132/EBBA H- 1 3 C _1 3 C -H DCH DCH' DHH Dec H- 13 C _1 2 C -H D H C D C H ' D D- 1 3 C _ 13 C -D H H H- 1 2 C _ 12 C C _ 12 -H C -D D H -331.94 -3.20 -335.14 -337.810.7 DHD -51.27 -0.33 -51.60 -52.110.7 H 1 3 1 2 5 -2412.9110.16 -265.7110.16 -337.8210.11 Dec e<? 3 -50.61 -2413.65 -1.44 -265.46 -3.32 -335.06 -380.5710.06 -41.1910.06 -7.9110.09 -444.0810.10 G/3 parameterst G/3 = 0.4116X10" erg G/3, = 0.7219X10" erg cm" G/3 = 0.0 G/3 = 0.0109X10" erg cm" rms error = 1.55 Hz* 5 -2363.03 -264.02 -331.75 -380.69 -41.06 -7.96 -444.19 D 2 -2413.26±0.13 -265.6910.13 -337.8110.27 -441.4610.20 -4.43 -0.07 -0.04 3.04 D D 1 -51.71 -2413.34 -1.69 -265.41 -3.44 -334.99 -3.01 -444.08 -376.26 -40.99 -7.92 -447.23 PCD DCD' H- -2361.63 -263.72 -331.55 -447.09 1 Results and Discussion / 73 Table 7 cont. Coupling Molecule Rigid Nonrigid Total Experiment 1132 H- 1 3 C = C-H 13 D H C DCH' DHH DCC H- C = 1 3 1 2 C-H D H C D CH' D H H D- 1 3 C = 1 3 C-D D D CD' D D C D D -5598.89 -625.59 -786.03 -1059.97 -173.20 -5772.09 -5.48 -631.07 -11.82 -797.85 9.00 -1050.97 -5770.87±0.18 -632.44±0.18 -804.59±0.21 -1043.40±0.42 -5600.90 -625.79 -786.30 170.32 -5771.22 -5.24 -631.03 -11.5 -797.80 -5770.72±0.40 -632.31±0.40 -805.13±0.27 -885.62 -96.48 -18.65 -1052.67 -17.58 -903.20 -0.36 -96.84 -0.16 -18.81 9.97 -1042.70 -902.28±0.10 -97.37 + 0.10 -20.14±0.16 -1042.99±0.16 Dec H _ i 2 C =1 2 C _ H -786.68 -11.18 -797.76 -804.6±0.3 -121.19 -1.30 -122.49 -124.210.3 D H H H- 1 2 C = 1 2 C-D DHD G/3 parameterst GB = 0.922X10" erg G/3, = -0.622X10" erg cm" G/3 = 0.0 G/3 = 0.1214X10" erg c m " rms error = 4.05 Hz* 1 3 eq 5 1 2 5 1 3 * Average of f f and C results t G is the anisotropy in the second rank tensor describing the mean field due to the liquid crystal. B is a second rank tensorial property of the solute. By, B and B are the derivatives of B with respect to the symmetric C-H stretch, C-C stretch and symmetric bend symmetry modes. t rms error between calcuated and experimental couplings. 1 1 3 2 3 Results and Discussion / 74 C-H stretch is thus the sum of the derivatives with respect to both stretching modes. The calculated dipolar couplings have been obtained using the force field of Strey and Mills [71] and the equilibrium geometry and definitions of the internal coordinates are also from this reference. For the rigid dipolar couplings we have . ^ <d a b W _ . > = d eq eq ab ab,*i sb.hl _ _ "2 2 m 2 n d v i ^kl ] 7 n n QQm=0 =0 <Q > m ^ m m m m ™m dd 2 i ^ „ + E (—— ) (3_Q a b ah ii fJ) 3Q n m n _ Qm-Qn=0 (45) <Q„Q„> <Q Q> m ++ •. •. •. n For a purety harmonic potential terms of order < Q >, < m Qm - » ^m » < e t c - are zero. The rigid couplings reported in table 7 have been calculated with only the harmonic terms (even powers of Q) included and the series has been truncated after the sixth order term. The non-rigid contribution to the couplings 0 has been calculated from expression (36) in the section on methane. The angular dependent parts of the expressions for the calculated classically from expression (40) acetylene the interaction energy is rigid and non-rigid couplings are (also in the section on methane). For sufficiently large that the expansion of the exponential in expression (40) can no longer be truncated after the linear term. Because of this, the integrals in expression (40) have been calculated numerically. As can be seen from table 7, quite good agreement between the total calculated couplings and experimental values is obtained. The rigid effect clearly dominates the observed couplings, but it is also apparent that non-rigid effects also make a considerabale contribution. The change in sign of the order parameter is explained very simply as a change in the sign of the anisotropy in the mean field, G. In 55wt per cent 1132 the order parameter is positive and the change in the order Results and Discussion / 75 parameter between the three liquid crystals is qualitatively consistent electric field gradient - molecular quadrupole moment mechanism. with the Thus, the ordering of acetylene can be explained as resulting from two mechanisms; (i) the electric field gradient molecular quadrupole moment mechanism and (ii) a second mechanism which makes a positive contribution to the order parameter. Although the fit obtained is quite good, the discrepancies lie well outside the error in the experimental dipolar couplings. It is also clear that the discrepancy between parameter theory and increases. carbon-carbon experiment One couplings due increases reason to for the as this the could anisotropy magnitude be in the the J of the order contribution to the coupling. the Thus, calculations were repeated with these couplings removed. These results are shown in table 8. As is obvious, better agreement with the experimental couplings is obtained. The values of the anisotropy in J , obtained from the difference between the calculated and experimental carbon-carbon dipolar couplings, are roughly the same in all 3 liquid crystals. However, the values obtained for the deuterated compound are consistently lower than those for the protonated species and all of the values are considerably larger than the theoretical values of 99.3 Hz [5, 72] and 85.0 Hz [5, 73]. A second possible reason for the discrepancies in the initial fit is the neglect of the anharmonic terms in the calculation of < d calculated using the anharmonic force field [14, £/ . > a D The linear term can be 16, 74] and for acetylene it is of the order of one per cent of the total coupling. Some initial fits with this term included were slightly worse than those shown in table 7. However, it is Results and Discussion TABLE 8 Acetylene in 1132, EBBA and 55wt per cent 1132/EBBA at 301K: Calculated and Experimental Dipolar Couplings and the Anisotropy in the Indirect C-C Coupling (Hz) Molecule Coupling Calculated Experiment EBBA H- 1 C 3 H- 1 C 3 _ 13 C - H _1 2 1386.41 1386.4310.04 DCH* 148.40 147.6510.04 DHH 188.78 189.3810.09 Dec 242.84 236.7710.05 1386.44 1386.4010.22 DCH C-H D H C 148.42 147.5710.16 D H 188.81 189.4010.22 DCD 218.74 218.7110.05 22.99 22.9710.05 D D 4.49 4.4010.11 Dec 248.33 239.0210.12 D H 188.83 189.410.2 29.06 29.210.2 D CH* H D- 1 C _1 C-D 3 3 D CD' D H- 1 C 2 _ 12 C-H H H - 1 C _1 C - D 2 2 DHD Gj3 parameters t GB = -0.2570X10erg G/3, = - 1 . 1 4 9 X 1 0 - erg cm" GB = 0.0 GB = - 0 . 1 6 9 6 X 1 0 " erg c m " rms error = 0.47 Hz* 1 3 eq 5 1 2 5 3 1 AJQC / 76 Results and Discussion / 77 Table 8 cont. Coupling Molecule Calculated Experiment AJrx 55wt per cent 1132/EBBA H- 1 3 C = 1 3 C-H -2412.96 -267.26 -336.64 -452.07 -2413.26±0.13 -265.69 + 0.13 -337.8110.27 -441.4610.20 -2413.37 -267.32 -336.72 -2412.9110.16 -265.7110.16 -337.8210.11 DCC -380.51 -41.32 -7.99 -449.50 -380.5710.06 -41.1910.06 -7.9110.09 -444.0810.10 D H -336.79 -337.810.7 DHD -51.87 -52.110.7 D H C Dor D H H Dec H- 1 3 C = 1 2 C-H D H C D CH* D H H D- 1 3 C =' C - D 3 DCD DcD 1 D D D HH 1 2 _ 1 2 C = C 1 s 2 C-H i 2 C _ H D G/3 parameters'!' GB = 0.4137X10" erg G/3, = 0.4396X10- erg c m " G/3 = 0.0 G/3 = -0.0578X10" erg c m " rms error = 0.92 Hz* 1 3 eq 5 1 2 5 3 1 305.53 156.47 Results and Discussion / 78 Table 8 cont. Molecule Coupling Calculated Experiment ** AJcc 1132 H- 1 3 C= 1 3 C-H H H c -5771.23 -635.97 -802.23 -1065.89 -5770.87±0.18 -632.4410.18 -804.5910.21 -1043.4010.42 D D . D H H -5770.49 -635.94 -802.19 -5770.7210.40 -632.3110.40 -805.1310.27 D D D < D D -902.72 -97.53 -18.90 -1056.76 -902.2810.10 -97.3710.10 -20.1410.16 -1042.9910.16 -802.16 -804.610.3 -123.18 -124.210.3 D D D D C C H H C H- 1 3 C= 1 2 C-H C H C H D- 1 3 C= 1 3 C-D C C D D D C C H- 1 H- 1 2 2 C= 1 C= 1 2 2 C-H D C-D D D H H H G/3 parameters t Gp = 0.928X10" erg G/3 , = -1.354X10erg cm" G/3 = 0.0 G/3 = -0.0548X10" erg cm" rms error = 2.11 Hz* 274.17 168.66 1 3 eq 5 1 2 5 1 3 *C-C coulpings not included in fit (see text). ** Average of H and C results *** AJcc I S the anisotropy in the indirect C-C coupling tensor, t G is the anisotropy in the second rank tensor describing the mean field due to the liquid crystal. /3 is a second rank tensorial property of the solute. /3,, /3 and /3 are the derivatives of /3 with respect to the symmetric C-H stretch, C-C stretch and symmetric bend symmetry modes. t rms error between calcuated and experimental couplings. 1 1 3 2 3 Results and Discussion / 79 clear that the contribution from the cubic term, which will partially cancel the linear term, must also be included. This work is not yet complete and it remains to be seen whether including these terms will result in better agreement with the experimental couplings. 5. Summary Although, the analysis presented in this section is not complete, it is clear that the vibration-rotation interaction couplings observed for acetylene. plaj's an important role in determining the When this interaction is included there is no need to invoke any specific interactions between solute and solvent. The negative order parameter in interactions [21], is experiences the Finally, there EBBA, accounted same is which external some evidence had for quite electric been explained naturally by field gradient in terms assuming of that specific acetylene as molecular hydrogen. of a contribution to the carbon-carbon dipolar couplings due to the anisotropy in the indirect spin-spin coupling. However, until all significant contributions to the calculated couplings account, this result remains somewhat speculative. have been taken into Results and Discussion / 80 D. SIZE A N D S H A P E E F F E C T S 1. Introduction In the preceding sections we have been able to make some significant conclusions about the nature of the interactions which lead to the orientation of hydrogen, methane and acetylene in nematic agreement with the observed solvents. For these dipolar couplings is small obtained solutes excellent by considering the solvent-solute interaction as being purely electrostatic and of second rank tensorial form and, in the cases of methane and acetylene, by taking into account the effects of the interaction between vibrations and rotations. For hydrogen we have also seen electric that field interaction the deuterium nuclei gradient between due this to their external in D liquid field 2 and HD experience crystal gradient environment. and the an external Moreover, the molecular quadrupole moment accounts for most of the orientational ordering of the hydrogens. For the deuterated methanes the intramolecular quadrupolar couplings depend upon two parameters: (i) the principal component of the electric field gradient tensor along the CD bond when the molecule is in its equilibrium geometry derivative of the electric field gradient with respect to the F 2 and (ii) the stretching mode. Both of these parameters should be molecular properties which depend only upon the electronic structure of the molecule. We have found that if the effect of the external electric field gradient, as measured from the hydrogen experiments, taken into account, the values of the parameters obtained by fitting to is the experimental couplings become liquid crystal independent as one would expect for a molecular property. This result suggests that deuterium nuclei in both hydrogen and methane experience the same electric field gradient. Results and Discussion / 81 For acetylene we have shown that although vibrations play an important role in determining the dipolar couplings, the interaction between the liquid crystal and the solute in its equilibrium geometry accounts for most of the observed orientational ordering. Moreover, a qualitative examination of the molecular order parameters in 1132 and EBBA suggests that the interaction between the external electric field gradient (as measured from the D 2 experiments) and the molecular quadrupole moment also plays an important role in determining the ordering of this solute. Finally, we 301.4 have shown that in a 55wt per cent mixture of ± 0.3K, the deuterons in D 2 1132/EBBA at and HD experience no external electric field gradient. It is important to point out, however, that both hydrogen and acetylene are appreciably ordered in this mixture and anisotropic couplings are observed for methane. This implies gradient and the that although molecular the quadrupole interaction between the moment is an electric field important orienting mechanism, it is not the only one. In this section we will use these results as a basis for studying a series of small solutes in an attempt to investigate the contributions to their ordering. We shall begin by making some assumptions based on the above results. shall assume that all solutes experience First, we the same electric field gradient as the deuterons in D - Thus, in the 55wt per cent 1132 mixture we assume that the 2 electric field contribute internal gradient to the molecular - molecular quadrupole ordering. Second, we motions may be shall done moment mechanism assume that separately the from the does not averaging over averaging over Results and Discussion / 82 molecular reorientation. This means that the dipolar coupling between two spins a and b with the solute may be written D = - ( * ^ T - ) <r "b> S , ab (46) b ab a where S =<lcos 0 2 ab abz - |> (^) is the order parameter for the a,b direction in the molecule. If we consider only the interaction between the solute and its average environment and neglect quantum mechanical effects the order parameters of the molecule are given by / i j ~ (3cos0 cos0 - (5^) exp(-U(n)/k T) dfl iZ 2f jZ B (48) e x p ( - U ( n ) / k T ) df} B where the integrations are over all orientations St=U{8^,8^2) and U(S2) is the mean potential describing the interaction between the solute and its environment. The problem, then, is to define U(fi). Results and Discussion / 83 2. Short Range Interaction Model c. The Interaction Potential In the limit of low solute concentration, the potential, U(fl), is determined by the intermolecular interactions between solute and solvent averaged over all motions of the molecules. In general, these interactions may be divided into two types: (i) long range interactions where the intermolecular distances are much larger than the molecular intermolecular dimensions distances and are (ii) shorter short than (or range interactions comparable to) where the the molecular dimensions. The mean potential may then be written as u(0)=u (n)+u (a). SR LR For the long range contribution to the potential, the liquid crystal provides an average electric electronic field, property of between the electric the field solute. gradient, Two etc. examples which of this interacts are the with some interaction permanent dipole moment or molecular polarizability of the solute and a mean electic field due to the liquid crystal. The interaction between the molecular quadrupole moment of the solute and the mean electric field gradient is also of this type. As we have seen this last mechanism plays a dominant role in the orientation of hydrogen. Hydrogen is an extremely small solute, however, and for important. larger solutes the short range interactions will become increasingly Results and Discussion / 84 b. The Short Range Potential The short range potential, which is a result of the coulombic repulsion between the electron clouds of the solvent and solute molecules, is effectively the energy required to displace the solvent molecules upon introduction of the solute. It is the angular dependence of this energy with respect to the liquid crystal director which will lead to orientation of the solute. We will assume required to displace the liquid crystal molecules director is negligible perpendicular to the averaging over the compared director. to the Then, to model the liquid crystal as in a direction parallel to the energy avoid relative positions of the that the energy required the solute for complexities and solvent displacements inherent molecules, in we an elastic tube aligned parallel to the external Z direction. A schematic representation of this is shown in figure 14. If the liquid crystal molecules are assumed to be rigid then the walls of the tube may be displaced in the X Y direction but will remain parallel to the Z direction. The displacement caused by the presence of the solute then leads to a restoring force on the solute which we assume to be proportional to the deformation of the tube in the X Y direction: F(Q) = - k / r(B,Q) dB = -kC(fi), (49) where r(B,Sl) is the vector in the XY plane from the origin to the tube surface, B is the angle between the laboratory X axis and r(B,Sl) (see figure 15), C(Q) is the circumference of the deformed tube, and k is the Hooke's law force constant Results and Discussion / 85 FIGURE 14 Elastic Tube Model for the Short Range Interactions T H 11 liquid crystal ,;, , / J- . M t::::H:::::::H::::*::::HHHn:H::::n:!:!::^ MnllIllllIIIIIIINIIIIflllllIIIIIIII=l|i^^iillllllllIl^pJtllll=llillflfi i " _ _ =_ The solute, represented by an array of van der Waals spheres, displaces the liquid crystal which in turn exerts a restoring force on the solute proportional to the displacement. The liquid crystal, shown here only in the plane of the paper for clarity, is displaced in such a way that the displaced volume is a tube whose walls are parallel to the laboratory fixed Z axis. fl is the angle between the molecule and liquid crystal axes. The notation fl is meant to describe the orientation of a molecule of any symmetry. For molecules with Cs or greater symmetry the single angle 8 is used to describe their orientation. v for the restoring force. The mean potential energy arising from the interaction of the solute with the elastic walls is then u^n) = -k / F(n) dc(Q) = | kc (n). 2 In general, the circumference of the deformed tube, C(O), will depend on how the interaction between the modeled. If we surface assume that of the molecule and the both surfaces are surface hard then of the tube C(fl) is is the Results and Discussion / 86 A. Definition of symbols used in expression (49) r is a vector from the origin to the surface of the projected molecule and (} is the angle between r and the X axis. B. Maximum circumference: The maximum circumference of the elastic tube representing the deformation of the liquid crystal solvent by a solute at some orientation is shown in bold face. C. Minimum circumference: The minimum circumference of the elastic tube representing the deformation of the liquid crystal solvent by a solute at some orientation is shown in bold face. Results and Discussion / 87 cirumference of the projection of the solute molecule onto the X Y plane as shown in figure 15. This circumference depends on how the solute molecule is described and how its outer array surface of nuclei with fixed is defined. positions Here, in the we model the molecule fixed solute as a rigid axis system. Each nucleus then has a sphere of radius r associated with it, where r is the van der Waals radius of the atom. The electron distribution around an atom in a molecule is, of course, not spherical and it could be represented by some other shape. For example, than the in aromatic systems the 7T orbitals extend a orbitals so that a better representation much further of the electron distribution would be an elipsoid or some complicated shape which more closely resembles the molecular orbitals. However, in the interest of keeping the model as simple as possible we shall describe the solute as a collection of overlapping spheres. The projection of the molecule in the X Y plane is thus an array of overlapping circles. From this array of circles we want to calculate the circumference of the projection of the outer surface of the molecule. This surface is also the surface of the elastic tube and represents inaccessible the boundary of the region of space which is to the solvent molecules. This region is, however, not easily defined and its circumference should be calculated somewhat differently for each solute. In order to avoid this, we approximate the circumference of the projection of the outer surface of the molecule by either the maximum (figure (figure 15b) or minimum 15c) circumference of the array of overlapping circles, both of which are defined in the same way for each solute molecule. A detailed description of the calculation of the maximum and minimum circumferences in terms of the angles 6 and <f> is given in appendix II. The definitions of the angles 8 and <j> are also given in the appendix. Results and Discussion / 88 c. Calculation of Order Parameters Having calculated C(J2) and thus U(£2), the order parameters of the solute may then be calculated by integrating numerically over all orientations of the solute. In terms of our model S^, for example, is given by: / / ( 3 cos 0 - l)exp(-kC (6\0)/2k T) sin0d0d0 z B 0—0 8=0 o ^zz TT _ 2 2 f TT/2 f " • (51) exp(-kC (0,0)/2k T) sin0d0d0 B 0=0 0=0 It is often assumed that the interaction potential may be expanded in terms of the Legendre polynomials and then truncated after the first term. In such a case the potential becomes independent of <f> for a solute with a Cg or higher symmetry axis. We have not made this assumption and symmetry requires only that the potential be invariant to the symmetry operations of the point groups of the solute and its environment. Thus our potential is independent of 0 only for solutes with a C^, axis. If we calculate the integrals in expression (51) using a trapezoidal rule and a step size of 7T/36, the values of the resulting order parameters are accurate to 4 figures. This model is appealing because it should be applicable to a large number of solutes and involves only one adjustable parameter, the force constant k, which should be same for all solutes. Results and Discussion / 89 3. Ordering of solutes in 55wt per cent 1132/EBBA a. Solutes with Cg or Higher Symmetry' In order to test the model described above, we will begin by measuring the order parameters of some simple solutes with Cg or higher symmetry in 55wt per cent 1132/EBBA where we assume that they experience no external electric field gradient due to the liquid crystal. Experimental dipolar couplings for these solutes dissolved crystals, EBBA have in 55wt per cent and 1132, 1132/EBBA and the two component liquid are given in table 9. Where necessary the couplings been obtained by fitting the experimental spectrum using the computer program LEQUOR [28]. All of these solutes have Cg or higher symmetry and it can be shown [1] that for such molecules there is only one independent order parameter, S , zz where z is the symmetry axis. Experimental values of S zz for the solutes dissolved in 55wt per cent 1132/EBBA have been obtained from the measured dipolar couplings using expression (46) and a structure from electron diffraction or microwave data. These results are shown in table 10. As a test of the model for the short range interactions the order parameters have also been calculated using expression (51). The results of this calculation are also given in table 10 and are plotted in figures 16a and 16b, against the experimental values. In figure 16a the minimum circumference has been, used and in figure 16b the maximum circumference. The points have been labeled in the same manner as the results in tables 9 and 10. The van der Waals radii have been taken from reference [75] and the values of k=5.0 and 4.0 dyne cm** for the minimum and maximum circumferences have been determined by a least squares fit to the lines of slope 1. As can be seen the correlations are quite good. This Results and Discussion / 90 TABLE 9 Spectral Parameters (Hz)* for Solutes with C3 or Higher Symmetry at 301.4K Solute Parameter 1132 55wt per cent 1132 EBBA a) 1,3,5-tribromobenzene DRH -213.56(3) -159.2(3) -111.90(9) b) c) 1,3,5-trichlorobenzene Benzene -213.23(3) -987.99(2) -191.73(2) -125.05(2) -156.60(3) -708.47(2) -137.15(3) -89.32(4) -148.12(3) -461.48(1) -88.46(1) -57.34(1) -2340.9(3) -382.5(6) -63.3(3) 1461.5(5) 1351.5(5) 203.69(9) 36.19(14) 5.77(13) -140.4(4) -137.22(14) 3059.3(5) 506.5(5) 83.1(5) -1943.6(5) -1817.2(5) 2853.10(14) -401.43(14) 2.395ppm 1656.75(3) -232.89(3) 1.319ppm 595.23(1) -82.47(1) 0.338ppm -4763.2(3) -4607.4(8) -506.3(3) -490.5(8) -834.6(7) -832.8(6) -644.6(3) -643.81(8) -2131.23(8) 1075.1(5) 1060.4(3) 114.9(5) 112.0(3) 180.2(12) 178.6(2) 145.8(4) 140.0(7) 2761.16(3) -219.31(1) 2.75(2) 1493.03(4) -118.74(4) 2.69(7) D H H ortho D eta para D m D d) Hydrogen D H H ODD e) Propyne*"" B HD D CH D 3 HCH 3 C. Shift f) Acetylene D H C D CH D C C D H H g) 2-butyne D HH' J HH' -234.73(8) -389.4(2) -298.56(7) 1200.06(3) -95.68(3) 2.77(5) Results and Discussion / 91 Table 9 cont. Parameter Solute h) i) 2,4-hexadiyne 2795.38(12) -74.36(8) 2.77(5) 1886.50(3) -334.44(2) 1227.34(2) -223.68(1) 878.36(1) -158.59(1) D H -837.97(6) -105.84(9) -132.50(5) -195.74(15) -646.08(4) -81.81(6) -102.48(3) -151.084(10) -545.56(8) -69.64(12) -87.12(9) -127.4(2) -571.06(3) -106.90(4) -76.26(4) -22.47(6) -3.14(8) 3.78(10) -583.42(3) -109.18(3) -77.78(4) -22.40(6) -3.07(8) 3.73(10) -812.73(2) -151.86(3) -108.09(4) -22.40(5) -3.03(6) 3.75(8) HH' J HH' D is trans c 1,3,5-trifluorobenzene EBBA 3824.06(11) -101.06(9) • 2.69(7) D D j) 55wt per cent 1132 4547.59(14) -120.51(9) 2.75(2) Dgem Allene 1132 D D HF,ortho HF,para D F F H k) Hexafluorobenzene Dortho Dnieta Dpara Jortho Jmeta Jpara 1) Bromomethane D H 988.41(9) m) Iodomethane D H 960.00(13) n) Nitrous Oxide DNN -63.21(10) -7641(10) o) Nitrogen**** NN -36.3 -48.7 H H D 326.59(13) 426.5(4) 223.02(9) * J couplings not reported here have been taken from the literature and the references are given in table 11. ** C. Shift is the proton chemical shift between CH and CH . resonance is to low field in all cases. 3 The CH 3 * * * 77ie first set of numbers refers to the proton spectrum and the second set to the 13C spectrum. **** Dip i Q ar C0U pH g n S from reference [134]. Results and Discussion / 92 FIGURE 16a Short Range Interaction Model (Minimum Circumference) S (exp) Experimental vs calculated order parameters in 55wt per cent 1132. The calculated order parameters have been fit by the least squares method to the solid lines of slope 1. The labeling of the points refers to the molecules indicated in tables 9 and 10. T = 301.4 K, k = 5.0 dyne/cm, F^z 0.0 esu, correlation coefficient = 0.985. = Results and Discussion / 93 FIGURE 16B Short Range Interaction Model (Maximum Circumference) -0.2 0.0 0.2 0.4 S (exp) Experimental vs calculated order parameters in 55wt per cent 1132. The calculated order parameters have been fit by the least squares method to the solid lines of slope 1. The labeling of the points refers to the molecules indicated in tables 9 and 10. T = 301.4 K, k = 4.0 dyne/cm, Fzz ~ 0.0 esu, correlation coefficient = 0.982. Results and Discussion / 94 T A B L E 10 Experimental and Calculated Order Parameters for Solutes with C 3 or Higher Symmetry Dissolved in 55wt per cent 1132/EBBA at 301.4K Solute S(exp)* minimum circumference a b c d e f g h i j k 1 1,3,5-tribromobenzene 1,3,5-trichlorobenzene Benzene Hydrogen H HD D Propyne Acetylene 2-butyne 2,4-hexadiyne Allene 1,3,5-trifluorobenzene Hexafluorobenzene Iodomethane 2 2 -0.215 -0.211 -0.181 -0.0007 -0.0008 -0.0009 0.160 0.089 0.140 0.359 0.129 -0.200 -0.219 0.022 -0.210 -0.196 -0.112 0.014 0.014 0.014 0.118 0.074 0.186 0.400 0.107 -0.155 -0.184 0.070 S(calc)** maximum circumference polarizability -0.259 -0.229 -0.175 0.012 0.012 0.012 0.132 0.074 0.187 0.207 0.202 -0.180 -0.197 0.106 -0.237 -0.225 -0.121 0.011 0.011 0.011 0.106 0.062 0.170 0.358 0.092 -0.177 -0.250 0.079 * Order parameters obtained from the experimental dipolar couplings given in table 9. ** Order parameters calculated using the short range interaction model with the minimum and maximum circumferences, and the polarizability electric field mechanism. supports the idea that in this mixture the ordering of solutes is dominated by their size and shape and that the external electric field gradient experienced by all solutes is small. The fact that the correlation is not affected greatly by using the maximum provide or minimum a reasonable circumference representation suggests of the surface that both approximations of the molecule within the limits of the model. As may be seen in table 10 the order parameter which is predicted the least Results and Discussion / 95 accurately and is in fact of the wrong sign, is that of hydrogen. We have seen, in the preceding sections that most of the orientational ordering of this solute is due to the electric field gradient - molecular quadrupole moment mechanism and that quantum mechanical effects must be taken into account in order to describe the isotope dependence of the order parameters. That the contributions to the ordering due to short range effects should be small for this molecule and that this contribution cannot be described using the same force constant as used to describe the orientation of the larger solutes, is not surprising since for small values of the circumference the repulsive forces between the solvent molecules themselves should allow the solute to tumble more freely. Moreover, it may be necessary range to consider quantum mechanical effects in the description of the short interactions since such effects do play a role in the long range interactions. In considering other possible mechanisms for the orientation of the solutes it is useful to consider briefly the many cases 16, [10, 12, nature 13, 17] of the molecular proposed properties mechanisms involved. In involve molecular properties which are approximately equal to the sum of contributions associated with each bond. It is important to note, however, that this is not true of the molecular quadrupole moment. usually depends The contribution associated in some way on the size of the atoms with a given bond involved so that the molecular property is a function of the size and shape of the molecule as a whole. Because also predict the of this, mechanisms order parameters involving such in the molecular properties should 55wt per cent 1132/EBBA mixture reasonably well. To investigate this possibility order parameters for the various Results and Discussion / 96 solutes have been calculated using as an orienting mechanism between the molecular polarizability of the solute (given in table the interaction 11) and the anisotropy of a mean square electric field due to the liquid crystal. The results of this calculation are shown in table 10 and are plotted in figure 17. The value of the anisotropy in the mean square elecric field has been obtained from a least squares fit of the calculated order parameters to the experimental values -i and is found to be 2.27X10 suggests that quadrupole if the r o o cm" . The correlation is again quite good. This iO contribution to - electric field gradient the solute mechanism ordering is from sufficiently the small, molecular then any molecular property which is related to the size and shape of the solute can be used to predict the order parameters. b. Solutes of Lower Symmetry Because solutes with C3 or higher symmetry have only one independent order parameter, it is not possible and variations Solutes of in lower the to distinguish between imperfections average symmetry, environment however, have experienced two or by more in the model different independent solutes. order parameters and because the order parameters must describe orientation in exactly the same environment, such solutes present a more rigorous test of the model. Moreover, deuteration of the liquid crystal allows us to monitor the experimental conditions by using the deuteron spectrum as an internal standard. A series of solutes with and &2h y B Tnme ^y r are listed in table 12 along with their experimental order parameters in 55wt per cent 1132/EBBA-d2. Order parameters for the solutes with Cg or higher symmetry have also been measured in this mixture and the results are also shown in table 12. The order parameters, Results and Discussion / 97 T A B L E 11 Solute Polarizabilities and References Solute Polarizability (cm 2 4) xx a-zz yy - a a a) 1,3,5-tribromobenzene b) 1,3,5-trichlorobenzene c) Benzene d) Hydrogen e) Propyne f) Acetylene g) 2-butyne h) 2,4-hexadiyne i) Allene j) 1,3,5-trifluorobenzene k) Hexafluorobenzene 22.39 18.88 11.73 0.71 4.54 2.78 5.29 7.78 4.47 11.56 12.18 22.39 18.88 11.73 0.71 4.54 2.78 5.29 7.78 4.47 11.56 12.18 13.34 11.36 6.54 1.01 7.57 4.53 9.51 12.42 9.02 6.19 6.19 1) Iodomethane m) Bromomethane 1) Acetone 2) Furan 3) Thiophene 4) Pyridine 5) Fluorobenzene 6) Chlorobenzene 7) T T F 6.47 4.81 7.15 7.75 9.51 11.94 11.88 14.52 17.62 6.47 4.81 5.16 5.36 7.27 6.21 6.32 7.42 4.96 8.92 6.74 7.05 7.48 11.02 11.02 11.37 15.84 25.74 f a References Structure J [76P [76]a [79] [81, 82] [79] [79] [79] [80, 85]d [97]b [87]b [87, 122]b [87]b [34]b [90] [93]a [93]a [90] [99]b [90] [104] [58] [84] [21] [86] [88] [91] [94] [94] [97] [100] [102] [83] 8) 2,6-difluoropyridine 9) Iodobenzene 10) 1,2-dichlorobenzene 11) 1,2-dicyanobenzene 12) 1,3-dichlorobenzene 13) 1,3-dinitrobenzene [105] [102] [107] [109] [110] [111] 14) 1,4-dichlorobenzene 15) 1,4-dibromobenzene [114] [115] [77] [78] [80] [10] [80] [80] [80] [80] [80] [80, 78] [78, 89] [80] [80] [92] [95] [96] [98] [101] [9 IF [83] [106]c [119]e [108]c [109] [108F [112, 113]c [115] [115] c a) Extrapolated to zero frequency using frequency dependence of benzene given in reference [79]. b) Extrapolated to zero frequency. c) The structure has been calculated from data in the given reference(s). d) The polarizability has been estimated using data from the given references. e) Regular hexagon assumed. f) Tetrathiofulvalene Results and Discussion / 98 FIGURE 17 Polarizability - Electric Field Mechanism S (exp) Experimental vs calculated order parameters in 55wt per cent 1132. T = 301.4 K, C<Eu > - <E >) = 1.59X105 statvolt cm'l, F = 0.0 esu. The calculated order parameters have been fit by the least squares method to the solid line of slope 1 (correlation coefficient = 0.947). The labeling of the points refers to the molecules indicated in tables 9 and 10. 2 Z X 1/2 zz Results and Discussion / 99 TABLE 12 Experimental and Calculated Solute Order Parameters in 55wt per cent 1132/EBBA-d2 at 301.4K Solute S(exp)£ S(calc) fit to all results S.Rb Acetone s yy Fur an Syy s Thiophene 22 xx yy $zz J D Pyridine J xx Syy Szz Fluorobenzene J xx D yy Szz Chlorobenzene ?xx ^yy Szz TTFi ?xx S 2,6-difluoropyridine zz ^xx yy Szz s separate fit to each solute Pol.c S.Rd Pole 0.0711 (2) -0.0720 (1) 0.0009 (1) 0.0754 0.0718 •0.0036 0.0366 0.0669 0.0303 0.0711 0.0696 0.0033 0.0408 0.0736 0.0340 0.0907 (1) •0.1365 (1) 0.0458 (2) 0.0636 •0.0894 0.0258 0.0469 0.0770 0.0302 0.0974 0.1336 0.0362 0.0841 0.1383 0.0520 0.0586 (1) •0.1490 (4) 0.0905 (5) 0.0547 •0.0998 0.0451 0.0057 0.0997 0.0940 0.0810 0.1474 0.0664 0.0061 •0.1235 0.1166 0.1093 (2) -0.1545 (4) 0.0452 (6) 0.0727 •0.1138 0.0412 0.1129 0.1632 0.0502 0.1016 0.1566 0.0550 0.1067 0.1555 0.0488 0.0507 (1) -0.1906 (1) 0.1399 (2) 0.0242 •0.1384 0.1142 0.0645 -0.1640 0.0995 0.0274 -0.1804 0.1530 0.0760 0.1947 0.1188 0.0070 (1) -0.2037 (2) 0.1967 (3) 0.0109 -0.1501 0.1610 0.0622 0.2196 0.1575 0.0201 0.1856 0.2057 0.0635 -0.2229 0.1594 -0.081 (3) 0.2514 (17) 0.3320 (17) 0.1513 -0.2570 0.4083 0.1299 -0.3520 0.4819 0.1129 •0.2189 0.3318 0.0499 •0.2828 0.3327 0.1529 (2) -0.1908 (2) 0.0379 (4) 0.1231 •0.1503 0.0272 0.1577 0.1868 0.0291 Results and Discussion / 100 Table 12 cont. Solute S(exp)a S(calc) fit to all results S.Rb Iodobenzene ^xx S zz 1,2-dichlorobenzene J xx Syy 1,2-dicyanobenzene xx J °yy S 1,3-dichlorobenzene 22 ^>xx Syy Szz 1,3-dinitrobenzene xx J *yy ^zz Pol.c -0.0594 0.1735 0.2329 -0.0642 0.1811 0.2453 0.0647 (1) -0.2248 (3) 0.1601 (2) 0.0254 0.1754 0.1500 0.0245 -0.2040 0.1795 0.063 (11) -0.247 (13) 0.184 (24) 0.0272 0.2100 0.1828 0.0254 -0.2255 0.2001 0.1958 (4) -0.2243 (3) 0.0285 (7) 0.2029 0.1870 0.0159 0.2092 0.1915 -0.0177 0.1961 (2) -0.2287 (2) 0.0326 (4) 0.2211 •0.2096 •0.0116 0.2090 0.1999 -0.0091 -0.0721 (1) -0.2467 (4) 0.3188 (5) 0.1046 0.2003 0.3049 -0.1123 -0.2105 0.3228 0.1420 0.2173 0.3594 -0.1420 -0.2170 0.3590 ^xx 1,4-dibromobenzeneg &XX yy Szz -0.1021 (1) 0.2546 (1) 0.3567 (2) Benzene Szz -0.1756 -0.1220 -0.1608 1,3,5-trifluorbenzene S'zz -0.2014 -0.1687 -0.1654 zz -0.2380 -0.2121 -0.2141 1,3,5-trichlorobenzene s J S.R.d -0.0243 (5) -0.2146 (6) 0.2389 (11) 1,4-dichlorobenzeneg "yy $zz separate fit to each solute Pole Results and Discussion / 101 Table 12 cont. Solute S(exp)a S(calc) fit to all results 1,3,5-tri- s bromobenzene hexafluoro- S benzene 2,4-hexadiyne Szz zz 2 2 separate fit to each solute S.Rb Pol.c -0.1886 -0.2268 -0.2434 -0.2280 -0.1990 -0.1804 0.3664 0.4439 0.1846 2-butyne 0.1943 0.2081 0.1670 Allene Szz 0.1302 0.1197 0.1302 Bromomethane Szz 0.0676 0.0664 0.0734 Iodomethane Szz 0.0650 0.0786 0.0942 Hydrogen -0.0008 0.0152 0.0109 Acetylene Szz s 0.1123 0.0825 0.0663 Propyne Szz 0.1559 0.1314 0.1177 zz S.Rd Pole a) Order parameters obtained from the experimental dipolar coupling constants. The bracketed numbers are errors. The axis system is chosen such that: z is the symmetry axis, y is perpendicular to the plane of the molecule, and x is perpendicular to y, and z. b) Order parameters calculated using a model for the short range repulsive interactions. The force constant governing the interaction has been obtained by a least squares fit of S x> Syy, and S for all solutes to the experimental values. X zz c) Order parameters calculated using the interaction between the molecular polarizability of the solute and the mean electric field due to the liquid crystal solvent. The anisotropy in the mean suare electric field has been obtained by a least squares fit of S , Syy, and S for all solutes to the experimental values. xx zz d) Order parameters calculated in the same manner as b) except that the force constant has been fit for each solute independently. e) Order parameters calculated in the same manner as c) except that the anisotropy in the mean electric field squared has been fit independently for each solute. f) For TTF z is along the C~C double bond joining the two rings, x is parallel to the vector joining the sulfur atoms within one ring, and y is perpendicular to x and z. g) For 1,4-dichlorobenzene and 1,4-dibromobenzene the z axis is chosen such that it passes through both halogen atoms. The y axis is perpendicular to the plane of the molecule and the x axis is perpendicular to y and z. Results and Discussion / 102 which are taken from references [39] and [116], have been calculated from the measured dipolar couplings using the computer program SHAPE [30] and a geometry from one of electron diffraction, microwave spectroscopy or NMR using liquid crystal solvents. References for the J couplings and geometries are given in table 11 and the definitions of the molecule fixed axes for the solutes are given in the footnotes to table 12. Order parameters for the solutes have been calculated using the short range interaction model and the minimum circumference of the projection of the molecule onto the XY plane (figure 15C). The value of the force constant, k, has been obtained by a least squares fit of the calculated order parameters Sxx , S,„, yj and S zz to the experimental values and is found to be 5.55 dyne cm"*. The order parameters calculated from this fit are listed in table 12 and are plotted against the experimental values in figure 18 (the solutes with Cg or higher symmetry are also included in this plot). As can be seen the fit is very good and the model is successful in explaining the results for all the solutes. These results may also help to distinguish among various mechanisms which correlate roughly as the size and shape of the solute. With this in mind order parameters for the solutes have been calculated using the polarizability-electric field interaction as the orienting mechanism. Unfortunately, the low symmetry of the solutes makes determination of their polarizability tensors more difficult because one needs to measure three properties of the solute which depend on the three independent elements of the polarizability tensor, a , a Q [99]. In some measurements recent zz publications [90, 93, 99] the authors of the Kerr and Cotton-Mouton effects with and a ^ — a^, have combined light scattering Results and Discussion / 103 FIGURE 18 Short Range Interaction Model S (exp) Experimental versus calculated order parameters for solutes in 55wt per cent 1132/EBBA-d2- The calculated order parameters, also shown in tables 10 and 11, have been obtained using the minimum circumference. For solutes with C2v or I>2h symmetry the values of Sxx> Syy, and S are plotted (circles). For solutes with C3 or higher symmetry only Szz is plotted (trianlges). The definitions of the axes, x, y, and z are given as footnotes to table 12. The force constant describing the short range interactions has been obtained from a least squares fit to the solid line of slope 1. T = 301.4K, k = 5.55 dyne cm-1, F^z 0.0 esu, correlation coefficient = 0.985. zz V = Results and Discussion / 104 measurements to obtain the elements of the polarizability tensors of six of the solutes with C and 2w s y m m e t ry used here. These solutes are labeled 1-6 in table 11. Their polarizabilities extrapolated to zero frequency are also given in table 11 along with calculated values for tetrathiofulvalene (TTF). These polarizabilities have been used to calculate order parameters for the solutes and the value of the mean square electric field has squares fit of S^, S^, and S •i 1.64X10 are o r 2 2 been obtained from a least to the experimental values, and is found to be o V cm . The results of this calculation are shown in table 12 and - plotted in figure 19. As can be seen, the agreement between the experimental and calculated values is quite good, although a slightly better fit is obtained with the short range interaction model (figure 18). As with the solutes of C3 or higher symmetry we are unable to clearly distinguish between these two mechanisms because both are functions of the size and shape of the solutes. In the above discussion we have tested our model by comparing order parameters for different solutes with one another. There is no reason to expect that all solutes will be governed by exactly the same mean field, as each may experience a different distribution of local environments. For example, it is well known that various tetrahedral solutes experience different local magnetic susceptibility anisotropies and therefore must he in different average environments [117, 135]. For the solutes of order parameters, S , xx S yy and D and S zz 2 n any two of the three non zero are independent of each other. Since all three order parameters must depend on the environment in the same way, they should be predicted exactly by the model if the form of the potential is correct. To investigate this, order parameters for these solutes have been calculated using Results and Discussion / 105 FIGURE 19 Polarizability - Electric Field Interaction 0.6 0.4 0.2 - O o ^ 0.0 -0.2 - -0.4 -0.4 -0.2 J 0.0 I I I 0.2 L 0.4 0.6 S (exp) Experimental versus calculated order parameters for solutes in 55wt per cent 1132IEBBA-d2' The calculated order parameters, also shown in tables 10 and 11, have been obtained using the interaction between the polarizability tensor of the solute, (given in table 12), and the mean electric field squared due to the liquid crystal. For solutes with C2v &2h symmetry the values of S x> Syy, and S are plotted (circles). For solutes with C3 or higher symmetry only S is plotted (triangles). The anisotropy in the mean electric field squared has been obtained from a least squares fit to the line of slope 1. T = 301.4K, (<E*> - <E >) = 1.35X10 & statvolt cm' , F z z 0.0 esu, correlation coefficient — 0.941. o r X V zz zz Z ± 1 = 112 Results and Discussion / 106 the size and shape model and fitting the force constant seperately for each solute. The order parameters obtained from this calculation are shown in table 12 and the force constants are shown in table 13. As can be seen the agreement between the experimental and calculated order parameters is excellent. This demonstrates very clearly that the scatter in figure 18 is primarily due to our assumption that each solute experiences exactly the same average environment and that the form of our potential is correct. Order parameters have also been calculated using the polarizability - electric field mechanism and fitting the anisotropy in the mean square electric field separately for each solute. These order parameters are also shown in table 12 and the values obtained for the anisotropy in the mean square electric field are shown in table 13. Although the agreement with experiment is good, in most cases it is not as good as obtained using the short range interaction model. In order to see more clearly how well the relationships between the order parameters for each solute are predicted we define an asymmetry parameter _ ll ~ S S (52) 22 such that for the experimental order parameters |Sgg|>|S22l |Snl- This definition > ensures that zero and 17 calculated from the experimental order parameters lies between one. Figures experimental values 20a and 20b show the comparison of 77 and those calculated from the between the separate fits to each solute for the short range and polarizability electric field interactions respectively. Although there are some deviations, the correlation shown in figure 20a is Results and Discussion / 107 T A B L E 13 Force Constants, Electric Fields, Liquid Crystal Deuteron Splittings and Liquid Crystal Order Matrix Elements (55wt per cent 1132/EBBA-d2 at 301.4K) 2 Solute k» (dyne cm-1) < E 2 U > ~ 1 > AI/QC (10 (kHz) cm" ) < E 2 DHD<* (Hz) S zz SQ.c.)e S^-S, 2 Acetone 5.37 2.26 16.91 738 0.610 0.016 Fur an 8.54 3.74 17.93 778 0.643 0.016 Thiophene 8.47 2.54 17.46 750 0.620 0.014 Pyridine 7.90 1.92 17.18 755 0.624 0.018 Fluorobenzene 7.48 2.46 17.65 765 0.633 0.016 Chlorobenzene 7.03 2.05 17.81 733 0.606 0.006 TTF 4.50 1.29 17.82 750 0.620 0.010 2,6-difluoropyridine 7.12 16.89 743 0.614 0.017 lodobenzene 5.83 17.58 733 0.606 0.009 1,2 -dichlorobenzene 6.64 17.68 750 0.620 0.010 1,2-dicy anobenzene 6.09 17.77 765 0.633 0.014 1,3-dichlorobenzene 5.72 17.75 750 0.620 0.011 1,3-dinitrobenzene 5.24 17.64 773 0.639 0.018 1,4-dichlorobenzene 5.87 17.60 753 0.623 0.013 1,4-dibromobenzene 5.55 17.84 765 0.633 0.014 Benzene 17.51 715 0.591 0.005 1,3,5 -trifluorobenzene 17.81 743 0.614 0.009 1,3,5 -trichlorobenzene 17.76 693 0.573 -0.003 1,3,5 -tribr omobenzene 17.59 753 0.623 0.013 Results and Discussion / 108 Table 13 cont. <E„ > 2 ka Solute <E*> S(l.c.)e (kHz) D D (Hz) Hexafluorobenzene 16.95 733 0.606 0.015 2,4-hexadiyne 17.76 693 0.573 -0.003 2-butyne 17.59 753 0.622 0.013 Allene 17.69 735 0.608 0.008 Bromomethane 17.47 730 0.604 0.009 Iodomethane 17.70 780 0.645 0.019 Acetylene 17.54 752 0.622 0.014 Propyne 17.65 705 0.583 0.001 (10l5v (dyne cm" ) 1 d H 2 SxX~Syy cm" ) 2 a) k is the force. constant which governs the short range interactions (see text). The values of k presented have been obtained by a least squares fit of the calculated values of S ^ , Syy, and S for the given solute to the experimental values. zz 2 „ 2 b) E | | > <Ej^ >is the anisotropy in the mean electric field squared due to the liquid crystal (see text). The values presented have been obtained by a least squares fit of the calculated values of S , Syy, and S of the given solute to the experimental values. < — xx zz c) A P Q IS the quadrupolar splitting of the deuterons ortho to the nitrogen in the EBBA-d.2 contained in the 55wt per cent 1132/EBBA-d2 mixture (see figure 5). d) DJJTJ is the dipolar coupling between the deuterons and the adjacent protons of the EBBA-d2 contained in the 55wt per cent 1132/EBBA-d2 mixture (see figure 5). e) S(l.c.) refers to the order parameters of the aniline ring of the contained in the 55wt per cent 1132IEBBA-d2 (see figure 5). EBBA-d2 Results and Discussion / 109 Experimental versus calculated values of the asymmetry parameter: 17 = (SJJ — S22) I £>33 (see text) for solutes with C2v T>2h symmetry in 55wt per cent 1132/EBBA-d2. o r A)Short Range Interaction Model The asymmetry parameters have been calculated using a model for the short range interactions. The force constant describing the interaction has been adjusted separately for each solute (see tables 11 and 13). The labeling of the points refers to table 11. T = 301.4K, Fzz = 0.0 esu, correlation coefficient = 0.7620. The solid line of slope 1 has been drawn to aid the eye. Results and Discussion / 110 FIGURE 20B Asymmetry Parameters for C2v and ®2h Solutes 1.2 03 0.8 / 05 Cd O / / 0.4 06 /CD2Q7 0.0 -0.4 A -0.4 i I 0.0 I I 0.4 I I 0.8 I I I 1.2 V (exp) Experimental versus calculated values of the asymmetry parameter: 17 = (SJJ — S22) I £>33 ( tiext) for solutes with C2v or D2h symmetry in 55wt per cent 1132/EBBA-d2. see BJPolarizability Model The asymmetry parameters have been calculated using the interaction between the polarizability tensor of the solute and the mean electric field squared due to the liquid crystal. The mean electric field squared has been adjusted separately for each solute (see tables 11 and 13). The labeling of the points refers to table 11. T = 301.4K, F22, 0.0 esu, correlation coefficient = —0.4795. The solid line of slope 1 has been drawn to aid the eye. = Results and Discussion / 111 resonably good for the short-range interaction model. By comparison the asymmetry parameters calculated using the polarizability - electric field mechanism (figure 20b) do not agree as well with experiment. While this may be a result of inaccuracies in the polarizabilities we have used, it is unlikely that this could 2 2 account for all of the discrepancy. Because the values of k and (<E|) > - <E^ >) have been adjusted separately for each solute the scatter in figures 20a and 20b arises entirely because of inaccuracies in the form of the potential. These results suggest that the polarizability - electric field mechanism cannot account for all of the orientation of the solutes, although it cannot be ruled out entirely. This is also consistent with the results of a recent study by Emsley et al. [9] on anthracene dissolved in several liquid crystals. They found that dispersion forces, which depend on the solute polarizability, could not account for the solute order parameters. However, it should be pointed out that they did not take the electric field gradient - molecular quadrupole moment mechanism into account. c. Liquid Crystal Spectrum as an Internal Standard We now turn our attention to information which may be obtained from the deuteron spectrum of the liquid crystal. In the preceding discussion we have seen that slightly different values for the mean field are required to explain the experimental order parameters for each solute. One reason for this may be the variations in experimental conditions between samples. Because the deuteron spectrum is strongly dependent on the ordering of the liquid crystal molecules, which in turn has a strong dependence on temperature [57] and solute concentration, we can use the spectrum to monitor the experimental conditions. The deuteron spectrum of EBBA-d2 in the nematic phase consists of four lines Results and Discussion / 112 which arise from the quadrupolar couplings of the deuterons and the dipolar couplings between each deuteron and the adjacent ring protons [118]. The dipolar couplings between the resolved deuterons and contribute to the and the other spins linewidth. In the in the molecule are not 55wt per cent 1132/EBBA-d2 mixture containing no solute at 301.4K the value obtained for the quadrupolar splitting, A V Q , is 17.94 kHz and for the dipolar coupling, DJJJJ, is 774 Hz. The contribution of J^^C—1 Hz) to the dipolar coupling has been ignored. We define the asymmetry parameter, rj^, for the electric field gradient tensor, q, along the CD bond as: <iaa - ^ c c 7?q = • (53) qbb where the a axis is in the plane of the aromatic ring and perpendicular to the CD bond, b is along the CD bond, and C is perpendicular to both a and b. A value of 0.04 [119] will has been obtained for rj from single crystal studies of benzene-dg and we shall use this value for the aromatic deuterons in EBBA-d . We 2 also assume a value of 185 kHz for the quadrupolar coupling constant, e qQ /h [118]. If we then take a value of 2.49X I O distance [80] and 120.0° 2 D x zz cm for the ortho HD for the CCD angle the following order parameters for the anline ring (see figure 21) are obtained: S_ = S = - 8 —0.312, S„, = x 0.634. FIGURE 21 Structure and Molecule Fixed Axes; EBBA-d2 yy —0.322, and Results and Discussion / 113 Values of S and zz - S for the EBBA-d2 aniline containing various solutes at a variety of temperatures 14 and 15b. In all cases — ring for samples are given in tables 13, is approximately zero so that the rigid aniline ring has almost axially symmetric ordering. The values of S xx — Syy shown in table 14 are in good agreement with the value of —0.013 obtained by Diehl and Tracey for the same liquid crystal [118]. It is interesting to note that although the values of S pure EBBA-d2 (table xx 14) — Syy are small they are consistently and consistently positive for the negative for 55wt per cent 1132/EBBA-d2 mixture (tables 13 and 15b). These order parameters may be used as an internal standard for comparing the experimental conditions between samples. However, because of the relatively large errors associated with D ^ we shall use the quadrupolar splitting AI>Q for such comparisons. For the purposes of making small corrections for the effects of finite dilution and temperature differences among the samples we shall assume (i) that F ^ for a given sample with a given liquid crystal quadrupolar splitting is the same as in the sample containing D2 with the same liquid crystal quadrupolar splitting, and (ii) that k for a sample with a given liquid crystal quadrupolar splitting is the same as splitting. for As a reference our reference sample we with will use the same results liquid for crystal quadrupolar the solutes furan and thiophene. The experimental order parameters for these two solutes as a function of temperatures in the 55wt per cent 1132/EBBA-d2 mixtures are reported in references [24, 39] and are shown in table 15a. The liquid crystal splittings are given in table 15b. The value of k for each of the two solutes at each temperature has been obtained by fitting the calculated order parameters to the Results and Discussion / 114 TABLE 14 EBBA-d2 AT 301.4K: Liquid Crystal Order Parameters", Deuteron Quadrupolar Splittings and Dipolar Couplings Sample a) b) c) d) e) f) g) h) i) j) k) 1) m) n) o) A ^Q** (kHz) DHD*** 21.30 20.92 21.29 21.18 21.28 21.25 21.32 21.46 21.33 21.22 20.58 21.24 21.22 21.48 21.04 805 830 823 753 785 810 820 835 813 823 775 833 820 814 834 1,3,5-tribromobenzene 1,3,5-trichlorobenzene Benzene Hydrogen Propyne Acetylene 2-butyne 2,4-hexadiyne Allene 1,3,5-trifluorobenzene Hexafiuorobenzene Bromome thane Iodomethane Furan Thiophene S SxX zz Syy (Hz) 0.666 0.686 0.680 0.623 0.649 0.670 0.678 0.690 0.672 0.680 0.641 0.689 0.678 0.673 0.690 -0.009 0.000 -0.005 -0.020 -0.014 -0.008 -0.006 -0.004 -0.008 -0.004 -0.010 -0.002 -0.005 -0.009 0.000 * The order parameters of EBBA-d2 refer to the aniline ring (see figure 5). * * A UQ is the quadrupolar splitting of the deuterons in EBBA-d2*** ^KD the dipolar coupling between the deuterons and adjacent protons in EBBA-d2I S experimental values. These results are shown in figure 22 as a function of A P Q . It is clear from the plot that the dependence of k on A P Q is roughly linear over the range that has been measured, and the slope is to a good approximation the the same for both solutes. It is also apparent that the small variation in the liquid crystal splittings shown in table 13 cannot account for the variation in k from one sample to another. These results suggests the solutes are governed by slightly different force constants, either because they sample different local environments, or because of inadequacies in our model, and not because of variations in the experimental conditions. Results and Discussion / 115 T A B L E 15A Furan and Thiophene in 55wt per cent 1132/EBBA-d : Experimental Solute Order Parameters 2 Solute Temp. (K) Furan Thiophene a Sxx Syy 304 306 308 310 312 314 320 325 0.0907(1) 0.0884(1) 0.0863(1) 0.0841(1) 0.0820(1) 0.0798(1) 0.0735(1) 0.0669(1) -0.1365(1) -0.1335(2) -0.1311(4) -0.1285(2) -0.1258(1) -0.1227(1) -0.1152(4) -0.1055(6) 0.0457(2) 0.0452(3) 0.0448(5) 0.0443(3) 0.0437(2) 0.0430(2) 0.0418(5) 0.0389(7) 304 306 308 310 312 314 320 325 0.0586(1) 0.0572(1) 0.0559(1) 0.0545(1) 0.0531(1) 0.0517(1) 0.0469(1) 0.0423(1) -0.1490(4) -0.1459(2) -0.1427(2) -0.1395(2) -0.1363(2) -0.1326(2) -0.1206(4) -0.1113(17) 0.0905(5) 0.0887(3) 0.0868(3) 0.0850(3) 0.0832(3) 0.0809(3) 0.0737(5) 0.0690(21) a) The temperature reported is the dial setting of the temperature control unit. A dial setting of 304K corresponds to an actual probe probe temperature of 301.4 ± 0.3K. The spectrum of D 2 in 55wt per cent 1132/EBBA-d 2 has been measured as a function of temperature and the results of these measurements are given in table 16. In figure 23 the electric field gradient is plotted as a function of A P Q . It is clear, that to a good approximation, the electric field gradient also varies linearly with the liquid crystal splitting. It is interesting to note, however, that the best line through the points does not pass through the origin. Since both the electric field gradient and the liquid crystal splitting must go to zero at the nematic to isotropic phase transition temperature, the plot cannot be linear over the entire nematic range of the liquid crystal solution. It is also clear from the figure that Results and Discussion / 116 TABLE 15B Furan and Thiophene in 55wt per cent 1132/EBBA-d2: Liquid Crystal Order Parameters, Deuteron Quadrupolar Splittings and Dipolar Coupling Constants Sample Temp.* ** HD (Hz) s Sx Syy- (K) **Q (kHz) Furan 304 306 308 310 312 314 320 325 17.93 17.58 17.33 17.02 16.74 16.46 15.49 14.44 778 732 715 698 722 698 663 628 0.643 0.605 0.591 0.577 0.597 0.577 0.548 0.519 0.016 0.008 0.007 0.006 0.014 0.011 0.012 0.013 Thiophene 304 306 308 310 312 314 320 325 17.46 17.12 16.81 16.57 16.25 15.83 14.68 13.57 750 715 732 715 698 663 610 575 0.620 0.591 0.605 0.591 0.577 0.548 0.504 0.475 0.014 0.009 0.016 0.014 0.013 0.008 0.007 0.009 D f zz X * The temperature reported is the dial setting of the temperature control unit. A dial setting of 304K corresponds to an actual probe probe temperature of 301.4 ± 0.3K. ** &VQ is the quadrupolar splitting of the deuterons ortho to the nitrogen in the EBBA-a2 contained in the 55wt per cent 1132/EBBA-d2 mixture (see figure 5). t Dfjp is the dipolar coupling between the deuterons and the adjacent protons of the EBBA-d2 contained in the 55wt per cent 1132/EBBA-d2 mixture. $ S^, S and S are elements of the order parameter tensor of the aniline ring of the EBBA-d2 contained in the 55wt per cent 1132IEBBA-d2 mixture. yy zz Results and Discussion / 117 FIGURE 22 Force Constant in 55wt per cent 1132/EBBA-d2 1 (dyne < s 8 Force constant M / 7 ty fi 13 1 3 1 1 15 1 1 17 A i / (kHz) Q 77ie /brce constant describing the short range interactions between solvent and solute versus the quadrupolar splitting of the deuterons of the EBBA-d2 contained in the 55wt per cent 1132IEBBA-d2- The variation in both quantities is obtained by varying the temperature. The value of the force constant is that obtained from a least squares fit of the calculated order parameters to the experimental values for furan (circles) and thiophene (triangles). The slope and the intercept of the best straight line through the points are; 0.6049 dyne cm~l kHz'l; —2.13 dyne cm-1. ZZ °-° > 301.4K, correlation coefficient = 0.9817. F = esu T = Results and Discussion / 118 FIGURE 23 Electric Field Gradients in 55wt per cent 1132/EBBA-d2 14 15 16 AI/Q 17 18 (kHz) The measured electric field gradient versus the quadrupolar splitting of the deuterons of the EBBA-d2 contained in the 55 wt per cent 1132/EBBA-d2. The electric field gradient is calculated [10, 31, 32] from the spectrum of D2 dissolved in 55wt per cent 1132/EBBA-d2. The variation in both quantities is obtained by varying the temperature. The slope and intercept of the best straight line through the points are 5.84X109 esu kHz-1 and -I.OIXIOU esu; correlation coefficient = 0.970. Results and Discussion / 119 T A B L E 16 D2 in 55wt per cent 1132/EBBA-d2 DHD* T* B** (K) (Hz) (Hz) (esuXlOll) (kHz) (Hz) 304 306 308 310 312 314 320 325 -125.26 -127.98 -130.25 -132.21 -134.08 -135.79 -138.82 -139.59 5.08 5.20 5.29 5.39 5.49 5.56 5.70 5.77 0.023 0.008 0.012 -0.019 -0.057 -0.058 -0.083 -0.144 17.71 17.43 17.24 16.98 16.80 16.57 15.68 14.90 760 720 715 740 760 720 695 663 D*** F Z Z * The temperature reported is the dial setting of the temperature control unit. A dial setting of 304K corresponds to an actual probe temperature of 301.4K ** B is the quadrupolar coupling of the deuterons in D2 *** D is the dipolar coupling between the deuterons in D2 tAvQ is the quadrupolar splitting of the deuterons ortho to the nitrogen in the EBBA-d2 contained in the 55wt per cent 1132/EBBA-d2 mixture (see figure 21). $ D H D is the dipolar coupling between the deuterons and the adjacent protons of the EBBA-d2 contained in the 55wt per cent 1132IEBBA-d2 mixture (see figure 21). over the range of liquid crystal splittings measured for the other samples at 301.4K (table 13) the electric field gradient experienced by D Thus, if all solutes experience the same mean electric 2 is very small. field gradient as D2, mechanisms involving this field gradient may be safely ignored for our samples. However, the possibility remains that different distribution of local environments field gradient. each solute experiences and hence a different a somewhat mean electric Results and Discussion / 120 d. Summary We will now briefly summarize the results for solutes in 55wt per cent 1132/EBBA before going on to discuss the results for other liquid crystals. First, we have shown that the order parameters of solutes dissolved in this mixture are quite strongly correlated with the size and shape of the solutes. Second, we have shown that a simple model for the short-range interactions gives somewhat better agreement with the experimental results than the interaction between the polarizability and a mean square electric field, due to the liquid crystal. Both of these mechanisms are, at least approximately, functions of the size and shape of the solutes. Finally, we have shown that the force constant necessary to describe the solute order parameters using the short-range interaction model is not exactly the same for all solutes, and that this is not a result of variations in the experimental conditions. 4. O r d e r i n g o f solutes i n other l i q u i d crystals Up to this point we have only considered the ordering of solutes in 55wt per cent 1132/EBBA where F z z is measured to be zero. We will now turn our attention to results in other liquid crystals. In these solvents F ^ is no longer zero so that if our model for the short range interactions is to find any general application, we must first be able to calculate the contribution to solute ordering from the electric field gradient - molecular quadrupole moment mechanism. Unfortunately, this requires a knowledge of the molecular quadrupole moments of the solutes which are not easily determined [56]. It is clear from the range of values which have been reported for solutes such as benzene [120-122] and acetylene [123-126], that the accuracy of the literature values is not always very Results and Discussion / 121 good. Values for the molecular quadrupole moments of 15 of the solutes we have studied are shown in table 17. We will use these values throughout to calculate the contribution to solute ordering from the electric field gradient - molecular quadrupole moment mechanism. We begin by considering the ordering of solutes in the two component liquid crystals, 1132 and EBBA, of the zero field gradient mixture. a. 1132 and EBBA Experimental order parameters for the solutes shown in table 17 have measured in 1132 and EBBA-d2 and these results are shown in table been 18. In order to calculate order parameters for these solutes we require both the force constant, k, and the electric field gradient, F z z , for the two solvents. We shall arbitrarily assume that the value of k is the same as found for the 55wt per cent 1132/EBBA-d2 mixture and we shall use the value of F z z obtained from on D2 in 1132 and EBBA-d2 [31, 39]. If these assumptions are correct and if the exponential in expression (48) can be expanded in a Taylor experiments series and truncated after the linear term then the difference between the order parameters in either mixture represents 1132 or EBBA-d2 and the 55wt per cent the contribution to ordering from 1132/EBBA-d 2 the molecular quadrupole moment electric field gradient mechanism. Calculated versus experimental values for {S(EBBA-d ) 2 - S(55wt per cent 1132)} and {S(1132) - S(55wt per cent 1132)} are shown in figures 24 and 25. Considering that there are no adjustable parameters involved and that differences the agreement between the experimental in order parameters are being plotted, and calculated values is really quite good. As can be seen, all three components of the order parameter tensor are Results and Discussion / 122 T A B L E 17 Molecular Quadrupole Moments Solute Molecular quadrupole moment (10-26 2) a e 1,3,5-trifluorobenzene 1,3,5-trichlorobenzene 1,3,5-tribromobenzene Hexafluorobenzene 2,4-hexadiyne 2-butyne Allene c 2 c H2 HD D2 Acetylene 5.35 0.653 0.648 0.642 5.5 ± 2.5 Qxx Qyy Qzz Thiophene Reference c m 0.94 -3.24 -4.80 9.50 2.59b . 1 5 2.38b 1.62C 2.60b 1.56 7.35 4.64 ± 1.1 Bromomethane Propyne Fur an u -7.80 ± 2.2 Benzene Iodomethane Hydrogend s 4.82 5.9 -6.1 Qxx Qyy 0.2 6.6 -8.3 Qzz 1-7 average of values from [120] and [121] [120] [128] [128] [122] estimate [39] estimate [39] estimate [39] upper limit [127] average of values from [129, 130] [130] [55] average of values from [123, 124] [131] [132] [132] a) For benzene, bromomethane, and acetylene the value presented is an average of the largest and smallest values reported in the literature. b) Estimated from solute order parameter in 1132 c) Estimated from solute order parameter in d) Average for ground vibrational state EBBA-d2 Results and Discussion / 123 T A B L E 18 Solute Order Parameters: 1132 and EBBA-d2 at 301.4K a) b) c) d) e) f) g) h) i) j) k) 1) m) n) Solute 1132 1,3,5-tribromobenzene 1,3,5-trichlorobenzene Benzene Hydrogen Propyne Acetylene 2-butyne 2,4-hexadiyne Allene 1,3,5-trifluorobenzene Hexafluorobenzene Bromomethane Iodomethane Furan -0.2883 -0.2878 -0.2519 0.0082 0.2542 0.1912 0.2592 0.4270 0.1979 -0.2590 -0.2144 0.0947 0.0915 0.1727 -0.2063 0.0336 0.1118 -0.1990 0.0871 Sxx Syy Szz o) Thiophene S>xx Syy Szz a) For the solutes with C3 V predicted quite well EBBA-d2 -0.1848 -0.2038 -0.1157 -0.0111 0.0553 -0.0585 0.1524 0.3144 0.0925 -0.1744 -0.3144 0.0368 0.0280 0.0017 -0.0565 0.0548 -0.0082 -0.0737 0.0819 or higher symmetry Szz is given. for furan (labeled as circles) and thiophene (labeled as triangles) in both liquid crystals. All of the points (with the exception of the ones corresponding to 1,3,5-trifluorobenzene and positive (table where Qmol is reported to be small 17) he in the correct quadrant. The considerable scatter in the points is most likely a result of our assumption of the transferability of the force constant and electric field gradient. The experimental errors associated with molecular quadrupole moments, which are reflected points, also contribute to the scatter. in the error bars on the For three of the solutes, 2,4-hexadiyne, 2-butyne and allene, no values for the molecular quadrupole moments have been reported in the literature. However, from the experimental values of {S(1132) — Results and Discussion / 124 FIGURE 24 Solutes in EBBA-d2: Short Range Interaction Model -0.2 0.0 0.2 S (exp) Experimental versus calculated values of {S(EBBA-d2) — S(55 wt per cent 1132/EBBA-d2)} for furan (circles), thiophene (triangles), and solutes with C3 or higher symmetry. S(EBBA-d2) is the order parameter of the solute in EBBA-d2 and S(55 wt per cent 1132IEBBA-d2) is the order parameter of the solute in 55wt per cent 1132IEBBA-d2- For the solutes with C3 or higher symmetry only S z plotted. The labeling of the points refers to table 10. The error bars are meant to represent the maximum possible error from Qmol the calculated values, and have been calculated as follows. Where more than one literature value for Q l is available, the error bar represents the order parameters calculated using the limits of the literature values. Where only one literature value is available, the accuracy of Qmol is taken equal to that of a similar molecule. For example the benzene value has been used for all benzene derivatives. 55wt per cent 1132/EBBA-d2: T = 301.4K, k = 5.55 dyne cm-1, F = 0.0 esu. EBBA-d2 : T = 301 AK, k = 5.55 dyne cm-1, Fzz —5.75X1QH esu. correlation coefficient = 0.858 V I S V Z I n mo z = z Results and Discussion / 125 FIGURE 25 Solutes in 1132: Short Range Interaction Model -0.2 0.0 0.2 S (exp) Experimental versus calculated values of (S(1132) — S(55wt per cent 1132/EBBA-d2)} for furan (circles), thiophene (triangles) and solutes with C3 or higher symmtery. S(1132) is the order parameter of a given solute in 1132 and S(55wt per cent 1132/EBBA-d2) is the order parameter of a given solute in 55wt per cent 1132/EBBA-d2- The calculated and experimental values have been obtained and are labeled in the same manner as for figure 24. 55wt per cent 1132/EBBA-d2: T = 301.4K, k = 5.55 dyne cm-1, F = 0.0 esu. 1132: T=301.4K, k = 5.55 dyne cm-1, F = 6.07X10H esu, correlation coefficient = 0.833. V z z z z Results and Discussion / 126 S(55wt per cent 1132)} and {S(EBBA-d ) - S(55wt per cent 1132)} estimates of 2 the molecular quadrupole moment of each solute can be made and are presented in table 17. There is little variation between the two values obtained for each solute and the signs of the molecular quadrupole moments the molecular structure. The value obtained for allene are consistent is also magnitude and sign with the upper limit reported in reference pleasing to note that the negative axes of acetylene and with consistent in [127]. It is also order parameters observed for the symmetry hydrogen in EBBA-d 2 are correctly predicted without invoking special mechanisms for these two solutes. These results demonstrate very clearly the importance of the electric field gradient - molecular quadrupole moment mechanism for the ordering of solutes in these two liquid crystals. b. 5CB-a,B-d As a final test of our model for the short-range interactions we shall consider the 4 order of solutes dissolved j3-d4-n-pentyl)-4-cyanobiphenyl in (5CB-a,B-d^). the nematic liquid Values of crystal in this liquid crystal at various temperatures (from reference [40]) are presented in table measured dipolar and quadrupolar coupling constants splittings of the a- and ^-deuterons of the defined as T/T^_j where is TJJ_J the temperature of the solution. The values of as observed in EBBA and 1132 [31, quadrupolar (also from reference splitting of the liquid to 2 19 with the and the quadrupolar The temperature T isotropic phase r is transition &re about 30 per cent as large 39] and are the same sign as in Measured order parameters for 9 solutes with C g in 5CB-a,B-d^ of D 5CB-a,B-d^. nematic 4-(a, y 1132. or higher symmetry dissolved [40]) are shown in table 20 along with the crystal deuterons. These experimental order Results and Discussion / 127 T A B L E 19 D 2 Bt T (K) T * 275 285 295 300 305 307 0.8929 0.9253 0.9578 0.9740 0.9903 0.9968 (Hz) D* (Hz) F Z (lOH esu) At> (kHz) (kHz) 901.41±0.18 709.19±0.33 524.85±0.23 424.93±0.18 299.37±0.20 219.45±0.20 -38.92+0.13 -30.43±0.15 -22.8010.10 -18.64+0.08 -13.00+0.08 -9.70±0.09 3.60±0.13 2.52±0.24 2.3010.16 2.1510.13 1.3010.13 1.2310.14 67.21 62.65 56.25 51.09 42.39 34.23 50.78 45.70 39.33 34.91 28.05 22.27 r * T = TIT^j r in 5CB a Z where 7JJI is the nematic isotropic phase transition temperature. t B is the quadrupolar coupling of the deuterons. $ D is the dipolar coupling between the two deuterons T A B L E 20 Order Parameters of Solutes in 5CB at 294K* Solute a) b) c) d) e) f) g) h) i) l,3,5tribromobenzene l,3,5trichlorobenzene benzene hydrogen propyne acetylene 2-butyne 2,4-hexadiyne allene S(exp) S(calc)** S(calc)*** (kHz) Aiv/J (kHz) -0.2160 -0.2124 -0.2170 50.05 34.21 -0.2139 -0.1920 -0.2116 53.61 37.21 -0.1352 0.0033 0.1668 0.1686 0.1624 0.3545 0.1160 -0.1353 0.0175 0.1475 0.1429 0.1774 0.3684 0.1105 -0.1158 0.0206 0.1618 0.1712 0.1738 0.4445 0.1209 44.03 56.25 53.17 58.11 48.04 55.68 52.20 29.50 39.33 36.62 41.18 32.61 39.03 36.43 * All samples were measured at T=294K except benzene which was measured at T=-296K. ** k = 4.22 dynecm-1 , F zz = 2.30X1011 esu for allsolutes. *** k and F for each solute scaled using Aj> (see text). k cm-1 d kf for the D2 sample at 394K = 5.15 dyne cm'l. zz an a 0 = —1.26 dyne Results and Discussion / 128 parameters may then interactions and taking be predicted into using account our the model interaction for the between short the range molecular quadrupole moments of the solutes and the electric field gradients given in table 19. For 2,4-hexadiyne, 2-butyne and allene, averages of the estimates of the molecular quadrupole moments obtained from the solute order parameters in 1132 and EBBA-d2 have been used. The value of k has again been obtained by a least squares fit and is found to be 4.26 dyne cm *. The results of the fit are - shown in table 20 and are plotted in figure 26. Excellent agreement between the calculated and experimental values is obtained. It is interesting to note that the values of the molecular quadrupole moments for 2,4-hexadiyne, 2-butyne and allene needed to explain the order parameters of these solutes in EBBA-d2 and 1132, also predict their orientation in 5CB-a,/~-d4. In calculating the solute order parameters we have assumed that the are infinitely however, dilute and at do not satisfy precisely these the same criteria. temperature. Differences in samples The experiments, temperature and concentration among the samples are apparent in the variation of liquid crystal deuteron quadrupolar splittings (A»» a and Ap^ in table 20). these splittings to correct the values of k and F Z 2 We will now use for these differences. For all spectra of the samples reported here, the value of Ai> ^ obtained for a given value of A»> a is not affected by the presence of the solute. This implies that the solutes do not significantly affect conformational averaging or the value of 7j (expression may be dilution. 52) of the solvent used to make molecules corrections for and means that either the of temperature effects Ai> a or Ap^ and finite Results and Discussion / 130 A plot of the calculated field gradient (table 19) versus the measured quadrupolar splitting of the a-deuterons of the 5CB-a,/3-d 4 (figure 27) shows a roughly linear correlation. The best straight line through the points has a slope of 6.62X10^ esu kHz"* and an /3-deuterons has intercept a slope of of —1.28X10** 7.82X10^ esu esu. (A similar curve kHz"* and an for intercept the of — 0.70X10** esu.) There is considerable scatter in the results, primarily because the dipolar and quadrupolar splittings enough to give accurate values of F calibration curve for F We k z z z z > of D 2 cannot be measured accurately This graph, however, can be used as a . obtain a similar curve for the force constant by assuming that the value of corresponding to a given liquid crystal splitting, Ai> , a solutes. The dependence of k on Ai> a is the same for all for some reference solute may then be used as a calibration curve for k for all solutes. We have shown in figure 22 that the force constant obtained for furan and thiophene in 55wt per cent 1132/EBBA-d2 depends roughly linearly on the liquid crystal deuteron quadrupolar splittings over a small temperature range. Thus, it seems reasonable to expect that a similar behaviour will be found for 5CB-a,/3-d . The force constant kj for 4 a given sample is then: i= k where &v l a is the k o+ (3k/3Ai» ), quadrupolar splitting 5CB-a,^-d^ for sample i. The values of k the dependence of k on A v use a (54) a due Q to the and bk/dAv a a- deuterons in the will be obtained from for a reference sample. For this purpose we will 1,3,5-trichlorobenzene. A plot of k versus A»> a for this sample is shown in figure 28. From this plot we obtain a value of 0.1141 dyne cm' * kHz"* for Results and Discussion / 129 FIGURE 26 Solutes in 5CB-a,/3-d4: Short Range Interaction Model -0.2 0.0 S 0.2 0.4 (exp) Experimental versus calculated order parameters in 5CB-a,B~d^. The points are labeled in the same manner as in table 20. The force constant describing the short range interactions has been obtained from a least squares fit to the solid line of slope 1. T = 294K; k = 4.22 dyne cm-1; F = 2.305X10H esu; correlation coefficient = 0.997 z z Results and Discussion / 131 FIGURE 27 Electric Field Gradients in 5CB-a,^-d4 9 on fc 2 30 40 50 Av a 60 70 (kHz) The measured electric field gradient versus the quadrupolar splitting of the a-deuterons of the 5CB-a, {}-d4. The electric field gradient is measured from the spectrum of D2 dissolved in the 5CB-a, $-d4. The variation in both quantities is obtained by varying the temperature. The slope and intercept of the best straight line through the points are 6.64X10$ esu kHz'l and —1.28X10U esu; correlation coefficient = 0.941. Results and Discussion / 132 FIGURE 2 8 Force constant in 5CB-a,^-d4 7 i 30 40 50 Av a 60 70 (kHz) The force constant describing the short range interactions between solvent and solute versus the quadrupolar splitting of the a-deuterons of the 5CB-a,0-d4. The value of the force constant is that needed to explain the orientation of l,3,5-trichiorobenzene;dkJdhVa 0.1141 dyne cm-lkHz-1; k = —1.26 dyne cm-1; correlation coefficient — 0.992. = 0 Results and Discussion / 133 9k/9FAi> gives a and —1.26 dyne c m -1 for k 9 k/9Ai>0 = 0.1329 dyne c m evident from figure 28, -1 Q (a similar curve using the kHz -1 /3-deuterons and k =-0.19 dyne cm ). -1 Q As is our assumption that the variation of k with A f a is linear is valid over the range of the experimental results. The values of k- and fo r each sample may now be obtained and order parameters for the other solutes may be calculated. The values obtained from this calculation are shown in table 20 and are plotted against the experimental order parameters in figure 29. Again excellent agreement with experiment is obtained. However, it is important to point out that unlike the order parameters presented in figure 26 only the value for 1,3,5-trichlorobenzene has been obtained from a fit to the experimental value. The excellent agreement for the other solutes indicates very clearly that in this liquid crystal the force constant and electric field gradient may be transferred from one solute to another provided that the variation in experimental conditions is taken into account. It is also interesting to note that 2,4-hexadiyne, 2-butyne the values of the molecular quadrupole moments of and allene, obtained from the order parameters of the solvents in EBBA-d2 and 1132, predict their ordering in SCB-a^-d^ very well. The order parameters calculated using the B -deuteron splitting to scale F ^ and k differ from those reported in table 20 by at most 3.1 per cent. The values obtained for k- for the D t 5.15 and 5.02 dyne c m -1 solution at 295K using the a and B deuterons are 2 respectively. Hence the calculated order parameters do not depend on which quadrupolar splitting is used to calculate F ^ and k. To further test the above scaling method we compare the experimental and Results and Discussion / 134 F I G U R E 29 Solutes in 5CB-a,/3-d4: Short Range Interaction Model -0.2 0.0 0.2 0.4 S (exp) Calculated versus experimental order parameters for solutes dissolved in 5CB-a, &-d4. The calculated order parameters have been obtained in the same manner as for figure 26 except that values of F z and k have been obtained from figures 28 and 29 using the quadrupolar splittings of the a-deuterons of the 5CB-a,(}-d4 (see text). The labeling of the points refers to table 20. correlation coefficient = 0.994. Z Results and Discussion / 135 FIGURE 30 2,4-hexadiyne in 5CB-a,/J-d4 0.6 0.4 o 00 0.2 0.0 0.0 0.2 0.4 0.6 S (exp) Calculated versus experimental order parameters for 2,4- hexadiyne in 5CB-a, B-d.4. The calculated order parameters have been obtained in the same manner as for figure 26 except that values of F and k have been obtained from figures 28 and 29 using the quadrupolar splittings of the a-deuterons of the 5CB-a, B-d4 (see text), correlation coefficient = 0.991. circles = 1 mole per cent, squares = 9 mole per cent. z z Results and Discussion / 136 calculated order parameters temperatures and F z z [133] obtained for two 2,4-hexadiyne samples at various with those calculated using the values of kj from figure 27 from figure 28. The results of this calculation are presented in figure 30. As can be seen fairly good agreement is obtained between the calculated and experimental results. However, the points lie on a straight line whose slope is slightly different from one. This could be a result of several factors. First, our assumption that all solutes entirely valid as we have experience exactly the same environment is not seen in 55wt per cent 1132/EBBA-d2. Second, the calculated values depend on the value of the molecular quadrupole moment which has been estimated from the order parameters of the solute in 1132 and EBBA. Given the number of assumptions agreement have values we have made in estimating this value, the in figure 30 is surprisingly good. Finally, the linear correlations we assumed from figures of A<v and LvQ a for 28 and 29 may not be valid over the range of these samples. making relatively small correction to k and F good. In addition, as Nonetheless, z z can be seen from figure for the purposes of our method would appear to be 30 parameters for both dilute and concentrated samples. the method predicts order IV. CONCLUSIONS The main point of this thesis has been an attempt to describe the ordering of small solutes in nematic liquid crystals. As we have seen our approach of beginning with very simple solutes and studying their spectra in some detail has been quite successful. The discovery, from the spectrum of deuterium [10, 15, 31, 32, 36], of the presence of an external electric field has lead to a better understanding of the ordering of molecular hydrogen and other solutes. It is now clear that the interaction between the molecular quadrupole moment and this field gradient accounts for a considerable portion of the observed ordering. This interaction does not account for all of the subtle effects observed for hydrogen where the field gradient is small. It is involved in the ordering of this solute. apparent that other mechanisms are For highly symmetric solutes such as methane the coupling of the internal motions and the reorientations of the solute also plays a significant role in determining the observed dipolar couplings. These effects are also seen for acetylene, although they account for a much smaller portion It of dependence non-rigid the of observed the ratios contributions to couplings. of the the is also clear dipolar couplings couplings and not in that the liquid acetylene because of crystal results changes in from the average structure as had previously been assumed [21]. We have also seen that the size and shape of a solute governs its ordering in 55wt per cent 1132/EBBA where the external field gradient is zero. This suggests that the short range interactions are the dominant factor in determining solute ordering and a very simple model for these interactions successfully predicts the order parameters of a wide range of solutes in this mixture. However, because solute properties such as the polarizability are also functions of the 137 size and shape it is difficult to Conclusions / 138 distinguish clearly between mechanisms involving such properties and the size and shape. The predicting the short range anisotropy interaction in the model is somewhat more order matrix of solutes with successful C2 V in and &2h symmetry. This result gives support to the idea that it is indeed the size and shape of a solute which determines its ordering. Finally we have shown that in liquid crystals where the field gradient is not zero that the combination of the short range interaction model and the electric field gradient quadrupole moment mechanism predicts the order parameters quite well. From this work, then, a very clear picture has been gained of the general form of the interactions involved in the orientation of solutes in nematic liquid crystals. This should provide a basis for a description of the system on a molecular level. The size and shape model may also provide a method of explaining the ordering of the liquid crystal molecules themselves. V. APPENDIX I Summary of Results for H , HD and D in Several Liquid Crystals 2 Temperature EBBA 298.0 298.0 298.0 320.0 320.0 320.0 1132 298.0 298.0 298.0 320.0 320.0 320.0 Phase V 298.0 298.0 298.0 310.0 1167 310.0 310.0 310.0 330.0 330.0 23.9wt per 310.0 49.3wt per 310.0 61.3wt per 303.0 66.1wt per 310.0 76.0wt per 310.0 SCB-a^-d, 310.0 Solute H HD D H HD D 2 B (Hz) 2 2 -1817.20±0.50 -1943.60±0.50 2 2 -1556.70 + 0.40 -1654.10±0.40 -9.279X10-3 -1.064X10-2 -1.001X10-2 -0.450X10-2 -655.8610.06 -706.9810.05 -572.7010.60 -534.3010.30 1141.2010.30 185.6810.06 30.6310.02 24.9810.20 152.3810.30 -3.967X10-3 -0.420X10-2 -0.450X10-2 -0.367X10-2 -0.345X10-2 2 -788.2610.43 33.6710.14 -0.494X10-2 2 -267.6310.55 11.1510.18 -0.164X10-2 2 -13.7510.11 0.4510.04 -0.007X10-2 2 116.0610.42 -5.3810.14 0.079X10-2 2 373.0710.20 -16.4610.60 0.242X10-2 2 524.8510.23 -22.8010.10 0.335X10- 2 H D HD D 2 2 2 H HD D D HD 1132/EBBA D 1132/EBBA D 1132/EBBA D 1132/EBBA D 1132/EBBA D 2 2 2 cent D 8.134X10-3 0.864X10*2 0.929X106.244X100.662X10-2 0.705X10*2 2670.5010.30 72.5010.30 443.2010.90 30.6310.02 1351.80±0.60 1461.4910.30 2 cent -1.063X10-2 -1.144X10-2 -1.220X10-2 -9.151X10-3 -0.983X10-2 -1.046X10-2 -1695.3010.90 -1588.1010.90 -706.9810.05 2 cent 3059.3010.50 506.5010.50 83.1010.20 2632.0010.30 434.9010.40 71.2110.10 1029.0010.10 1106.8011.40 2 cent S(exp) -2340.9010.30 -382.5010.60 -63.2710.30 -1796.1010.30 -292.8011.60 -48.0010.80 H HD D H HD D cent D (Hz) 139 2 3 2 Appendix I / 140 Summary of Hydrogen Results cont. Temperature Solute B (Hz) 50wt per cent 1132/EBBA (sample 1) H 293.7 295.4 H 298.9 H 300.6 H 302.4 H 304.1 H 305.9 H 307.6 H 309.4 H 311.1 H H 312.9 314.6 H 316.4 H 318.2 H HD 290.1 -251.93±0.07 293.7 HD -254.20±0.07 295.4 HD -255.85±0.07 297.2 HD -257.96±0.07 HD 298.9 -259.21±0.07 HD -260.27±0.21 301.6 HD 303.3 -260.08±0.21 HD 305.0 -260.07±0.21 HD 306.7 -259.9810.21 D 290.1 -251.9310.07 D 293.7 -254.2010.07 295.4 -255.8510.07 D 297.2 -257.9610.07 D D 298.9 -259.2110.07 D -260.2710.21 301.6 303.3 -260.0810.21 D D -260.0710.21 305.0 -259.9810.21 306.7 D 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 D (Hz) S(exp) 448.7310.47 455.4810.47 457.8410.47 459.3610.47 459.7210.47 460.5410.47 459.9910.47 458.8210.47 456.9810.47 454.5510.47 451.4610.47 447.6710.47 443.4010.47 438.5110.47 65.3310.14 65.8710.14 66.4210.14 66.9610.14 67.3910.14 67.7110.14 68.0410.14 67.8110.14 67.7110.14 10.5010.02 10.6110.02 10.6610.02 10.7810.02 10.8110.02 10.9510.07 10.8610.07 10.8910.07 10.8910.07 -1.559X L0-3 -1.582X 10-3 -1.591X 10-3 -1.596X 10-3 -1.598X LO-3 -1.600X 10-3 -1.599X 10-3 -1.595X 10-3 -1.588X 10-3 -1.580X 10-3 -1.569X 10-3 -1.556X 10-3 -1.541X 10-3 -1.524X 10-3 -1.476X 10-3 -1.488X 10-3 -1.501 X 10-3 -1.513X 10-3 -1.523 X 10-3 -1.530X 10-3 -1.537X 10-3 -1.532X 10-3 -1.530X 10-3 -1.541X 10-3 -1.557X 10-3 -1.565X 10-3 -1.582X 10-3 -1.587X 10-3 -1.607 X 10-3 -1.594X 10-3 -1.599X 10-3 -1.599X 10-3 Appendix I / 141 Summary of Hydrogen Results cont. Temperature Solute B (Hz) 50wt per cent 1132/EBBA (sample 2) 293.6 H 295.4 H 298.9 H 300.6 H 302.4 H 304.1 H 305.9 H 307.6 H 309.4 H 311.1 H 312.9 H 314.6 H 316.4 H 318.2 H 293.7 HD -282.9110.35 295.4 HD -284.9410.49 298.9 HD -287.6910.30 300.7 HD -287.4410.35 302.4 HD -288.1410.35 304.2 HD 305.9 HD 307.7 HD 309.4 HD 311.2 HD 312.9 HD 314.7 HD 316.4 HD 318.2 HD 293.7 D -294.1715.00 295.4 D -298.3810.69 D 298.9 -301.3910.35 300.7 D -301.0410.35 302.4 D -302.3010.35 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 D (Hz) S(exp) 448.7310.47 455.4810.47 457.8410.47 459.3610.47 459.7210.47 460.5410.47 459.9910.47 458.8210.47 456.9810.47 454.5510.47 451.4610.47 447.6710.47 443.4010.47 438.5110.47 77.7012.00 78.8210.71 78.0610.71 78.7710.71 78.6210.71 79.6210.71 78.3910.71 78.0410.71 77.4710.71 77.1910.71 76.4210.71 76.0410.71 74.9510.71 74.1310.71 11.5311.65 12.2010.23 12.6010.20 12.5810.20 13.0910.12 -1.559X 10-3 -1.582X 10-3 -1.591X 10-3 -1.596X 10-3 -1.598X 10-3 -1.600X 10-3 -1.599X 10-3 -1.595X 10-3 -1.588X 10-3 -1.580X 10-3 -1.569X 10-3 -1.556X 10-3 -1.541X 10-3 -1.524X L0-3 -1.755X 10-3 -1.781X 10-3 -1.764X 10-3 -1.780X 10-3 -1.776X 10-3 -1.799X 10-3 -1.771X 10-3 -1.764X 10-3 -1.751X 10-3 -1.745X 10-3 -1.727X 10-3 -1.719X 10-3 -1.694X 10-3 -1.676X 10-3 -1.692 X 10-3 -1.791X 10-3 -1.850X 10-3 -1.847X 10-3 -1.922X 10-3 Appendix I / Summary of Hydrogen Results cont. Temperature 55wt per cent Solute H H 304.1 H H 309.4 311.1 312.9 D (Hz) S(exp) 2 200.3710.09 2 203.6910.09 - 6 . 9 6 3 X ]LO-4 -7.078X]LO-4 1132/EBBA 300.6 302.4 305.9 307.6 B (Hz) H H H 2 207.0610.09 209.6410.09 2 212.1710.09 2 214.5010.23 2 2 216.4310.23 217.8310.23 -7.455X]LO-4 - 7 . 5 2 3 X ]LO-4 - 7 . 5 7 2 X ]LO-4 218.3510.23 - 7 . 5 9 0 X ]LO-4 35.7810:14 36.2810.14 36.7310.14 - 8 . 0 8 4 X ]LO-4 - 8 . 1 9 8 X ]LO-4 - 8 . 3 0 0 X ]LO-4 37.1510.14 - 8 . 3 9 5 X ]LO-4 - 8 . 4 7 5 X ]LO-4 - 8 . 5 4 1 X ]LO-4 H H HD - 1 3 5 . 3 9 + 0.14 HD HD -137.22±0.14 - 1 3 8 . 6 8 + 0.14 305.9 307.7 309.4 HD HD - 1 3 9 . 9 4 + 0.14 -141.00±0.14 37.5010.14 HD -141.92±0.14 37.7910.14 311.2 HD -143.53±0.14 38.0210.14 312.9 314.7 300.7 302.4 304.2 HD HD D 38.2010.26 38.1510.14 2 -142.66±0.26 - 1 4 2 . 6 4 + 0.14 -138.02±0.29 D 2 -140.40±0.43 314.6 300.7 302.4 304.2 305.9 307.7 309.4 311.2 312.9 314.7 2 2 D D D D D D D 2 2 2 2 2 2 2 -142.10±0.30 -143.62±0.35 -145.10±0.14 -146.12±0.14 - 1 4 6 . 8 3 + 0.14 -146.8810.14 -146.83±0.14 - 7 . 1 9 6 X ]LO-4 - 7 . 2 8 6 X ]LO-4 -7.374X]LO-4 5.5910.14 5.7710.13 5.8810.10 5.9910.12 6.1210.05 6.1510.05 6.1610.05 6.1310.05 6.0810.05 - 8 . 5 9 3 X ]LO-4 - 8.634 X ]LO-4 - 8 . 6 2 3 X ]LO-4 - 8 . 2 0 6 X ]LO-4 - 8 . 4 7 0 X ]LO-4 - 8 . 6 3 2 X ]LO-4 - 8 . 7 9 4 X ]LO-4 - 8 . 9 8 6 X ]LO-4 - 9 . 0 3 0 X ]LO-4 - 9 . 0 4 5 X ]LO-4 - 9 . 0 0 2 X ]LO-4 - 8 . 9 2 9 X ]LO-4 142 Appendix I / 143 Summary of Hydrogen Results cont. Temperature Solute 60wt per cent 1132/EBBA 291.3 H 293.2 H 294.2 H 295.2 H 296.2 H 297.1 H 298.1 H 299.1 H 300.0 H 301.0 H 302.0 H 303.0 H 304.9 H 291.3 HD 293.3 HD 294.3 HD 295.2 HD 296.2 HD 297.2 HD 298.1 HD 299.1 HD 300.1 HD 301.1 HD 302.0 HD 303.0 HD 293.3 D 302.0 D 303.0 D B (Hz) 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 4.4210.21 1.5510.14 -3.4010.11 -3.9410.11 -5.5810.11 -8.1510.11 -9.8610.14 -11.5110.14 -14.2810.14 -16.1310.18 -18.2110.11 -20.6310.11 9.4610.11 -12.1010.11 -13.7510.11 D (Hz) S(exp) -42.6110.23 -34.2510.12 -29.7510.12 -27.4210.24 -24.0010.19 -19.3910.47 -16.0110.21 -13.5910.14 -7.3410.19 -2.7610.38 0.6010.23 3.9610.14 10.2510.16 -3.8210.28 -2.7910.14 -1.4410.11 -1.3010.11 -0.9310.11 -0.0410.11 0.4610.14 0.9110.14 1.6210.14 2.1010.18 2.8010.11 3.4910.11 -0.7210.04 0.2110.04 0.4510.04 1.480X10-4 1.190X10-4 1.034X10-4 0.953X10-4 0.834X10-4 0.674X10-4 0.556X10-4 0.472X10-4 0.255X10-4 0.096X10-4 -0.021X10-4 -0.138X10-4 -0.356X10-4 8.629X10-5 6.303X10-5 3.253X10-5 2.937X10-5 2.101X10-5 0.090X10-5 -1.039X10-5 -2.056X10-5 -3.660X10-5 -4.745X10-5 -6.327X10-5 -7.886X10-5 1.057X10-4 -0.308X10-4 -0.661X10-4 Appendix I / 144 Summary of Hydrogen Results cont. Temperature Solute 70wt per cent 1132/EBBA HD 293.7 295.4 HD 297.2 HD HD 305.9 HD 307.7 309.4 HD HD 311.2 HD 312.9 D 293.7 295.4 D D 297.2 D 305.9 307.7 D 309.4 D 311.2 D 312.9 D 2 2 2 2 2 2 2 2 The results for 1132, tables 2, 3 and 9. EBBA, B (Hz) 363.4810.07 352.0510.07 339.8610.07 280.8210.07 270.1210.07 258.3610.07 246.6710.07 235.8510.07 363.4810.07 352.0510.07 339.8610.07 280.8210.07 270.1210.07 258.3610.07 246.6710.07 235.8510.07 D (Hz) S(exp) -52.0710.14 -49.3710.14 -46.5010.14 -31.2610.14 -28.8910.14 -26.1410.14 -23.4410.14 -19.7510.14 -16.0010.02 -15.4810.02 -14.9810.02 -12.4010.02 -11.9710.02 -11.4510.02 -10.9610.02 -10.4410.02 1.176X10-3 1.115X10-3 1.051X10-3 0.706X10-3 0.653X10-3 0.591X10-3 0.530X10-3 0.446X10-3 2.348X10-3 2.272X10-3 2.199X10-3 1.821X10-3 1.757X10-3 1.681X10-3 1.609X10-3 1.533X10-3 Phase V and 1167 are from reference [10]. See also VI. APPENDIX II A. C A L C U L A T I O N OF MAXIMUM A N D MINIMUM CIRCUMFERENCES In this appendix we will describe the calculation of the minimum and maximum circumferences of the projection of a solute onto the XY plane. These circumferences are used in section D of chapter III in the model for the short range interactions. 1. Projection of solute In order to project the molecule onto the X Y plane we define the transformation from space fixed coordinates to molecule fixed coordinates in terms of two rotations of the axis system. The first rotation is in a clockwise direction about the space fixed Y axis through an angle 8. This is then followed by a clockwise rotation through an angle <f> about the resulting z axis. In terms of the angles 8 and 4> the coordinates of the projection of the i nucleus of the solute onto t n the X Y plane are given by: The Xj- = *j- cosfl cos<f> — yj cos0 sin$ + z± sin0 (55) Yj = x- s'm<j> + yj- cos0 . (56) projection of the molecule in the X Y plane is thus an array of overlapping circles centred at the points (Xj,Yj). The projected circle associated with the fa nucleus is given by: (X - Xj) 2 + (Y - 145 Yj) 2 = r 2 f . (57) Appendix II / 146 The maximum and minimum associated with the j values of X and Y on the projected circle nucleus are given by: Xmaxj = X ; + rj Xminj = Xj - r ; - r (59) £ Ymaxj = Yj + r Yminj = Y (58) f (60) ; (61) If Xmaxy and Ymax^ are the largest of the various Xmaxj and Ymaxj, and Xminl and Ymin^, are the smallest of the various Xmin-i and Ymin;i then the m • 7 points (Xmaxy,Yy), (Xmin^Yj), (X^Ymax^), and (X ,Ymin ) m lie m on the outer surface of the array of projected circles. 2. Maximum circumference The maximum circumference may then be calculated from the eqn.(57) by starting at (X-,YmaxO and moving around circle j points satisfying in small steps until the intersection with an adjacent circle, k, is found. Then starting at the intersection point we move around circle k in the same manner until the next intersection point is found and so on until the point (Xy,YmaxO is found again. The sum of the distances between adjacent points will be equal to the circumference in the limit of a small enough step. 3. Minimum circumference For the minimum circumference we begin by drawing line segments between all pairs of nuclei such that each involved and does not intersect line the segment is tangential internuclear vector at to the two circles a point which lies between the two nuclei. We then define the angle 6 between the Y axis and a Appendix II / 147 given line segment. Starting with the nucleus m which has the smallest Ymin we take the nuclei with Xmax- >^ Xmax m and calculate 6 for the line segments between each of these nuclei and nucleus m. We label the nucleus whose line segment has the largest 6 as m'. This line segment is then taken as a section of the projection of the outer surface of the molecule. Its contribution to the circumference is the distance between the intersection points of the line segment with the two circles plus the arc length of circle m between ( X , Ymin ) and m m the intersection point of the line segment with circle m. The procedure is then repeated starting at nucleus m' and we continue in this manner until the nucleus with the maximum Xmax is found. The contribution to the circumference between the other three pairs of outer points may be calculated in the manner to give the total circumference. same VII. APPENDIX III A. SOURCE C O D E FOR T H E C A L C U L A T I O N OF EXPRESSION (51) PROGRAM SIZE(INPUT, ^ OUTPUT); ************************************************** This program calculates the order parameters a set of molecules based on their size and shape. The order parameter i s the c a l s s i c a l average weighted using the potential calculated by the function POTEN which uses the minimum circumference of the projection of the molecule i n the XY plane. **************************************************^ CONST MAXATOM = 20; MAXMOL = 2 0 ; K = 1.38044E-16; PI = 3.141592645; TYPE ATOMRECORD = RECORD P, Q, R, RADIUS, X, Y: REAL; RIJ: ARRAY [ l . . MAXATOM] OF REAL; END; MOLRECQRD = RECORD NAME: PACKED ARRAY [ l . . 30] OF CHAR; NATOM, XMAXNUC, YMAXNUC, XMINNUC, YMINNUC: INTEGER; FKON, EFG, QXX, QYY, QZZ: REAL; DIFF, SEXP, SXX, SYY, SZZ: ARRAY [ l . . 3] OF REAL; ATOM: ARRAY [ l . . MAXATOM] OF ATOMRECORD; END; MOLARRAY = ARRAY [ l . . MAXMOL] OF MOLRECORD; VAR MOLECULE: MOLARRAY; MOLPT, NMOL: INTEGER; FCON, FZZ, T: REAL; LXTAL: PACKED ARRAY [1 OUTCODE: (REED, RITE); CIRCFILE: FILE OF REAL; FUNCTION DASIN(VAR FUNCTION DATAN2(VAR . . 30] OF CHAR; X:REAL):REAL; FORTRAN; X, Y:REAL):REAL; FORTRAN; PROCEDURE READINPUT; ! 148 Appendix III / 149 VAR I, MOLPT, ATOMPT: INTEGER; CODE: CHAR; BEGIN PAGE(OUTPUT); WRITELNC INPUT '); WRITELN (' + '); WRITELN; READLN(CODE); CASE CODE OF 'R' : OUTCODE := REED; •W : OUTCODE := RITE; END(*CASE*); FOR I := 1 TO 30 DO LXTAL [ I ] := '; I := 0; WHILE NOT EOLN(INPUT) DO BEGIN I := I + 1; READ(LXTAL[I]); END (*WHILE*); READLN(NMOL, T ) ; WRITELN(' 'LIQUID CRYSTAL: ' , LXTAL); WRITELNC ' , 'TEMPERATURE : ' , T : 7: 1); WRITELN; WRITELN(' MOLECULE','COORDINATES':26,'RADIUS':19, 'QUADRUPOLE MOMENT *: 26); WRITELN('(XE-8 CM)':34,'(XE-8 CM)':22, '(XE-26 ESU CM2)':23); WRITELN('X' : 23, 'Y' : 8, ' Z ' : 8, 'XX': 30, ' Y Y ' : 6, ' Z Z ' : 8 ) ; WRITELN(' + ' , ' ' :42, ' ' :42); FOR MOLPT := 1 TO NMOL DO" WITH MOLECULE[MOLPT] DO BEGIN FOR I := 1 TO 30 DO NAME[I] := ' ' ; I : 0• WHILE NOT EOLN(INPUT) DO BEGIN I := I + 1; READ(NAME[I]) END (*WHILE*); READLN(FKON, EFG); WRITELNC ' , 'FIELD GRADIENT:', EFG: 12); WRITELNC ' , "FORCE CONSTANT:', FKON: 12); READ(SEXP[1J,SEXP[2],SEXP[3], NATOM, QXX, QYY, QZZ); WRITELN; WRITELNC ' , NAME, QXX: 37: 2, QYY: 6: 2, QZZ: 6: 2); QXX := QXX * 1.0E-26; QYY := QYY * 1.0E-26; QZZ := QZZ * 1.0E-26; FOR ATOMPT := 1 TO NATOM DO WITH ATOM[ATOMPT] DO = Appendix III / 150 BEGIN READLN(P, Q, R, RADIUS); WRITELNC ' , P: 24: 4, Q: 8: R: 8: 4, RADIUS:12: END (*WITH*); END (*WITH*); END (*READINPUT*); 4, 2); PROCEDURE PROJECTNUCLEKTHETA, PHI: REAL); ^************************************************** This procedure projects the coordinates of the solute onto the XY plane as a function of i t s orientation which are given by the angles 8 and t/>. ******************************************** VAR ATOMPT: INTEGER; BEGIN WITH MOLECULE[MOLPT] DO BEGIN FOR ATOMPT := 1 TO NATOM DO WITH ATOM[ATOMPT ] DO BEGIN Y := P * SIN(PHI) + Q * COS(PHI); X := P * COS(THETA) * COS(PHI) Q * COS(THETA) * SIN(PHI) + R * SIN(THETA); END (*WITH*); END (*WITH*); END (*PROJECTNUCLEI*) ; PROCEDURE OUTERNUCLEI; ^******************************************* This procedure finds the points with the maximum and mimimum X and Y values for the projection of a given molecule i n the XY plane ******************************************* ^ VAR ATOMPT: INTEGER; YMAX, YMIN, XMAX, XMIN: REAL; BEGIN WITH MOLECULE[MOLPT] DO BEGIN YMAX := 0.0; YMIN := 0.0; XMAX := 0.0; XMIN := 0.0; FOR ATOMPT := 1 TO NATOM DO WITH ATOM[ATOMPT] DO Appendix III / 151 END BEGIN IF (Y + RADIUS - YMAX) > BEGIN YMAX := Y + RADIUS; YMAXNUC := ATOMPT; END (*IF*); IF (Y - RADIUS - YMIN) < BEGIN YMIN := Y - RADIUS; YMINNUC := ATOMPT END (*IF*); IF (X + RADIUS - XMAX) > BEGIN XMAX := X + RADIUS; XMAXNUC := ATOMPT; END (*IF*); IF (X - RADIUS - XMIN) < BEGIN XMIN := X - RADIUS; XMINNUC := ATOMPT END (*IF*); END (*WITH*); END (*WITH*); (*OUTERNUCLEI*); 1.0E-5 THEN - l.OE-5 THEN l.OE-5 THEN - l.OE-5 THEN PROCEDURE ROTATEAXES; ^***************************************************** This procedure rotates the space fixed axis system through 90 degrees about Z so that the procedure GETCIRCUM need only calculate the section of the circumference from YMIN to XMAX after each of four 90 degree r o t a t i o n s . ******************************************************\ VAR ATOMPT, SAVINT: INTEGER; XOLD: REAL; BEGIN WITH MOLECULE[MOLPT] DO BEGIN FOR ATOMPT := 1 TO NATOM DO WITH ATOM [ATOMPT] DO BEGIN XOLD := X; X := Y; Y := - XOLD; END (*WITH*); SAVINT := XMAXNUC; XMAXNUC := YMAXNUC; YMAXNUC := XMINNUC; XMINNUC := YMINNUC; YMINNUC := SAVINT; END (*WITH*); END (*ROTATEAXES*); Appendix III / 152 PROCEDURE GETRIJ; (*************************************************** This procedure calculates the distances between pairs of n u c l e i for the projection of the solute i n the XY plane. These distances are used i n calculating the circumference. ******************************************** VAR XI, YI: REAL; I, J : INTEGER; BEGIN WITH MOLECULE[MOLPT] DO BEGIN FOR I := 1 TO NATOM DO BEGIN WITH AT0M[I] DO BEGIN XI := X; YI := Y END (*WITH*); FOR J := 1 TO NATOM DO WITH ATOM[j] DO RIJ[I] := SQRT((XI - X) * (XI - X) + (YI - Y) * (YI - Y ) ) ; END (*FOR*); END (*WITH*); END (*GETRIJ*); PROCEDURE GETCIRCUM(VAR SEG: REAL); ^ *************************************************** This procedure calculates the section of the circumference between YMIN and XMAX ****************************************************\ VAR ATOMPT, If J: XI, YI, RI, XMAXJ, INTEGER; XMAX, XMAXI, ALPHAJ, ALPHA, GAMMA, DELTAJ, DELTA, COSGAM, SINGAM, BETA: REAL; BEGIN CASE OUTCODE OF Appendix III / 153 REED: READ(CIRCFILE,SEG); RITE: BEGIN WITH MOLECULE[MOLPT] DO BEGIN J := 0; SEG := 0.0; BETA := 0.0; I := YMINNUC; WITH ATOM[XMAXNUC] DO XMAX := X + RADIUS; WHILE I O XMAXNUC DO BEGIN DELTA := 0.0; WITH AT0M[I] DO BEGIN XI := X; YI := Y; RI := RADIUS; XMAXI := X + RADIUS; END (*WITH*); IF ABS (XMAX - XMAXI) < 1E-5 THEN BEGIN WITH ATOM[XMAXNUC] DO IF RIJ[I] O 0.0 THEN BEGIN ALPHA := (RADIUS - RI) / R I J [ l ] ; ALPHA := DASIN(ALPHA); SINGAM := (X - XI) / R l J f l l ; COSGAM := (Y - YI) / R I J [ I J ; GAMMA := DATAN2(SINGAM,COSGAM); DELTA := ALPHAJ + GAMMA; J := XMAXNUC; END (*IF*) ELSE J := XMAXNUC; END (*IF*) ELSE BEGIN FOR ATOMPT := 1 TO NATOM DO WITH ATOM[ATOMPT] DO BEGIN XMAXJ := X + RADIUS; IF ((XMAXJ - XMAXI) > - 1.0E-5) AND (ATOMPTOI) AND (RIJ[I] o 0.0) THEN BEGIN ALPHAJ :=(RADIUS-RI)/RIJ[I]; ALPHAJ := DASIN(ALPHAJ); SINGAM := (X - XI) / R I J [ l ] ; COSGAM := (Y - YI) / R I J [ l ] ; GAMMA:= DATAN2(SINGAM,COSGAM); DELTAJ := ALPHAJ + GAMMA; IF (DELTAJ - DELTA) > -1.0E-5 THEN BEGIN DELTA := DELTAJ; ALPHA := ALPHAJ; Appendix III / 154 J := ATOMPT; END (*IF*); END (*IF*); END (*WITH*); END (*ELSE*); WITH AT0M[J] DO BEGIN SEG := SEG + RI * ((PI / 2.0 - DELTA) - BETA); BETA := PI / 2.0 - DELTA; SEG := SEG + R I J [ l ] * COS(ALPHA); END (*WITH*); I 5= J * END (*WHILE*); SEG := SEG + ATOM[ I ] .RADIUS * (PI / 2.0 - BETA); WRITE(CIRCFILE,SEG); END (*WITH*); END(*RITE*); END (*CASE*); END (*GETCIRCUM*); FUNCTION POTEN(THETA, PHI: REAL): REAL; ^******************************************** This function calculates the p o t e n t i a l energy of the molecule as a function of the angles 8 and 4>. The p o t e n t i l a has the form U = ~i K C Where C i s the circumference of the projection of the molecule i n XY plane 2 *******************************************\ VAR C, SEG, PXX, PYY, PZZ: REAL; J : INTEGER; BEGIN (*P0TEN*) WITH MOLECULE[MOLPT] DO BEGIN C : 0.0" IF OUTCODE = RITE THEN BEGIN PROJECTNUCLEI(THETA, PHI) ; GETRIJ; OUTERNUCLEI; END(*IF*); GETCIRCUM(SEG); C := C + SEG; FOR J := 1 TO 3 DO BEGIN IF OUTCODE = RITE THEN = Appendix III / ROTATEAXES; GETCIRCUM(SEG); C := C + SEG; END (*FOR*); PXX := (3.0 * SIN(THETA) * SIN(THETA) * COS(PHI) * COS(PHI) - 1.0) / 2.0; PYY : = (3.0 * SIN(THETA) * SIN (THETA) * SIN(PHI) * SIN(PHI) - 1.0) / 2.0; PZZ := (3.0 * COS (THETA) * COS (THETA) - 1.0) / 2.0; U := (FCON * C * C - FZZ * (QXX * PXX + QYY * PYY + PZZ * QZZ)) / (2.0 * K * T ) ; U := EXP(- U) * SIN(THETA); END (*WITH*); POTEN := U; END (*POTEN*); PROCEDURE INTEGRAL(FUNCTION F(X, Y: REAL): REAL; L, U, L I , U l : REAL); /************************************************** This procedure calculates the order parameters of the molecule by integrating the function F(X,Y) from X=L to X=U and Y=L1 to Y=U1. INTEGRAL i s c a l l e d with: F(X,Y) = POTEN(THETA,PHI) ******************************************** VAR Fl, F2, F3, F4, INC,INC1, THETA, THETA1, SUM, PHI1, PHI: REAL; I, Ji NLOOPS1, NLOOPS: INTEGER; BEGIN WITH MOLECULE [MOLPT] DO BEGIN FOR I := 1 TO 3 DO BEGIN SZZ[I] := 0.0; SYY[I] := 0.0; SXX[I] := 0.0; END (*FOR*); SUM := 0.0; F2 := 0.0; NLOOPS1 :=36; NLOOPS := 36; INC := (U - L) / NLOOPS; INC1 := ( U l - L l ) /NLOOPS1; 155 Appendix III / 156 THETA := L ; PHI := L I ; F4 := F(THETA, PHI); FOR I := 1 TO NLOOPS1 DO BEGIN F3 := F4; PHI1 := PHI; PHI := I * INC1 + L I ; F4 := F(THETA, PHI); F2 := F2 + (F3 + F4) * INC1 / 2.0; SXX[2] := SXX[2] + (COS(PHI) * COS(PHI) * F4 + COS(PHIl) * COS(PHIl) * F3) * INC1 / 2.0; SYY[2] := SYY[2] + (SIN(PHI) * SIN(PHI) * F4 + SIN(PHIl) * SIN(PHIl) * F3) * INC1 / 2.0; END (*FOR*); SZZ[2] := F2 * (1 / 2) * (3 * COS(THETA) * COS(THETA) - i ) ; SYY [2] := SYY[2] * (3 / 2) * SIN(THETA) * SIN(THETA); SXX [2] := SXX[2J * (3 / 2) * SIN(THETA) * SIN(THETA); FOR I = 1 TO NLOOPS DO BEGIN • F2; SZZ 1' := SZZ 2' SYY 1' := SYY '2' SXX ' l ' := SXX '2' THETA1 := THETA; THETA := INC * I + L ; F2 := 0.0; SZZ[2] := 0.0; SYY[2l := 0.0; SXX[2] := 0.0; PHI := L I ; F4 := F( THETA, PHI); FOR J := 1 TO NLOOPS 1 DO BEGIN F3 := F4; PHI1 := PHI; PHI := J * INC1+ L I ; F4 := F(THETA, PHI); F2 := F2 + (F3 + F4) * INC1 / 2.0; SXX[2] := SXX[2] + (COS(PHI) COS(PHI) * F4 + COS(PHIl) * COS(PHIl) * F3) * INC1 / 2.0; SYY[2] := SYY[2] + (SIN(PHI) * SIN(PHI) * F4 + SIN(PHIl) * SIN(PHIl) * F3) * INC1 / 2.0; END (*FOR*); SZZ[2] := F2 * (1 / 2) * (3 * COS (THETA) COS(THETA) - 1); SYY [21 := SYY SYY[2] (3 / 2) * SIN(THETA) * SIN(THETA); SXX[2J SXX [2] := SXX (3 / 2) * SIN(THETA) * SIN(THETA); SUM : = SUM + (FI F2) * INC / 2, := SXX[3] SXX (SXX[1'] + SXX 2' ) * INC / 2.0; := SYY" SYY (SYY[1" + SYY '2' ) * INC / 2.0; := SZZ ] SZZ (SZZ[1 J + SZZ '2' ) * INC / 2.0; END (*FOR*); SXX[3] := SXX[3] / SUM - 1 / 2; Appendix III / END SYY[3] := SYY[3] / SZZ[3] := SZZ[3] / WRITELN; WRITELN('SXX = ' : WRITELN('SYY = ' : WRITELNCSZZ = ' : END (*WITH*); (*INTEGRAL*); SUM SUM; 20, 20, 20, - 1 / 2 ; SXX[3l); SYY[3J); SZZ[3]); PROCEDURE PRINTOUT; ^ ************************************* This procedure writes the output **************************************^ VAR If MOLPT: INTEGER; BEGIN PAGE(OUTPUT); WRITELN(' OUTPUT '); WRITELN (' + ') ; WRITELN; WRITELN(' ' , 'LIQUID CRYSTAL: ' , LXTAL); WRITELN(' ' , "FIELD GRADIENT:', FZZ: 12); WRITELNC ' , 'FORCE CONSTANT:', FCON: 12); WRITELNC *, 'TEMPERATURE : ' , T: 7: 1); WRITELN; WRITELN(' ' , 'MOLECULE': 8, 'S(CALC)': 33, 'S(EXP)': 8, 'PERCENT ERROR': 17); WRITELN; FOR MOLPT := 1 TO NMOL DO WITH MOLECULE[MOLPT] DO BEGIN WRITEC NAME); WRITE('SXX',SXX[3]:8:4,SEXP[l]:8:4,DIFF[l]:14:3); WRITELN; WRITELN('SYY*:34,SYY[3]:8:4,SEXP[2]:8:4,DIFF[2]:14:3); WRITELN('SZZ':34,SZZ[3]:8:4,SEXP[3]:8:4,DIFF[3]:14:3); WRITELN; END (*WITH*); END (*PRINTOUT*); (* This i s the main body of the program*) BEGIN RESET(INPUT,'UNIT=SCARDS'); REWRITE(OUTPUT,'UNIT=SPRINT,NOCC*); READINPUT; CASE OUTCODE OF REED: RESET(CIRCFILE,'UNIT=1,BLOCK'); RITE: REWRITE(CIRCFILE,'UNIT=1,BLOCK'); END(*CASE*); PAGE(OUTPUT); FOR MOLPT := 1 TO NMOL DO WITH MOLECULE [MOLPT] DO (*SIZE*) 157 Appendix III / BEGIN FZZ := EFG; FCON := FKON; WRITEC ',NAME); WRITELN; INTEGRAL(POTEN, 0.0, PI / 2.0 ,0,PI); SEXP DIFF [ 1 ] = ABS((SXX 3 ] / SEXP - SEXP ?.) / SEXP = ABS((SYY DIFF SEXP = ABS((SZZ DIFF 11} / SEXP END (*WITH*); PRINTOUT; CLOSE(CIRCFILE); END (*SIZE*). I]: 100.0; 100.0; 100.0; 158 VIII. BIBLIOGRAPHY 1. A. Saupe, Z. Naturforsch. 1 9 A , 161 (1964) 2. P. Diehl and C.L. Khetrapal, NMR Basic Principles and Progress, (1969) 3. J.W. Emsley and J.C. Lindon, "N.M.R. Spectroscopy Solvents" Permagon Press, Oxford (1975) 4. C.L. Khetrapal, A.C. Kunwar, Adv. Liq. Cryst. 6, 173 (1983) 5. J . Lounila and J . Jokisaari, Progr. NMR Spectrosc. 6. D. Catalano, C. Forte, C.A. Veracini and C.Zannoni, Israel J . Chem. 23, 283 (1983) 7. T. Vaananen, J . Jokisaari and J . Lounila, J . Mag. Reson. 54, 436 (1983) 8. J . Jokisaari, T. Vaananen and J . Lounila, Mol. Phys. 4 5 , 141 (1982) 9. J.W. Emsley, R. Hashim, G.R. Luckhurst and G.N. Shilstone, 1, 437 (1986) 10. E . E . Burnell, C.A. de Lange and J . G . Snijders, Phys. Rev. A 25, 2339 (1982) 11. J . M . Anderson, J . Mag. Reson., 4 , 231 (1971) 12. E.T. Samulski, Ferroelectrics 30, 83 (1980) 13. J.C. Robertson, (1971) 14. J . G . Snijders, C.A. de Lange and E . E . Burnell, J . Chem. Phys. 7 7 , 5386 (1982) 15. J . G . Snijders, C.A. de Lange and E . E . Burnell, J . Chem. Phys. 7 9 , 2964 (1983) 16. J . G . Snijders, C.A. de Lange and E . E . Burnell, Israel J . Chem. 23, 269 C T . Yim and D.F.R. 1, 1 Using Liquid Crystal 1 5 , 249 (1982) Liq. Cryst., Gilson, Can. J . Chem. 4 9 , 2345 (1983) 17. A. Saupe, Mol. Cryst. 1, 527 (1966) 18. J . Nehring and A. Saupe, Mol. Cryst. Liq. Cryst. 8, 403 (1969) 19. J.W. Emsley, R. Hashim, G.R. Luckhurst, G. Rumbles and F.R. Viloria, Mol. Phys. 4 9 , 1321 (1983) 159 Bibliography / 160 20. R. Hashim, G.R. Luckhurst and S. Romano, Mol. Phys. 56, 1217 (1985) 21. P. Diehl, S. Sykora, W. Niederberger 14, 260 (1974) 22. E . E . Burnell, A . J . van der Est, G.N. Patey, Snijders, Bull. Mag. Res. (in press) (1987) 23. G.R. Luckhurst in "Nuclear Magnetic Resonance of Liquid Crystals" J.W. Emsley ed. (D. Reidel Publ. Co. 1983) 24. F. Reinitzer, Monatsh. Chem. 9, 421 (1888) 25. O. Lehmann, Z. Phys. Chem. 4, 462 (1889) 26. G. Friedel, Ann. de Chem. 18, 273 (1922) 27. G.W. Gray and P.A. Winsor, "Liquid Crystals and Plastic Crystals", 1, 1 (1974) 28. P. Diehl, H. Kellerhals Progress", 6, 1, (1972) 29. J . Lounila and P. Diehl, J . Mag. Reson. 56, 254 (1984) 30. P. Diehl, P.M. Henrichs, W. Niederberger, Mol. Phys. 20, 139 (1971) 31. G.N. Patey, E . E . Burnell, J . G . Snijders and C.A. de Lange, Chem. Phys. Lett. 99, 271 (1983) 32. P.B. Barker, A . J . van der Est, E . E . Burnell, G.N. Patey, C.A. de Lange and J . G . Snijders, Chem. Phys. Lett. 107, 426 (1984) 33. E . E . Burnell and C.A. de Lange, J . Chem. Phys. 76, 3474 (1982) 34. J . Lounila, Mol. Phys. 58, 897 (1986) 35. D.L. Gray and A.G. Robiette, Mol. Phys. Rev. 37, 1901 (1979) 36. A . J . van der Est, P.B. Barker, E . E . Burnell, C.A. de Lange Snijders, Mol. Phys. 56 161, (1985) 37. E . E . Burnell, C.A. de Lange, Chem. Phys. Lett., 76, 268 (1980) 38. A . J . van der Est, M.Y. Kok and E . E . Burnell, Mol. Phys. (1987) and E . E . Burnell, J . Mag. Reson. and E . Lustig, C.A. de "NMR Basic Lange, J.G. Principles and and J . G . 60, 397 Bibliography / 161 39. M.Y. Kok, A . J . van der Est and E . E . Burnell, Liq. Cryst. (submitted) (1987) 40. A. Weaver, A . J . van der Est, J.C.T. Rendell, G.L. Hoatson, and E . E . Burnell, Liq. Cryst. (in press) (1987) 41. G.S. Bates, E . E . Burnell, G.L. Hoatson, Weaver, Chem. Phys. Lett. 134, 161 (1987) 42. J . L . Charlton, B. Kadai and H.K. Lai, Can. J . Chem., 52, 3577 (1974). 43. P. Keller and L. Liebert, Solid State Phys. Suppl. 1 4 , 19 (1978) 44. Kodak puplication no. JJ-14, Eastman Liquid Crystal Products (1973) 45. J.W. Mellor and G.D. Parkes, "Mellor's Longman's, Green and Co. (1939) 46. W. Pauli, Phys. Rev. 58, 716 (1940) 47. N.F. Ramsey, Phys. Rev. 91, 303 (1953) 48. N.F. Ramsey, "Molecular Beams" Oxford University Press, London (1956) 49. R.F. Code and N.F. Ramsey, Phys. Rev. A4, 1945 (1971) 50. C A . de Lange, J . G . Snijders and E . E . Burnell, in "Nuclear Magnetic Resonance of Liquid Crystals" J.W. Emsley (ed.) D. Reidel Publishing Co. (1985) 51. A.D. Buckingham, J . Chem. Phys. 36, 3096 (1962) 52. B.P. Stoicheff, Can. J . Phys. 35, 730 (1957) 53. G.F. Newell, Phys. Rev. 78, 711 (1950) 54. R.V. Reid, Jr. and M.L. Vaida, Phys. Rev. A 7, 1841 (1973) 55. J . L . Hunt, J.D. Poll and L. Wolniewicz, Can. J . Phys. 62, 1719 (1984) 56. A.D. Buckingham, "Electric Moments of Molecules" in "Physical Chemistry: An Advanced Treatise" H . Eyring, D. Henderson and W. Jost eds. Vol. IV, pg 349 Academic Press (1970) 57. W. Maier and A. Saupe, Z. Naturforsch. 14 A, 882 (1959) 58. P. Diehl, (1979) H . Bosiger P. Palffy-Muhoray Modern and H . Zimmermann, G.S. Bates Inorganic J . Mag. Reson. and A. Chemistry" 33, 113 Bibliography / 162 59. R. Ader and A. Loewenstein, Mol. Phys. 24, 455 (1972) 60. R. Ader and A. Loewenstein, Mol. Phys. 30, 199 (1976) 61. A. Loewenstein, Chem. Phys. Lett. 38, 543 (1976) 62. A. Amanzi, P.L. Barili, Lett 44, 110 (1976) 63. A. Frey and R.R. Ernst, Chem. Phys. Lett. 49, 75 (1977) 64. P. Diehl, A.C. Kunwar and H. Zimmermann, J . Mag. Reson. (1978) 65. I.Y. Wei and C.S. Johnson, J . Mag. Reson. 23, 259 (1976) 66. J . G . Snijders, W. van der Meer, E . J . Baerends and C.A. de Lange, J . Chem. Phys., 79, 2970 (1983) 67. M . Gussoni, Adv. in Infrared and Raman Spec. 6, 61 (1979) 68. C. Eckart, Phys. Rev. 47, 552 (1935) 69. B. Crawford, J . Chem. Phys. 20, 977 (1952) 70. M . Gussoni and S. Abbate, J . Mol. Spectrosc. 62, 53 (1976) 71. G. Strey and L M . Mills, J . Mol. Spectrosc. 39, 103 (1976) 72. A.D. Buckingham and I. Love, J . Mag. Reson. 2, 338 (1970) 73. H . Nakatsuji, H. Kato, I. Morishima and T. Yonezawa, Chem. Phys. Lett. 4, 607 (1970) 74. L M . Mills, J . Phys. Chem. 80, 1187 (1976) 75. A. Bondi, J . Phys. Chem. 68, 441, (1964) 76. G.L.D. Ritchie and J . Vrbancich, Austr. J . Chem. 35, 869 (1982) 77. Structure assumed the same as has been taken from [78] 78. A . Allmennigen, I. Hargittai, J . Brunvoll, A. Domenicano J . Molec. Struct. 116, 119 (1984) 79. G.R. Alms, A.K. Burnham and W.H. Flygare, J . Chem. Phys. 63, 3321 (1975) P. Chidichimo and C.A. Veracini, Chem. Phys. 30, 621 1,3,5-trichlorobenzene; the C-Br distance and S. Samdal, Bibliography / 163 80. M.D. Harmony, V . M . Laurie, R.L. Kuczowski, R.H. Schwendeman, D.A. Ramsay, F . J . Lovas, W.J. Lafferty and A.G. Maki, J . Chem. Phys. Ref. Data 8, 619 (1979) 81. K.B. MacAdam and N.F. Ramsey, Phys. Rev. A6, 898 (1972) 82. G.A. Victor and A. Dalgarno, J . Chem. Phys. 50, 2535 (1969) 83. T.C. Wong, (1977) 84. G. Englert, A. Saupe and J.P. Weber, Z. Naturforsch. A 23, 152 (1968) 85. G. Nagarajan, Acta Phys. Pol. 30, 743 (1966) 86. E . Sackmann, J . Chem. Phys. 51, 2984 (1969) 87. M.P. Bogaard, A.D. Buckingham, R.K. Pierens and A . H . White, J . Chem. Soc, Farad. Trans. I, 74, 3008 (1978) 88. C T . Yim and D.F.R. Gilson, Can. J . Chem. 46, 2783 (1968) 89. A. Allmennigen, O. Bastiasen, R. Siep and H . Siep, Acta. Chem. Scand. 18, 2115 (1964) 90. A.K. Burnham, T.D. Gierke, J . Chem. Phys. 73, 4822 (1980) 91. M.S. Gopinathan and P.T Narasimhan, J . Magn. Res. 6, 147 (1972) 92. H.J.M. Bowen, J . Donohue, D.G. Jenkin, O. Kennard, P.J. Wheatly and D.H. Whiffin, "Tables of Interatomic Distances and Configuration in Molecules and Ions", L . E . Sutton, D.G. Jenkin, A.D. Mitchell and L . C . Cross Eds. Special Publication No. 11, Chem. Soc. London (1958) 93. G.R. Dennis, I.R. Gentle, G.L.D. Ritchie and C.G. Andrieu, J . Chem. Soc. Farad. Trans. H, 79, 539 (1983) 94. R.C. Long Jr., S.L. Baughcum and J . H . Goldstein, J . Magn. Res. 7, 253 (1972) 95. B. Bak, D. Christensen, W.B. Dixon, L. Hansen-Nygaard, J.Rastrup-Andersen and M . Schottlander, J . Mol. Spec. 9, 124 (1962) 96. B. Bak, D. Christensen, W.B. Dixon and Rastrup-Andersen, J . Mol. Spec. 7, 58 (1961) 97. E . E . Burnell and C A . de Lange, Mol. Phys. 16, 95 (1969) E . E . Burnell and L . Weiler, Chem. L. Phys. Lett. 50, 243 Hansen-Nygaard, J. Bibliography / 164 98. B. Bak, L. Nygaard-Hansen 361 (1958) and J . Rastrup-Andersen, J . Mol. Spec. 2, 99. P.B. Lukins and G.L.D. Ritchie, J . Phys. Chem. 89, 1312 (1985) 100. J . Gerritsen and C. MacLean, Recueil, 91, 1393 (1972) 101. L . Nygaard, J . Bojesen, T. Pederson and J . Rastrup-Andersen, Struct. 2, 209 (1968) 102. H.B. Evans Jr., A.R. Tarpley and J . H . Goldstein, 2552 (1968) 103. F . Michel, H . Nery and G. Roussy, Compt. Rendu, 278B, 203 (1974) 104. R.M. Metzger and C H . Rhee, Mol. Cryst. Liq. Cryst. 85, 81 (1982) 105. W.A. Thomas and G.E. Griffin, Org. Magn. Res. 2, 503 (1970) 106. O.L. Stiefvater, Z. Naturforsch. 30A, 1765 (1975) 107. P. Diehl and C.L. Khetrapal, Mol. Phys. 15, 201 (1968) 108. M . Onda and I. Yamaguchi, J . Mol. Struct. 34, 1 (1976) 109. P. Diehl, J . Amrein, H . Bbsiger (1982) 110. J . M . Read Jr., R.W. Crecely, R.S. Butler Goldstein, Tet. Lett. 10, 1215 (1968) 111. J . Martin and B.P. Dailey, J . Chem. Phys. 37, 2594 (1962). 112. O.G. Batyukhnova, N.I. Sadova, L . V . Vilkov and Yu. A. Pankrushev, J . Mol. Struct. 97, 153 (1983) 113. J . Trotter and C.S. Williston, Acta Cryst. 21, 285 (1966) 114. E . E . Burnell, P. Diehl (1974) 115. E . E . Burnell and M.A. Sweeny, Can. J . Chem. 52, 3565 (1974) 116. M.Y. Kok, M.Sc. Thesis, University of British Columbia (1986) 117. A.D. Buckingham, E . E . Burnell and C.A. de Lange, J . Amer. Chem. Soc. 90, (1968) 118. P. Diehl and A.S. Tracey, Mol. Phys. 30, 1917 (1975) J . Phys. J . Mol. Chem. 72, and F. Moia, Org. Mag. Res. 18, 20 and W. Niederberger, and J . E . Loemker, J . H . Can. J . Chem. 52, 151 Bibliography / 165 119. R.G. Barnes and J.W. Bloom, J. Chem. Phys. 57, 3082 (1972) 120. J. Vrbancich and G.L.D. Ritchie, J. Chem. Soc. Farad. Trans. II, 76, 648 (1980) 121. R.L. Shoemaker and W.H. Flygare, J. Chem. Phys. 51, 2988 (1969) 122. M.R. Battaglia, A.D. Buckingham 78, 421 (1981) 123. D.E. Stogryn and A.P. Stogryn, Mol. Phys. 11, 371 (1966) 124. S.G. Kukolich, W.G. Read, J.A. Shea and E.J. Campbell, J. Amer. Chem. Soc. 105, 6423 (1983) 125. R.D. Amos and J.H. Williams, Chem. Phys. Lett. 66, 471 (1979) 126. W. Ravenek and F.M.M. Geurts, Chem. Phys. 90, 73 (1984) 127. W.C. Pringle, (1983) 128. G.L.D. Ritchie, Chem. Phys. Lett. 93, 410 (1982) 129. C L . Norris, E.F. Pearson and W.H. Flygare, J. Chem. Phys. 60, 1758 (1974) 130. P.G. Del Conde, P.S. Bagus Acta. 45, 121 (1977) 131. R.L. Shoemaker and W.H. Flygare, J. Amer. Chem. Soc. 91, 5417 (1969) 132. D.H. Sutter and W.H. Flygare, J. Amer. Chem. Soc. 91, 4063 (1969) 133. G.S. Bates, E.E. Burnell, G.L. Hoatson, P. Palffy-Muhoray Weaver, Chem. Phys. Lett, (in press) (1987) 134. E.E. Burnell and CA. de Lange, Chem. Phys. Lett. 106, 413 (1984) 135. E.E. Burnell and CA. de Lange, Chem. Phys. Lett.(in press) (1987) and J.H. Williams, Chem. Phys. Lett. S.M. Jacobs and D.H. Rosenblatt, Mol. Phys. 50, 205 and C.W. Bauschlicher Jr., Theor. Chim. and A.
- Library Home /
- Search Collections /
- Open Collections /
- Browse Collections /
- UBC Theses and Dissertations /
- NMR of solutes in nematic liquid crystals : an investigation...
Open Collections
UBC Theses and Dissertations
Featured Collection
UBC Theses and Dissertations
NMR of solutes in nematic liquid crystals : an investigation of the mechanisms of orientational ordering Van der Est, Arthur James 1987
pdf
Page Metadata
Item Metadata
Title | NMR of solutes in nematic liquid crystals : an investigation of the mechanisms of orientational ordering |
Creator |
Van der Est, Arthur James |
Publisher | University of British Columbia |
Date Issued | 1987 |
Description | Dipolar and quadrupolar couplings measured from NMR experiments, and order parameters obtained from these couplings, are reported for a number of small solutes dissolved in several nematic liquid crystals. These results are discussed in terms of the solute-solvent interactions. It has been shown that the interaction between the solute molecular quadrupole moment and a mean external electric field gradient due to the liquid crystal accounts for most but not all of the ordering of molecular hydrogen. The remaining contribution to ordering is discussed in terms of possible mechanisms. The anisotropic couplings observed for methane and acetylene are discussed in terms of a model which takes into account the interaction between the vibrations and rotations of the solute. Excellent agreement between the observed and calculated dipolar couplings is obtained. Evidence is given that these two solutes experience the same field gradient as molecular hydrogen. In a mixture of 55wt per cent 1132 (Merck ZLI 1132) and 45wt per cent EBBA (N-(4-ethoxybenzylidene)-4'-n-butylaniline at 301.4K the deuterons in D₂ experience no external electric field gradient. The order parameters of a series of solutes in this mixture are calculated in terms of a simple model for the short range interactions which depend on the size and shape of the solute. These calculated order parameters are in very good agreement with the experimental results. In liquid crystals where the field gradient is not zero it is shown that the combination of the short range interaction model and the electric field gradient - molecular quadrupole moment mechanism predicts the order parameters very well. |
Subject |
Molecular dynamics Nuclear magnetic resonance spectroscopy Liquid crystals |
Genre |
Thesis/Dissertation |
Type |
Text |
Language | eng |
Date Available | 2010-08-20 |
Provider | Vancouver : University of British Columbia Library |
Rights | For non-commercial purposes only, such as research, private study and education. Additional conditions apply, see Terms of Use https://open.library.ubc.ca/terms_of_use. |
DOI | 10.14288/1.0060386 |
URI | http://hdl.handle.net/2429/27557 |
Degree |
Doctor of Philosophy - PhD |
Program |
Chemistry |
Affiliation |
Science, Faculty of Chemistry, Department of |
Degree Grantor | University of British Columbia |
Campus |
UBCV |
Scholarly Level | Graduate |
Aggregated Source Repository | DSpace |
Download
- Media
- 831-UBC_1987_A1 V36.pdf [ 8.26MB ]
- Metadata
- JSON: 831-1.0060386.json
- JSON-LD: 831-1.0060386-ld.json
- RDF/XML (Pretty): 831-1.0060386-rdf.xml
- RDF/JSON: 831-1.0060386-rdf.json
- Turtle: 831-1.0060386-turtle.txt
- N-Triples: 831-1.0060386-rdf-ntriples.txt
- Original Record: 831-1.0060386-source.json
- Full Text
- 831-1.0060386-fulltext.txt
- Citation
- 831-1.0060386.ris
Full Text
Cite
Citation Scheme:
Usage Statistics
Share
Embed
Customize your widget with the following options, then copy and paste the code below into the HTML
of your page to embed this item in your website.
<div id="ubcOpenCollectionsWidgetDisplay">
<script id="ubcOpenCollectionsWidget"
src="{[{embed.src}]}"
data-item="{[{embed.item}]}"
data-collection="{[{embed.collection}]}"
data-metadata="{[{embed.showMetadata}]}"
data-width="{[{embed.width}]}"
async >
</script>
</div>
Our image viewer uses the IIIF 2.0 standard.
To load this item in other compatible viewers, use this url:
http://iiif.library.ubc.ca/presentation/dsp.831.1-0060386/manifest