Monte Carlo simulation of miscibility of polymer blends with repulsive interactions: effect of chain structure

The effects of the chain structure and the intramolecular interaction energy of an A/B copolymer on the miscibility of the binary blends of the copolymer and homopolymer C have been studied by means of a Monte Carlo simulation. In the system, the interactions between segments A, B and C are more repulsive than those between themselves. In order to study the effect of the chain structure of the A/B copolymer on the miscibility, the alternating, random and block copolymers were introduced in the simulations, respectively. The simulation results show that the miscibility of the binary blends strongly depends on the intramolecular interaction energy ((ε) over bar (AB)) between segments A and B within the A/B copolymers. The higher the repulsive interaction energy, the more miscible the A/B copolymer and homopolymer C are. For the diblock copolymer/homopolymer blends, they tend to form micro phase domains. However, the phase domains become so small that the blend can be considered as a homogeneous phase for the alternating copolymer/ homopolymer blends. Furthermore, the investigation of the average end-to-end distance ((h) over bar) in different systems indicates that the copolymer chains tend to coil with the decrease Of (ε) over bar (AB) whereas the (h) over bar of the homopolymer chains depends on the chain structure of the copolymers.

Site symmetry dependence of repulsive interactions between chemisorbed oxygen atoms on Pt{100}-(1x1)

Ab initio total energy calculations using density functional theory with the generalized gradient approximation have been performed for the chemisorption of oxygen atoms on a Pt{100}-(1 x 1) slab. Binding energies for the adsorption of oxygen on different high-symmetry sites are presented. The bridge site is the most stable at a coverage of 0.5 ML, followed by the fourfold hollow site. The atop site is the least stable. This finding is rationalized by analyzing the ''local structures'' formed upon oxygen chemisorption. The binding energies and heats of adsorption at different oxygen coverages show that pairwise repulsive interactions are considerably stronger between oxygen atoms occupying fourfold sites than those occupying bridge sites. Analysis of the partial charge densities associated with Bloch states demonstrates that the O-Pt bond is considerably more localized at the bridge site. These effects cause a sharp drop in the heats of adsorption for oxygen on hollow sites when the coverage is increased from 0.25 to 0.5 ML. Mixing between oxygen p orbitals and Pt d orbitals can be observed over the whole metal d-band energy range.

Competition between local potentials and attractive particle-particle interactions in superlattices

Naturally occuring or man-made systems displaying periodic spatial modulations of their properties on a nanoscale constitute superlattices. Such modulated structures are important both as prototypes of simple nanotechnological devices and as particular examples of emerging spatial inhomogeneity in interacting many-electron systems. Here we investigate the effect different types of modulation of the system parameters have on the ground-state energy and the charge-density distribution of the system. The superlattices are described by the inhomogeneous attractive Hubbard model, and the calculations are performed by density-functional and density-matrix renormalization group techniques. We find that modulations in local electric potentials are much more effective in shaping the system's properties than modulations in the attractive on-site interaction. This is the same conclusion we previously [M.F. Silva, N.A. Lima, A.L. Malvezzi, K. Capelle, Phys. Rev. B 71 (2005) 125130.] obtained for repulsive interactions, suggesting that it is not an artifact of a specific state, but a general property of modulated structures. (c) 2007 Elsevier Ltd. All rights reserved.

Protein interactions studied by SAXS:effect of ionic strength and protein concentration for BSA in aqueous solutions

We have studied a series of samples of bovine serum albumin (BSA) solutions with protein concentration, c, ranging from 2 to 500 mg/mL and ionic strength, I, from 0 to 2 M by small-angle X-ray scattering (SAXS). The scattering intensity distribution was compared to simulations using an oblate ellipsoid form factor with radii of 17 x 42 x 42 A, combined with either a screened Coulomb, repulsive structure factor, S-SC(q), or an attractive square-well structure factor, S-SW(q). At pH = 7, BSA is negatively charged. At low ionic strength, I <0.3 M, the total interaction exhibits a decrease of the repulsive interaction when compared to the salt-free solution, as the net surface charge is screened, and the data can be fitted by assuming an ellipsoid form factor and screened Coulomb interaction. At moderate ionic strength (0.3-0.5 M), the interaction is rather weak, and a hard-sphere structure factor has been used to simulate the data with a higher volume fraction. Upon further increase of the ionic strength (I >= 1.0 M), the overall interaction potential was dominated by an additional attractive potential, and the data could be successfully fitted by an ellipsoid form factor and a square-well potential model. The fit parameters, well depth and well width, indicate that the attractive potential caused by a high salt concentration is weak and long-ranged. Although the long-range, attractive potential dominated the protein interaction, no gelation or precipitation was observed in any of the samples. This is explained by the increase of a short-range, repulsive interaction between protein molecules by forming a hydration layer with increasing salt concentration. The competition between long-range, attractive and short-range, repulsive interactions accounted for the stability of concentrated BSA solution at high ionic strength.

Interaction of water with the fluorine-covered anatase TiO2(001) surface

The interaction of water with the fluorine-covered (001) surface of anatase titanium dioxide (TiO2) has been studied within the framework of density functional theory (DFT). Our results show that water dissociation is unfavorable due to repulsive interactions between surface fluorine and oxygen. We also found that the reaction of hydrofluoric acid with a surface hydroxyl group to form a surface Ti–F bond is exothermic, while the removal of fluorine from the surface needs additional energy of about half an eV. Therefore, water molecules are predicted to remain intact at the interface with the F-terminated anatase (001).

The structures of the ammonium halides

A comparison with the alkali halides suggests that all the ammonium halides should occur in the NaCl centre-of-mass structure. Experimentally, at room temperature and atmospheric pressure, only NH4I crystallizes in this structure, while NH4F is found in the ZnO structure, and NH4C1 and NH4Br occur in the CsCl structure. We show that a distributed charge on the NH4+ ion can explain these structures. Taking charges of + 0.2e on each of the five atoms in NH4+, as suggested by other studies, we have recomputed the Madelung energy in the cases of interest. A full ionic theory including electrostatic, van der Waals and repulsive interactions then explains the centre-of-mass structures of all the four ammonium halides. The thermal and pressure transitions are also explained reasonably well. The calculated phase diagram of NH4F compares well with experiment. Barring the poorly understood NH4F(II) phase, which is beyond the scope of this work, the other features are in qualitative agreement. In particular, the theory correctly predicts a pressure transition at room temperature from the ZnO structure directly to the CsCl structure without an intermediate NaCl phase. A feature of our approach is that we do not need to invoke hydrogen bonding in NH4F.

The structures of the ammonium halides

A comparison with the alkali halides suggests that all the ammonium halides should occur in the NaCl centre-of-mass structure. Experimentally, at room temperature and atmospheric pressure, only NH,1 crystallizes in this structure, while NH,F is found in the ZnO structure, and NH&I and NH,Br occur in the CsCl structure. We show that a distributed charge on the NH,+ ion can explain these structures. Taking charges of + 0.2e on each of the five atoms in NH,+, as suggested by other studies, we have recomputed the Madelung energy in the cases of interest. A full ionic theory including electrostatic, van der Waals and repulsive interactions then explains the centre-of-mass structures of all the four ammonium halides. The thermal and pressure transitions are also explained reasonably well. The calculated phase diagram of NH,F compares well with experiment. Barring the poorly understood NH,F(II) phase, which is beyond the scope of this work, the other features are in qualitative agreement. In particular, the theory correctly predicts a pressure transition at room temperature from the ZnO structure directly to the CsCl structure without an intermediate NaCl phase. A feature of our approach is that we do not need to invoke hydrogen bonding in NH,F.

Enhanced charge and spin currents in the one-dimensional disordered mesoscopic Hubbard ring

We consider a one-dimensional mesoscopic Hubbard ring with and without disorder and compute charge and spin stiffness as a measure of the permanent currents. For finite disorder we identify critical disorder strength beyond which the charge currents in a system with repulsive interactions are larger than those for a free system. The spin currents in the disordered repulsive Hubbard model are enhanced only for small U, where the magnetic state of the system corresponds to a charge-density wave pinned to the impurities. For large U, the state of the system corresponds to localized isolated spins and the spin currents are found to be suppressed. For the attractive Hubbard model we find that the charge currents are always suppressed compared to the free system at all length scales.

Persistent currents in a one-dimensional ring for a disordered Hubbard model

We consider a one-dimensional Hubbard model in the presence of disorder. We compute the charge stiffness for a mesoscopic ring as a function of the size L, which is a measure of the persistent currents. We find that for finite disorder the persistent currents of the system with repulsive interactions are larger than those of the system with attractive interactions. This counterintuitive result is due to the fact that local-density fluctuations are reduced in the presence of repulsive interactions.

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 translational and rotational mobility and orientational preference of ethane in one-dimensional channels

Structural and dynamical properties of ethane in one-dimensional channels of AlPO4-5 and carbon nanotube have been investigated at dilute concentration with the help of molecular dynamics simulation. Density distributions and orientational structure of ethane have been analyzed. Repulsive interactions seem to play an important role when ethane is located in the narrow part of the AlPO4-5 channel. In AlPO4-5, parallel orientation is predominant over perpendicular orientation except when ethane is located in the broader part of the channel. Unlike in the case of single-file diffusion, our results in carbon nanotube show that at dilute concentrations the mean squared displacement, mu(2)(t) approximate to t(alpha), alpha = 1.8. The autocorrelation function for the z-component of angular velocity of ethane in space-fixed frame of reference shows a pronounced negative correlation. This is attributed to the restriction in the movement of ethane along the x- and y- directions. It is seen that the ratio of reorientational correlation times does not follow the Debye model for confined ethane but it is closer to the predictions of the Debye model for bulk ethane.

An improved method for calculating activities from distribution equilibria

The method of Gibbs-Duhem integration suggested by Speiser et al. has been modified to derive activities from distribution equilibria. It is shown that, in general, the activities of components in melts with a common anion can be calculated, without using their standard Gibbs energies of formation, from eqUilibrium ratios and the knowledge of activities in the metal phase. Moreover, if systems are so chosen that the concentration of one element in the metal phase lies in the Henry's law region (less than 1 %), information on activities in the metal phase is not required. Conversely, activities of elements in an alloy can be readily calculated from equilibrium distribution ratios alone, if the salt phase in equilibrium contains very small amounts of one element. Application of the method is illustrated using distribution ratios from the literature on AgCI-CuCI, AgBr-CuBr, and CuDo.5 -PbD systems. The results indicate that covalent bonding and van der Waals repulsive interactions in certain types of fused salt melts can significantly affect the thermodynamic properties of mixing.

Nanoparticle induced miscibility in LCST polymer blends: critically assessing the enthalpic and entropic effects

The use of copolymer and polymer blends widened the possibility of creating materials with multilayered architectures. Hierarchical polymer systems with a wide array of micro and nanostructures are generated by thermally induced phase separation (TIPS) in partially miscible polymer blends. Various parameters like the interaction between the polymers, concentration, solvent/non-solvent ratio, and quenching temperature have to be optimized to obtain these micro/nanophase structures. Alternatively, the addition of nanoparticles is another strategy to design materials with desired hetero-phase structures. The dynamics of the polymer nanocomposite depends on the statistical ordering of polymers around the nanoparticle, which is dependent on the shape of the nanoparticle. The entropic loss due to deformation of polymer chains, like the repulsive interactions due to coiling and the attractive interactions in the case of swelling has been highlighted in this perspective article. The dissipative particle dynamics has been discussed and is correlated with the molecular dynamics simulation in the case of polymer blends. The Cahn Hillard Cook model on variedly shaped immobile fillers has shown difference in the propagation of the composition wave. The nanoparticle shape has a contributing effect on the polymer particle interaction, which can change the miscibility window in the case of these phase separating polymer blends. Quantitative information on the effect of spherical particles on the demixing temperature is well established and further modified to explain the percolation of rod shaped particles in the polymer blends. These models correlate well with the experimental observations in context to the dynamics induced by the nanoparticle in the demixing behavior of the polymer blend. The miscibility of the LCST polymer blend depends on the enthalpic factors like the specific interaction between the components, and the solubility product and the entropic losses occurring due to the formation of any favorable interactions. Hence, it is essential to assess the entropic and enthalpic interactions induced by the nanoparticles independently. The addition of nanoparticles creates heterogeneity in the polymer phase it is localized. This can be observed as an alteration in the relaxation behavior of the polymer. This changes the demixing behavior and the interaction parameter between the polymers. The compositional changes induced due to the incorporation of nanoparticles are also attributed as a reason for the altered demixing temperature. The particle shape anisotropy causes a direction dependent depletion, which changes the phase behavior of the blend. The polymer-grafted nanoparticles with varying grafting density show tremendous variation in the miscibility of the blend. The stretching of the polymer chains grafted on the nanoparticles causes an entropy penalty in the polymer blend. A comparative study on the different shaped particles is not available up to date for understanding these aspects. Hence, we have juxtaposed the various computational studies on nanoparticle dynamics, the shape effect of NPs on homopolymers and also the cases of various polymer blends without nanoparticles to sketch a complete picture on the effect of various particles on the miscibility of LCST blends.

Quantitative relationships between molecular structures, environmental temperatures and octanol-air partition coefficients of PCDD/Fs

A concise quantitative model that incorporates information on both environmental temperature M and molecular structures, for logarithm of octanol-air partition coefficient (K-OA) to base 10 (logK(OA)) of PCDDs, was developed. Partial least squares (PLS) analysis together with 14 quantum chemical descriptors were used to develop the quantitative relationships between structures, environmental temperatures and properties (QRSETP) model. It has been validated that the obtained QRSETP model can be used to predict logK(OA) of other PCDDs. Molecular size, environmental temperature (T), q(+) (the most positive net atomic charge on hydrogen or chlorine atoms in PCDD molecules) and E-LUMO (the energy of the lowest unoccupied molecular orbital) are main factors governing logK(OA) of PCDD/Fs under study. The intermolecular dispersive interactions and thus the size of the molecules play a leading role in governing logK(OA). The more chlorines in PCDD molecules, the greater the logK(OA) values. Increasing E-LUMO values of the molecules leads to decreasing logK(OA) values, implying possible intermolecular interactions between the molecules under study and octanol molecules. Greater q(+) values results in greater intermolecular electrostatic repulsive interactions between PCDD and octanol molecules and smaller logK(OA) values. (C) 2002 Elsevier Science B.V. All rights reserved.