Theory of solvation in polar nematics

Vitaly Kapko,Dmitry V. Matyushov

DOI: https://doi.org/10.1063/1.2178318

2005-12-15

Abstract:We develop a linear response theory of solvation of ionic and dipolar solutes in anisotropic, axially symmetric polar solvents. The theory is applied to solvation in polar nematic liquid crystals. The formal theory constructs the solvation response function from projections of the solvent dipolar susceptibility on rotational invariants. These projections are obtained from Monte Carlo simulations of a fluid of dipolar spherocylinders which can exist both in the isotropic and nematic phase. Based on the properties of the solvent susceptibility from simulations and the formal solution, we have obtained a formula for the solvation free energy which incorporates experimentally available properties of nematics and the length of correlation between the dipoles in the liquid crystal. Illustrative calculations are presented for the Stokes shift and Stokes shift correlation function of coumarin-153 in 4-n-pentyl-4'-cyanobiphenyl (5CB) and 4,4-n-heptyl-cyanopiphenyl (7CB) solvents as a function of temperature in both the nematic and isotropic phase.

Chemical Physics,Computational Physics

What problem does this paper attempt to address?

The main problem that this paper attempts to solve is to develop a linear response theory to describe the solvation process of ions and dipolar solutes in anisotropic, axially symmetric polar solvents. Specifically, the focus of the research is on the solvation phenomena in polar nematic liquid crystals. The following are the key points of the paper:

1. **Problem Background**:
- The solvation problem is a long - standing and active research topic in the field of physical chemistry.
- Calculating the solvation free energy is particularly complex because it is affected by multiple factors, including short - range cavity formation energy, medium - range dispersion and induction forces (non - polar solvation), and long - range electrostatic interactions (polar solvation).

2. **Research Objectives**:
- Through the linear response theory, construct a solvation response function, which is based on the projection of solvent dipole susceptibility onto rotation invariants.
- Provide a quantitative framework to analyze steady - state and time - resolved optical spectra and make several experimentally verifiable predictions.

3. **Main Results**:
- Derive the solvation free energy formula, combining the experimental properties of nematic liquid crystals and the correlation length between dipoles.
- Discover that the anisotropy of solvation depends on the dielectric anisotropy of the solvent: in solvents with positive dielectric anisotropy, the solvation anisotropy is negative, and vice versa.
- The solvation free energy is discontinuous at the nematic - isotropic phase transition point, and this discontinuity is affected by the size of the solute and the microscopic dipole correlation in the liquid - crystal solvent.

4. **Specific Applications**:
- The paper provides illustrative calculations of the Stokes shift and its correlation function of fluorescein - 153 in 4 - n - pentyl - 4' - cyanobiphenyl (5CB) and 4,4 - n - heptyl - cyanobiphenyl (7CB) solvents when the temperature changes in the nematic and isotropic phases.

### Mathematical Expressions

To understand these concepts more clearly, the following are the Markdown - format representations of some key formulas:

- Solvation response function:
\[
P(\mathbf{r}_1)=\int\chi(\mathbf{r}_1,\mathbf{r}_2)\cdot\mathbf{E}_0(\mathbf{r}_2)d\mathbf{r}_2
\]

- Solvation free energy:
\[
\mu =-\frac{1}{2}\int P(\mathbf{r}_1)\cdot\mathbf{E}_0(\mathbf{r}_1)d\mathbf{r}_1
\]

- Chemical potential in k - space:
\[
\mu =-\frac{1}{2}\int\frac{d\mathbf{k}_1d\mathbf{k}_2}{(2\pi)^6}\tilde{\mathbf{E}}_0(\mathbf{k}_1)\cdot\tilde{\chi}(\mathbf{k}_1,\mathbf{k}_2)\cdot\tilde{\mathbf{E}}_0(-\mathbf{k}_2)
\]

- Dipolar solvation field:
\[
F^{(d)}_\alpha=\frac{m_{0,\alpha}}{R_1^3R_\alpha}
\]

- Anisotropic solvation free energy:
\[
\mu^{(d)}_c =-\frac{m_0^2}{2R_1^3}\left[R_x+(R_z - R_x)\cos^2\theta_0\right]
\]

Through these formulas and theoretical frameworks, the authors are able to deeply explore and explain the complexity of solvation behavior in polar nematic liquid crystals and its changes depending on temperature and solute shape.