Continuous configuration time-dependent self-consistent field method for polyatomic quantum dynamical problems

A new continuous configuration time-dependent self-consistent field method has been developed to study polyatomic dynamical problems by using the discrete variable representation for the reaction system, and applied to a reaction system coupled to a bath. The method is very efficient because the equations involved are as simple as those in the traditional single configuration approach, and can account for the correlations between the reaction system and bath modes rather well. (C) American Institute of Physics.

Self-consistent microscopic treatment of the effects of self-motion of the probe on ionic and dipolar solvation dynamics

A simple but self-consistent microscopic theory for the time dependent solvation energy of both ions and dipoles is presented which includes, for the first time, the details of the self-motion of the probe on its own solvation dynamics. The theory leads to several interesting predictions. The most important of them is that, for dipolar solvation, both the rotational and the translational motions of the dipolar solute probe can significantly accelerate the rate of solvation. In addition, the rotational self-motion of the solute can also give rise to an additional mechanism of nonexponentiality in solvation time correlation functions in otherwise slow liquids. A comparison between the present theoretical predictions and the recent experimental studies of Maroncelli et al. on solvation dynamics of aniline in l-propanol seems to indicate that the said experiments have missed the initial solvent response up to about 45 ps. After mapping the experimental results on the redefined time scale, the theoretical results can explain the experimental results for solvation of aniline in 1-propanol very well. For ionic solvation, the translational motion is significant for light solutes only. For example, for Li+ in water, translational motion speeds up the solvation by about 20%. The present theory demonstrates that in dipolar solvation the partial quenching of the self-motion due to the presence of specific solute-solvent interactions (such as H-bonding) may lead to a much slower solvation than that when the self-motion is present. This point has been discussed. In addition, we present the theoretical results for solvation of aniline in propylene carbonate, Here, the solvation is predicted to be complete within 15-20 ps.

Effect of Solvent Molecular Size on the Self-Assembly of Amphiphilic Diblock Copolymer in Selective Solvent

Real-space self-consistent field theory (SCFT) is employed to study the effect of solvent molecular size on the self-assembly of amphiphilic diblock copolymer in selective solvent. The phase diagrams in wide ranges of interaction parameters and solvent molecular size were obtained in present study. The results indicate that the solvent molecular size is a key factor that determines the self-assembly of amphiphilic diblock copolymer. The self-assembled morphology changes from circle-like micelle to line-like micelle, then to loop-like micelle by decreasing the solvent molecular size in a wide range of solvent selectivity. We analyze and discuss this change in terms of the solvent solubility and the entropy contribution.

Finite-stretching corrections to the Milner-Witten-Cates theory for polymer brushes

This paper investigates finite-stretching corrections to the classical Milner-Witten-Cates theory for semi-dilute polymer brushes in a good solvent. The dominant correction to the free energy originates from an entropic repulsion caused by the impenetrability of the grafting surface, which produces a depletion of segments extending a distance $\mu \propto L^{-1}$ from the substrate, where $L$ is the classical brush height. The next most important correction is associated with the translational entropy of the chain ends, which creates the well-known tail where a small population of chains extend beyond the classical brush height by a distance $\xi \propto L^{-1/3}$. The validity of these corrections is confirmed by quantitative comparison with numerical self-consistent field theory.

Converting the nanodomains of a diblock-copolymer thin film from spheres to cylinders with an external electric field

We investigate the ability of an applied electric field to convert the morphology of a diblock-copolymer thin film from a monolayer of spherical domains embedded in the matrix to cylindrical domains that penetrate through the matrix. As expected, the applied field increases the relative stability of cylindrical domains, while simultaneously reducing the energy barrier that impedes the transition to cylinders. The effectiveness of the field is enhanced by a large dielectric contrast between the two block-copolymer components, particularly when the low-dielectric contrast component forms the matrix. Furthermore, the energy barrier is minimized by selecting sphere-forming diblock copolymers that are as compositionally symmetric as possible. Our calculations, which are the most quantitatively reliable to date, are performed using a numerically precise spectral algorithm based on self-consistent-field theory supplemented with an exact treatment for linear dielectric materials.

Electric field alignment in thin films of cylinder-forming diblock copolymer

We investigate thin films of cylinder-forming diblock copolymer confined between electrically charged parallel plates, using self-consistent-field theory ( SCFT) combined with an exact treatment for linear dielectric materials. Our study focuses on the competition between the surface interactions, which tend to orient cylinder domains parallel to the plates, and the electric field, which favors a perpendicular orientation. The effect of the electric field on the relative stability of the competing morphologies is demonstrated with equilibrium phase diagrams, calculated with the aid of a weak-field approximation. As hoped, modest electric fields are shown to have a significant stabilizing effect on perpendicular cylinders, particularly for thicker films. Our improved SCFT-based treatment removes most of the approximations implemented by previous approaches, thereby managing to resolve outstanding qualitative inconsistencies among different approximation schemes.

Undulation instability in block-copolymer lamellae subjected to a perpendicular electric field

We examine the stability of lamellar stacks in the presence of an electric field, E-0, applied normal to the lamellae. Calculations are performed with self-consistent field theory (SCFT) supplemented by an exact treatment of the electrostatic energy for linear dielectric materials. The calculations identify a critical electric field, E-0*, beyond which the lamellar stack becomes unstable with respect to undulations. This E-0* rapidly decreases towards zero as the number of lamellae in the stack diverges. Our quantitative predictions for E-0* are consistent with previous experimental measurements by Xu and co-workers.

On the ultrafast charge migration and subsequent charge directed reactivity in Cl center dot center dot center dot N halogen-bonded clusters following vertical ionization

In this article, we have presented ultrafast charge transfer dynamics through halogen bonds following vertical ionization of representative halogen bonded clusters. Subsequent hole directed reactivity of the radical cations of halogen bonded clusters is also discussed. Furthermore, we have examined effect of the halogen bond strength on the electron-electron correlation-and relaxation-driven charge migration in halogen bonded complexes. For this study, we have selected A-Cl (A represents F, OH, CN, NH2, CF3, and COOH substituents) molecules paired with NH3 (referred as ACl:NH3 complex): these complexes exhibit halogen bonds. To the best of our knowledge, this is the first report on purely electron correlation-and relaxation-driven ultrafast (attosecond) charge migration dynamics through halogen bonds. Both density functional theory and complete active space self-consistent field theory with 6-31+G(d, p) basis set are employed for this work. Upon vertical ionization of NCCl center dot center dot center dot NH3 complex, the hole is predicted to migrate from the NH3-end to the ClCN-end of the NCCl center dot center dot center dot NH3 complex in approximately 0.5 fs on the D-0 cationic surface. This hole migration leads to structural rearrangement of the halogen bonded complex, yielding hydrogen bonding interaction stronger than the halogen bonding interaction on the same cationic surface. Other halogen bonded complexes, such as H2NCl:NH3, F3CCl:NH3, and HOOCCl:NH3, exhibit similar charge migration following vertical ionization. On the contrary, FCl:NH3 and HOCl:NH3 complexes do not exhibit any charge migration following vertical ionization to the D-0 cation state, pointing to interesting halogen bond strength-dependent charge migration. (C) 2015 AIP Publishing LLC.

Lattice Boltzmann study of hydrodynamic effects in lamellar ordering process of two-dimensional quenched block copolymers

By incorporating self-consistent field theory with lattice Boltzmann method, a model for polymer melts is proposed. Compared with models based on Ginzburg-Landau free energy, our model does not employ phenomenological free energies to describe systems and can consider the chain topological details of polymers. We use this model to study the effects of hydrodynamic interactions on the dynamics of microphase separation for block copolymers. In the early stage of phase separation, an exponential growth predicted by Cahn-Hilliard treatment is found. Simulation results also show that the effect of hydrodynamic interactions can be neglected in the early stage.

The phase behavior of comblike block copolymer A(m+1)B(m)/homopolymer A mixtures

The phase behaviors of comblike block copolymer A(m+1)B(m)/homopolymer A mixtures are studied by using the random phase approximation method and real-space self-consistent field theory. From the spinodals of macrophase separation and microphase separation, we can find that the number of graft and the length of the homopolymer A have great effects on the phase behavior of the blend. For a given composition of comblike block copolymer, increasing the number of graft does not change the macrophase separation spinodal curve but decreases the microphase separation region. The addition of a small quantity of long-chain homopolymer A increases the microphase separation of comblike block copolymer/homopolymer A mixture.

Defect evolution and hydrodynamic effects in lamellar ordering process of two-dimensional quenched block copolymers

The effects of hydrodynamic interactions on the lamellar ordering process for two-dimensional quenched block copolymers in the presence of extended defects and the topological defect evolutions in lamellar ordering process are numerically investigated by means of a model based on lattice Boltzmann method and self-consistent field theory. By observing the evolution of the average size of domains, it is found that the domain growth is faster with stronger hydrodynamic effects. The morphological patterns formed also appear different. To study the defect evolution, a defect density is defined and is used to explore the defect evolutions in lamellar ordering process. Our simulation results show that the hydrodynamics effects can reduce the density of defects. With our model, the relations between the Flory-Huggins interaction parameter chi, the length of the polymer chains N, and the defect evolutions are studied.

STUDY OF COARSENING OF MICRODOMAINS IN BLOCK COPOLYMERS BY LATTICE BOLTZMANN METHODS

In order to understand the coarsening of microdomains in symmetric diblock copolymers at the late stage, a model for block copolymers is proposed. By incorporating the self consistent field theory with the free energy approach Lattice Boltzmann model, hydrodynamic interactions can be considered. Compared with models based on Ginzburg-Landau free energy, this model does not employ phenomenological free energies to describe systems. The model is verified by comparing the simulation results obtained using this method with those of a dynamical version of the self consistent mean field theory. After that,the growth exponents of the characteristic domain size for symmetric block copolymers at late stage are studied. It is found that the viscosity of the system affects the growth exponents greatly, although the growth exponents are all less than 1/3 Furthermore, the relations between the growth exponent, the interaction parameter and the chain length are studied.

Nucleation in A/B/AB blends: Interplay between microphase assembly and macrophase separation

We study the interplay between microphase assembly and macrophase separation in A/B/AB ternary polymer blends by examining the free energy of localized fluctuation structures (micelles or droplets), with emphasis on the thermodynamic relationship between swollen micelles (microemulsion) and the macrophase-separated state, using self-consistent field theory and an extended capillary model. Upon introducing homopolymer B into a micelle-forming binary polymer blend A/AB, micelles can be swollen by B. A small amount of component B (below the A-rich binodal of macrophase coexistence) will not affect the stability of the swollen micelles. A large excess of homopolymer, B, will induce a microemulsion failure and lead to a macrophase separation.

Comparing the morphology and phase diagram of H-shaped ABC block copolymers and linear ABC block copolymers

By using a combinatorial screening method based on the self-consistent field theory (SCFT) for polymers, we have investigated the morphology of H-shaped ABC block copolymers (A(2)BC(2)) and compared them with those of the linear ABC block copolymers. By changing the ratios of the volume fractions of two A arms and two C arms, one can obtain block copolymers with different architectures ranging from linear block copolymer to H-shaped block copolymer. By systematically varying the volume fractions of block A, B, and C, the triangle phase diagrams of the H-shaped ABC block copolymer with equal interactions among the three species are constructed. In this study, we find four different morphologies ( lamellar phase ( LAM), hexagonal lattice phase ( HEX), core-shell hexagonal lattice phase (CSH), and two interpenetrating tetragonal lattice (TET2)). Furthermore, the order-order transitions driven by architectural change are discussed.

Automatic generation of active coordinates for quantum dynamics calculations: application to the dynamics of benzene photochemistry

A new practical method to generate a subspace of active coordinates for quantum dynamics calculations is presented. These reduced coordinates are obtained as the normal modes of an analytical quadratic representation of the energy difference between excited and ground states within the complete active space self-consistent field method. At the Franck-Condon point, the largest negative eigenvalues of this Hessian correspond to the photoactive modes: those that reduce the energy difference and lead to the conical intersection; eigenvalues close to 0 correspond to bath modes, while modes with large positive eigenvalues are photoinactive vibrations, which increase the energy difference. The efficacy of quantum dynamics run in the subspace of the photoactive modes is illustrated with the photochemistry of benzene, where theoretical simulations are designed to assist optimal control experiments