Influence of substrate interaction and confinement on electric field induced transition in symmetric block copolymer thin films
(in press, Physical Review E)
Abstract
In the present work, we study morphologies arising due to competing substrate interaction, electric field and confinement effects on a symmetric diblock copolymer. We employ a coarse grained nonlocal CahnHilliard phenomenological model taking into account the appropriate contributions of substrate interaction and electrostatic field. The proposed model couples the OhtaKawasaki functional with Maxwell equation of electrostatics, thus alleviating the need for any approximate solution used in previous studies. We calculate the phase diagram in electric fieldsubstrate strength space for different film thicknesses. In addition to identifying the presence of parallel, perpendicular and mixed lamellae phases similar to analytical calculations, we also find a region in the phase diagram where hybrid morphologies (combination of two phases) coexist. These hybrid morphologies arise either solely due to substrate affinity and confinement or are induced due to the applied electric field. The dependence of the critical fields for transition between the various phases on substrate strength, film thickness and dielectric contrast is discussed. Some preliminary 3D results are also presented to corroborate the presence of hybrid morphologies.
pacs:
I Introduction
Selfassembly of block copolymers has been an actively pursued field of study because of its wide technological implications Bates (1991); Bates and Fredrickson (1990). Depending upon the volume fraction of the components and segregation regime, block copolymers exhibit a range of periodic morphologies such as lamellae, gyroids, cylinders, spheres etc.Bates (1991); Bates and Fredrickson (1990) As we intend to study symmetric diblock copolymers which order into lamellar structure, the rest of the paper will focus only on this particular morphology.
In the absence of any external field, a symmetric diblock copolymer forms domains of lamellar morphology with various degrees of alignment. Practical applications, however, require complete alignment of the microphase separated domains. In general this can be achieved by application of external fields, primary among them being substrate field ThurnAlbrecht et al. (2000); Heier et al. (1999), shear field Albalak and Thomas (1993); Balsara et al. (1994) and electric field Amundson et al. (1991); ThurnAlbrecht et al. (2002); Böker et al. (2002). Since the microphase separation generally takes place on a substrate, it is impractical to assume that the experiments are devoid of surface effects. Typically, interfacial energy difference between the two blocks in contact with a substrate (i.e , where and are the interfacial energies between a monomer A or B and substrate S) can cause surface induced ordering resulting in parallel arrangement of the domains with respect to the surface Fredrickson (1987); Shull (1992); Brown and Chakrabarti (1994). If the copolymer system is confined between two rigid substrate walls, two different interaction cases may be considered, (a) the walls are symmetric i.e. both walls attract the same monomer and (b) antisymmetric walls i.e. both walls attract different monomers Turner (1992); Walton et al. (1994); Kikuchi and Binder (1994); Brown and Chakrabarti (1995); Wang et al. (2000). As shown in Fig. 1 (a) and (b), the system either forms integral or halfintegral number of lamellae, where is the equilibrium lamellar spacing, depending on whether the wall is symmetric or antisymmetric respectively Turner (1992); Walton et al. (1994); Kikuchi and Binder (1994); Brown and Chakrabarti (1995); Wang et al. (2000). If the film thickness is incommensurate with the lamellar spacing, the copolymers are said to be in a frustrated state Lambooy et al. (1994). Frustration can be accomodated by deviation of equilibrium spacing from to which is typical in thicker films where the deviation can span over a large number of lamellae Walton et al. (1994); Lambooy et al. (1994). On the contrary, frustration in thin films is relieved by a change of configuration from parallel to perpendicular lamellae Turner (1992); Walton et al. (1994); Kikuchi and Binder (1994); Brown and Chakrabarti (1995); Wang et al. (2000) as shown in Fig. 1(c). In the present work, we restrict our study to the case of symmetric walls.
Electric field, on the other hand is gaining popularity in guiding selfassembly of block copolymers because of the ease with which it can be applied, especially in thin films Xu et al. (2003, 2004). In contrast to stable parallel configuration in presence of substrate interaction, presence of electric field makes the perpendicular arrangement more stable. The reason for this can be rationalised in terms of dielectric permittivity mismatch between the two monomer components Landau et al. (1984); Pereira and Williams (1999); Tsori (2009). In electrostatics, the electric field isolines prefer to pass through the regions of higher permittivity Landau et al. (1984); Pereira and Williams (1999). Since this is readily achieved in a perpendicular state, the system minimizes its free energy by adopting perpendicular configuration.
A lot of factors can aid or deter the evolving phase morphology in presence of electric field such as the segmental interaction or segregation ( is the Flory interaction parameter which scales inversely with temperature and is the degree of polymerization), the substrate interaction strength and the confinement which can act as geometric barrier. The copolymers can be easily modulated at low segregation owing to weak segmental interaction between the two blocks Tsori and Andelman (2002); Matsen (2006a). The substrate interaction can hinder the alignment process by electric field as then the applied field has to overcome the interfacial interaction Xu et al. (2003). Confinement (film thickness) on the other hand can aid or hinder the alignment process by electric field depending upon its relation with the natural lamella period Tsori and Andelman (2002); Ashok et al. (2001).
The existing literature on parallel to perpendicular lamella transition can be grouped into three categories depending upon the methodology of study: analytical Pereira and Williams (1999); Ashok et al. (2001); Tsori and Andelman (2002), experimental Xu et al. (2003, 2004) and numerical Kyrylyuk et al. (2003); Lyakhova et al. (2006); Lin et al. (2005); Matsen (2006a). While most of the earlier theoretical work in this regard focus on strong segregation limit (SSL), Pereira and Williams (1999); Ashok et al. (2001); Tsori and Andelman (2002) literature on weak segregation limit (WSL) Tsori and Andelman (2002) is relatively sparse. A common finding in the SSL is the presence of three stable configurations, namely, parallel, perpendicular and mixed with respect to the substrate. The mixed structure contains parallel lamellae closer to the substrate and perpendicular lamellae in the middle, depending upon substrate interaction strength and electric field Pereira and Williams (1999); Tsori and Andelman (2002). If , the electric field causes a direct transition from parallel to perpendicular state Pereira and Williams (1999); Tsori and Andelman (2002). However, if, , two critical fields establish in thicker films Tsori and Andelman (2002). At lower critical field, the parallel arrangement changes to mixed morphology, while, on increasing electric field results in perpendicular lamellae Tsori and Andelman (2002). In WSL, the mixed morphology is not predicted and the transition from parallel to perpendicular is direct Tsori and Andelman (2002).
Experimental studies Xu et al. (2003, 2004) focusing on intermediate to strong segregation regimes have reported some observations, which do not completely converge with analytical findings. It is reported that a completely perpendicular arrangement in presence of electric field can only occur when the interfacial energy difference between the two blocks with the substrate is balanced Xu et al. (2003) i.e. . Any mismatch results in a mixed morphology Xu et al. (2003). The discrepancy between the analytical and experimental observations are attributed to the pathway dependence of the alignment process in the presence of an electric field Xu et al. (2003). Analytical studies on the other hand focus only on the final equilibrium configuration. In reality, however, the microdomains can get kinetically trapped in metastable states and a high activation energy may be needed to reach the stable state.
Thus, numerical studies can provide an efficient bridge between the analytical calculation and experiments. Dynamic density functional theory (DDFT) Kyrylyuk et al. (2003); Kyrylyuk and Fraaije (2006); Lyakhova et al. (2006); Xu et al. (2005), selfconsistentfield theory (SCFT) Lin et al. (2005); Matsen (2006a) and cell dynamics simulations (CDS) Pinna and Zvelindovsky (2008) have previously been employed to study electric field induced alignment.
The focus of most numerical studies have either been on mechanism of alignment of ordered domains Böker et al. (2003); Zvelindovsky and Sevink (2003); Schmidt et al. (2005); Pinna et al. (2009); Ruppel et al. (2013) or on orderorder transition from one morphology to the other in bulk samples Ly et al. (2007, 2013); Li et al. (2013); Schmidt et al. (2010); Pester et al. (2015). The effect of substrate interaction and confinement have not been investigated thoroughly which could significantly alter the phase morphologies and critical electric field for transition. Matsen studied the stability of monolayer Matsen (2005) and multilayer Matsen (2006b) lamellae films but did not consider substrate affinity. Lyakhova et al. Lyakhova et al. (2006) and Kyrylyuk et al. Kyrylyuk and Fraaije (2006) did consider substrate interaction, but the effect of film thickness was not considered. Moreover, the presence of any intermediate or mixed phases were not observed. Thus, the primary objective of the present work is to systematically investigate the effect of substrate interaction and confinement on the resulting morphologies. We adopt a coarsegrained CahnHilliard approach to study the selfassembly of symmetric diblock copolymers under competing substrate interaction, electric field and confinement.
Secondly, one of the serious shortcoming of most previous works involve an approximate solution of the electrostatic field. The perturbed solution of the Maxwell equation is based on the assumption of weak fractional variation of the dielectric constant Amundson et al. (1993, 1994). Though appropriate in the proximity of the orderdisorder transition(ODT) temperature, the results can be significantly marred as the segregation increases. Such an assumption is relaxed in the present study by coupling the OhtaKawasaki functional to Maxwell equation to calculate the electrostatic field distribution.
Ii Theoretical Model and Numerical Methods
ii.1 Theoretical Model
The diffuse interface approach to model block copolymers follows the OhtaKawasaki free energy functional Ohta and Kawasaki (1986), which includes a longrange interaction term in addition to short range interaction terms (bulk + interfacial energy) in the CahnHilliard model Cahn and Hilliard (1958) to account for chain interactions. The bulk part of the free energy functional of a diblock copolymer can be written as Bahiana and Oono (1990); Liu and Goldenfeld (1989); Chakrabarti et al. (1989); Choksi, Peletier and Williams (2009),
where the terms in the first integral constitute the shortrange interactions, while the second integral consists of the longrange interaction. is an order parameter that denotes the local concentration difference between the two components, . is the bulk free energy taken to be of the form to account for two stable phases () below the critical temperature. is the gradient energy coefficient which penalises gradients in the order parameter. is a numerical parameter that determines the extent of microphase separation and scales as , where N is the number of segments. G is the Green’s function having the property = . The parameters and are related to the polymer architecture through the relations Ohta and Kawasaki (1986); Choksi, Peletier and Williams (2009),
(2)  
(3) 
where is the Kuhn statistical length or the average monomer space size and is the relative molecular weight which is a measure of the length of A monomer chain compared to the whole macromolecule. The segregation is determined by and as Choksi, Peletier and Williams (2009),
(4) 
The free energy of the domain surface in presence of attracting walls is written as, Schmidt and Binder (1985); Tsori and Andelman (2001)
(5) 
The above expression results from a Taylor series expansion of bare surface energy Schmidt and Binder (1985). The terms and have special physical interpretation. denotes the surface chemical potential difference. A positive value expresses preferential attraction of B component and viceversa Schmidt and Binder (1985); Tsori and Andelman (2001). The term takes into account the deviation of Flory parameter at the surface Tsori and Andelman (2001)
To account for two confining walls at two ends in the current study, we rewrite the above expression in terms of a function as,
(6) 
where, and are the interaction strengths of the wall at and respectively. In the present study the term is set to zero, i.e. we neglect any deviation of interaction from the bulk. This specific choice of surface potential results in short range interaction. The electrostatic contribution to the free energy functional can be obtained as Landau et al. (1984); Amundson et al. (1993),
(7) 
where, is the permittivity of free space, is the volume occupied by one polymer chain and is the dielectric permittivity which is taken to be phase dependent. is the space dependent potential due to the applied voltage. We apply a linear interpolation of the permittivity between the two phases assuming the polymer to behave as a linear dielectric material by,
(8) 
The assumption of linear dielectric behaviour has been previously employed in SCFT calculations Lin et al. (2005); Matsen (2006a). Thus the total free energy functional in units of can be finally written as,
(9) 
Substituting the derived expressions we get,
(10)  
In absence of the electrostatic free energy contribution, the model reduces to the one from Brown et al. Brown and Chakrabarti (1995), that was used to study surface induced ordering of block copolymers. As evident, the free energy functional is dependent upon , and their spatial derivatives. The minimization of the functional, then, requires the evaluation of the variation with respect to each individual variable i.e. and Garcia, Bishop, and Carter (2004); Matsen (2006). The kinetic evolution of the conserved order parameter follows the dynamics of Model B framework Hohenberg and Halperin (1977),
(11) 
where the chemical potential is defined as,
(12) 
Hence the kinetic equation can be expressed as,
(13) 
Additionally, the variation of the functional with respect to , assuming the electric field to be in local equilibrium, leads to Garcia, Bishop, and Carter (2004),
(14) 
The above expression is nothing but the Maxwell equation which provides the spatial distribution of .
We remark that thermal fluctuations can be accounted by adding a noise term in Eq. 11, where, is the strength of the noise, the reciprocal of which is roughly equal to the quench depth and is noise distribution following the fluctuationdissipation theorem . It is worth mentioning that thermal fluctuations are important as far as the orderdisorder transitions in weak segregation regime are concerned. However, the equilibrium morphology which is the focus of the present work is not influenced upon incorporation of stochastic noise. The present claim is corroborated by Ref. Brown and Chakrabarti (1995), where the effect of confining surfaces was studied and the final morphology was found to be independent of noise, both in WSL as well as SSL. Therefore, we assert that the neglect of stochastic noise in the present work is reasonably well justified.
ii.2 Numerical methods and parameters
We nondimensionalize all the quantities selecting characteristic energy scale , length scale and time , The free energies are rescaled using , by the lattice cell size and the time by . The terms and is nondimensionalized as, and . The dielectric constants are rendered nondimensional using and the electric field by . Using representative values of and , we obtain .
The simulation geometry chosen for the present study is shown in Fig. 2. The setup consists of two rigid surfaces at and confining the copolymer film across which a constant voltage is applied. Dirichlet boundary conditions are applied for voltage at and with and while Neumann boundary condition is applied at and . Therefore, electric field is aligned along ydirection. The confining substrates, also attract one of the copolymers which is controlled by the numerical parameter as mentioned before. The appropriate boundary condition to account for attracting substrates translate into, Cahn (1977); Puri and Binder (1992)
(15) 
Additionally, no mass transport is allowed through the rigid surface by applying a noflux boundary condition at the surfaces, Puri and Binder (1992); Brown and Chakrabarti (1994)
(16) 
Periodic boundary condition is applied for in xdirection. We solve Eq. (13) using an explicit finite difference method where the spatial derivatives are discretized using central difference which is second order accurate in space and temporal discretization using first order Euler technique. The laplace equation in Eq. (14) is solved iteratively using SuccessiveOverRelaxation (SOR) method. The initial guess for is tailored by providing a linear initial profile in ydirection (corresponding to constant electric field, since, ) to facilitate faster convergence. The various nondimensional model parameters are selected as , , , , , , . We remark, that the present results are not influenced by the choice of grid resolution ( and ). To this end, we are able to replicate our numerical results (at ) with finer grid spacing (). In order to scale up the timestep width which scales as for CahnHilliard equation, we conveniently resort to a larger grid spacing. The values of the permittivity closely resembles a PSPMMA copolymer system Amundson et al. (1994); Ashok et al. (2001), though other values have also been used in the literature Lyakhova et al. (2006). The values of and correspond to a segregation of . The surface interaction strength is varied as 0.1, 0.5, 1.0 and 1.5. Since we intend to study symmetric walls i.e. both the surfaces attract the same monomer, we select . Moreover, we are interested to study the effect of confinement. Simulations are carried out for different box sizes in ydirection, . The natural lamellar spacing is around grid points, so that the selected film thicknesses allows us to study systems with , and . Another implication of the above mentioned values of substrate interaction and film thickness is that the surface induced ordering length is greater than the film thickness. In other words, this implies that in the absence of electric field, lamellae, parallel to the substrate, span across the entire film. The surface induced ordering length for the smallest substrate affinity of is around . The magnitude of electric field is tuned by changing the value of applied voltage, and by normalising it with the box size , i.e., to maintain the same electric field for different box sizes. The box size in xdirection is kept fixed as in all the simulations. The initial microstructure is generated by assigning a computational noise between about the average composition () corresponding to a disordered state. The system is then allowed to evolve in presence of electrostatic field and attracting substrates.
To gain insights during the microstructure evolution process we define two parameters, average density profile along y direction Brown and Chakrabarti (1994); Yan and Xie (2007) and degree of alignment Hori et al. (2007); Zhang and Qiao (2013) as,
(17)  
(18) 
where and are the fourier space wave vectors in x and y direction and = . is the magnitude of intensity of the fourier power spectrum defined as . denotes the system size and the terms in the angular bracket imply the product of and its complex conjugate in fourier space. In cases, where the alignment is parallel to the substrates i.e parallel lamellae along the y direction form, , as there is no relevant periodicity along this direction. As a result, the value of degree of alignment parameter is . In the opposite case, the alignment is perpendicular to the surface, and the value of degree of alignment parameter is . Thus, a value of implies parallel lamellae, while, a value of implies a perpendicular lamellae. A parallel lamellae arrangement in y direction is characterized by an oscillatory average density profile, whereas a flat profile about corresponds to a perpendicular lamellar arrangement.
Iii Results
iii.1 Effect of electric field and surfaces
We first study some typical morphologies arising due to the interplay of substrate interaction, confinement and electric field followed by an evaluation of the resulting phase diagram. The result is categorized into three different regimes depending on the film thickness .
iii.1.1 Thicker films with
The combined effect of the substrate and electric field for model parameters , and is presented in Fig. 3. The phase separation initiates from the surface leading to the formation of parallel lamellae. However, at t = 80, the effect of electric field sets in, leading to undulations which ultimately break up the inner layers into smaller domains. Subsequently, the smaller domains coalesce and get aligned in the direction of the electric field. This phenomenon proceeds outwards layer by layer resulting in a perpendicular lamellar arrangement due to energetic consideration.
The average density profile in Fig. 3(e) during early stages (corresponding to t = 10) is oscillatory near the surface, due to the formation of enriched and depleted layers. With time the oscillatory profile develops throughout the film thickness. At t = 180, the innermost oscillation dies out and is replaced by a flatter profile which highlights the destruction of parallel structure at the center of the film. Much later, (corresponding to t = 5000) the density profile becomes flat in the bulk of the film. However the average value shows small enrichment layers at the immediate vicinity of the surfaces even though the microstructure at the final timestep appears to be completely perpendicular. The kinetics of alignment is presented in Fig. 3(f). The value of is 1 during early stages corresponding to parallel ordering along the surface. There is a smooth temporal transition from 1 to a value closer to +0.9 which depicts the formation of perpendicular lamellae.
The influence of increasing the magnitude of substrate interaction strength and electric field for the same film thickness is presented in Fig. 4. The mechanism that leads to the formation of perpendicular lamellae is essentially the same as earlier, i.e local lamellae disruption and coalescence. However, as a result of greater substrate interaction, the parallel lamellae near the surfaces i.e. the wetting layer never break resulting in alternate enriched and depleted layers at the boundaries. In literature such a morphology is termed as mixed Pereira and Williams (1999); Tsori and Andelman (2002); Xu et al. (2003). Though the surface induced ordering length is greater than the film thickness of our study, the effect of surface is predominant closer to the walls and fades as we move away from the walls. In other words this implies that the effect of surface is non uniform over the whole ordering length. As a result, when the electric field drives the domain alignment perpendicular to the surface, above a threshold interaction strength , and below a threshold electric field , the substrate interaction dominates near the walls resulting in a few parallel layers. Meanwhile the effect of electric field is predominant at the center (away from the wall) and is able to induce a change in configuration in this region.
The degree of alignment achieved in the direction of electric field is around 75 in Fig. 4(e). Due to higher magnitude of electric field, faster kinetics is observed as can be seen by either comparing the microstructures in Fig. 3 and Fig. 4 or by comparing the slope of in Fig. 3(e) and Fig. 4(e) during the transition period. The latter has a steep transition region as compared to the smoother transition region in the former case.
Finally we discuss the consequences of further increasing the electric field to whilst keeping the other two parameters and unaltered. The results are shown in Fig. 5. An interesting phenomenon to observe is the mechanism of alignment by the electric field. In contrast to the previous cases, the parallel ordering never goes beyond two layers. Instead, electric field is sufficiently high to orient the composition fluctuations in the nonphase separated region, leading to the appearance of perpendicular lamellae at the middle of the film. Subsequently, the parallel layers near the walls also collapse and the rearrangement of perpendicular lamellae proceeds by defect annihilation mechanism Amundson et al. (1994). The average density profile in Fig. 5(e) shows enrichment layers at the walls at all times but the oscillatory profile, characteristic of the parallel lamellae configuration never develops at the film center. The transition regime of alignment kinetics is abrupt as compared to the earlier cases. The value of saturates to a value of +0.8 which constitutes to 90 alignment in the direction of the applied field.
iii.1.2 Films with
The next study focuses on configurations when the film thickness is comparable to the bulk lamellar spacing. In particular we discuss the case of film thickness of which is approximately equal to . For two cases with (Fig. 6(a)) and (Fig. 6(b)) are presented. The electric field is kept constant at . For low substrate interaction strengths, a parallel arrangement is found to be stable as in Fig. 6(a). However, with increasing substrate interaction strength we see, circular domains emerging at the center of the film as in Fig. 6(b). To differentiate the effect of electric field and substrate interaction responsible for this phenomenon, we consider the microstructural pattern in the absence of electric field. The resulting evolution (not shown) is similar, comprising of an inner layer of circular domains. Therefore, it can be inferred that the effect is solely driven by the substrate confinement independently of the applied electric field. The results can also be interpreted in terms of interference of composition waves. For higher , the film thickness, , is close to half integral of equilibrium lamellar spacing. A destructive interference takes place at the center due to the composition waves emanating from the walls and an inner lamellae cannot be maintained. As a result circular domains start appearing in the middle. The breakup of the inner lamellae into circular domains is similar to the formation of holes as described by the meanfield theory of Shull Shull (1992). The same phenomenon is absent at low , presumably because the destructive interference at the film center is not sufficiently strong. An average density plot including the transition at for the case of lower is presented in Fig.6(d). It can be verified that the innermost layer is thinner () than the next two adjacent layer on both sides() and a parallel arrangement can still be maintained with the innermost layer being in a compressed state.
The average density profile presented in Fig. 6e) shows the asymmetry in the roots of . The average value at the center is slightly B rich, clearly signifying the absence of any parallel lamellae structure. The slight asymmetry in the average density points indicates, either the evolution of a perpendicular phase or an in plane asymmetric phase. The microstructure in Fig. 6b) indicates the second possibility where circular domains coexist at the center simultaneouly with the wetting layers. The presence of such hybrid structures (combination of two different phases) have previously been reported in cylinder forming systems at similar film thickness Huinink et al. (2000). However their transition in electric field has not been reported previously. Although a 3D simulation is desirable to adequately address the issue, we can certainly make some predictions from the current 2D study. The average density profile for at different field strengths is presented in Fig. 7. All the plots correspond to . At low electric field strengths, , the plot is similar to Fig. 6d), comprising of parallel lamellae. With a slight increase in the electric field strength, the value of at the center of the film shifts towards zero (slightly B rich). The density profile is similar to Fig. 6e), and denotes the appearance of a structure other than lamellae. In this case it is not the effect of substrate and confinement alone that causes this transition, but the presence of electric field does play its part. With further increase of electric field and , the surface enrichment decreases. At the same time, the average value at the center of the film shifts to positive values and the profile tends to get flatter. Though it is clear that perpendicular phases now span, atleast in the middle of the film, the exact nature is very difficult to determine precisely in 2D simulation. We speculate that the perpendicular phases are either cylindrical structures (for ) or perpendicular lamellae (for ).
iii.1.3 Films with
Next, we decrease the film thickness below the equilibrium lamellar spacing, to . The results corresponding to model parameters and are shown in Fig. 8. The early stage microstructure corresponding to t = 10 is a superposition of parallel and perpendicular lamellae. With time (t= 30 and 70), the system evolves through a metastable antisymmetric configuration and transforms into a perpendicular state. This observation is also corroborated by the average density profile. The result is consistent with the findings of Walton et al. Walton et al. (1994) who argue in favor of a transient antisymmetric arrangement in symmetric thin films during the formation of vertical configuration.
Even with the smallest electric field strength (used in the course of this study), a perpendicular arrangement is seen to be stable. Even in the absence of the electric field, a stable perpendicular arrangement establishes implying that the geometrical confinement predominates over the electric field. Interestingly, reaches the value of +1 in this case. This observation points towards the following important fact : When the arrangement is guided by substrate confinement, deviation from perfect perpendicular morphology is negligible. However, significant deviation in perpendicularity is observed when the ordering is achieved due to the application of electric field.
For , parallel arrangement is found to be stable in the absence and at low strength of electric fields. At higher electric field strengths, the microstructure is a superposition of parallel and perpendicular lamellae as can be seen in Fig. 9. Because of the small film thickness, the substrate interaction is predominant and electric field is not able to completely eradicate the previous surface ordering phenomenon. Correspondingly saturates to a value of 0 which is midway between parallel and perpendicular configuration.
iii.2 Phase diagram
We construct a phase diagram to summarize the influence of electric field , surface interaction strength and film thickness on the evolving morphologies by classifying them into parallel, perpendicular and mixed category. We designate the morphology mixed only when there exists at least one completely parallel layer ( since the layers closer to the substrate are onehalf of the inner layers) . Classification based on such criteria will allow us to compare the resulting phase diagram with the analytical calculations Tsori and Andelman (2002); Pereira and Williams (1999). Any other combination of phases is denoted as hybrid structure. The resulting configuration stability diagram is shown in Fig. 10. The following points can be appreciated,

For a given film thickness, the magnitude of applied electric field to induce a perpendicular arrangement increases with increasing magnitude of substrate interaction strength. Any deviation from this generality occurs only for very thin films e.g. (Fig. 10(e)) and closer to halfintegral lamellar thickness i.e (in our case , Fig. 10(c)).

At low strength of substrate interaction () and electric fields, parallel arrangement is found to be stable. However beyond a certain critical value of the electric field e.g. corresponding to , the configuration oscillates between perpendicular and parallel configuration (Fig. 10 (a)(e)). The present findings accentuate the previous analytical results Tsori and Andelman (2002). In general : As decreases, the effect of substrate interaction becomes more prominent. Therefore a higher magnitude of electric field is required to induce a transition from parallel to perpendicular configuration. However, for incommensurable with the bulk lamellar spacing (say ) i.e. halfway between integral lamellar spacings, the free energy of the parallel configuration is maximum and hence a lower electric field can induce a perpendicular transition.

For unstrained films e.g. , the critical electric field required for a parallel to perpendicular transition scales as Tsori and Andelman (2002). Though the numerical calculations are carried out at discrete values of electric field, we can still verify this behavior: for h = 0.1 and lies between and (Fig. 10 (a)). Taking for , for can be calculated to be which lies between the values of and in Fig. 10(b). Similar trend can be verified for in Fig. 10(a) and (b).

For higher substrate interaction strengths, i.e , the range of stability of parallel configuration increases with decreasing film thickness (Fig. 10(a), (b), (d), (e)). For , the parallel configuration is more stable.

Mixed morphologies are stable only for thicker films () and for higher substrate interaction strength, (Fig. 10(a), (b)). Interestingly, the electric field required for the transition of parallel to mixed morphology, i.e. the first critical field is independent (or at most weakly dependent due to the discretized nature of the phase diagram) of substrate interaction strength (Fig. 10(a) and (b)). However, the second critical field i.e. the field required to convert mixed to perpendicular morphology is dependent strongly upon the substrate interaction strength (approximately linearly).

Next, we consider the variation of the critical fields with film thickness for unstrained films () for fixed substrate interaction strengths of (Fig. 10 (a) and (b)). Clearly, both the critical fields depend upon the film thickness. The dependence of the first critical field, though, is stronger than the second critical field.

Interesting morphologies arise at where the film thickness is incommensurable with the lamella period. At low substrate affinities a usual parallel to perpendicular lamellae transition is observed but at higher affinities a wide range of hybrid structures results. The exact nature, though, is not clear in the present study.
iii.3 Role of dielectric contrast
The effect of increasing the dielectric contrast i.e. on the final configuration, will now be discussed. Increasing the contrast between the two blocks implies that the material is more responsive to an applied electric field and, as a consequence, if is large, the resulting phase diagram is governed by the relative mismatch. In the studies so far, we consider ( and ), which are close to the values reported by Amundson et al. Amundson et al. (1994) ( and ) for PSPMMA copolymer. If we now increase the permittivity difference by selecting and as reported in Ref. Lyakhova et al. (2006); Tsori and Andelman (2002); Matsen (2006b), for the same copolymer system, an increase in electrostatic free energy contribution is expected for the same magnitude of applied field.
We restrict the discussion to unstrained films (). The resulting phase diagram is presented in Fig. 11 and the corresponding values of the critical electric fields are drastically lowered. According to the analytical calculations Tsori and Andelman (2002); Pereira and Williams (1999), this decrease is proportional to . Considering the critical electric field to be for and and (Fig. 10(a)), the critical field on increasing the dielectric contrast to according to above equality (on holding and constant) yields a critical value of around which complies well with Fig. 11(a). A similar behaviour is retrieved for as well.
The nature of the phase diagram changes dramatically and the region of mixed morphology in the phase diagram is diminished. With enhanced dielectric contrast, the dependence of the critical fields on the substrate interaction strength for a given film thickness becomes rather weak. This is contrary to the behavior at low dielectric contrast where the second critical field (mixed to perpendicular) displayed a strong dependency on substrate interaction strength. However, the dependency of the critical fields on film thickness for a given substrate interaction strengths is similar to that at low dielectric contrast i.e. both fields depend on the film thickness with the dependency of the second critical field being higher than the first.
Iv Discussions and Conclusions
We conclude the paper with a critical assessment of the results from the present study to the experiments, analytical and SCFT calculations. The most relevant experimental study is due to Xu et al. Xu et al. (2003, 2004). They studied thin films of varying thickness () Xu et al. (2004). For film thicknesses less than , substrate interaction was found to be dominant resulting in parallel arrangement, even on application of electric field of . The substrate induced ordering length in their study was about and the segmental interaction or the segregation was . The segregation in the present work is relatively weaker than their experiments () and ordering length is around (with lowest substrate interaction strength). The critical field for lowest substrate interaction from our study is of the order of and for and respectively for a dielectric contrast of . The second value is actually quiet high and may well exceed the dielectric breakdown of the material and in such cases only parallel ordering would be exhibited. We remark that in recent experiments, electric field of the order of have been applied Crossland et al. (2010). However, if we examine the critical field for enhanced dielectric contrast of (the same as in experiment), the values are around and for and respectively, which are well within experimental range. Given a higher segregation in experiments, we can expect the critical fields to be higher than that in the present study. Infact using SCFT Matsen Matsen (2006b) calculated the critical field to be around for the same experimental conditions of Xu et al. ( and ) for , though substrate interaction was not explicitly considered. Moreover if we compare the critical field for from their work, the critical field would turn about to be roughly for dielectric contrast of and (Fig 5(a) in Ref. Matsen (2006b)).
We additionally compare our findings to that of Lin et al. Lin and Schick (2006) who studied sphere to cylinder transition. Using SCFT calculations, a complete phase diagram was calculated. Their segregation also corresponds to . The present work is thus complementary to their study. For weak substrate interaction, the critical field (maximum value) calculated by them is around and for and respectively (Fig. 7 in ref. Lin and Schick (2006)). Considering that sphere to cylinder transition generally takes place at field strength lower than parallel to perpendicular lamella transition, our values of and are quite in agreement with SCFT calculations. A similar comparison can also be made at higher substrate strengths.
The calculated phase diagrams are similar in spirit to that by Lyakhova et al. Lyakhova et al. (2006). In both studies, the phase diagram is obtained from dynamic microstructure evolution rather than static calculations Tsori and Andelman (2002); Lin and Schick (2006). The authors investigated parallel to perpendicular transition of lamellar morphology using dynamic SCFT coupled to perturbed solution of Maxwell equation for thin films of and segregation of . Mixed morphologies were however not observed in that study. In the present study, the system was allowed to evolve from a disordered state under combined electric and substrate field, whilst in the study of Lyakhova et al., electric field was applied to well developed microstructures. Possible difference can arise because of the initial level of ordering. Our results are qualitatively similar to the results of Lin et al. Lin and Schick (2006) who observed the presence of mixed phases in film thickness as low as in cylinder forming systems at a similar segregation.
We next compare our phase diagram to the analytical calculation of Tsori et al. Tsori and Andelman (2002). The authors computed the phase diagrams both in weak and strong segregation regime. The results presented correspond to an intermediate regime. Though our results are closer to WSL, the phase diagram presented in Fig . 10 is similar to the analytical calculation of Tsori et. al in SSL (Fig. 8 and 9 in ref. Tsori and Andelman (2002)). We remark that the two phase diagrams correspond to two different segregation regimes and are based on different assumptions. Phase diagram calculated by Tsori et al. corresponds to SSL and is based on the assumption of finite surface ordering length and high dielectric contrast ( as compared to our ). In WSL, they assumed the ordering length to be greater than the film thickness and mixed morphology was not considered. In the present study, mixed morphology is observed inspite of the ordering length being lower than the film thickness. We believe that the consideration of finite surface ordering length in context of WSL (and in the occurence of mixed morphology in particular) might have been an over assumption. Moreover, only a single intermediate phase i.e. a mixed lamellae morphology was considered. The results of the present study, however, indicates the presence of other intermediate or hybrid structures, thus altering the phase diagram significantly from analytical theories.
We conclude by briefly discussing the influence of wall interaction characteristics on the equilibrium morphologies. In the present work, we have restricted the study to symmetric substrate interaction. However, in principle two additional cases are possible. The substrate can be (i) antisymmetric i.e. both walls attract different monomers with same strength or (ii) asymmetric i.e. both walls have preference towards same or different monomers, but possess different interaction strength. We believe that the consideration of antisymmetric case may not lead to any new geometries other than the ones reported here. Only the region of their respective occurrence in the phase diagram might change, given that film thickness corresponding to integral number of lamellae spacing would then be the frustrated state and halfintegral, the natural state. On the contrary, asymmetric interaction can potentially engender an additional type of mixed morphology, where the system adopts a parallel configuration on one side (where interaction strength is comparatively higher) and perpendicular configuration on the other (where interaction strength is weaker). However, it remains to be seen if the competing electric field can stabilize other morphologies, for instance, cylindrical that has been reported in the present study.
To summarize, we have studied the morphology evolution of a symmetric diblock copolymer under competing substrate interaction and electric field using a coupled OhtaKawasaki functional and Maxwell equation. By solving the full Maxwell equation, we do not assume weak dielectric inhomogeneity, making the model equally applicable irrespective of segregation. A good agreement with the analytical and SCFT calculations, amply demonstrates the predictive capability of the proposed model. A distinct advantage of coarse graining is the accessability to large scale simulation, especially in three dimension. Moreover, in 3D the nature of mixed/hybrid morphologies are well defined and such a simple classification (parallel, perpendicular and mixed lamellae) may not be suffice. Our 2D results does point out that in the incommensurate films in the regime , this interplay of substrate, confinement and electric field leads to rich hybrid structures and even the occurence of a parallel lamellae to perpendicular cylinder transition. Infact in recent experiment on gyroid forming copolymers by Crossland et al. Crossland et al. (2010) a large number of coexisting morphologies were observed at low dielectric contrast. We present some preliminary 3D results in Fig. 12 and 13 which corroborates the findings of our present 2D study. The exact nature of these hybrid structures and their transition in electric field will be communicated shortly.
Acknowledgements.
The work has been supported financially by the ministry of the state BadenWuerttemberg (Mittelbau programme).References
 Bates (1991) F. S. Bates, Science 251, 898 (1991).
 Bates and Fredrickson (1990) F. S. Bates and G. H. Fredrickson, Annu. Rev. Phys. Chem. 41, 525 (1990).
 ThurnAlbrecht et al. (2000) T. ThurnAlbrecht, J. Schotter, G. Kästle, N. Emley, T. Shibauchi, L. KrusinElbaum, K. Guarini, C. Black, M. Tuominen, and T. P. Russell, Science 290, 2126 (2000).
 Heier et al. (1999) J. Heier, J. Genzer, E. J. Kramer, F. S. Bates, S. Walheim, and G. Krausch, J. Chem. Phys. 111, 11101 (1999).
 Albalak and Thomas (1993) R. J. Albalak and E. L. Thomas, J. Polym. Sci., Part B: Polym. Phys. 31, 37 (1993).
 Balsara et al. (1994) N. Balsara, B. Hammouda, P. Kesani, S. Jonnalagadda, and G. Straty, Macromolecules 27, 2566 (1994).
 Amundson et al. (1991) K. Amundson, E. Helfand, D. D. Davis, X. Quan, S. S. Patel, and S. D. Smith, Macromolecules 24, 6546 (1991).
 ThurnAlbrecht et al. (2002) T. ThurnAlbrecht, J. DeRouchey, T. P. Russell, and R. Kolb, Macromolecules 35, 8106 (2002).
 Böker et al. (2002) A. Böker, A. Knoll, H. Elbs, V. Abetz, A. H. Müller, and G. Krausch, Macromolecules 35, 1319 (2002).
 Fredrickson (1987) G. H. Fredrickson, Macromolecules 20, 2535 (1987).
 Shull (1992) K. R. Shull, Macromolecules 25, 2122 (1992).
 Brown and Chakrabarti (1994) G. Brown and A. Chakrabarti, J. Chem. Phys. 101, 3310 (1994).
 Turner (1992) M. S. Turner, Phys. Rev. Lett. 69, 1788 (1992).
 Walton et al. (1994) D. Walton, G. J. Kellogg, A. M. Mayes, P. Lambooy, and T. P. Russell, Macromolecules 27, 6225 (1994).
 Kikuchi and Binder (1994) M. Kikuchi and K. Binder, J. Chem. Phys. 101, 3367 (1994).
 Brown and Chakrabarti (1995) G. Brown and A. Chakrabarti, J. Chem. Phys. 102, 1440 (1995).
 Wang et al. (2000) Q. Wang, Q. Yan, P. F. Nealey, and J. J. de Pablo, J. Chem. Phys. 112, 450 (2000).
 Lambooy et al. (1994) P. Lambooy, T. P. Russell, G. J. Kellogg, A. M. Mayes, P. D. Gallagher, and S. K. Satija, Phys. Rev. Lett. 72, 2899 (1994).
 Xu et al. (2003) T. Xu, C. J. Hawker, and T. P. Russell, Macromolecules 36, 6178 (2003).
 Xu et al. (2004) T. Xu, Y. Zhu, S. P. Gido, and T. P. Russell, Macromolecules 37, 2625 (2004).
 Landau et al. (1984) L. D. Landau, J. Bell, M. Kearsley, L. Pitaevskii, E. Lifshitz, and J. Sykes, Electrodynamics of continuous media, Vol. 8 (elsevier, 1984).
 Pereira and Williams (1999) G. Pereira and D. Williams, Macromolecules 32, 8115 (1999).
 Tsori (2009) Y. Tsori, Rev. Mod. Phys. 81, 1471 (2009).
 Tsori and Andelman (2002) Y. Tsori and D. Andelman, Macromolecules 35, 5161 (2002).
 Matsen (2006a) M. W. Matsen, Macromolecules 39, 5512 (2006a).
 Ashok et al. (2001) B. Ashok, M. Muthukumar, and T. P. Russell, J. Chem. Phys. 115, 1559 (2001).
 Kyrylyuk et al. (2003) A. Kyrylyuk, G. J. A. Sevink, A. V. Zvelindovsky, and J. Fraaije, Macromol. Theor. Simul. 12, 508 (2003).
 Lyakhova et al. (2006) K. Lyakhova, A. V. Zvelindovsky, and G. J. A. Sevink, Macromolecules 39, 3024 (2006).
 Lin et al. (2005) C.Y. Lin, M. Schick, and D. Andelman, Macromolecules 38, 5766 (2005).
 Kyrylyuk and Fraaije (2006) A. V. Kyrylyuk and J. G. H. Fraaije, J. Chem. Phys. 125, 164716 (2006).
 Xu et al. (2005) T. Xu, A. V. Zvelindovsky, G. J. A. Sevink, K. Lyakhova, H. Jinnai, and T. P. Russell, Macromolecules 38, 10788 (2005).
 Pinna and Zvelindovsky (2008) M. Pinna and A. V. Zvelindovsky, Soft Matter 4, 316 (2008).
 Böker et al. (2003) A. Böker, H. Elbs, H. Hänsel, A. Knoll, S. Ludwigs, H. Zettl, A. V. Zvelindovsky, G. J. A. Sevink, V. Urban, V. Abetz, et al., Macromolecules 36, 8078 (2003).
 Zvelindovsky and Sevink (2003) A. V. Zvelindovsky and G. J. A. Sevink, Phys. Rev. Lett. 90, 049601 (2003).
 Schmidt et al. (2005) K. Schmidt, A. Böker, H. Zettl, F. Schubert, H. Hänsel, F. Fischer, T. M. Weiss, V. Abetz, A. V. Zvelindovsky, G. J. A. Sevink, et al., Langmuir 21, 11974 (2005).
 Pinna et al. (2009) M. Pinna, L. Schreier, and A. V. Zvelindovsky, Soft Matter 5, 970 (2009).
 Ruppel et al. (2013) M. Ruppel, C. W. Pester, K. M. Langner, G. J. A. Sevink, H. G. Schoberth, K. Schmidt, V. S. Urban, J. W. Mays, and A. BoÌker, ACS nano 7, 3854 (2013).
 Ly et al. (2007) D. Q. Ly, T. Honda, T. Kawakatsu, and A. V. Zvelindovsky, Macromolecules 40, 2928 (2007).
 Ly et al. (2013) D. Ly, M. Pinna, T. Honda, T. Kawakatsu, and A. V. Zvelindovsky, J. Chem. Phys. 138, 074904 (2013).
 Li et al. (2013) S. Li, Y. Jiang, Y. Ji, and X. Wang, Polymer 54, 6636 (2013).
 Schmidt et al. (2010) K. Schmidt, C. W. Pester, H. G. Schoberth, H. Zettl, K. A. Schindler, and A. Boker, Macromolecules 43, 4268 (2010).
 Pester et al. (2015) C. W. Pester, K. Schmidt, M. Ruppel, H. G. Schoberth, and A. BoÌker, Macromolecules (2015).
 Matsen (2005) M. W. Matsen, Phys. Rev. Lett. 95, 258302 (2005).
 Matsen (2006b) M. W. Matsen, Soft Matter 2, 1048 (2006b).
 Amundson et al. (1993) K. Amundson, E. Helfand, X. Quan, and S. D. Smith, Macromolecules 26, 2698 (1993).
 Amundson et al. (1994) K. Amundson, E. Helfand, X. Quan, S. D. Hudson, and S. D. Smith, Macromolecules 27, 6559 (1994).
 Ohta and Kawasaki (1986) T. Ohta and K. Kawasaki, Macromolecules 19, 2621 (1986).
 Cahn and Hilliard (1958) J. W. Cahn and J. E. Hilliard, J. Chem. Phys. 28, 258 (1958).
 Bahiana and Oono (1990) M. Bahiana and Y. Oono, Phys. Rev. A 41, 6763 (1990).
 Liu and Goldenfeld (1989) F. Liu and N. Goldenfeld, Phys. Rev. A 39, 4805 (1989).
 Chakrabarti et al. (1989) A. Chakrabarti, R. Toral, and J. D. Gunton, Phys. Rev. Lett. 63, 2661 (1989).
 Choksi, Peletier and Williams (2009) R. Choksi, M. A. Peletier, and J. F. Williams, SIAM J. Appl. Math 69, 1712 (2009).
 Schmidt and Binder (1985) I. Schmidt and K. Binder, J. Phys. (Paris) 46, 1631 (1985).
 Tsori and Andelman (2001) Y. Tsori and D. Andelman, Eur. Phys. J. E 5, 605 (2001).
 Garcia, Bishop, and Carter (2004) R. E. Garcia , C. M. Bishop and W. C. Carter, Acta Mater. 52, 11 (2004).
 Matsen (2006) M. W. Matsen, J. Chem. Phys. 124, 074906 (2006).
 Hohenberg and Halperin (1977) P. C. Hohenberg and B. I. Halperin, Rev. Mod. Phys. 49, 435 (1977).
 Cahn (1977) J. W. Cahn, J. Chem. Phys. 66, 3667 (1977).
 Puri and Binder (1992) S. Puri and K. Binder, Phys. Rev. A 46, R4487 (1992).
 Yan and Xie (2007) L.T. Yan and X.M. Xie, J. Chem. Phys. 126, 064908 (2007).
 Hori et al. (2007) H. Hori, O. Urakawa, O. Yano, and Q. TranCongMiyata, Macromolecules 40, 389 (2007).
 Zhang and Qiao (2013) G. Zhang and G. Qiao, J. Chem. Phys. 139, 134903 (2013).
 Huinink et al. (2000) H. Huinink, J. BrokkenZijp, M. Van Dijk, and G. J. A. Sevink, J. Chem. Phys. 112, 2452 (2000).
 Crossland et al. (2010) E. J. Crossland, S. Ludwigs, M. A. Hillmyer, and U. Steiner, Soft Matter 6, 670 (2010).
 Lin and Schick (2006) C.Y. Lin and M. Schick, J. Chem. Phys. 125, 034902 (2006).