A novel self-consistent-field lattice model for block copolymers

We develop a self-consistent-field lattice model for block copolymers and propose a novel and general method to solve the self-consistent-field equations. The approach involves describing the polymer chains in a lattice and employing a two-stage relaxation procedure to evolve a system as rapidly as possible to a free-energy minimum. In order to test the validity of this approach, we use the method to study the microphases of rod-coil diblock copolymers. In addition to the lamellar and cylindrical morphologies, micellar, perforated lamellar, gyroid, and zigzag structures have been identified without any prior assumption of the microphase symmetry. Furthermore, this approach can also give the possible orientation of the rods in different structures.

Mean field lattice model for adsorption isotherms in zeolite NaA

Using a lattice model for adsorption in microporous materials, pure component adsorption isotherms are obtained within a mean field approximation for methane at 300 K and xenon at 300 and 360 K in zeolite NaA. It is argued that the increased repulsive adsorbate-adsorbate interactions at high coverages must play an important role in determining the adsorption behavior. Therefore, this feature is incorporated through a "coverage-dependent interaction'' model, which introduces a free, adjustable parameter. Another important feature, the site volume reduction, has been treated in two ways: a van der Waal model and a 1D hard-rod theory [van Tassel et al., AIChE J. 40, 925 (1994)]; we have also generalized the latter to include all possible adsorbate overlap scenarios. In particular, the 1D hard-rod model, with our coverage-dependent interaction model, is shown to be in best quantitative agreement with the previous grand canonical Monte Carlo isotherms. The expressions for the isosteric heats of adsorption indicate that attractive and repulsive adsorbate-adsorbate interactions increase and decrease the heats of adsorption, respectively. It is concluded that within the mean field approximation, our simple model for repulsive interactions and the 1D hard-rod model for site volume reduction are able to capture most of the important features of adsorption in confined regions. (C) 1999 American Institute of Physics. [S0021-9606(99)70515-5].

Study of self-assembly of symmetric coil-rod-coil ABA-type triblock copolymers by self-consistent field lattice method

The self-assembly of symmetric coil-rod-coil ABA-type triblock copolymer melts is studied by applying self-consistent field lattice techniques in a three-dimensional space. The self-assembled ordered structures differ significantly with the variation of the volume fraction of the rod component, which include lamellar, wave lamellar, gyroid, perforated lamellar, cylindrical, and spherical-like phases. To understand the physical essence of these phases and the regimes of occurrence, we construct the phase diagram, which matches qualitatively with the existing experimental results. Compared with the coil-rod AB diblock copolymer, our results revealed that the interfacial grafting density of the separating rod and coil segments shows important influence on the self-assembly behaviors of symmetric coil-rod-coil ABA triblock copolymer melts. We found that the order-disorder transition point changes from f(rod)=0.5 for AB diblock copolymers to f(rod)=0.6 for ABA triblock copolymers. Our results also show that the spherical-like and cylindrical phases occupy most of the region in the phase diagram, and the lamellar phase is found stable only at the high volume fraction of the rod.

Toda and Volterra lattice equations from discrete symmetries of KP hierarchies

The discrete models of the Toda and Volterra chains are being constructed out of the continuum two-boson KP hierarchies. The main tool is the discrete symmetry preserving the Hamiltonian structure of the continuum models. The two-boson currents of KP hierarchy are being associated with sites of the corresponding chain by successive actions of discrete symmetry.

Temperature and pressure dependence of phase separation of trans-decahydronaphthalene/polystyrene solution

The cloud-point temperatures (T-c1's) of ti-ans-decahydronaphthalene (TD)/polystyrene (PS, M-w = 270 kg/mol) solutions were determined by fight scattering measurements over a range of temperatures (1-16 degreesC), pressures (100-900 bar), and compositions (4.2-21.6 vol% polymer). The system phase separates upon cooling and the T-c1 was found to increase with the rising pressure for the constant composition. In the absence of special effects this finding indicates positive excess volumes. The special attention was paid to the demixing temperatures as a function of the pressure for the different polymer solutions and the plots in the T-volume fraction plane and P-volume fraction plane. The cloud-point curves of polymer solutions under changing pressures were observed for different compositions, demonstrates that the TD/PS system exhibits UCST (phase separation upon cooling) behavior. With this data the phase diagrams under pressure were calculated applying the Sanchez-Lacombe (SL) lattice fluid theory. Furthermore, the cause of phase separation, i.e., the influence of Flory-Huggins (FH) interaction parameter under pressure was investigated.

Thermodynamics of phase behavior in PEO/P(EO-b-DMS) homopolymer and block co-oligomer mixtures under pressure

The cloud-point temperatures (T-cl's) of poly(ethylene oxide) (PEO) and poly(ethylene oxide)-block-polydimethylsiloxane (P(EO-b-DMS)) homopolymer and block-oligomer mixtures were determined by turbidity measurements over a range of temperatures (105 to 130degrees), pressures (1 to 800 bar), and compositions (10-40 wt.-% PEO). The system phase separates upon cooling and T-cl was found to decrease with an increase in pressure for a constant composition. In the absence of special effects, this finding indicates negative excess volumes. Special attention was paid to the demixing temperatures as a function of the pressure for the different polymer mixtures and the plots in the T-phi plane (where phi signifies volume fractions). The cloud-point curves of the polymer mixture under pressures were observed for different compositions. The Sanchez-Lacombe (SL) lattice fluid theory was used to calculate the spinodals, the binodals, the Flory-Huggins (FH) interaction parameter, the enthalphy of mixing, and the volume changes of mixing. The calculated results show that modified P(EO-b-DMS) scaling parameters with the new combining rules can describe the thermodynamics of the PEO/P(EO-b-DMS) system well with the SL theory.

Effects of pressure and molecular weight on the miscibility of polystyrene and cyclohexane

With the aid of Sanchez-Lacombe lattice fluid theory (SLLFT), the phase diagrams were calculated for the system cyclohexane (CH)/polystyrene (PS) with different molecular weights at different pressures. The experimental data is in reasonable agreement with SLLFT calculations. The total Gibbs interaction energy, g*(12) for different molecular weights PS at different pressures was expressed, by means of a universal relationship, as g(12)* =f(12)* + (P - P-0) nu*(12) demixing curves were then calculated at fixed (near critical) compositions of CH and PS systems for different molecular weights. The pressures of optimum miscibility obtained from the Gibbs interaction energy are close to those measured by Wolf and coworkers. Furthermore, a reasonable explanation was given for the earlier observation of Saeki et al., i.e., the phase separation temperatures of the present system increase with the increase of pressure for the low molecular weight of the polymer whereas they decrease for the higher molecular weight polymers. The effects of molecular weight, pressure, temperature and composition on the Flory Huggins interaction parameter can be described by a general equation resulting from fitting the interaction parameters by means of Sanchez-Lacombe lattice fluid theory.

Correlation equations for average deposition rate coefficients of nanoparticles in a cylindrical pore

Nanoparticle deposition behavior observed at the Darcy scale represents an average of the processes occurring at the pore scale. Hence, the effect of various pore-scale parameters on nanoparticle deposition can be understood by studying nanoparticle transport at pore scale and upscaling the results to the Darcy scale. In this work, correlation equations for the deposition rate coefficients of nanoparticles in a cylindrical pore are developed as a function of nine pore-scale parameters: the pore radius, nanoparticle radius, mean flow velocity, solution ionic strength, viscosity, temperature, solution dielectric constant, and nanoparticle and collector surface potentials. Based on dominant processes, the pore space is divided into three different regions, namely, bulk, diffusion, and potential regions. Advection-diffusion equations for nanoparticle transport are prescribed for the bulk and diffusion regions, while the interaction between the diffusion and potential regions is included as a boundary condition. This interaction is modeled as a first-order reversible kinetic adsorption. The expressions for the mass transfer rate coefficients between the diffusion and the potential regions are derived in terms of the interaction energy profile. Among other effects, we account for nanoparticle-collector interaction forces on nanoparticle deposition. The resulting equations are solved numerically for a range of values of pore-scale parameters. The nanoparticle concentration profile obtained for the cylindrical pore is averaged over a moving averaging volume within the pore in order to get the 1-D concentration field. The latter is fitted to the 1-D advection-dispersion equation with an equilibrium or kinetic adsorption model to determine the values of the average deposition rate coefficients. In this study, pore-scale simulations are performed for three values of Peclet number, Pe = 0.05, 5, and 50. We find that under unfavorable conditions, the nanoparticle deposition at pore scale is best described by an equilibrium model at low Peclet numbers (Pe = 0.05) and by a kinetic model at high Peclet numbers (Pe = 50). But, at an intermediate Pe (e.g., near Pe = 5), both equilibrium and kinetic models fit the 1-D concentration field. Correlation equations for the pore-averaged nanoparticle deposition rate coefficients under unfavorable conditions are derived by performing a multiple-linear regression analysis between the estimated deposition rate coefficients for a single pore and various pore-scale parameters. The correlation equations, which follow a power law relation with nine pore-scale parameters, are found to be consistent with the column-scale and pore-scale experimental results, and qualitatively agree with the colloid filtration theory. These equations can be incorporated into pore network models to study the effect of pore-scale parameters on nanoparticle deposition at larger length scales such as Darcy scale.

Field Equilibrium Finite Elements for Structural Analysis

Many finite elements used in structural analysis possess deficiencies like shear locking, incompressibility locking, poor stress predictions within the element domain, violent stress oscillation, poor convergence etc. An approach that can probably overcome many of these problems would be to consider elements in which the assumed displacement functions satisfy the equations of stress field equilibrium. In this method, the finite element will not only have nodal equilibrium of forces, but also have inner stress field equilibrium. The displacement interpolation functions inside each individual element are truncated polynomial solutions of differential equations. Such elements are likely to give better solutions than the existing elements.In this thesis, a new family of finite elements in which the assumed displacement function satisfies the differential equations of stress field equilibrium is proposed. A general procedure for constructing the displacement functions and use of these functions in the generation of elemental stiffness matrices has been developed. The approach to develop field equilibrium elements is quite general and various elements to analyse different types of structures can be formulated from corresponding stress field equilibrium equations. Using this procedure, a nine node quadrilateral element SFCNQ for plane stress analysis, a sixteen node solid element SFCSS for three dimensional stress analysis and a four node quadrilateral element SFCFP for plate bending problems have been formulated.For implementing these elements, computer programs based on modular concepts have been developed. Numerical investigations on the performance of these elements have been carried out through standard test problems for validation purpose. Comparisons involving theoretical closed form solutions as well as results obtained with existing finite elements have also been made. It is found that the new elements perform well in all the situations considered. Solutions in all the cases converge correctly to the exact values. In many cases, convergence is faster when compared with other existing finite elements. The behaviour of field consistent elements would definitely generate a great deal of interest amongst the users of the finite elements.

Tokamak magnetic field lines described by simple maps

The magnetic field line structure in a tokamak can be obtained by direct numerical integration of the field line equations. However, this is a lengthy procedure and the analysis of the solution may be very time-consuming. Otherwise we can use simple two-dimensional, area-preserving maps, obtained either by approximations of the magnetic field line equations, or from dynamical considerations. These maps can be quickly iterated, furnishing solutions that mirror the ones obtained from direct numerical integration, and which are useful when long-term studies of field line behavior are necessary (e.g. in diffusion calculations). In this work we focus on a set of simple tokamak maps for which these advantages are specially pronounced.

Free expansion of fermionic dark solitons in a boson-fermion mixture

We use a time-dependent dynamical mean-field-hydrodynamic model to study the formation of fermionic dark solitons in a trapped degenerate Fermi gas mixed with a Bose-Einstein condensate in a harmonic as well as a periodic optical-lattice potential. The dark soliton with a 'notch' in the probability density with a zero at the minimum is simulated numerically as a nonlinear continuation of the first vibrational excitation of the linear mean-field-hydrodynamic equations, as suggested recently for pure bosons. We study the free expansion of these dark solitons as well as the consequent increase in the size of their central notch and discuss the possibility of experimental observation of the notch after free expansion.

Discrete breathers in a two-dimensional Fermi-Pasta-Ulam lattice

Using asymptotic methods, we investigate whether discrete breathers are supported by a two-dimensional Fermi-Pasta-Ulam lattice. A scalar (one-component) two-dimensional Fermi-Pasta-Ulam lattice is shown to model the charge stored within an electrical transmission lattice. A third-order multiple-scale analysis in the semi-discrete limit fails, since at this order, the lattice equations reduce to the (2+1)-dimensional cubic nonlinear Schrödinger (NLS) equation which does not support stable soliton solutions for the breather envelope. We therefore extend the analysis to higher order and find a generalised $(2+1)$-dimensional NLS equation which incorporates higher order dispersive and nonlinear terms as perturbations. We find an ellipticity criterion for the wave numbers of the carrier wave. Numerical simulations suggest that both stationary and moving breathers are supported by the system. Calculations of the energy show the expected threshold behaviour whereby the energy of breathers does {\em not} go to zero with the amplitude; we find that the energy threshold is maximised by stationary breathers, and becomes arbitrarily small as the boundary of the domain of ellipticity is approached.

Topology of toroidal helical fields in non-circular cross-sectional tokamaks

The ordinary differential magnetic field line equations are solved numerically; the tokamak magnetic structure is studied on Hefei Tokamak-7 Upgrade (HT-7U) when the equilibrium field with a monotonic q-profile is perturbed by a helical magnetic field. We find that a single mode (m, n) helical perturbation can cause the formation of islands on rational surfaces with q = m/n and q = (m +/- 1, +/- 2, +/- 3,...)/n due to the toroidicity and plasma shape (i.e. elongation and triangularity), while there are many undestroyed magnetic surfaces called Kolmogorov-Arnold-Moser (KAM) barriers on irrational surfaces. The islands on the same rational surface do not have the same size. When the ratio between the perturbing magnetic field B-r(r) and the toroidal magnetic field amplitude B(phi)0 is large enough, the magnetic island chains on different rational surfaces will overlap and chaotic orbits appear in the overlapping area, and the magnetic field becomes stochastic. It is remarkable that the stochastic layer appears first in the plasma edge region.

Self-Assembly of Rod-Terminally Tethered Three-Armed Star-Shaped Coil Block Copolymer: Investigation of the Presence of the Branching in the Coil to the Self-Assembled Behavior

Self-assembled behavior of rod-terminally tethered three-armed star-shaped coil block copolymer melts was studied by applying self-consistent-field lattice techniques in three-dimensional (3D) space. Similar to rod-coil diblock copolymers, five morphologies were observed, i.e., lamellar, perforated lamellar, gyroidlike, cylindrical and sphericallike structures, while the distribution of the morphologies in the phase diagram was dramatically changed with respect to that Of rod-coil diblock copolymers.