Self-consistent field theory simulations of polymers on arbitrary domains

We introduce a framework for simulating the mesoscale self-assembly of block copolymers in arbitrary confined geometries subject to Neumann boundary conditions. We employ a hybrid finite difference/volume approach to discretize the mean-field equations on an irregular domain represented implicitly by a level-set function. The numerical treatment of the Neumann boundary conditions is sharp, i.e. it avoids an artificial smearing in the irregular domain boundary. This strategy enables the study of self-assembly in confined domains and enables the computation of physically meaningful quantities at the domain interface. In addition, we employ adaptive grids encoded with Quad-/Oc-trees in parallel to automatically refine the grid where the statistical fields vary rapidly as well as at the boundary of the confined domain. This approach results in a significant reduction in the number of degrees of freedom and makes the simulations in arbitrary domains using effective boundary conditions computationally efficient in terms of both speed and memory requirement. Finally, in the case of regular periodic domains, where pseudo-spectral approaches are superior to finite differences in terms of CPU time and accuracy, we use the adaptive strategy to store chain propagators, reducing the memory footprint without loss of accuracy in computed physical observables.

[1]  Matsen.  Stabilizing new morphologies by blending homopolymer with block copolymer. , 1995, Physical review letters.

[2]  Kris T. Delaney,et al.  Shape optimization for DSA , 2016, SPIE Advanced Lithography.

[3]  Koji Asakawa,et al.  2.5-inch disk patterned media prepared by an artificially assisted self-assembling method , 2002 .

[4]  Henry I. Smith,et al.  Fabrication of nanostructures with long-range order using block copolymer lithography , 2002 .

[5]  Steven G. Johnson,et al.  The Design and Implementation of FFTW3 , 2005, Proceedings of the IEEE.

[6]  G. Kalosakas,et al.  Improved numerical algorithm for exploring block copolymer mesophases , 2002 .

[7]  S. P. Lloyd,et al.  Least squares quantization in PCM , 1982, IEEE Trans. Inf. Theory.

[8]  G. Strang.  On the Construction and Comparison of Difference Schemes , 1968 .

[9]  G. Fredrickson,et al.  Microdomain Ordering in Laterally Confined Block Copolymer Thin Films , 2007 .

[10]  Frédéric Gibou,et al.  Robust second-order accurate discretizations of the multi-dimensional Heaviside and Dirac delta functions , 2008, J. Comput. Phys..

[11]  P. Gennes,et al.  Effects of polymer solutions on colloid stability , 1979 .

[12]  Len Czuba.  Application of Plastics in Medical Devices and Equipment , 2014, Handbook of Polymer Applications in Medicine and Medical Devices.

[13]  Pingwen Zhang,et al.  Boundary effects in confined copolymer system and compressible SCFT model , 2014, J. Comput. Appl. Math..

[14]  E. Helfand.  Theory of inhomogeneous polymers: Fundamentals of the Gaussian random‐walk model , 1975 .

[15]  E. Helfand,et al.  Block Copolymer Theory. III. Statistical Mechanics of the Microdomain Structure , 1975 .

[16]  M. Matsen.  Phase Behavior of Block Copolymer/Homopolymer Blends , 1995 .

[17]  Kris T. Delaney,et al.  Thermodynamic and kinetic aspects of defectivity in directed self‐assembly of cylinder‐forming diblock copolymers in laterally confining thin channels , 2014 .

[18]  Kris T. Delaney,et al.  Field-theoretic simulations of directed self-assembly in cylindrical confinement: placement and rectification aspects , 2014, Advanced Lithography.

[19]  Laurence W. McKeen,et al.  3 – Plastics Used in Medical Devices , 2014 .

[20]  S. Osher,et al.  A level set approach for computing solutions to incompressible two-phase flow , 1994 .

[21]  C. M. Bates,et al.  Multiblock Polymers: Panacea or Pandora’s Box? , 2012, Science.

[22]  F. Gibou.  A fast hybrid k-means level set algorithm for segmentation , 2005 .

[23]  Frédéric Gibou,et al.  A level set approach for diffusion and Stefan-type problems with Robin boundary conditions on quadtree/octree adaptive Cartesian grids , 2013, J. Comput. Phys..

[24]  J. Sethian,et al.  Fronts propagating with curvature-dependent speed: algorithms based on Hamilton-Jacobi formulations , 1988 .

[25]  E. Kramer,et al.  SCFT Simulations of Thin Film Blends of Block Copolymer and Homopolymer Laterally Confined in a Square Well , 2009 .

[26]  Héctor D. Ceniceros,et al.  Spectral collocation methods for polymer brushes. , 2011, The Journal of chemical physics.

[27]  G. Fredrickson.  The equilibrium theory of inhomogeneous polymers , 2005 .

[28]  R. Segalman.  Patterning with block copolymer thin films , 2005 .

[29]  Steven G. Johnson,et al.  The Fastest Fourier Transform in the West , 1997 .

[30]  Carsten Burstedde,et al.  p4est: Scalable Algorithms for Parallel Adaptive Mesh Refinement on Forests of Octrees , 2011, SIAM J. Sci. Comput..

[31]  J. Strain.  Tree Methods for Moving Interfaces , 1999 .

[32]  Frank S. Bates,et al.  Origins of Complex Self-Assembly in Block Copolymers , 1996 .

[33]  M. Matsen.  The standard Gaussian model for block copolymer melts , 2002 .

[34]  Frédéric Gibou,et al.  A second order accurate level set method on non-graded adaptive cartesian grids , 2007, J. Comput. Phys..

[35]  Mohammad Mirzadeh,et al.  Parallel level-set methods on adaptive tree-based grids , 2016, J. Comput. Phys..

[36]  Frédéric Gibou,et al.  Efficient symmetric discretization for the Poisson, heat and Stefan-type problems with Robin boundary conditions , 2010, J. Comput. Phys..

[37]  Frédéric Gibou,et al.  Geometric integration over irregular domains with application to level-set methods , 2007, J. Comput. Phys..

[38]  Frédéric Gibou,et al.  A supra-convergent finite difference scheme for the variable coefficient Poisson equation on non-graded grids , 2006, J. Comput. Phys..

[39]  J. Tukey,et al.  An algorithm for the machine calculation of complex Fourier series , 1965 .