Surface-Induced Nanostructures and Phase Diagrams of ABC Linear Triblock Copolymers under Spherical Confinement: A Self-Consistent Field Theory Simulation

We investigate the nanostructures and phase diagrams of ABC linear triblock copolymers confined in spherical cavities by using real-space self-consistent field theory. Various 3D morphologies, such as spherical concentric lamellae, dumbbell-like cylinder, and rotational structures, are identified in the phase diagrams, which are constructed on the basis of the diameters of spherical cavities and the interaction between the polymers and preferential surfaces. We designate specific monomer-monomer interactions and block compositions, with which the polymers spontaneously form a cylindrical morphology in bulk, and firstly study morphology transformation with a neutral surface when a confining radius progressively increases. We then focus on phase morphologies under the preferential surfaces and consolidate them into phase diagrams. The spherical radius and the degree of preferential interactions can obviously induce the formation of a cylindrical morphology. Theoretical results correspond to an amount of recent experimental observations to a high degree and contribute to synthesising functional materials.


Introduction
Block copolymers are extensively circulated for their ability to self-assemble and form ordered morphologies at a nanoscale, thereby providing a development platform for fabrication techniques of multifunctional materials. Spontaneous periodic structures possess various potential technological applications in several areas, such as drug delivery, nanolithography, and photonic crystals [1][2][3][4]. Elementary operations of diversifying the formation of block copolymers contain the modifications of molecular architectures, the degree of polymerisation, the volume fractions of components, and the interactions among monomers [5,6], with which the progress of functional materials is insufficient to satisfy the development demands. Various external techniques have been developed in theoretical and experimental studies because the physical morphologies and properties of polymer nanoparticles are not only subject to the internal controlling parameters of polymers but also influenced by external fields and space conditions. For example, hybrid nanomaterials that incorporate inorganic nanoparticles into self-assembled block copolymers under spherical nanoparticle confinement may be used in electronic and optics devices, which can combine the inherent properties of the two units and generate desirable properties [7]. Our work aims to reveal the mechanisms of phase separation and recombination into various nanostructures at the level of molecules. copolymers confined in incompressible spherical geometries. Each linear triblock copolymer consists of three distinctive chemical compositions. The polymerization index of each chain is specified as N, which is set to 400 in our system. The monomeric volumes of the three chemical compositions are assumed to be identical, v A = v B = v C , and the statistical segment length of a polymer is denoted as a. For convenience, in subsequent equations and diagrams, the lengths and volumes are all expressed in units of mean-square gyration radius, R g = a √ N/6. A, B and C monomers per chain are expressed as N A = f A N, N B = f B N, and N C = (1 − f A − f B )N, correspondingly. Fredrickson et al. proposed a masking technique to realise thin-film confinement effects [39]; here, we utilise the technique to realise the spherical confinement effects, which are particularly appropriate for pseudo-spectral numerical methods with periodic boundary conditions. The key to realising geometric constraint is introducing an interfacial density function φ W ( → r ), and a modified incompressibility constraint.
where φ i ( → r ) represents the volume fraction of species i. To enforce the polymer segments into the spherical space, we model a spherical cavity using a cosine function for φ W ( → r ), where r is the distance from the centre of the spherical cavity, R is the radius, and T is the thickness of the spherical wall. In the piecewise function, φ W (r) = 1 defines the impenetrable hard wall; moreover, 0 < φ W (r) < 1 represents the blending regions of the surface and copolymer melts, and φ W (r) = 0 represents the regions of full copolymer melts. The confining geometry with the distribution of Equation (2) controlled by the spherical radius R can be implemented given a relatively small interface region represented by a small T. However, an extremely sharp interface indicates a high spatial resolution and reduces the convergence speed of iteration computational equations. Thus, we select T = 0.45R g , which is appropriate for our confinement model. The essence of free energies is the competition between entropy, −TS, and internal energy, U. In the presence of the interactions between polymers and walls, a sum of the quadratic term must be incorporated into the Helmholtz free energy per polymer chain in a temperature T.
where k B is the Boltzmann constant; χ ij is the Flory-Huggins parameter between species i and j, i, j ∈ (A, B, C); χ iW denotes the dimensionless intensity of the interaction between species and preferential surfaces. The local incompressibility constraint in the function is ensured by a Lagrange multiplier, ξ( → r ). Q is the single-chain partition function expressed as in which the statistical weight function q( → r , s) denotes a segment of a chain with the free end (s = 0) of A block having its sth segment at → r . The value of contour variable ranges from s = 0 to 1, the terminal of C block. Therefore, the local segment density can be obtained as and The propagator q( → r , s) satisfies the following differential modified diffusion equations (MDE): and with the initial condition q( is the chemical potential field produced by surrounding chains. According to the mean-field approximation, minimising free energy with respect to monomer densities, ∂F/∂φ i = 0, leads to the following expressions for chemical potentials, The effective volume V occupied by copolymers is defined as under the spherical confinement, and the effective diameters listed in this study are calculated in accordance with the formula, V = 1 6 πd 3 eff . Real-space combinatorial screening method, originally proposed by Drolet and Fredrickson [40], is the general principle for solving the self-consistent Equations (4)- (9). A relax process is assumed by an iteration step, starting with an initially guessed chemical potential field ω i of species i. In accordance with the second-order operator-split method improved by Rasmussen et al. [41,42], the 'forward' q( → r , s) and 'backward' q + ( → r , s) end-integrated propagators can be exactly implemented numerically. Here, we use the linear relation of ω new i and ω new i , to realise the constant update, where ∆t = 0.1 is effective in solving SCFT calculations in our work. The relaxation steps are repeated until the relative difference among the free energies of the sequential iteration steps is reduced to a satisfactory precision, which is less than 10 −5 . The implementation of this process is derived from a minimisation procedure of free energy; this procedure is based on the commensuration between the size of a simulation box and the periodicity of equilibrium morphologies [42,43].
In the present study, the polymers are designed into a rectangular box of K x × K y × K z = 64 × 64 × 64 microlattices with a periodic boundary condition in bulk space, and we adjust the simulation box lengths in two directions to reach stable structures with the lowest free energies. The polymer chain parameters, f A : f B : f C = 0.2 : 0.6 : 0.2 and χ AB N = χ BC N = χ AC N = 35, are designated throughout the study, and we discover the microphase separation structure at the parameters in bulk by optimising the domain sizes. We adjust the grid spacing in two directions, x-and y-axes, to search for the stable structure, where ∆x and ∆y vary from 0.178 R g to 0.192 R g . Figure 1 illustrates the free energy potential with an equal increment ∆ = 0.0992R g in ∆x and ∆y. A distinct depression area and a rock bottom exist on the free energy surface, which is surrounded by metastable structures. The dimensionless lowest free energy, F/nk B T = 8.578047, locates at the simulation boundary lengths of l x = l y = 11.80R g . For a spherical confinement, the periodic boundary condition is unnecessary because the spherical diameters directly affect the self-assembled conformations, and the grid spacing is set to ∆x = ∆y = ∆z = 0.15R g . For the phase behavior with variable spherical surfaces, the diameters and strengths of the surface field are incrementally varied to investigate the significant phases of ABC linear triblock copolymers. We select the most stable phase from all structures produced by proposing different initial guesses of the mean fields by comparing the free energies.

Results
The phase space of the ABC linear triblock copolymer system within the nanosphere is spanned by four physical parameters, namely, , , , and , where i = A, B and C. In this study, we concentrate on the influence of confinement, including a subset of the parameter space, and , and designate specific parameters of the copolymer model. We consider symmetric ABC linear polymers, = = = 35, with which the system is already in an intermediateseparated regime, to facilitate the theoretical results to match the experiment [44,45]. Afterwards, we firstly analyse the phase structure of a symmetric ABC triblock copolymer with the volume fractions of : : = 0.2: 0.6: 0.2 in bulk (i.e., without the spherical confinement). On the basis of this foundation structure, we explore the effects of a neutral spherical confinement on a 1D phase diagram by gradually increasing the effective diameter from 2.60 to 8.60 . Subsequently, the effects of attraction and repulsion between the confinement walls and the block copolymers on the morphology are described on the − plane. For an intelligible measure of preferential adsorption strength, we define a strength proportion in our phase diagrams, = ⁄ . To demonstrate the inner structures clearly, the 3D structures are described in red, green and blue, which illustrate the A-, B-and C-rich regions in a space, respectively, even the sectional plots.

Morphologies in Bulk and Phase Diagram under Neutral Spherical Confinement
In this section, we investigate the effects of the confinement spatial scale on the cylinder-formed ABC linear triblock copolymers with a neutral spherical surface. According to the optimising result in Figure 1, we obtain the two interpenetrating tetragonally arranged cylinder phase (TET2) exhibited in Figure 2a, which has been reported in theoretical and experimental studies [45][46][47][48][49]. Plot (b) displays that A-and C-blocks are gathered in the alternate circular regions, whilst outside space is covered with B-block. The results indicate that the periods in two mutually perpendicular orientations, L1 and L2, are nearly equal, and the optimal structure is generated when the simulation box is approximately twice the size of the phase period. Perpendicular to the columnar axis, we plot the distributions of the A-terminals, = 0, C-terminals, = 1, and the AB junction points, s = , in Figure 2c-e, respectively, using red, blue, and black. The plots in (c) and (d) demonstrate that the ends of the A and C blocks all gather near the central regions of the respective rich domain, and the distribution of conjunctions forms a diffused, perfect circle between the A-and B-rich domains, as presented in Figure 2e. In accordance with the segment distributions and orientational ordering, a

Results
The phase space of the ABC linear triblock copolymer system within the nanosphere is spanned by four physical parameters, namely, f i , χ ij N, χ iW N, and d eff , where i = A, B and C. In this study, we concentrate on the influence of confinement, including a subset of the parameter space, χ iW N and d eff , and designate specific parameters of the copolymer model. We consider symmetric ABC linear polymers, χ AB N = χ BC N = χ AC N = 35, with which the system is already in an intermediate-separated regime, to facilitate the theoretical results to match the experiment [44,45]. Afterwards, we firstly analyse the phase structure of a symmetric ABC triblock copolymer with the volume fractions of f A : f B : f C = 0.2 : 0.6 : 0.2 in bulk (i.e., without the spherical confinement). On the basis of this foundation structure, we explore the effects of a neutral spherical confinement on a 1D phase diagram by gradually increasing the effective diameter d eff from 2.60 R g to 8.60 R g . Subsequently, the effects of attraction and repulsion between the confinement walls and the block copolymers on the morphology are described on the χ iW − d eff plane. For an intelligible measure of preferential adsorption strength, we define a strength proportion λ iW in our phase diagrams, λ iW = χ iW /χ ij . To demonstrate the inner structures clearly, the 3D structures are described in red, green and blue, which illustrate the A-, B-and C-rich regions in a space, respectively, even the sectional plots.

Morphologies in Bulk and Phase Diagram under Neutral Spherical Confinement
In this section, we investigate the effects of the confinement spatial scale on the cylinder-formed ABC linear triblock copolymers with a neutral spherical surface. According to the optimising result in Figure 1, we obtain the two interpenetrating tetragonally arranged cylinder phase (TET 2 ) exhibited in Figure 2a, which has been reported in theoretical and experimental studies [45][46][47][48][49]. Plot (b) displays that A-and C-blocks are gathered in the alternate circular regions, whilst outside space is covered with B-block. The results indicate that the periods in two mutually perpendicular orientations, L 1 and L 2 , are nearly equal, and the optimal structure is generated when the simulation box is approximately twice the size of the phase period. Perpendicular to the columnar axis, we plot the distributions of the A-terminals, s = 0, C-terminals, s = 1, and the AB junction points, s = f A , in Figure 2c-e, respectively, Polymers 2018, 10, 1276 6 of 19 using red, blue, and black. The plots in (c) and (d) demonstrate that the ends of the A and C blocks all gather near the central regions of the respective rich domain, and the distribution of conjunctions forms a diffused, perfect circle between the A-and B-rich domains, as presented in Figure 2e. In accordance with the segment distributions and orientational ordering, a schematic is presented in Figure 2f for the column-captured domains in the x-y plane, in which B middle block has to bridge them on the dotted lines. schematic is presented in Figure 2f for the column-captured domains in the x-y plane, in which B middle block has to bridge them on the dotted lines. We model a spherical cavity in the centre of the simulation space to study the influence of confined environment on the self-assembled nanostructure and ensure that the polymer molecules cannot escape from the hard wall. The main features of the morphologies of ABC linear triblock copolymers confined in spherical cavities of uniform variable diameters, in which the cavity surfaces are neutral to the confined polymers, are investigated. Six structures are observed under the neutral spherical confinement illustrated in Figure 3a with a proportionable size. These structures are the partially wrapped oblate spheroid (WOS), spherical concentric lamellae (SCL), dumbbell-like cylinder (DLC), semiring-triangle (ST) structure, quadruple texture (QT) structure and irregular ringlike (IR) structure. The structure images listed in each table grid are the overall structure, profile and configuration snapshots of A-, B-and C-blocks, correspondingly.
In the first field depicted in Figure 3a, the WOS structure observed in a small restricted space contains two combined parts, namely, a hollow hemisphere formed by B block and a single solid ellipsoid formed by C block. The hollow hemisphere tightly wraps the solid ellipsoid. On the interface between the explicit regular structure and the confinement wall, even inside the wall, monomers of A block and a fraction of B block are mixed well because of the compression of confinement and the We model a spherical cavity in the centre of the simulation space to study the influence of confined environment on the self-assembled nanostructure and ensure that the polymer molecules cannot escape from the hard wall. The main features of the morphologies of ABC linear triblock copolymers confined in spherical cavities of uniform variable diameters, in which the cavity surfaces are neutral to the confined polymers, are investigated. Six structures are observed under the neutral spherical confinement illustrated in Figure 3a with a proportionable size. These structures are the partially wrapped oblate spheroid (WOS), spherical concentric lamellae (SCL), dumbbell-like cylinder (DLC), semiring-triangle (ST) structure, quadruple texture (QT) structure and irregular ring-like (IR) structure. The structure images listed in each table grid are the overall structure, profile and configuration snapshots of A-, B-and C-blocks, correspondingly.
In the first field depicted in Figure 3a, the WOS structure observed in a small restricted space contains two combined parts, namely, a hollow hemisphere formed by B block and a single solid ellipsoid formed by C block. The hollow hemisphere tightly wraps the solid ellipsoid. On the interface between the explicit regular structure and the confinement wall, even inside the wall, monomers of A block and a fraction of B block are mixed well because of the compression of confinement and the cosine function model for the wall. An interchange of A-and B-rich regions is permissible in view of the complete folded symmetry of the polymer chains, which are the same as the phases shown subsequently. Higuchi et al. reported a similar structure of diblock copolymers under a spherical 3D confinement [50][51][52]. Given that the volume of the confinement space increases and monomers become unbridled, the middle component completely encases a balled-up structure, whilst the other minority component remains hidden. The SCL structure described above is also observed in several other copolymer systems, such as diblock copolymers in spherical nanopores [53]. Onion-like layering systems, each layer being rich in one of the three monomer components, are presented in the study of ABC star triblock copolymers confined in a spherical cavity [54]. The three monomer components all begin to mutually separate, and the rotational symmetric DLC structure occurs when the enclosed space enlarges. An exchangeable minority segment forms a dumbbell standing in the centre of the structure through the entire structure. The other segment discontinuously appears at the equator of the nanoparticle, where the majority middle domain resembles a large cylindrical shell. Subsequently, two small disks of A-monomers are elongated to a semiring shape. The C-monomers form a triangular column domain centred in a B-rich matrix in the ST structure, where the B middle domain resembles a boat anchor. With a continuous increase in a confined volume, the A-monomers are distorted to interconnect and form a coronal structure which resembles a triangular pyramid with arc edges. The C-monomers form a tetragonal column domain centred in a B-rich matrix, in which the latter provides a background for the minority structures to interweave. Overall, the QT structure consists of four identical regions with a threefold rotational symmetry. A normal corona-shell-core structure observed in A 2 B 8 C 8 terpolymer system under rigid spherical confinement is similar to our results [55]. When the spherical diameter d eff approximates the edge length of the square simulation box K x × ∆x, an IR structure occurs rather than analogous cylinder structures in bulk. the complete folded symmetry of the polymer chains, which are the same as the phases shown subsequently. Higuchi et al. reported a similar structure of diblock copolymers under a spherical 3D confinement [50][51][52]. Given that the volume of the confinement space increases and monomers become unbridled, the middle component completely encases a balled-up structure, whilst the other minority component remains hidden. The SCL structure described above is also observed in several other copolymer systems, such as diblock copolymers in spherical nanopores [53]. Onion-like layering systems, each layer being rich in one of the three monomer components, are presented in the study of ABC star triblock copolymers confined in a spherical cavity [54]. The three monomer components all begin to mutually separate, and the rotational symmetric DLC structure occurs when the enclosed space enlarges. An exchangeable minority segment forms a dumbbell standing in the centre of the structure through the entire structure. The other segment discontinuously appears at the equator of the nanoparticle, where the majority middle domain resembles a large cylindrical shell.
Subsequently, two small disks of A-monomers are elongated to a semiring shape. The C-monomers form a triangular column domain centred in a B-rich matrix in the ST structure, where the B middle domain resembles a boat anchor. With a continuous increase in a confined volume, the A-monomers are distorted to interconnect and form a coronal structure which resembles a triangular pyramid with arc edges. The C-monomers form a tetragonal column domain centred in a B-rich matrix, in which the latter provides a background for the minority structures to interweave. Overall, the QT structure consists of four identical regions with a threefold rotational symmetry. A normal corona-shell-core structure observed in A B C terpolymer system under rigid spherical confinement is similar to our results [55]. When the spherical diameter approximates the edge length of the square simulation box × ∆x, an IR structure occurs rather than analogous cylinder structures in bulk. The 1D phase diagram for varying from 2.60 to 8.60 is demonstrated in Figure 3b. However, regardless of whether increases or decreases outside the range, the structures at both ends survive intact, and the data are inconsequential to our work. In calculating the phase diagram, the increment in is set to 0.12 near the phase boundary and is 0.30 in other cases. In The 1D phase diagram for varying d eff from 2.60 R g to 8.60 R g is demonstrated in Figure 3b. However, regardless of whether d eff increases or decreases outside the range, the structures at both ends survive intact, and the data are inconsequential to our work. In calculating the phase diagram, the increment in d eff is set to 0.12 R g near the phase boundary and is 0.30 R g in other cases. In Figure 3b, with the increase in effective diameter of a spherical cavity, the ordered morphologies described above appear successively. The SCL structure occupies the maximum span, and the DLC structure is briefly present in the diagram. When d eff is greater than one-half of the edge length K x × ∆x/2, the ordered morphologies change actively and are complex. We depict the dependence of the free energy F/nk B T as a function of d eff /R g in Figure 4 to demonstrate that the structures illustrated in Figure 3 transform. The three free energy branches, which correspond to the DLC, ST and QT structures, can be viewed here. The crossing points determine the first-order phase transition among the three studied phases, which precisely testify the phase boundaries exhibited in Figure 3b. A current conclusion of the spherical confinement system is that cylinder-forming linear triblock copolymers assemble into intricate and symmetric structures that are widely divergent with the structures in bulk because the hard wall thwarts the periodicity and applies isotropic compression. Experimentally, there have been some methods to obtain spherical nanoparticles. Zhu and Chen et al. have reported that the polymeric nanorods obtained in anodic aluminium oxide template can transform to spheres upon thermal annealing, which is driven by Rayleigh instability [56,57]. A self-organized precipitation method used to prepare block copolymer nanoparticles can control the particle diameter from several tens of nanometers to several micrometers by changing the preparation conditions, such as concentration of the polymers and evaporation speed of solvent [51]. structure is briefly present in the diagram. When is greater than one-half of the edge length × ∆x 2 ⁄ , the ordered morphologies change actively and are complex. We depict the dependence of the free energy ⁄ as a function of ⁄ in Figure 4 to demonstrate that the structures illustrated in Figure 3 transform. The three free energy branches, which correspond to the DLC, ST and QT structures, can be viewed here. The crossing points determine the first-order phase transition among the three studied phases, which precisely testify the phase boundaries exhibited in Figure 3b. A current conclusion of the spherical confinement system is that cylinder-forming linear triblock copolymers assemble into intricate and symmetric structures that are widely divergent with the structures in bulk because the hard wall thwarts the periodicity and applies isotropic compression. Experimentally, there have been some methods to obtain spherical nanoparticles. Zhu and Chen et al. have reported that the polymeric nanorods obtained in anodic aluminium oxide template can transform to spheres upon thermal annealing, which is driven by Rayleigh instability [56,57]. A selforganized precipitation method used to prepare block copolymer nanoparticles can control the particle diameter from several tens of nanometers to several micrometers by changing the preparation conditions, such as concentration of the polymers and evaporation speed of solvent [51]. 5

Morphologies and Phase Diagram in Spherical Cavity with AC-Attractive Surface
In addition to the confinement degree, a short-range interaction between the surfaces and copolymers plays a crucial role in the frustrated packing of compared block copolymers, which has been verified in multiple systems [39,54,[58][59][60][61]. Here, we focus on the effects of the confinement degree and preferential surface. On the basis of the 1D research route with the variable , we expand the dimensions of the phase diagram by assuming an alterable interaction between the surfaces and copolymers.
which represents that A and C components attach to the wall, and the B component is repelled from the surface. The strengths of the three polymer-surface interactions, λ , are constantly kept synchronised in the following study. However, a high-precision description of the morphologies in the 2D parameter space is formidable. Considering this case, − , is increased from 0.1 to 0.8 in an increment of 0.1 and from 2.6 to 8.6 in an increment of 0.6 , thereby producing 88 phase points in the phase

Morphologies and Phase Diagram in Spherical Cavity with AC-Attractive Surface
In addition to the confinement degree, a short-range interaction between the surfaces and copolymers plays a crucial role in the frustrated packing of compared block copolymers, which has been verified in multiple systems [39,54,[58][59][60][61]. Here, we focus on the effects of the confinement degree and preferential surface. On the basis of the 1D research route with the variable d eff , we expand the dimensions of the phase diagram by assuming an alterable interaction between the surfaces and copolymers.
which represents that A and C components attach to the wall, and the B component is repelled from the surface. The strengths of the three polymer-surface interactions, λ, are constantly kept synchronised in the following study. However, a high-precision description of the morphologies in the 2D parameter space is formidable. Considering this case, −λ ACW , is increased from 0.1 to 0.8 in an increment of 0.1 and d eff from 2.6 R g to 8.6 R g in an increment of 0.6 R g , thereby producing 88 phase points in the phase diagram. Figure 5a displays six structures with their internal conformations found in the abovementioned parameter array; these structures are the monochromic sphere (MS), SCL structure, Frisbee-ring (FR) structure, discrete lattice-lamellae (DLL), and two types of symmetrical Frisbee (SF-I and SF-II). The MS, SCL and FR structures exhibited in the left column in Figure 5a present the same exterior morphology, similar to homogeneous, which can be distinguished by the profile images. The MS structure closely approximates a homogeneous sphere, except for the mixing region inside the wall, which is removed from the structure images. We call the mixing region a disordered wetting layer, which has been mentioned by Chen and Liang under a thin film system [37]. The A-, B-and C-blocks gather in the wetting layer in contrast to their results given the combined effect of the compression and the attraction from the surfaces. The SCL structure survives from the AC-attractive surface field because of its strong symmetrical stability, thereby indicating that the block chain is insensitive to this polymer-wall interaction. The internal conformation of the FR structure transforms into a layered architecture, which is linked to a relatively large space, as presented in the superstratum in Figure 5b. A minority block forms a ring on the centre plane, whereas the other block forms two Frisbees, symmetrical to each other about the midplane. The FR profile image shows that the concaves of the Frisbees face outward to the wall.
The DLL structure contains one spherical core in the centre and possesses the mirror symmetry. This phenomenon can be considered an evolution from the SCL to the DLL structure with an increase in confinement diameter d eff under a weak wall-polymer interaction, as illustrated in the column of −λ ACW = 0.1 in Figure 5b. Several A-monomers break away from the wall and are discretely embedded in the B-block matrix considering the relaxed restriction and weak preferential potential. Moreover, two types of SF structures are observed in the last two columns in Figure 5a. The SF type contains two Frisbees formed by two minority components at the mirror symmetric locations. The profile image clearly demonstrates that the concaves of the Frisbees all face outward to the interior. In contrast to SF-II, a small A-rich disk and a B-rich disk are embedded on the outer layer of SF-I. Overall, A-and C-rich regions appear alternately along the symmetric axis. In particular, these regions have the same formation mechanism and internal interactions but with a large space in the SF-I model. The spherical surface releases a few minority components. In contrast to the case under the neutral surface, these frustrated structures are generally succinct and of notable regularity. Figure 5b plots the location of these morphologies represented by the shrunken profile images. In Figure 5b, the MS structures are all located at low spherical diameters. The area must be a narrow space that can only accommodate the majority of the middle component and ensure that the middle block bridges on both sides of the wall. The SCL structure occupies the middle area of the phase diagram. A niche for the SCL with any interface interaction λ exists. With a fixed d eff , such as d eff = 4.4 R g , the SCL phase translates to the MS phase with the increase in −λ ACW ; with a fixed −λ ACW , the SCL phase also translates to the MS phase with the decrease in d eff . In this particular parameter environment, the space confinement and surface adsorption field have the same effect on the morphologies. The DLL and FR structures, which correspond to the QT and IR structures in the neutral system, are stabilised only in the upper-left corner of the phase diagram, divided by a red dashed line, of a parameter space (−λ ACW = 0.1). A reasonable explanation is that the minority components adhere to the spherical cavity walls and then form several layered structures with the fixed linear topological constraint when the adsorption interactions between the walls and blocks are applied. The combination of strong AC-attractive fields and large spherical space leads to the two homologous SF structures that show a high stability in the phase diagram. The phase transition between FR and SF phases are subtle; therefore, the transition process from FR to SF-I is obtained along the phase path = 8.0 as an example. Figure 6a displays the free energies as functions of the AC-attractive strength − . To highlight the cross, we subtract the free energies from a linear function.
The two free energy curves are satisfied at approximately − = 0.115, which denotes an order-order transition (OOT). The free energy is higher in the SF-I phase at the − < 0.115 region than in the FR phase, thereby indicating that it is a metastable phase. However, the SF-I phase stabilises at the − > 0.115 region, where its free energy is low. The intersection also suggests that the transition between the FR and SF-I phases is a first-order phase transition. We divide the free energy into two parts, namely, the internal energy ⁄ and entropic energy − ⁄ , to discuss their contribution separately in Figure 6b,c. From a global perspective, ⁄ and − ⁄ The phase transition between FR and SF phases are subtle; therefore, the transition process from FR to SF-I is obtained along the phase path d eff = 8.0R g as an example. Figure 6a displays the free energies as functions of the AC-attractive strength −λ ACW . To highlight the cross, we subtract the free energies from a linear function.
The two free energy curves are satisfied at approximately −λ ACW = 0.115, which denotes an order-order transition (OOT). The free energy is higher in the SF-I phase at the −λ ACW < 0.115 region than in the FR phase, thereby indicating that it is a metastable phase. However, the SF-I phase stabilises at the −λ ACW > 0.115 region, where its free energy is low. The intersection also suggests that the transition between the FR and SF-I phases is a first-order phase transition. We divide the free energy into two parts, namely, the internal energy U/nk B T and entropic energy −S/nk B , to discuss their contribution separately in Figure 6b,c. From a global perspective, U/nk B T and −S/nk B decrease with −λ ACW for both phases. The increase in surface adsorption strength gradually reduces the internal energies and the disordered degree of the system. A sudden drop in the entropic energy of the SF-I phase from −λ ACW = 0.11 to 0.12 occurs, thus decreasing the free energy. Therefore, the entropic energy is the principal factor behind this phase transition. decrease with − for both phases. The increase in surface adsorption strength gradually reduces the internal energies and the disordered degree of the system. A sudden drop in the entropic energy of the SF-I phase from − = 0.11 to 0.12 occurs, thus decreasing the free energy. Therefore, the entropic energy is the principal factor behind this phase transition.

Morphologies and Phase Diagram in Spherical Cavity with B-Attractive surface
In this subsection, we examine the inverse adsorption potential, the B-attractive surfaces, to continue the theme: which represents that the spherical cavity prefers B-monomers and repels A-and C-monomers. Two minority components are no longer constrained but freely interact with the middle component inside the cavity; thus, the structures observed here can be diversified and intricate. On the basis of the

Morphologies and Phase Diagram in Spherical Cavity with B-Attractive surface
In this subsection, we examine the inverse adsorption potential, the B-attractive surfaces, to continue the theme: which represents that the spherical cavity prefers B-monomers and repels A-and C-monomers. Two minority components are no longer constrained but freely interact with the middle component inside the cavity; thus, the structures observed here can be diversified and intricate. On the basis of the structural features and symmetries observed in the 3D patterns, we divide the morphologies into four groups displayed in Figure 7. structural features and symmetries observed in the 3D patterns, we divide the morphologies into four groups displayed in Figure 7. (a): Two (or three)-coloured sphere (TCS). Four structures belong to TCS group displayed in Figure 7a. The TCS-I structure is a combination of an A-rich hemisphere and a C-rich hemisphere. In contrast to the MS structure, TCS-I is covered by a B-filled wetting layer. The TCS-II structure replicates the appearance of TCS-I with a small B-rich core that exists in the sphere centre. The core (a): Two (or three)-coloured sphere (TCS). Four structures belong to TCS group displayed in Figure 7a. The TCS-I structure is a combination of an A-rich hemisphere and a C-rich hemisphere. In contrast to the MS structure, TCS-I is covered by a B-filled wetting layer. The TCS-II structure replicates the appearance of TCS-I with a small B-rich core that exists in the sphere centre. The core region enlarges and transforms into a disk that separates the A-and C-rich regions in the TCS-III structure. The three TCS structures are different in the B-rich region; however, these structures are homologous and follow the same formation mechanism. Li and Shi et al. observed Janus nanoparticles formed by A 8 B 8 C 8 triblock copolymer confined in small nanopores by using a simulated annealing technique [62]. The Janus nanoparticles have echoes of the formation mechanism from our TCS structures and are considered in single-lamellar structures [63]. TCS-IV contains three disks, of which the C-rich disk is located in the centre region and the A-rich disks on the two sides.
(b): Disk ring (DR). In Figure 7b, from the outward appearance of DR-I, it has the same layered construction as that of TCS-IV. However, C-monomers form a ring instead of a disk and encircle a B-rich cylinder. In a DR-II structure, The B-rich cylinder expands, and several monomers slip the leash from the C-rich ring. In the DR-III structure, the core region contains basically B monomers, the C-rich region forms two rings and the A-rich region forms a large ring accompanied by two small disks. On the basis of DR-III, a new core located inside the B-rich region is formed by A-monomers with the profile map, that is, DR-IV.
(c): Sphere-ring (SR) structure. In Figure 7c, SR-I with axial symmetry includes two identical block domains, in which two disks and two rings formed by A-and C-monomers embed each other. B-monomers gather into a central spheroidal region. When A and C block domains are separated from each other and no longer pack tightly together given environmental change, B-monomers spill out into the interstices, and a minority block forms a solid sphere inside the B-rich region, with which the structure transforms to SR-II. In the SR-III structure, B-monomers form a spherical core, and the Aand C-rich regions form two rings accompanied by a small disk. Compared with SR-I, a four-layer structure turns into a six-layer structure of SR-III, in which each of the two minority blocks forms one more parallel ring. The uniform SR structure can also be assembled by ABC star triblock copolymers with f A : f B : f C = 0.3 : 0.3 : 0.4 in a preferential spherical cavity [54]. An MC study has reported a similar stacked lamella with alternate A and B domains of ABA triblock copolymers under soft confinement [64]. Obviously layered structures are common in a spherical space. SR-IV has the same layered system which contains two disks, that is, two A-rich rings and two C-rich rings with SR-III. The only structural difference is that an additional A-or C-formed sphere occupies the central position of a B-formed core. SR-II and SR-IV structures with newly formed cores are novel discoveries in our study.
(d): Rotational structure (RS), treble texture (TT) structure and QT structure. In Figure 7d, RS-I includes three uniform discrete patches mounted on one large triquetrum core with a threefold rotational symmetry. RS-II and RS-III have the same geometrical framework with a fourfold rotational symmetry. C-block forms a hollow shell with an approximate shape of a cube, in which the central region is filled with B-monomers. A-block forms six patches uniformly inlaying in the surface of the C-formed cube. In addition, an A-rich core exists in RS-III but none in RS-II. In a study of cylinder-forming diblock copolymers, Liang and Shi et al. also found perforated cubic structures [65]. In the TT structure, three arcs formed by a C-block connect the two poles and incise the surface into three identical domains. The A-rich region forms a triangular domain centred in a B-rich matrix with a threefold rotational symmetry. After double checking, we still observe the QT structure as a stable phase under the B-attractive surfaces. In a previous study, virus-like nanoparticles with three and four patches, which correspond to TT and QT structures, have been reported [62]. The overall patchy morphology is largely influenced by the volume fraction of the middle block in linear triblock copolymers. Figure 8 plots the location of the 17 morphologies which are different from those under the opposite attractive surface field. Among these morphologies, only the QT structure is not destroyed by the B-attractive surface field but is stable in a small phase region, where d eff = 7.4 R g and −λ BW < 0.3. For the convenience of observation, the morphologies are represented by the shrunken profile images in the phase diagram. A division by different morphology regions is implemented. Homologous structures, such as TCS-I, II and III, belong to the same area. These structures are presented in the bottom-left region of the phase diagram, in which the spherical volume and adsorption strength are small. With the increase in the strength of the wall preference, TCS-IV is observed in a relatively wide range of −λ BW ≈ 0.4 − 0.8 region; with the increase in the diameter of the spherical cavity, DR-I and DR-II are observed in a narrow range of −λ BW ≈ 0.1 − 0.2 region. A TT phase point appears between QT and DR-II phase under severe existing conditions, −λ BW = 0.1 and d eff = 6.8 R g . The two columnar structures, namely, TT and QT, are written in water in our phase space. With the further increase in the spherical volume, QT transforms into SR-II located on the top-left region. SR-I, III, IV, RS-II and III gather in the middle region of the phase diagram. Among these morphologies, SR-I requires a relatively higher confinement degree than others, and SR-III and RS-II remain stable when −λ BW continuously fluctuates upward. However, SR-I cannot bear the powerful adsorption force and has its structure reorganised. A growing −λ BW with d eff ≈ 3.8 − 5.0 R g leads to the introduction of RS-I and DR-III. DR-IV is a product of large spherical cavity, in which d eff ≥ 8.0 R g . Not every DR-IV structure observed in our work seems as regular as the illustration in Figure 7b. From an overall perspective, the DR-IV structures obtained when −λ BW ≥ 0.6 have an alternative five-layer structure; however, the rings are distorted to varying degrees because the large confinement space and strong surface field are comparable to the interaction among blocks. This phase diagram covers numerous complex dynamic processes, which demonstrate that the morphologies of linear triblock copolymers are sensitive to the AC-repulsive surfaces. left region. SR-I, III, IV, RS-II and III gather in the middle region of the phase diagram. Among these morphologies, SR-I requires a relatively higher confinement degree than others, and SR-III and RS-II remain stable when − continuously fluctuates upward. However, SR-I cannot bear the powerful adsorption force and has its structure reorganised. A growing − with ≈ 3.8 − 5.0 leads to the introduction of RS-I and DR-III. DR-IV is a product of large spherical cavity, in which ≥ 8.0 . Not every DR-IV structure observed in our work seems as regular as the illustration in Figure   7b. From an overall perspective, the DR-IV structures obtained when − ≥ 0.6 have an alternative five-layer structure; however, the rings are distorted to varying degrees because the large confinement space and strong surface field are comparable to the interaction among blocks. This phase diagram covers numerous complex dynamic processes, which demonstrate that the morphologies of linear triblock copolymers are sensitive to the AC-repulsive surfaces. We study the transition from RS-II to SR-III along the phase path = 6.8 as an example to expose the dynamic process. Figure 9 demonstrates the internal energy ⁄ and entropic energy − ⁄ as functions of the B-attractive strength − . In Figure 9a, the internal energy of RS-II decreases progressively, but the internal energy curve of SR-III has an inflection point, which provides favourable RS-II at a large surface field. By contrast, RS-II has a favourable entropic contributionat a small surface field of approximately − < 0.463, which dominates the chain configuration, as exhibited in Figure 9b. The preferential surfaces play an important role in the phase behaviour reflected by the variations in internal and entropy energies. We study the transition from RS-II to SR-III along the phase path d eff = 6.8 R g as an example to expose the dynamic process. Figure 9 demonstrates the internal energy U/nk B T and entropic energy −S/nk B as functions of the B-attractive strength −λ BW . In Figure 9a, the internal energy U of RS-II decreases progressively, but the internal energy curve of SR-III has an inflection point, which provides favourable RS-II U at a large surface field. By contrast, RS-II has a favourable entropic contribution −S at a small surface field of approximately −λ BW < 0.463, which dominates the chain configuration, as exhibited in Figure 9b. The preferential surfaces play an important role in the phase behaviour reflected by the variations in internal and entropy energies.  , that contributes to the free energy for RS-II and SR-III phases.

Discussion
We studied the self-assembly of bulk cylinder-forming ABC triblock copolymers under a spherical confinement by using a real-space SCFT. We detailedly discussed the spontaneous morphologies and phase diagrams on the degree of spherical confinement and the strength of the surface preference characterised by parameters and , correspondingly. The bulk tetragonally arranged a cylinder structure, and the chain conformations were described by analysing its segment distributions. A rigid spherical cavity with a finite thickness was considered, and we positioned the cylinder-forming copolymers to observe its recombination in a neutral surface. We then supplemented varying AC-attractive and B-repulsive potentials on the shield wall and discovered several novel structures obtained in − space. Furthermore, we aimed to find several distinctive phase behaviours in an opposite interfacial adsorption potential; in particular, a series of complex and correlative structures with spatial symmetries was demonstrated.
For copolymers under neutral confinement, six atypical phases, namely, WOS, SCL, DLC, ST, QT and IR, were considered stable phases exhibited sequentially in the 1D phase diagram, associated with a changing diameter of a spherical cavity. Systems located at the largest middle range of a length scale in the phase diagram assembled into the SCL structure, thus showing a strong symmetry. The complexity of the internal structure increased with the augmentation of the confinement volume and the frequent transformation of phase structure. On account of the masking technique, the wall thickness might affect the morphologies to a certain extent in our work. On the AC-attractive surface, five new phases, namely, MS, FR, DLL, SF-I and SF-II, were found; these phases had no relation with the cylinder structure. Among these phases, the SF-type phase presented a relatively strong stability in a large spherical space. In addition, the SCL phase remained and filled a large quota in Figure 5b together with the MS phase.
On the AC-repulsive surface, various symmetrical morphologies were demonstrated; among these morphologies, the layered structures were the mainstream. The homologous structures, such as SR-I and SR-II, were derived from the sensitivity of a polymer reconstruction for short-chainrepulsive potential. Several structures, such as DR-III and SR-III, were also confirmed in other polymer systems and external environment. We analysed the free energies and internal and entropic

Discussion
We studied the self-assembly of bulk cylinder-forming ABC triblock copolymers under a spherical confinement by using a real-space SCFT. We detailedly discussed the spontaneous morphologies and phase diagrams on the degree of spherical confinement and the strength of the surface preference characterised by parameters d eff and λ iW , correspondingly. The bulk tetragonally arranged a cylinder structure, and the chain conformations were described by analysing its segment distributions. A rigid spherical cavity with a finite thickness was considered, and we positioned the cylinder-forming copolymers to observe its recombination in a neutral surface. We then supplemented varying AC-attractive and B-repulsive potentials on the shield wall and discovered several novel structures obtained in λ − d eff space. Furthermore, we aimed to find several distinctive phase behaviours in an opposite interfacial adsorption potential; in particular, a series of complex and correlative structures with spatial symmetries was demonstrated.
For copolymers under neutral confinement, six atypical phases, namely, WOS, SCL, DLC, ST, QT and IR, were considered stable phases exhibited sequentially in the 1D phase diagram, associated with a changing diameter of a spherical cavity. Systems located at the largest middle range of a length scale in the phase diagram assembled into the SCL structure, thus showing a strong symmetry. The complexity of the internal structure increased with the augmentation of the confinement volume and the frequent transformation of phase structure. On account of the masking technique, the wall thickness might affect the morphologies to a certain extent in our work. On the AC-attractive surface, five new phases, namely, MS, FR, DLL, SF-I and SF-II, were found; these phases had no relation with the cylinder structure. Among these phases, the SF-type phase presented a relatively strong stability in a large spherical space. In addition, the SCL phase remained and filled a large quota in Figure 5b together with the MS phase.
On the AC-repulsive surface, various symmetrical morphologies were demonstrated; among these morphologies, the layered structures were the mainstream. The homologous structures, such as SR-I and SR-II, were derived from the sensitivity of a polymer reconstruction for short-chain-repulsive potential. Several structures, such as DR-III and SR-III, were also confirmed in other polymer systems and external environment. We analysed the free energies and internal and entropic energies of adjacent phases and verified the first-order phase transitions. The influence of preferential surface on entropy energy mainly optimised the phase transitions. Our calculation results would contribute to fabricating advanced polymer materials and could be used for reference in subsequent research.