Interface structures in ionic liquid crystals

Ionic liquid crystals (ILCs) are anisotropic mesogenic molecules which additionally carry charges. This combination gives rise to a complex interplay of the underlying (anisotropic) contributions to the pair interactions. It promises interesting and distinctive structural and orientational properties to arise in systems of ILCs, combining properties of liquid crystals and ionic liquids. While previous theoretical studies have focused on the phase behavior of ILCs and the structure of the respective bulk phases, in the present study we provide new results, obtained within density functional theory, concerning (planar) free interfaces between an isotropic liquid $L$ and two types of smectic-A phases ($S_A$ or $S_{AW}$). We discuss the structural and orientational properties of these interfaces in terms of the packing fraction profile $\eta(\vec{r})$ and the orientational order parameter profile $S_2(\vec{r})$ concerning the tilt angle $\alpha$ between the (bulk) smectic layer normal and the interface normal. The asymptotic decay of $\eta(\vec{r})$ and of $S_2(\vec{r})$ towards their values in the isotropic bulk is discussed, too.


I. INTRODUCTION
Ionic liquid crystals (ILCs) are pure ionic systems, solely composed of cations (+) and anions (−).Moreover, at least one of the ion species is characterized by a highly anisotropic molecular shape [1].This anisotropic shape is typically due to long alkyl-chains which are attached to charged moieties.Although the alkyl-chains exhibit a rather strong flexibility, due to microphase segregation of the charged parts and of the alkyl-chains, liquid-crystalline phases are indeed observable among ILCs [1][2][3].In the past decades various types of ILCs have been synthesized [1,3].Different combinations of, e.g., (charged) imidazolium rings and alkyl-chains allow one to tune not only the length of the ionic mesogenes but also the location of their charges, i.e., the intra-molecular charge distribution.Thereby one is able to promote distinctive properties of ILCs, for instance, a high thermal and high electrochemical stability, which might be beneficial for technological applications [1,[3][4][5][6].(We note, that here the term "mesogene" refers to any kind of molecule which gives rise to the formation of mesophases, irrespective of the underlying microscopic mechanism.Accordingly, the aforementioned anisotropic molecules, which form mesophases via microphase segregation, are considered to be mesogenes.) A specific example of an ILC system, which has been studied, e.g., in Refs.[7,8], is composed of cations with long alkyl-chains attached (1-dodecyl-3methylimidazolium) and significantly smaller anions (iodide).For such an ILC system, one observes a liquid crystalline structure, in particular the smectic-A phase S A .(The S A phase is characterized by layers of particles which are well aligned with the layer normal and the layer spacing is of the size of the particle length.)The layer structure of the large cations leads to a locally increased concentration of anions in between the layers of cations [7].Thereby, the nanostructure of the cations gives rise to "pathways" for the anions, which increase the conductivity measurable in the direction parallel to the layers.Therefore this particular type of an ILC system is a promising candidate for technological applications, e.g., as electrolyte in dye-sensitized solar cells (DSSCs) [7,9].
While the complexity of the underlying interactions gives rise to these interesting properties of ILCs, it is at the same time very challenging to study these systems within theory or simulations.Previous theoretical studies [10,11] of ILC systems have been able to reduce this complexity by considering a simplified description of ILC systems, which incorporates, however, the generic properties of ILCs.They rely on an effective one-species description in which one of the ion species (referred to as counterions) is not accounted for explicitly, but is incorporated as a continuous background, giving rise to the screening of the coions.On the contrary, the coions are modeled as ellipsoidal particles.Thus, the anisotropic molecular shape, which gives rise to the formation of mesophases, and the (screened) electrostatic interaction are both incorporated by this approach.Of course, this is a simplified representation of any realistic ionic liquid crystalline system.However, it allows one to study the interplay of the two key features, i.e., an anisotropic molecular shape and the presence of charges, which are omnipresent in ILC systems.Yet it should be noted, that ILC systems exhibiting a significant difference in size of the cations and of the anions (e.g., the aforementioned example of 1-dodecyl-3-methylimidazolium) might be can-didates which come closest to the present theoretical representation of ILCs, as the size difference rationalizes in parts the idea of structureless point-like counterions.
As a first step, such a model allows one to study the phase behavior of ILC systems and thereby to gain insight about how molecular properties, e.g., the aspectratio or the charge distribution of the molecules, affect the phase behavior of such types of ILCs.A comprehensive understanding of the relation between the underlying molecular properties and the resulting phase behavior is inter alia, necessary for a systematic synthesis of ILCs, which should meet specific material properties.Furthermore, theoretical guidance is beneficial for finding and exploring novel materials properties which might occur in ILC systems.For instance, in Ref. [11] a new smectic-A structure (S AW ) has been observed, which exhibits an alternating layer structure.In between layers of elongated particles, which prefer to be oriented parallel to the layer normal, like in the ordinary S A phase, one observes secondary layers in which the particles prefer to be oriented perpendicular to it.Due to this alternating structure the layer spacing of this new S AW phase is significantly w ider compared to the ordinary S A phase.The S AW structure is stabilized by charges which are located at the tips of the molecules.This shows in an exemplary way how the combination of liquid-crystalline behavior and electrostatics can lead to an interesting and novel phenomenology.
The aim of the present investigation is to extend the analysis by studying spatially inhomogeneous systems of ILCs.This is done by investigating how the structural and orientational properties of ILC systems are affected by the presence of a free interface between coexisting bulk states.Both smectic-A phases, S A and S AW , observed in Ref. [11] can be in coexistence with the isotropic liquid phase L. This is of intrinsic interest, because it allows one to investigate interfaces which interpolate between a structured and orientationally ordered (i.e., smectic) phase and an isotropic, homogeneous, and thus structureless, fluid phase.In particular, the transition in the structural and in the orientational order allows one to study the interplay of both properties while they build up at the interface.Although there are theoretical analyses [12][13][14][15][16][17][18] concerning related types of free interfaces, in these studies the constituent particles are plain liquid crystals without any charges.On the other hand, there is a vast number of theoretical studies on ionic fluids.The thermodynamic behavior [19][20][21][22] as well as the structure [23][24][25][26][27] of these types of fluids, in which longranged Coulomb interactions are present, have been intensively studied.However, ionic systems are often analyzed assuming a simple geometry of the particles, such as a spherical shape of the particles like in the restricted primitive model [28][29][30].In this regard, the present study attempts to analyze the aforementioned type of interface between an isotropic and a smectic phase by accounting for an anisotropic particle shape combined with the presence of charges.
Moreover, different orientations between the interface normal and the smectic layer normal are possible.In this context, an interesting question addresses the equilibrium tilt angle between the interface and the smectic layer normal.This angle may provide insight into nucleation and growth phenomena which are affected by the dependence of the interfacial tension on the orientation of the considered structure [31,32].
The present study is structured as follows: In Sec.II the model and the employed density functional theory approach are presented.Our results for the interfaces between the isotropic liquid L and the considered smectic-A phases S A or S AW are discussed in Sec.III.Finally, in Sec.IV we summarize the results and draw our conclusions.

II. MODEL AND METHODS
This section presents in detail the molecular model of ILCs as employed here.In particular, we discuss the intermolecular pair potential, which can be applied to a wide range of ionic and liquid crystalline materials due to its flexibility provided by a large set of parameters.
This model is studied by (classical) density functional theory (DFT), which will be applied to spatially inhomogeneous systems, in particular free interfaces formed between coexisting bulk phases.The methodological and technical details of the present DFT approach are described in Sec.II B.

A. Molecular model and pair potential
We consider a coarse-grained description of the ILC molecules as rigid prolate ellipsoids of length-to-breadth ratio L/R ≥ 1 (see Fig. 1).Thus, the orientation of a molecule is fully described by the direction ω ω ω(φ, ϑ) of its long axis, where φ and ϑ denote the azimuthal and polar angle, respectively.
The two-body interaction potential consists of a hard core repulsive and an additional contribution U GB + U es beyond the contact distance Rσ, the sum of which can be attractive or repulsive: where r r r 12 := r r r 2 −r r r 1 denotes the center-to-center distance vector between the two particles labeled as 1 and 2, and ω ω ω i , i = 1, 2, are their orientations with |ω ω ω i | = 1.The contact distance Rσ(r r r12 , ω ω ω 1 , ω ω ω 2 ) depends on the orientations of both particles and on the direction of the centerto-center distance vector, which is expressed by the unit vector r r r12 := r r r 12 /|r r r 12 |.In Eq. ( 1), we have subdivided the contributions beyond the contact distance |r r r Cross-sectional view of two ILC molecules in the plane spanned by the orientations ω ω ωi, i = 1, 2, of their long axis.The particles are treated as rigid prolate ellipsoids, characterized by their length-to-breadth ratio L/R ≥ 1.Their orientations are fully described by the direction of their long axis ω ω ωi, i = 1, 2; r r r12 is the center-to-center distance vector.The charges of the ILC molecules (blue dots) are located on the long axis at a distance D from their geometrical center.The counterions are not modeled explicitly, but they are implicitly accounted for in terms of a background, giving rise to the screening of the charges of the ILC molecules.
The second contribution U es (r r r 12 , ω ω ω 1 , ω ω ω 2 ) in Eq. ( 1) is the electrostatic repulsion of ILC molecules.Within the scope of the present study, the counterions are not modeled explicitly.They will be considered to be much smaller in size than the ILC molecules such that they can be treated as a continuous background.On the level of linear response, this background gives rise to the screening of the pure Coulomb potential between two charged sites on a length scale given by the Debye screening length λ D such, that the effective electrostatic interaction of the ILC molecules is given by The charges q are located symmetrically on the long axis of the ILC molecules at a distance D from the geometrical center of the particles (compare Fig. 1).The prefactor γ = q 2 /(4πε) of dimension [energy] × [length] characterizes the electrostatic energy scale, where ε denotes the permittivity.In principle, the Debye screening length is a function of temperature T and of the number density ̺ c of the counter ions.Thus, it depends on the thermodynamic state of the fluid.However, in the present model λ D is taken to be a constant parameter.In order to compare results, obtained within this model, with data from actual physical systems, one could measure the value of the Debye screening length experimentally and tune the model parameter λ D accordingly.In Fig. 2 we illustrate the full pair potential (Eq. ( 1)) beyond the contact distance for certain choices of the parameters.The two top panels, (a) and (b), show the pure Gay-Berne potential (uncharged liquid crystals), which is predominantly attractive in the space outside the overlap volume (salmon-colored area).The shape of this over-lap volume changes by varying the particle orientations as well as by changing the length-to-breadth ratio L/R.However, these dependences are not apparent from Fig. 2, because L/R = 4 and the particle orientations ω ω ω i are kept fixed for all panels.In panel (b) the anisotropy parameter ε R /ε L = 4 is chosen to be two times larger than for panel (a) (ε R /ε L = 2).Thus, the ratio of the well depth at the tails and at the sides is increased.The two bottom panels, (c) and (d), show the same choices for the Gay-Berne parameters as for panel (a), but the electrostatic repulsion of the charged groups on the molecules, illustrated by blue dots, is included (γ/(Rε 0 ) = 0.25).In panel (c) the loci of the two charges of the particles coincide at their centers (i.e., D/R = 0) while in panel (d) they are located near the tips (D/R = 1.8).For both cases with charge, the effective interaction range is significantly increased compared with the uncharged case and is governed by the Debye screening length, chosen as It is worth mentioning, that the present model cannot be considered as a quantitatively valid description of any realistic ionic liquid crystal system.A screened electrostatic pair interaction of the Yukawa form (Eq. ( 5)) is the extreme case of the effective pair potential between ions in a (dilute) electrolyte at high temperatures.Nonetheless, for the purpose of the present theoretical study, which is concerned with the basic microscopic mechanisms and the generic molecular properties present in ILC systems, the employed model is appropriate as it incorporates the following key properties of ILCs: First, a sufficiently anisotropic shape (prolate) of the particles, i.e., they can be considered as (calamitic) mesogenes.In this context, an assessment of the bulk phase behavior, depending on the length-to-breadth ratio of the particles, is provided in Ref. [11].In particular, the relevance of a sufficiently anisotropic shape (i.e., L/R > 2) for observing genuine smectic phases is discussed.Second, the ionic properties of ILCs are incorporated such that they reflect the main feature of ionic fluids, i.e., the effective interaction of the ionic compounds via a screened electrostatic pair interaction.Although the chosen functional form given by Eq. ( 5) cannot be considered as a quantitatively reliable representation, it still accounts for the fact that the actual ion-ion pair interaction in an ionic fluid is indeed short-ranged, rather than long-ranged, as it is the case for the bare Coulomb interaction.
In conclusion, Eq. ( 5) is characterized by an effective interaction strength γ/R (which will be numerically expressed as the relative interaction strength γ/(Rε 0 ) compared to the interaction strength ε 0 of the Gay-Berne potential), an effective interaction range λ D , and an effective location D of the charge sites inside the coions.In order to represent specific ILC molecules by a particular set of parameters of the present model, one would tune the independent model parameters, i.e., L/R, ε R /ε L , γ/(Rε 0 ), D/R, and λ D /R such that the resulting total pair potential U (r r r 12 , ω ω ω 1 , ω ω ω 2 )/ε 0 (compare Eq. (1) and Fig. 2) resembles (qualitatively) the actual pair poten-tial of the considered ILC molecules.In this regard, it is worth mentioning that in principle comparisons of our effective theory with particle simulations can be made, related to the study by Saielli et al. [35], who performed molecular dynamics (MD) simulations for a mixture of (ellipsoidal) Gay-Berne and (spherical) Lennard-Jones particles.Additionally, both species carry charges and therefore resemble cations and anions, respectively.Our ad-hoc pair potential (Eq. ( 1)) of the coions can be compared with the effective interaction, which can be determined as the logarithm of the particle-particle distribution function of the elongated cations in the MD simulations.
We note, that the choices L/R = 4 and ε R /ε L = 2, which are used throughout our analysis, are comparable to those used in previous studies (see, e.g., Refs.[10,11,35,36]) for similar kinds of particles.While these values of the Gay-Berne parameters give rise to the formation of smectic phases, the occurrence of nematic phases is typically observed for much larger values of L/R and ε R /ε L [37].

B. Density functional theory
The degrees of freedom of the particles (compare Sec. II A) are fully described by the positions r r r of their centers and the orientations ω ω ω of their long axes.Thus, within density functional theory, an appropriate variational grand potential functional βΩ[̺] of position-and orientation-dependent number density profiles ̺(r r r, ω ω ω) has to be found; the equilibrium density profile minimizes the functional.The grand potential functional for uniaxial particles, in the absence of external fields, can generically be expressed as where the integration domains V and S denote the system volume and the full solid angle, respectively.The first term in Eq. ( 7) is the purely entropic free energy contribution of non-interacting uniaxial particles, where β = 1/(k B T ) denotes the inverse thermal energy, µ the chemical potential, and Λ the thermal de Broglie wavelength.
The last term is the excess free energy βF [̺] in units of k B T , which incorporates the effects of the inter-particle interactions.Minimizing Eq. ( 7) leads to the Euler-Lagrange equation, which implicitly determines the equilibrium density profile ̺(r r r, ω ω ω): where FIG. 3. Sketch of the interface structure under consideration.Consider a planar interface, illustrated by the horizontal red line, between the isotropic bulk liquid L, imposed as the boundary condition at z → −∞, and the smectic-A phase SA (or SAW ), imposed as the boundary condition at z → +∞.Thus, the interface normal (red vertical arrow) points into the z-direction.At the top, the tails of four layers of particles of the (ordinary) SA phase are visible, which are well aligned with the smectic layer normal n n n := sin(α) x x x + cos(α) ẑ ẑ ẑ.In the bulk SA phase, the system is periodic in the direction of the smectic layer normal n n n with periodicity d which is a multiple of the smectic layer spacing.For the SA phase d turns out to be two times the distance between neighboring smectic layers (see Sec. II B below Eq. ( 15)).Thus, for a given tilt angle α between the interface normal and the smectic layer normal n n n, the system is periodic in x-direction with periodicity dx = d/ sin(α).Note that the interface structure is translationally invariant in the y-direction for all angles 0 ≤ α ≤ π/2.For α = 0 the system exhibits in addition translational invariance in the x-direction.
is the one-particle direct correlation function.It is fully determined by the excess free energy functional βF [̺].
The excess free energy functional is the characterizing quantity of the underlying many-body problem.However, in general it is not known exactly so that one has to adopt appropriate approximations of it.Following the approach of our previous study [11] concerning the bulk phase behavior of ILCs, in the spirit of Ref. [38] a weighted density expression for βF [̺] is considered: where βψ(r r r, ω ω ω, [̺]) denotes the effective one-particle potential.It is a functional of the so-called projected density ̺(r r r, ω ω ω): where P 2 (y) = (3y 2 − 1)/2 is the Legendre polynomial of degree 2. We point out that ̺(r r r, ω ω ω) represents an expansion of the density profile ̺(r r r, ω ω ω) in terms of a secondorder Fourier, and a second-order Legendre series, respectively.Thus, the coefficients Q i (r r r) are the corresponding expansion coefficients, which will be defined below.It is worth mentioning, that although the projected density ̺(r r r, ω ω ω) might take negative values, this does not imply an unphysical behavior as the actual density ̺(r r r, ω ω ω) is determined from the Euler-Lagrange equation (i.e., Eq. ( 8)) and thus is strictly positive.The following three types of bulk phases can be studied within this particular framework [11]: First, isotropic liquids with Q 0 = const 0 and Q i = 0 for i > 0. Second, nematic liquids with While for isotropic and nematic liquids the system is translationally invariant in all spatial directions, in the case of smectic-A phases the system is periodic in the direction of the smectic layer normal n n n with periodicity d, which is a multiple of the smectic layer spacing.For smectic-A phases the director is parallel to the smectic layer normal n n n and therefore the occurrence of rotationally symmetric distributions of the orientations ω ω ω around n n n, incorporated by the dependence on ω ω ω•n n n in Eq. ( 11), are plausible.We note that odd Fourier-modes in the projected density ̺(r r r, ω ω ω) vanish for bulk smectic-A phases, if the coordinate system is chosen such that the origin is located at the center of one of the smectic layers due to the mirror symmetry of smectic layers around their center.This is a direct consequence of the underlying point symmetry of the particles considered here (see Fig. 1).
Considering additional terms, corresponding to the odd modes in the second-order Fourier expansion of the density ̺(r r r, ω ω ω), would only give rise to a shift of the location of the bulk smectic layers.Although for systems with interfaces the odd modes in general do not vanish, here we neglect these contributions completely.The implications of additionally considering the odd terms (up to second order) are discussed in Appendix A. Both approaches are weighted-density-like approximations of the exact free energy functional.A priori, it is not obvious which one leads to better results, because considering more terms of the Fourier series leads only to a more accurate representation of ̺(r r r, ω ω ω) by the projected density ̺(r r r, ω ω ω).However, this does not imply that the resulting free energy functional βF [̺] is closer to its exact form, because independent of the choice for ̺(r r r, ω ω ω) it relies on the Parsons-Lee approach for its reference part and on the so-called modified mean-field approximation for the excess part (see below).Nevertheless, our approach of considering in Eq. ( 11) only the even modes up to second order captures the three types of bulk phases L, N , as well as S A , S AW , which are relevant for the present study in the same way as the full second-order Fourier expansion.
The effective one-particle potential βψ(r r r, ω ω ω) consists of two contributions.The first one incorporates the hardcore interactions via the well-studied Parsons-Lee functional [39,40], where f M (r r r − r r r ′ , ω ω ω, ω ω ω ′ ) is the Mayer f-function [41] of the hard core pair interaction potential and J (Q 0 ) modifies the corresponding original Onsager free energy functional (i.e., the second-order virial approximation) such that the Carnahan-Starling equation of state [40] is reproduced for spheres, i.e., for L = R [42,43]: where η 0 (Q 0 ) = Q 0 LR 2 π/6 (for Q 0 see Eq. ( 11)) denotes the mean packing fraction within one (bulk) smectic layer.It is proportional to the coefficient Q 0 , which gives the mean density within such a smectic layer (see below).The original Onsager functional is recovered in Eq. ( 12) by replacing The second contribution to the effective one-particle potential βψ[̺] takes into account the interactions beyond the contact distance (see the case |r r r 12 | ≥ Rσ in Eq. ( 1)) within the modified mean-field approximation [44], which is a variant of the extended random phase approximation (ERPA) [45]: The present study is devoted to the analysis of free interfaces which are formed between coexisting bulk phases.In particular, the planar interfaces between the isotropic liquid L and the two different types of smectic-A phases (S A or S AW , see Sec.III) will be considered, for which the interface normal is expected to be parallel to the z-direction (see Fig. 3).Due to the isotropy of the liquid phase L, the direction of the smectic layer normal n n n(α) := sin(α) x x x + cos(α) ẑ ẑ ẑ (15) can be chosen to lay in the x-z-plane.Its orientation is fully determined by the tilt angle α.For α = 0 the smectic layer normal n n n = ẑ ẑ ẑ points into the z-direction, like the interface normal, while for α = π/2 it points into the x-direction and thus it is perpendicular to the interface normal.The interfacial systems considered here are translationally invariant in the y-direction and show a periodic structure in the x-direction with a periodicity d x = d/ sin(α) (compare Fig. 3) where d is a multiple of the smectic layer spacing.(We note that the value of d is determined by the corresponding bulk density distribution which minimizes the grand potential functional, i.e., maximizes the bulk pressure (see Sec. 2.2.2 in Ref. [11]).It turns out that, for the S AW phase d equals the smectic layer spacing, while for the S A phase it equals two times the layer spacing, because for the S A phase one obtains bulk solutions 11), (16), and ( 17)).Thus the periodicity d along the layer normal n n n is twice the smectic layer spacing, i.e., the distance between neighboring layers.)For α = 0, d x diverges and the system is translationally invariant in the x-direction, too.
As mentioned above, the coefficients Q i (r r r) in Eq. ( 11) arise from expanding ̺(r r r, ω ω ω) in a second-order Legendreand Fourier-series [11]: with where T (r r r − r r r ′ ) is a cut-off function which defines the integration domain Due to the translational invariance in y-direction the extension of the integration domain V d can be chosen arbitrarily in the y-direction.
Due to the periodicity of ̺(r r r, ω ω ω) in the x-direction, this choice of the integration domain V d leads to coefficients Q i (z) (Eq.( 16)) which depend only on z, i.e., on the coordinate parallel to the interface normal.
For 0 < α ≤ π/2 one could also consider an integration domain which has a non-vanishing extent in z-direction.However, such a choice has at least two disadvantages: First, unlike d x , which corresponds to the periodicity of the system in x-direction, for 0 < α ≤ π/2 there is no obvious choice for the extent of V d parallel to the interface normal.Additionally, there is no unique choice for the geometrical shape of the integration domains; besides using a rectangular form, one could also use any other (two-dimensional) geometrical object as integration domain V d .In this sense the slice of length d x perpendicular to the interface normal is a simple but also consistent choice.Second, this choice renders the evaluation numerically less demanding, because it requires only a one-dimensional integration (exploiting the translational invariance in y-direction), instead of evaluating a twodimensional integral.We note, that an infinite extent of the integration domain parallel to the interface normal leads to coefficients Q i which are independent of the position r r r and therefore cannot be used to obtain interface profiles.
If α = 0, i.e., the smectic layer normal n n n = ẑ ẑ ẑ is parallel to the interface normal, d x diverges and the the system is translationally invariant in x-and y-direction.In this case, the integration domain V d has an extent of length d in z-direction, i.e., T (r r r −r r r ′ ) = Θ(d/2−|z −z ′ |) with arbitrary extent in the lateral dimensions x and y.As before, the coefficients Q i (z) depend only on the z-coordinate.It is worth mentioning, that for all tilt angles 0 ≤ α ≤ π/2 the correct (constant) bulk values of the coefficients Q i are recovered, although for 0 < α < π/2 the orientation of the integration domain V d (recall that V d is a slice of width d x in x-direction for all α ∈ (0, π/2]) changes with respect to the direction of the smectic layer normal n n n(α).However, because the integration domain covers a full period d x in x-direction, it gives the same values for the coefficients Q i in the bulk phases, as for evaluating the coefficients Q i with an integration domain parallel to the smectic layer normal n n n, which is the case for α = 0 and π/2.
Equation ( 8) has been solved numerically (utilizing a Picard scheme with retardation) by using Eq. ( 19) as well as the (constant) bulk values of the coefficients The structural properties and the orientational order at the free interface are analyzed in terms of the interface profiles of the packing fraction with the number density n(r r r) := S d 2 ω ̺(r r r, ω ω ω), and in terms of the orientational order parameter f (r r r, ω ω ω) := ̺(r r r, ω ω ω)/n(r r r) describes the orientational distribution.

C. Gibbs dividing surface
The position z η of the interface is determined by the density profile ̺(r r r, ω ω ω) for which we have adopted the notion of the Gibbs dividing surface [41]: where η 0 (r r r) = Q 0 (r r r) LR 2 π/6 is the mean packing fraction at position r r r.The quantities η 0,L = η(z → −∞) and η 0,SA = η(z → ∞) are the bulk values of η 0 (r r r) in the isotropic liquid phase L and the smectic-A phase S A (or S AW ), respectively.The interface position z η in Eq. ( 22) corresponds to the location of a step-like profile such that the number of particles in excess and in deficit of the bulk values is the same on both sides of the interface.Taking the derivative of the left-hand side h η (z) of Eq. ( 22) with respect to z leads to h ′ η (z = z η ) := η 0,SA − η 0,L which is a constant.Therefore h η (z) = (η 0,SA − η 0,L )z + h η (0) is a linear function and one has to evaluate h η (0) only once in order to obtain using Eq. ( 22), i.e., h η (z η ) = 0.While z η can be interpreted as the location of the transition in the structure from the isotropic liquid L to the smectic-A phase S A (or S AW ), replacing η by S 2 in Eq. ( 22) defines a position which corresponds to the transition in the orientational order from one phase to the other.Note, that instead of using the mean packing fraction η 0 or the mean orientational order parameter S 20 in Eq. ( 22), for determining the interface positions, in principle, one could also use the profiles η(r r r) and S 2 (r r r) directly.However, the disadvantage of this latter approach is that in the smectic-A bulk phase S A (or S AW ) the profiles η(r r r) and S 2 (r r r) are still functions of the position r r r (via the projection r r r • n n n onto the layer normal n n n).Typically, this prevents the use of the latter generalized Eqs. ( 23) and ( 24) for determining z η and z S2 .Instead, one has to solve Eq. ( 22) numerically, which requires many iterations depending on the desired accuracy.
Nevertheless, in the particular case α = π/2 the interface normal and the smectic layer normal are perpendicular.Due to the translational invariance of the smectic phases perpendicular to their layer normal, here the density profile η(z → ∞) and the orientational order parameter profile S 2 (z → ∞) do not depend on z for z → ∞ in the smectic bulk, but they depend only on the x-coordinate.Thus, for α = π/2 one can define interface contours zη (x) and zS2 (x), analogously to z η and z S2 : where m ∈ {η, S 2 }.

D. Interfacial tension
The interfacial tension Γ is a measure of the excess amount of work needed to form an interface between coexisting bulk phases [41].Accordingly, it can be calculated by determining the increase in the grand potential βΩ[̺] of the interface system in excess of the bulk grand potential βΩ 0 := −βpV which is given by the bulk pressure p (see Eq. ( 26) in Ref. [11]) times the system volume V: where A is the cross-sectional area of the system in lateral directions to the interface normal.Hence, Γ * (α) has the dimension 1/area.The pressure p coex := p(T coex , µ coex , d) at coexistence (T, µ) = (T coex , µ coex ) is the same in the isotropic liquid L and the smectic-A phase S A or S AW with the equilibrium layer spacing d.The equilibrium tilt angle α eq minimizes the interfacial tension Γ * (α = α eq ) (see Sec. III D).

III. RESULTS
In this section we present results for free interfaces formed between the isotropic liquid L and the smectic-A phase S A or S AW .The discussion focuses on two kinds of ionic liquid crystals (ILCs) which are described by the pair interaction potential U (r r r 12 , ω ω ω 1 , ω ω ω 2 ) (Eq. ( 1)), introduced in Sec.II A: First, ILCs with charges in the center, i.e., D = 0 (see Fig. 1 and Eq. ( 5)), and second, ILCs with charges at the tips, i.e., D/R = 1.8.In particular the structural and orientational properties of the interface are discussed in terms of the packing fraction profile η(r r r) and the orientational order parameter profile S 2 (r r r) for various relative orientations between the interface normal and the smectic layer normal, i.e., for different tilt angles α (see Fig. 3).All results presented here have been obtained via the density functional approach described in Sec.II B.
A. Interface normal parallel to the smectic layer normal (α = 0) First, we consider the case that the interface normal is parallel to the normal of the smectic layers, i.e., α = 0 (see Fig. 3).Both point into the z-direction and due to translational invariance in the x-and y-directions, the packing fraction η(z) and the orientational order parameter S 2 (z) are functions solely of the spatial coordinate z.For the case of an ionic liquid crystal with L/R = 4, ε R /ε L = 2, γ/(Rε 0 ) = 0.045, λ D /R = 5, and D = 0, i.e., the charges are localized in the center of the molecule, the bulk phase behavior is shown in the T * -η 0 -phase diagrams of Fig. 4(a) where T * = kT /ε 0 and η 0 = Q 0 LR 2 π/6 are the reduced temperature and the mean packing fraction, respectively.Within the considered temperature range T * ∈ [0.9, 1.65] solely a firstorder phase transition from the isotropic liquid phase L to the ordinary smectic-A phase S A occurs.The S A phase is characterized by a layer structure with a smectic layer spacing d/R ≈ 4.3, which is comparable to the particle length L/R = 4. Within the smectic layers the particles are well aligned with the smectic layer normal n n n.The blue lines in Fig. 4 L/R = 4 comparable with the particle length L. The particles in the layers are well aligned with the layer normal n n n.In panel (b), i.e., for D/R = 1.8 (the charges being located at the tips of the molecules), another smectic-A structure, referred to as the SAW phase can be observed at low reduced temperatures T * .The SAW phase exhibits an alternating structure, consisting of primary layers of particles being parallel to the layer normal and secondary layers in which the particles prefer to be perpendicular to it.This leads to an increased layer spacing d/R ≥ 7.5.The black dotted line in panel (b) marks the triple point at T * ≈ 1.23 for which the isotropic liquid L, the ordinary smectic-A phase SA, and the SAW phase are in three-phase-coexistence.A detailed description of the structural properties of the smectic-A phases SA and SAW , including illustrations of their microstructure, are provided in Ref. [11].The black dots (•) in panel (a), respectively the red dots (•) in (b), mark the coexisting bulk states at the reduced temperature T * = 1.3, respectively 0.9, imposed as boundary conditions for the free interfaces shown in Figs. 5 and 7-9.
and the light blue area in between the coexistence lines represents the two-phase region.
The L-S A -interface is shown in Fig. 5 for T * = 1.3.In the phase diagram in Fig. 4(a) the corresponding two coexisting bulk states are marked by black dots (•).Panels (a) and (b) show the packing fraction profile η(z) along the interface normal and the orientational order parameter profile S 2 (x), respectively.The black dashed vertical line in panel (a) marks the position z η of the Gibbs dividing surface, which is defined by Eq. ( 23).Correspondingly, the black dashed vertical line in panel (b) marks the position z S2 (Eq.( 24)).Apparently, the two interface positions z η and z S2 , which are related to the interfacial transition in the structure and in the orientational order, respectively, differ from each other.In  4(a).The tilt angle is α = 0, i.e., the smectic layer normal n n n = ẑ ẑ ẑ is parallel to the interface normal (see Fig. 3).For z/R > 0 the last layers of the SA phase are visible, in which the particles are still well aligned with the z-axis, indicated by large values of the orientational order parameter S2(z/R) > 0.8 within these layers.For z/R < 0 the layer structure of the density dies out rapidly and the orientational order vanishes as well.Ultimately, the isotropic bulk limit will be approached for z → −∞.However, already for z/R < −10 the profiles have de facto reached their bulk limits in the isotropic liquid L. The black dashed lines refer to the interface positions zη and zS 2 , respectively, calculated via Eqs.( 23) and (24).The difference (zη − zS 2 )/R ≈ 2.45 − 1.66 = 0.79 between the two interface positions is considerably smaller than the smectic layer spacing d/R ≈ 4.28 L/R = 4. Therefore the orientational order of the SA phase vanishes within the last smectic layer while approaching the isotropic liquid L. centers, i.e., D = 0, while the green curve shows data points for D/R = 1.8.The blue curve corresponds to a system of ordinary (uncharged) liquid crystals (OLCs) described by L/R = 4, ε R /ε L = 2, and γ/(Rε 0 ) = 0.The phase diagram for OLCs is not shown here; it is presented in Fig. 4(a) of Ref. [11].Within the considered temperature ranges, in all three cases the differences are at most as large as the length of the particle diameter R, which in turn is much smaller than the smectic layer spacing d/R ≈ 4.3 which is comparable to the particle length L, because the particles within the smectic layers are well aligned with the z-direction, indicated by S 2 (z) > 0.8 in the centers of the smectic layers.Thus, the small size of the differences shows that in these cases the transition in the orientational order and in the fluid structure go along with each other.As soon as the smectic layer structure FIG. 6.
The difference (zη − zS 2 )/R between the Gibbs dividing surface position zη (Eq.( 23)), and the surface position zS 2 (Eq.( 24)), which corresponds to the transition of the orientational order at the interface, are shown for three cases.First, an ordinary (uncharged) liquid crystal (OLC; blue curve); second, ILCs with charges in their center, i.e., D = 0 (violet curve); and, third, ILCs with charges at the tips, i.e., D/R = 1.8 (green curve).Here, the smectic layer normal n n n = ẑ ẑ ẑ is parallel to the interface normal, i.e., α = 0.In all cases studied, the differences (zη − zS 2 )/R are smaller than the smectic layer spacing d L, which for the SA phase is comparable to the particle length L/R = 4. Thus, the loss of orientational order occurs within the last smectic layer before approaching the isotropic liquid L. The inset shows data for the L-SAW -interface, which are accessible for D/R = 1.8 at sufficiently low temperatures T * .Although the difference (zη − zS 2 )/R is enlarged for 0.7 < T * ≤ 0.9, it is still considerably smaller than the layer spacing d/R ≈ 7.5 and decreases rapidly upon decreasing the temperature T * .Hence, for α = 0, the orientational order of the smectic-A phase, either SA or SAW , vanishes directly with the disappearance of the layer structure at the interface.
dies out, the orientational order vanishes as well.
While for ILCs with charges in their center, within the considered temperature range, only L-S A -coexistence is observable (see Fig. 4(a)).For ILCs with the charges at the tips, such as in the case L/R = 4, ε R /ε L = 2, γ/(Rε 0 ) = 0.045, λ D /R = 5, and D/R = 1.8, the bulk phase behavior changes significantly at low temperatures, i.e., for T * < 1.23.The bulk phase diagram in Fig. 4(b) shows that in this case the distinct smectic-A phase S AW occurs for intermediate mean packing fraction η 0 .The S AW phase is characterized by an alternating layer structure of smectic layers with a majority of particles being oriented parallel to the smectic layer normal n n n and a minority of particles localized in secondary layers which prefer orientations perpendicular to the smectic layer normal.Due to this alternating layer structure the smectic layer spacing d/R ≈ 7.5 is increased for the S AW phase.A detailed discussion of the structural and orientational properties of this new and peculiar smectic-A phase, in particular concerning the bulk density and the orientational order parameters profiles, is given in Ref. [11].
In Fig. 7 the L-S AW -interface profiles η(z) and S 2 (z) are shown for α = 0 and T * = 0.9.In the phase diagram in Fig. 4(b) the corresponding coexisting bulk states are marked by red dots (•).On the right hand side of Fig. 7 the alternating layer structure of the bulk S AW phase is evident.In the main layers the majority of the particles (η(z) > 2) has orientations parallel to the z-axis (S 2 (z) > 0.8) and in the secondary layers, formed by less of them (η(z) ≈ 0.6), the particles prefer orientations perpendicular to the z-axis (S 2 (z) < 0).For the L-S AWinterface the difference (z η − z S2 )/R ≈ 2.6 of the two interface positions is increased compared to the L-S Ainterface (see Fig. 6), because the smectic layer spacing d/R ≥ 7.5 in the S AW phase is enlarged, too.As before, the orientational order directly vanishes with the disappearance of the layer structure.Furthermore, the inset in Fig. 6 shows that (z η −z S2 )/R decreases upon lowering the temperature.Thus the difference z η − z S2 becomes smaller relative to the layer spacing d, such that the direct vanishing of the orientational order associated with the disappearance of the layer structure is observable for the whole temperature range considered here.For z → −∞ the isotropic bulk liquid L and for z → ∞ the bulk of the SA phase is approached.The decaying red stripes at the upper part of these plots show the tails of the smectic layers located at x/R ≈ 0, ±d/R, ±2d/R where d/R ≈ 4.28 is the smectic layer spacing.The black dashed lines mark the interface positions zη and zS 2 calculated via Eqs.( 23) and ( 24), while the white dotted lines mark the interface contours zη(x) and zS 2 (x) calculated via Eq.( 25).The difference (zη−zS 2 )/R ≈ 0.58−(−1.51)= 2.09 is larger than the particle diameter R, which is the relevant geometrical property of the particles at this interface, because for α = π/2 the particles in the SA layers are well aligned with the x-axis and therefore they are oriented perpendicular to the direction of the interface normal.The orientational order of the smectic-A phase persists up to a few particle diameters into the liquid phase, unlike the case α = 0, in which the disappearance of the layer structure causes a direct vanishing of the orientational order within the last layer (see Figs. 5-7).
B. Interface normal perpendicular to the smectic layer normal (α = π/2) For α = π/2 the interface normal and the smectic layer normal are perpendicular to each other.The smectic layer normal points into the x-direction and the interface normal into the z-direction (see Fig. 3).The associated L-S A -interface at T * = 1.3 for an ILC system with the charges concentrated at the center, described by the parameter set L/R = 4, ε R /ε L = 2, γ/(Rε 0 ) = 0.045, λ D /R = 5, and D = 0, is shown in Fig. 8.The corresponding bulk phases are given by the state points marked by black dots (•) in the phase diagram in Fig. 4 23) and ( 24), while the white dotted lines mark the interface contours zη(x) and zS 2 (x), which have been calculated via Eq.( 25).The differences cated by large values of the orientational order parameter S 2 (x, z) > 0.8 in the layers.
The black dashed lines in Fig. 8 show the interface positions z η and z S2 calculated from Eqs. ( 23) and ( 24), while the white dotted lines show the interface contours zη (x) and zS2 (x) obtained from Eq. ( 25).The contour lines zη (x) and zS2 (x) at the centers of the tails of the smectic layers, e.g., at x/R ≈ 2.14, are very close to z η and z S2 , respectively.This suggests that the two distinct definitions of the interface positions, i.e., using either Eqs.( 23) and (24) or Eq. ( 25), are consistent with each other, because the majority of the particles in the smectic phase are located close to the centers of the smectic layers.In Fig. 8(a) the packing fraction interface contour zη (x) exhibits discontinuities for lateral positions x at which the smectic bulk packing fraction η SA (x) := η(x, z → ∞) takes the same value η L = η(x, z → −∞) as in the isotropic liquid L, i.e., η SA (x) = η L .Thus, the numeri-cal calculation of the Gibbs dividing surface via Eq.( 25) leads to a divergence due to the vanishing denominator.This can be considered as an artifact, which, however, occurs only at the particular lateral positions x.Nevertheless, the benefit of considering zη (x) and zS2 (x) as interface positions is their dependence on the lateral coordinate x.In particular, for the case of the L-S AWinterface it is necessary to consider zη (x) and zS2 (x) in order to study the interface at the main layers and at the secondary layers separately (see below).
Interestingly, if the layer normal and the interface normal are perpendicular, one observes a significant difference (z η −z S2 )/R ≈ 0.72−(−1.76)= 2.48 between the interface position z η , corresponding to the structural transition, and z S2 corresponding to the transition in the orientational order between the coexisting phases.Hence, the alignment of the particles with the x-axis persists a few particle diameters deeper into the liquid phase L than the layer structure of the S A phase is maintainedunlike in the case α = 0, i.e., in which the smectic layer normal is parallel to the interface normal, for which the orientational order directly vanishes when the smectic layers disappear (see Sec. III A).We note, that the vanishing of the orientational order significantly after (upon approaching the interface from the orientational ordered phase) the structural transition associated with the density profile, has already been observed previously [18] in the case of the interface between an isotropic liquid and a plastic-triangular crystal (PTC).
For the type of ILCs with the charges at the tips, at low temperatures the new wide smectic-A phase S AW can be observed (see Fig. 4(b)).It is characterized by an alternating structure of layers in which the particles are predominantly parallel to the layer normal n n n = x x x (like in the S A phase) and layers of particles which are preferentially perpendicular to the layer normal.The free interface formed between the isotropic liquid L and the S AW phase for T * = 0.9 and α = π/2 is shown in Fig. 9.The red regions in Fig. 9(a) show the layers of particles (at x = 0 and x/R ≈ ±d/R = ±7.5)being parallel to the layer normal, while in between (at x/R ≈ ±d/(2R) = ±3.75) in light blue color the secondary layers are visible.The dark blue color at x/R ≈ ±d/(2R) = ±3.75 in panel (b) shows that the orientational order parameter S 2 (x, z) is negative at the location of the secondary layers, because there the particles are preferentially perpendicular to the layer normal.The interface at the parallel layers behaves very much like the L-S A interface, as can be inferred from the (white) interface contours zη (x/R = 0, ±7.5)/R ≈ 0.81 and zS2 (x/R = 0, ±7.5)/R ≈ −1.83 which show that the orientational ordering of the S AW phase persists into the liquid phase L for a few particle diameters.This is also apparent from the interface positions z η /R ≈ 1.0 and z S2 /R ≈ −2.3, depicted by the black dashed lines in Fig. 9. Conversely, at lateral positions x/R ≈ d/(2R) = ±3.75associated with the centers of the intermediate layers, it turns out that the orientational order undergoes the transition before the layer structure vanishes if one approaches the interface from the S AW side (z S2 (x/R = ±3.75)/R≈ 3.39 and zη (x/R = ±3.75)/R≈ −0.34; in order to guide the eye the magenta dots (•) in Fig. 9 mark these positions).This behavior is opposite to the above one and is presumably related to the fact, that the secondary layers consist of particles being preferentially perpendicular to the layer normal; unlike the particles in the main layers of the S AW phase or the particles in the S A layers, these particles do not align with the layer normal n n n = x x x.Instead they are avoiding an orientation parallel to it.While the transition across the L-S A interface -from alignment with the layer normal towards an isotropic orientational distribution -results in an increase of the effective particle diameter in the y-and z-direction, for the secondary S AW layers the effective diameter is decreased from the S AW phase towards the isotropic liquid L. In Fig. 9 there are discontinuities in the (white) interface contour lines zη (x) and zS2 (x), as in Fig. 8.These discontinuities occur at lateral positions x at which the packing fraction η(x, z → ±∞) or the orientational order parameter S 2 (x, z → ±∞) take the same value in the isotropic bulk, i.e., for z → −∞, as in the S AW bulk, i.e., for z → ∞.

C. Asymptotic behavior
In this section we discuss how the interface profiles of the packing fraction η(r r r) and the orientational order parameter S 2 (r r r) attain their respective values η L and S 2,L in the bulk liquid L. In Fig. 10  .In order to elucidate the view angle on these 3-dimensional logarithmic plots, the interface profiles η(x, z) and S 2 (x, z) are shown in addition as contour plots (see Fig. 8) at the base of the respective plot.
Interestingly, while for D = 0 the periodic structure of the profiles η(x, z) and S 2 (x, z) in x-direction is clearly apparent also in the decays ln |η(x, z) − η L | and ln |S 2 (x, z) − S 2,L | far away from the L-S A -interface (z/R < −20 in Figs.10(a) and (b)), for D/R = 1.8 (panels (c) and (d)) the decays vary only little as function of x.This distinct behavior can be a signature of the respective molecular charge distributions, because if the charges are localized at the centers of the molecules, due to the layer structure in the S A phase the charges are also localized at the centers of the smectic layers, while for D/R = 1.8 the charges are less localized along the lateral direction x.Close to the interface (z/R > −20) the structure is very similar in both cases and, as will be discussed later, it is the hard-core repulsion which is the dominant contribution here.
Turning the view parallel to the x-axis, one obtains Accordingly, the decay of the envelope is determined by the Debye screening length λ D /R = 5, highlighted by the orange lines in Fig. 11.It is worth mentioning that a DFT study [46] of the asymptotic behavior of the liquidvapor interface has yielded, unlike the present findings, a decay length l b larger than the Debye screening length λ D for a hard sphere system with additional Yukawa in-teraction.While in the present study the Yukawa potential is purely repulsive, in Ref. [46] using an attractive Yukawa potential is indispensable, because a sufficiently strong attraction is needed for liquid-vapor coexistence to occur.Interestingly, the asymptotic behavior of the orientational order parameter at far distances, i.e., for z/R < −60, differs from the electrostatic decay and another regime (highlighted by blue lines in Fig. 11) with a larger decay length ξ GB /R ≈ 10 sets in.This longer-ranged decay is due to the Gay-Berne interaction U GB which is verified by calculating the interface profile for an ordinary liquid crystal (OLC) without charges (compare the insets of panels (a) and (b) of Fig. 11).For the OLC, at far distances, i.e., z/R < −30, the same large decay length ξ GB /R ≈ 10 is observed.However, the amplitudes of the decay of the packing fraction and of the orientational order parameter differ significantly.(The blue line in panel (a) intersects the ordinate at ln |η − η L | ≈ −25, whereas the blue line in (b) intersects the ordinate at ln |S 2 − S 2,L | ≈ −20.)For D = 0, it turns out that for the orientational order parameter the crossover from the electrostatic decay towards the Gay-Berne decay occurs at z/R ≈ −67 (this position is marked by the red arrow in Fig. 11(b)), whereas for the case D/R = 1.8 the crossover occurs at z/R ≈ −45 (see the red arrow in Fig. 11(d)).Ultimately, the larger Gay-Berne decay length ξ GB /R ≈ 10 will also become apparent in the decay profile of the packing fraction.However, due to the smaller amplitude of the Gay-Berne decay of the density compared with the decay of the orientational order parameter (compare the insets in Figs.11(a) and (b)), in the present case the crossover occurs further away from the interface (in Fig. 11(a) the intersection of the orange line and the blue line is located at z/R ≈ −121 (not visible) and in Fig. 11(c) at z/R ≈ −97 (also not visible)).However, at very far distances z/R < −80, the magnitudes ln |η − η L | −25 are very small and cannot be resolved numerically.For this reason, in Figs.11(a) and (c) crossovers from the electrostatic regime to the Gay-Berne regime are not shown.
We note that, although the Gay-Berne potential U GB decays algebraicly ∝ (r 12 /R) −6 (see Eq. ( 2)), here the Gay-Berne decay is exponential, because solving the Euler-Lagrange equation in Eq. ( 8) requires the evaluation of the ERPA contribution βψ ERPA of the effective one particle potential βψ (see Eqs. ( 14) and ( 19)).The numerical calculation of this integral (which extends over the whole volume V of the system) requires a truncation in terms of a cut-off distance of the integral which leads to an exponential decay of this contribution, instead of a power law decay ∝ (z/R) −3 [46][47][48], as it is expected for the full Gay-Berne potential U GB .(The exponent 3 arises because the asymptotic behavior of an interfacial density profile, generated by long-ranged forces, varies proportional to the corresponding (total) potential, which acts on a test particle at a distance z from the interface and which is due to the pair interaction between the particles in one of the two coexisting phases (which are separated by the considered interface) and the test particle.Thus, via an integration of the Gay-Berne pair interaction, which decays ∝ (r 12 /R) −6 , over a half-space, one obtains the corresponding total potential decaying ∝ (z/R) −3 [47][48][49].) For z/R → −∞ the algebraic decay of the Gay-Berne interaction potential always dominates the exponential decay due to the screened electrostatic interaction, independent of the relative strength of the electrostatic and the Gay-Berne interaction potential.A variation of their relative strength γ/(Rε 0 ) would only lead to a shift of the location of the corresponding crossovers in the density and the order parameter profiles (see the red arrows in Fig. 11) caused by altering the amplitudes of the respective decays of the two interactions.
Close to the interface, i.e., for −20 < z/R < −5, in the insets of Fig. 11 one can observe an exponential decay with a decay length ξ PL /R ≈ 1.9 (depicted by the black lines) which arises from the pure hard-core Parsons-Lee contribution βψ PL .Thus ξ PL can be identified as the isotropic-liquid bulk correlation length of the pure hard-core system.Interestingly, while the hard-core correlation length ξ PL is observable in OLCs -within both the η and the S 2 profiles (at distances z/R ∈ [−20, −5] the respective decays closely follow the black lines which depict the hard-core decay in the insets of Figs.11(a) and (b)) -, for ILCs this decay is visible only within the S 2 profile.Only for the S 2 profile the amplitude of the hard-core decay is large enough, such that the hard-core correlation length ξ PL is observable before the electrostatic decay becomes dominant.The insets in Figs.11(c) and (d) show the interface profiles calculated for the pure hard-core system (βψ := βψ PL ) in order to verify that the decay close to the interface, i.e., for −20 < z/R < −5, is governed by the hard-core interaction.
Finally, it is worth mentioning that for all cases shown in Fig. 11, the structural and orientational properties close to the interface, i.e., for z/R > −10, agree very well.Thus, it is the hard-core interaction which determines the structural and orientational properties close to the interface, while the electrostatic and the Gay-Berne contributions dominate further away from the interface.At intermediate distances electrostatics dominates the decay of the interface profiles whereas far away from the interface ultimately the attractive Gay-Berne interaction dominates.Furthermore, the positions of the crossovers between these regimes are distinct for the packing fraction profile and the orientational order parameter profile.

D. Tilted interfaces
In this section we discuss the dependence of the structural and orientational properties of the liquid-smecticinterface on the tilt angle α.In Fig. 12 the L-S A -FIG. 12.The L-SA interface profiles η(x, z) (Eq.( 20)) and S2(x, z) (Eq.( 21)) for α = π/4 and T * = 1.3 are shown.
Here, an ILC with charges localized at its center is considered (L/R = 4, εR/εL = 2, γ/(Rε0) = 0.045, λD/R = 5, and D = 0).For z → −∞ the isotropic liquid bulk L is approached and for z → ∞ the bulk of the SA phase is attained, i.e., the interface normal is parallel to the z-axis.The red stripes at the top of the contour plots show the tails of the smectic layers.The black dashed lines mark the interface positions zη/R ≈ 5.56 and zS 2 /R ≈ 2.79 calculated via Eqs.( 23) and (24).Similar to the case α = π/2 (see Fig. 8), to a certain extent the orientational order persists into the liquid phase L.
interface profiles η(x, z) and S 2 (x, z) are shown for the reduced temperature T * = 1.3 (see the black dots (•) in Fig. 4(a)) and α = π/4.Here, we consider the case of ILCs with the charges localized in the center (L/R = 4, ε R /ε L = 2, γ/(Rε 0 ) = 0.045, λ D /R = 5, and D = 0).Like in the case α = π/2, (see Sec. III B) i.e., the interface normal and the smectic layer normal n n n = x x x are perpendicular, a persisting orientational order can be observed at the interface: The structural transition occurs at z η /R ≈ 5.56, whereas the transition in the orientational order between the two phases takes place at z S2 /R ≈ 2.79 which is a few diameters deeper in the isotropic liquid.
In Fig. 13 the interfacial tension Γ * (α) given by Eq. ( 26) and the distance z η − z S2 between the interface positions associated with the mean packing fraction η 0 (x) and the mean orientational order parameter S 20 (x) are shown as function of the tilt angle α.In Fig. 13(a) the case of the L-S A -interface for ILCs with the charges at their center is considered for T * = 1.Both the interfacial tension Γ * (α) (black dots, •) and the distance z η − z S2 (orange dots, •) exhibit a global minimum at α = 0 and a second, local minimum at α = π/2.Thus, the equilibrium tilt angle α eq = 0 corresponds to the configuration in which the interface normal and the smectic layer normal n n n = ẑ ẑ ẑ are parallel, whereas the corresponding perpendicular orientation α = π/2 is metastable.This increase in the interfacial tension Γ * below α = π/2 suggests that the configuration, in which the interface normal and the layer normal are orthogonal, should be observable without resorting to any external stabilizing field which could be provided, e.g., by a suit- ably structured substrate.This metastability of the tilt angle α = π/2 can be checked also via computer simulations.Interestingly, the increase of the interfacial tension below α = π/2 is accompanied by an increase in the distance z η − z S2 , suggesting that maintaining the local orientational order in the isotropic liquid beyond the smectic layers costs free energy.Consistently, in the case α eq = 0, for which the orientational order vanishes directly with the disappearance of the smectic layers, the cost in free energy is lowest.Apparently, for α = 0 the interfacial tension Γ * (α = 0) ≈ 0.006 is significantly smaller than for all other angles α shown in Fig. 13(a).For technical reasons we did not study small tilt angles α > 0 and FIG.14. Same as Fig. 12.Here, the L-SAW interface profiles η(x, z) (Eq.( 20)) and S2(x, z) (Eq.( 21)) are shown for α = π/3 and T * = 0.9.To this end, an ionic liquid crystal with charges at the tips is considered (L/R = 4, εR/εL = 2, γ/(Rε0) = 0.045, λD/R = 5, and D/R = 1.8).For z → −∞ the isotropic liquid bulk L is approached and for z → ∞ the bulk of the SAW phase, i.e., the interface normal is parallel to the z-axis.The transition in the structure occurs at zη/R ≈ 2.28 and the transition in the orientational order does so at zS 2 /R ≈ −0.68.
hence cannot comment on the functional form of Γ * (α) for 0 < α < π/6 in the case D/R = 0 or for 0 < α < π/4 in the case D/R = 1.8.This is indicated by connecting the data points at α = 0 and π/6 by dashed lines.(For the same reason, in (b) the data points at α = 0 and π/4 are connected by dashed lines.)It has been pointed out in Sec.II B, that due to the crossover at the tilt angle α = 0 from a periodic system towards one which is translationally invariant in lateral direction x, the integration domain V d for evaluating the coefficients Q i (r r r) (see Eq. ( 16)) is not continuously evolving at α = 0.For α > 0 it is a slice of length d x = d/ sin(α) in xdirection, while for α = 0 it is the subsystem of length d in z-direction at position r r r. (For α = 0 the extent in x-and y-direction is arbitrary due to the translational invariance in lateral direction.)In order to describe a continuous variation of the interfacial tension Γ * (α) for all tilt angles α ∈ [0, π/2], one thus needs to consider a different approach, which does not rely on a projected density and thereby on the direction of the bulk smectic layer normal n n n throughout the whole interface structure.Nonetheless, our above approach still allows one to compare the interfacial tension Γ * (α) for the extreme cases α = 0 and π/2, thus predicting which one of the two is preferred.Furthermore, our approach provides an understanding of the local increase in Γ * (α) below α = π/2, as one observes an increasing distance z η − z S2 between the transition in the structural and the orientational order at the interface.
Figure 13(b) shows data for the L-S AW -interface at T * = 0.9 for ILCs with charges located at the tips.Around α = π/2 the interfacial tension (black squares, ) is a rather flat function of α taking values around Γ * ≈ 0.07.The slight variations in Γ * for α ∈ [π/4, π/2] might be caused by the numerical evaluation of Eq. ( 8) which has to be done separately for each tilt angle α.Consistently, the distance z η − z S2 (orange squares, ) does not vary much as function of the tilt angle α.As above, the equilibrium tilt angle α eq = 0 corresponds to the configuration in which the interface normal and the smectic layer normal n n n = ẑ ẑ ẑ are parallel.Finally, in Fig. 14, we show the contour plot of the L-S AW -interface for α = π/3 and T * = 0.9 for an ILC system with D/R = 1.8, illustrating the structure of this type of interface.

IV. SUMMARY AND CONCLUSIONS
Free interfaces in systems composed of ionic liquid crystals (ILCs) have been studied within density functional theory (see Sec. II B).In particular, the discussion has been focused on two kinds of ionic liquid crystals: first, ILCs with the charges localized at the center of the molecules, i.e., D = 0 (see Figs. 1 and 2), and, second, ILCs with the charges at the tips of the molecules, i.e., D/R = 1.8.All other model parameters, i.e., L/R = 4, ε R /ε L = 2, γ/(Rε 0 ) = 0.045, λ D /R = 5, are identical in both cases.Therefore the two kinds differ solely by the charge distribution within the molecules.
For D = 0 coexistence between the isotropic liquid L and the ordinary smectic-A phase S A can be observed at a sufficiently large mean packing fraction η 0 (see Fig. 4(a)).The S A phase is characterized by a layered structure in the direction of the smectic layer normal n n n with a smectic layer spacing d ≈ L comparable to the particle length L. Within the smectic layers the particles are well aligned with the smectic layer normal.The phase behavior of ILCs is altered by varying the molecular charge distribution, as can be inferred from comparing the case D = 0 (i.e., charges at the center) and D/R = 1.8 (i.e., charges at the tips, see Fig. 4(b)).At sufficiently low temperatures a new smectic-A phase has been observed, which is referred to as the S AW phase [11].The S AW phase shows an alternating structure of layers with the majority of the particles being oriented parallel to the smectic layer normal n n n and the minority of the particles localized in secondary layers which prefer orientations perpendicular to n n n.Due to the alternating layer structure, the smectic layer spacing d/R ≈ 7.5 in the S AW phase is increased compared with the spacing in the S A phase.
For a parallel orientation of the smectic layer normal n n n = ẑ ẑ ẑ and the L-S A -interface normal, i.e., for α = 0 (see Fig. 3), it turns out that the interface locations z η and z S2 , associated with the transition in the structural and in the orientational order, respectively, are very close to each other (see Fig. 5).In fact, Fig. 6 shows that for the whole temperature range considered here, the difference z η − z S2 < d in the two interface positions is smaller than the smectic layer spacing d.Hence, for α = 0 the orientational order vanishes within the last smectic layer at the L-S A -interface.Concerning the interface positions, Fig. 6 demonstrates that ILCs with D/R = 1.8 and ordinary (uncharged) liquid crystals with L/R = 4 and ε R /ε L = 2 exhibit qualitatively the same results.Considering the L-S AW -interface (see Fig. 7) one observes an increase in z η − z S2 , but it remains significantly smaller than the smectic layer spacing d/R ≈ 7.5.Thus, for α = 0 it turns out that the loss of orientational order coincides with the disappearance of the layer structure of the respective smectic-A phase at the interface towards the isotropic liquid.This holds for all parameter values studied here.
Interestingly, for α = π/2, i.e., changing the relative orientation of the smectic layer normal n n n = x x x and the interface normal such that they are perpendicular to each other, leads to qualitative changes in the interfacial properties: a periodic structure of the interface in lateral direction x can be observed, which is a direct consequence of the periodicity in the bulk smectic-A phase with the smectic layer spacing d (see Figs. 3, 8, and 9).For the L-S A -interface (see Fig. 8) one observes considerable differences (z η −z S2 )/R 2 between the interface positions.Thus, the (nearly) parallel orientations of particles in the S A layers persists a few particle diameters R into the liquid phase L, unlike the case α = 0, for which the orientational order vanishes directly with the breakdown of the S A layer structure at the interface, i.e., within the last smectic layer.Due to the periodicity in (lateral) xdirection, in the case α = π/2 one indeed observes a qualitative change in the structure of the L-S AW -interface compared to the L-S A -interface.While at the tails of the S AW main layers the interface also features an orientational order which continues further into the liquid phase L than the layer structure ((z η (x) − zS2 (x))/R ≈ 2.6).
For the secondary layers it is the layer structure that persists deeper into the L phase than the orientational order ((z η (x) − zS2 (x))/R ≈ −3.73).The opposite behavior at the main, respectively secondary, layers is presumably driven by the orientational properties of the respective kinds of layers: in the main layers the particles are well aligned with the smectic layer normal n n n = x x x and therefore show an effective diameter in the y-z-plane which is comparable to the particle diameter R.However, in the secondary layers (here with S 2 (x, z) < 0) the particles avoid orientations parallel to the x-axis, giving rise to an considerably larger effective radius.Upon approaching the liquid phase L, this effective radius increases for the main layers of the S AW phase, whereas it decreases for the secondary layers.
In Sec.III C the asymptotic behavior of the interface profiles has been studied.In particular, in Figs. 10 and 11 the L-S A -interface for α = π/2 has been considered for the two ILC systems with D/R = 0 and 1.8.For D = 0, i.e., with the charges being localized at the center, the periodic structure of the interface is apparent from the quantities ln |η(x, z) − η L | and ln |S 2 (x, z) − S 2,L |, showing the logarithmic deviations of the profiles η(x, z) and S 2 (x, z) from their respective liquid bulk values η L FIG. 15.Same as Fig. 14.Here, the L-SAW interface profiles η(x, z) (Eq.( 20)) and S2(x, z) (Eq.( 21)) are calculated for α = π/3 and T * = 0.9 by using the projected density containing odd Fourier-modes up to second order.The profiles are qualitatively equivalent to those obtained without using the odd modes in the projected density ̺(r r r, ω ω ω) (see Eq. ( 11)).In agreement with the results shown in Fig. 14 one observes an orientational order (within the main layers of the SAW phase) persisting up to a few particle diameter R into the liquid phase ((zη − zS 2 )/R = 0.24 − (−2.23) = 2.47).
second-order ones within the projected density ̺(r r r, ω ω ω).Including these terms, ̺(r r r, ω ω ω) takes the following modified form: (A3) The coefficients Q i with i = 6, • • • , 9 vanish for the considered bulk phases, because smectic-A phases exhibit mirror-symmetry with respect to the layer center.
In general, at interfaces they do not vanish.In order to compare the corresponding interface profiles η(r r r) =  20) and (21), respectively; S is the full solid angle) obtained from solving the Euler-Lagrange equation, i.e., Eq. ( 8), by using the projected density without the odd terms given by Eq. ( 11) and the projected density containing these terms, i.e., by using Eq.(A2), the case α = π/3 (see Eq. ( 15)) and the L-S AW -interface shown in Fig. 14 are considered again.In Fig. 15 the two respective profiles are shown by using Eq.(A1): for both η(r r r) and S 2 (r r r) there are no qualitative differences compared with Fig. 14.The interface positions z η /R ≈ 0.24 and z S2 /R ≈ −2.23 are shifted in z-direction compared to the results shown in Fig. 14.But their distance (z η − z S2 )/R ≈ 2.47 is comparable to the previous results ((z η − z S2 )/R ≈ 2.96 in Fig. 14).Hence, in qualitative agreement with the results shown in Fig. 14 one observes a persisting orientational order (within the main layers of the S AW phase) up to a few particle diameters R into the liquid phase.

FIG. 2 .
FIG. 2. Contour-plots of the pair potential U for |r r r12| ≥ Rσ in the x-z-plane for four cases of particles with fixed length-tobreadth ratio L/R = 4 and fixed orientations.In each panel the centers of both particles lie in the plane y = 0.In order to illustrate the orientations of the ellipsoids, they have been included in the plots at contact with relative direction r r r12 = x x x.The set of points at contact in the x-z-plane is illustrated by the black curve, and the centers of the particles are shown by small black dots.Panel (a): uncharged liquid crystal with εR/εL = 2. Panel (b): uncharged liquid crystal with εR/εL = 4.With this choice the anisotropy of the potential is increased slightly.Panel (c): ILC with εR/εL = 2, D/R = 0, λD/R = 5, γ/(Rε0) = 0.25.Panel (d): ILC with εR/εL = 2, D/R = 1.8, λD/R = 5, γ/(Rε0) = 0.25.In (c) and (d) the loci of the charges are indicated as blue dots.The salmon-colored area is the excluded volume for given orientations of the two particles.

2
the considered interfaces between the isotropic liquid L and the smectic-A phases S A or S AW exhibit periodic structures in the x-direction with periodicity d x = d/ sin(α).Here, V d is a slice of length d x in x-direction with a vanishing extension in z-direction centered at position r r r, i.e., T (r r r − r r r ′ ) = Θ(d x /2 − |x − x ′ |)δ(z − z ′ ) where Θ(x) and δ(x) are the Heaviside step function and the Dirac delta function, respectively.The index d of the integration domain V d indicates that V d corresponds to a region which is specified by the periodicity d.

FIG. 4 .
FIG.4.Bulk phase diagrams for (a) ionic liquid crystals with L/R = 4, εR/εL = 2, γ/(Rε0) = 0.045, λD/R = 5, and D = 0 and (b) with L/R = 4, εR/εL = 2, γ/(Rε0) = 0.045, λD/R = 5, and D/R = 1.8.For D = 0, i.e., the charges being concentrated in the center of the molecules, solely a first-order phase transition from the isotropic liquid phase L to the ordinary smectic-A phase SA occurs at sufficiently high mean packing fractions η0.The ordinary smectic-A phase SA is characterized by a layer structure with smectic layer spacing d/R ≈ 4.3L/R = 4 comparable with the particle length L. The particles in the layers are well aligned with the layer normal n n n.In panel (b), i.e., for D/R = 1.8 (the charges being located at the tips of the molecules), another smectic-A structure, referred to as the SAW phase can be observed at low reduced temperatures T * .The SAW phase exhibits an alternating structure, consisting of primary layers of particles being parallel to the layer normal and secondary layers in which the particles prefer to be perpendicular to it.This leads to an increased layer spacing d/R ≥ 7.5.The black dotted line in panel (b) marks the triple point at T * ≈ 1.23 for which the isotropic liquid L, the ordinary smectic-A phase SA, and the SAW phase are in three-phase-coexistence.A detailed description of the structural properties of the smectic-A phases SA and SAW , including illustrations of their microstructure, are provided in Ref.[11].The black dots (•) in panel (a), respectively the red dots (•) in (b), mark the coexisting bulk states at the reduced temperature T * = 1.3, respectively 0.9, imposed as boundary conditions for the free interfaces shown in Figs.5 and 7-9.

Fig. 6 ,FIG. 5 .
FIG. 5.The L-SA-interface profile of the packing fraction η(z), panel (a), and the orientational order parameter S2(z), panel (b), are shown for an ionic liquid crystal with L/R = 4, εR/εL = 2, γ/(Rε0) = 0.045, λD/R = 5, and D = 0, i.e., the charges are concentrated in the center of the molecules.The free interface between the isotropic liquid L (imposed as boundary condition for z → −∞) and the ordinary smectic-A phase SA (i.e., z → ∞) is considered for the reduced temperature T * = 1.3.The corresponding coexisting bulk states are marked by the black dots (•) in the phase diagram in Fig.4(a).The tilt angle is α = 0, i.e., the smectic layer normal n n n = ẑ ẑ ẑ is parallel to the interface normal (see Fig.3).For z/R > 0 the last layers of the SA phase are visible, in which the particles are still well aligned with the z-axis, indicated by large values of the orientational order parameter S2(z/R) > 0.8 within these layers.For z/R < 0 the layer structure of the density dies out rapidly and the orientational order vanishes as well.Ultimately, the isotropic bulk limit will be approached for z → −∞.However, already for z/R < −10 the profiles have de facto reached their bulk limits in the isotropic liquid L. The black dashed lines refer to the interface positions zη and zS 2 , respectively, calculated via Eqs.(23) and(24).The difference (zη − zS 2 )/R ≈ 2.45 − 1.66 = 0.79 between the two interface positions is considerably smaller than the smectic layer spacing d/R ≈ 4.28 L/R = 4. Therefore the orientational order of the SA phase vanishes within the last smectic layer while approaching the isotropic liquid L.

FIG. 7 .
FIG. 7.For α = 0, the L-SAW -interface profiles η(z) and S2(z) are shown for ILCs with charges at the tips (L/R = 4, εR/εL = 2, γ/(Rε0) = 0.045, λD/R = 5, and D/R = 1.8)at the reduced temperature T * = 0.9 (see the red dots (•) in Fig.4(b)).For z → −∞ the isotropic liquid bulk L is approached whereas for z → ∞ the SAW bulk is attained.The difference (zη − zS 2 )/R ≈ 6.31 − 3.72 = 2.59 between the two interface positions is larger than the one of the L-SAinterface (compare Figs.5 and 6) but it is still smaller than the smectic layer spacing d/R = 7.5.Therefore the orientational order of the SAW phase also vanishes within the range of the last smectic layer at the interface.

FIG. 8 .
FIG. 8.The L-SA-interface profiles η(x, z), panel (a), and S2(x, z), panel (b), are shown for T * = 1.3 (see the black dots (•) in Fig. 4(a)) and α = π/2.Accordingly, the smectic layer normal n n n = x x x and the interface normal (parallel to the z-axis) are perpendicular.Here, ILCs with charges at the center are considered, described by the parameter set L/R = 4, εR/εL = 2, γ/(Rε0) = 0.045, λD/R = 5, and D = 0.For z → −∞ the isotropic bulk liquid L and for z → ∞ the bulk of the SA phase is approached.The decaying red stripes at the upper part of these plots show the tails of the smectic layers located at x/R ≈ 0, ±d/R, ±2d/R where d/R ≈ 4.28 is the smectic layer spacing.The black dashed lines mark the interface positions zη and zS 2 calculated via Eqs.(23) and (24), while the white dotted lines mark the interface contours zη(x) and zS 2 (x) calculated via Eq.(25).The difference (zη−zS 2 )/R ≈ 0.58−(−1.51)= 2.09 is larger than the particle diameter R, which is the relevant geometrical property of the particles at this interface, because for α = π/2 the particles in the SA layers are well aligned with the x-axis and therefore they are oriented perpendicular to the direction of the interface normal.The orientational order of the smectic-A phase persists up to a few particle diameters into the liquid phase, unlike the case α = 0, in which the disappearance of the layer structure causes a direct vanishing of the orientational order within the last layer (see Figs.5-7).

FIG. 10 .
FIG. 10.L-SA interface profiles of η(x, z) and S2(x, z) for T * = 10 and α = π/2.Accordingly, the smectic layer normal n n n = x x x and the interface normal (parallel to the z-axis) are perpendicular.Panels (a) and (b) show the logarithmic deviations ln |η(x, z) − ηL| and ln |S2(x, z) − S2,L| of the packing fraction and the orientational order parameter from their bulk values in the isotropic liquid L for an ILC with the charges concentrated at the center of the molecule, i.e., for D = 0. Panels (c) and (d) show ln |η(x, z)−ηL| and ln |S2(x, z)−S2,L| for an ILC with the charges at the tips, i.e., for D/R = 1.8.Note that on the base of each plot the interface profiles η(x, z) and S2(x, z) are shown in order to elucidate the viewing angle on the interface.The local height of the manifold above the base corresponds to the given color code.Interestingly, for D = 0 the periodic structure is still apparent even far away from the L-SA-interface, i.e., z/R < −20, unlike the case D/R = 1.8, for which the profiles are rather flat in lateral direction x.This can be related to the strong localization of charges at the centers of the smectic layers for D = 0, pronouncing the periodic structure, while for D/R = 1.8 the charge sites are spread and less localized along the x-direction.

FIG. 11 .
FIG. 11.The same quantities as shown in Fig.10.Panels (a) and (b) correspond to the case D = 0 presenting ln |η(x, z) − ηL| and ln |S2(x, z) − S2,L|, respectively, whereas panels (c) and (d) correspond to the case D/R = 1.8.However, here the direction of view is parallel to the x-axis so that the manifold from Fig.10is projected onto the plane spanned by the vertical axis and the z axis.Away from the interface, i.e., for z/R < −20, the decay length for ln |η(x, z) − ηL| can be identified as the Debye screening length λD/R = 5 for both cases (a) D = 0 and (c) D/R = 1.8.From the inset in panel (a), which shows ln |η(x, z) − ηL| for the corresponding (uncharged) ordinary liquid crystal with L/R = 4 and εR/εL = 2, it is apparent that the contributions due to the Gay-Berne potential (the asymptotics of which is indicated by the blue line) and due to the hard-core interaction (the asymptotics of which is depicted by the black line) are much weaker than the (screened) electrostatic contribution and do not play a role within the range of ln |η(x, z) − ηL| considered here.(In order to guide the eye, the blue and black lines are also shown in the main plots.Apparently, in (a) and (c) the blue and black lines are far below the respective profiles.)However, for ln |S2(x, z) − S2,L|, i.e., for panels (b) and (d), one observes crossovers -indicated by the intersection of the orange and blue lines at z/R ≈ −67 in (b) and z/R ≈ −45 in (d) (compare the red arrows in the respective plots) -from the electrostatic regime towards the decay governed by the Gay-Berne contribution with decay length ξGB/R ≈ 10.Such crossovers occur within the considered range z/R ∈ [−80, 0], because for the orientational order parameter the amplitude of the decay, due to the Gay-Berne interaction, is larger than for the packing fraction (compare the intersections of the blue lines with the ordinates in panels (a) and (b)).Due to the hard-core interaction, for z/R > −20 the decay length ξPL/R ≈ 1.9 (Parsons-Lee, black lines) is visible for the ordinary liquid crystal in the insets of (a) and (b) as well as for ln |S2(x, z) − S2,L| of the two considered ILCs.(Due to the small amplitudes of the hard-core contributions to ln |η(x, z) − ηL|, for the ILC considered here, this decay has not been observed.)In order to confirm, that the decay length ξPL/R ≈ 1.9 is indeed due to the hard-core interaction, the insets of the panels (c) and (d) show ln |η(x, z) − ηL| and ln |S2(x, z)−S2,L| of the pure hard-core system (βψ := βψPL).Interestingly, ln |η(x, z) − ηL| and ln |S2(x, z) − S2,L| behave very similarly close to the interface, i.e., z/R > −10, for all three kinds of systems studied here.This suggests that the structure and the orientational properties close to the interface are governed by the hard-core interaction which enters into the present DFT approach (see Secs.II A and II B).

4 FIG. 13 .
FIG.13.The (reduced) interfacial tension Γ * (α) (Eq.(26), black line) and the distance zη − zS 2 between the transition in the structural and the orientational order (orange line) as function of the tilt angle α.In panel (a) the L-SA interface at T * = 1 is considered for ILCs with their charges localized at the center (L/R = 4, εR/εL = 2, γ/(Rε0) = 0.045, λD/R = 5, and D = 0 (see Fig.4(a)).There are two minima: the global minimum at the equilibrium tilt angle αeq = 0 (i.e., interface normal and smectic layer normal are parallel) and a local minimum at α = π/2 which shows that the orthogonal orientation of the smectic layer normal and the interface normal is a metastable configuration.The increase of the interfacial tension below α = π/2 is accompanied by an increase of the distance zη − zS 2 .This suggests that maintaining to a certain extent the local orientational order in the isotropic liquid beyond the smectic layers costs free energy.For technical reasons we did not study small tilt angles α > 0. Hence we cannot comment on the functional form of Γ * (α) for 0 < α < π/6 in the case D/R = 0 or for 0 < α < π/4 in the case D/R = 1.8.This is indicated by connecting the data points at α = 0 and π/6 by dashed lines in (a) (see the discussion in the main text of Sec.III D).In panel (b) the L-SAW interface, which is accessible for ILCs with their charges at the tips (L/R = 4, εR/εL = 2, γ/(Rε0) = 0.045, λD/R = 5, and D/R = 1.8), is considered for T * = 0.9 (see Fig.4(b)).Also in this case the equilibrium tilt angle αeq = 0 corresponds to the parallel orientation of the interface normal and the layer normal.Below α = π/2, as function of α the interfacial tension is rather flat, taking the value Γ * ≈ 0.07.Thus, for the L-SAW interface the perpendicular orientation of the interface normal and of the smectic layer normal corresponds to a labile configuration.(Analogously to panel (a), the data points at α = 0 and π/4 in (b) are connected by a dashed line.)We note that Γ * (α) is symmetric around α = π/2, due to the mirror-symmetry of the particles.