Selective substitution of alkali cations in mixed alkali glass by solid-state electrochemistry

Electrosubstitution of alkali cations in mixed-alkali glass containing both Na₂O and K₂O for other monovalent metal cations (M⁺=Li⁺, Ag⁺, and Cs⁺) was investigated using a solid-state electrochemical method. The fundamental electrolysis system consists of anode/M⁺-conducting microelectrode/glass/Na-β″-Al₂O₃/cathode, where M⁺ is substituted for the alkali metal ions in the glass under an applied electric field. Li⁺ ions attacked only Na⁺ sites, and Ag⁺ ions replaced Na⁺ sites more readily than K⁺. In contrast, Cs⁺ ions simultaneously substituted for both Na⁺ and K⁺ sites. The substitution behavior appears to depend on the difference in ionic conductivity between K⁺ and Na⁺ and the radius of the dopant. This mechanism was discussed qualitatively.     

Shear viscosity of molten alkali halides from equilibrium and nonequilibrium molecular-dynamics simulations

The shear viscosity of molten NaCl and KCl was calculated through equilibrium (EMD) and nonequilibrium molecular-dynamics (NEMD) simulations in the canonical (N,V,T) ensemble. Two rigid-ion potentials were investigated, namely, the Born-Mayer-Huggins-Tosi-Fumi potential and the Michielsen-Woerlee-Graaf-Ketelaar potential with the parameters proposed by Ladd. The NEMD simulations were performed using the SLLOD equations of motion [D. J. Evans and G. P. Morriss, Phys. Rev. A 30, 1528 (1984)] with a Gaussian isokinetic thermostat and the results are compared with those obtained from Green-Kubo EMD (N,V,T) simulations and experimental shear viscosity data. The NEMD zero strain rate shear viscosity, eta(0), was obtained by fitting a simplified Carreau-type equation and by application of mode-coupling theory, i.e., a eta-gamma(1/2) linear relationship. The values obtained from the first method are found to be significantly lower than those predicted by the second. The agreement between the EMD and NEMD results with experimental data is satisfactory for the two potentials investigated. The ion-ion radial distribution functions obtained with the two rigid-ion potentials for both molten salts are discussed in terms of the differences between the two models.     

Dynamic light scattering in mixed alkali metaphosphate glass forming liquids

We report the first ever photon correlation spectroscopy performed on single alkali and mixed alkali metaphosphate glasses at refractory temperatures above the glass transition. We find not only a significant decrease in the glass transition temperature but also a decrease in fragility for the mixed alkali composition as compared with the single akali glasses. We argue that structural relaxation in these polymeric oxide glasses is largely controlled by the cross linking cations and that the changes in fragility that we observed are a reflection of changes in the cooperativity of structural relaxation wrought by the substantial decrease in the ion mobility that accompanies the mixing of alkali ions.     

Vibrational spectra and structure of mixed alkali borotungstate glasses: Evidence of mixed alkali effect

Mixed alkali borotungstate glasses with xLi₂O–(30 − x)Na₂O–10WO₃–60B₂O₃ (0 ≤ x ≤ 30) composition were prepared by melt quench technique. FT-IR and Raman spectroscopic studies were employed to investigate the structure of all the prepared glasses. Acting as complementary techniques, both IR and Raman measurements revealed that the network structure of the present glasses mainly based on BO₃ and BO₄ units placed in different structural groups. Raman spectra confirm the IR results regarding the presence of tungsten ions mainly as WO₆ groups. In the present work, the mixed alkali effect (MAE) has been investigated in the above glass system using FTIR and Raman studies.     

Structural properties of a calcium aluminosilicate glass from molecular-dynamics simulations: a finite size effects study

We study a calcium aluminosilicate glass of composition (SiO(2))(0.67)-(Al(2)O(3))(0.12)-(CaO)(0.21) by means of molecular-dynamics simulations, using a potential made of two-body and three-body interactions. In order to prepare small samples that can subsequently be studied by first principles, the finite size effects on the liquid dynamics and on the glass structural properties are investigated. We find that finite size effects affect the Si-O-Si and Si-O-Al angular distributions, the first peaks of the Si-O, Al-O, and Ca-O pair correlation functions, the Ca coordination, and the oxygen atoms' environment in the smallest system (100 atoms). We give evidence that these finite size effects can be directly attributed to the use of three-body interactions.     

The mechanism of vesicle fusion as revealed by molecular dynamics simulations

We describe molecular dynamics simulations elucidating the molecular details of the process of fusion for small lipid vesicles. The simulations are based on a coarse grained (CG) lipid model that accurately represents the lamellar state of a variety of phospholipids and enables us to observe intermediate stages during fusion at near atomic detail. Simulations were conducted on a variety of systems containing common phospholipids such as phosphatidylcholine (PC), phosphatidylethanolamine (PE), lysoPC, and mixtures of the above. The fusion intermediates found are in general agreement with the stalk-pore mechanism. Transient pores sometimes form adjacent to the stalk, however, resulting in the mixing of lipids from the outer and inner monolayers. The speed of stalk formation and the opening of the fusion pore can be modulated by altering the lipid composition in qualitative agreement with experimental observations.     

Hydrogen bonding in supercritical tert-butanol assessed by vibrational spectroscopies and molecular-dynamics simulations

We have investigated the state of aggregation in supercritical tert-butanol (T = 523 K,0.05 < rho < 0.4 g cm(-3)) by means of vibrational spectroscopies (infrared and Raman) and molecular-dynamics (MD) simulations. A quantitative band shape analysis of the spectra associated with the OH stretching mode of tert-butanol has been done using activities computed by ab initio calculations on small clusters. This allows us to determine the degree of hydrogen bonding and populations of oligomers. These latter quantities have been derived from MD simulations and very consistent results are found with experiments. These results show that hydrogen bond still exist in supercritical tert-butanol and that the fluid mainly consists of oligomers smaller than tetramers.     

Calculation of the surface tension of oxygen using molecular-dynamics simulations

The surface tension of oxygen at the liquid-vapor interface is calculated for the temperature range of 60-90 K using molecular-dynamics simulations and is shown to be within 1.0% error of experimental values for most of the temperatures studied. The potential used here is the same as in our previous study on liquid oxygen alone [S. D. Bembenek and B. M. Rice, J. Chem. Phys. 113, 2354 (2000)] and was optimized with an innovative statistical-mechanical method. The potential does not use a long-range cutoff nor a tail correction, which are usually considered necessary to obtain accurate values for the surface tension. We reason that the accuracy in surface tension is directly related to our parametrization method for the potential.     

Structural change in Li and Na aluminophosphate glasses: evidence of a "structural mixed alkali effect"

The short- and long-range structure of a series of single and mixed aluminophosphate glasses with the general composition [xNa(2)O (46 - x)Li(2)O], [yAl(2)O(3) (54 - y)P(2)O(5)] is analyzed using (31)P and (27)Al magic-angle spinning (MAS) NMR as well as small-angle X-ray scattering. These series of glasses allow analyzing both the effect of alumina incorporation in these glasses, for small alumina content (y = 0, 4, 8), and the structural changes associated with the so-called mixed alkali effect (x = 0, 11.5, 23, 34.5, 46). Our results indicate that aluminum is mainly octahedrally coordinated in these glasses and that there is most likely some segregation of the Al(OP)(6) species. In the pure phosphate glasses, we observe a "classical" continuous variation of the structural properties with the relative alkali content, but in the aluminophosphate, both local and long-range structural results reveal for the first time some nonlinear change as a function of the relative alkali content.     

Langevin stabilization of molecular-dynamics simulations of polymers by means of quasisymplectic algorithms

Algorithms for the numerical integration of Langevin equations are compared in detail from the point of view of their accuracy, numerical efficiency, and stability to assess them as potential candidates for molecular-dynamics simulations of polymeric systems. Some algorithms are symplectic in the deterministic frictionless limit and prove to stabilize long time-step integrators. They are tested against other popular algorithms. The optimal algorithm depends on the main goal: accuracy or efficiency. The former depends on the observable of interest. A recently developed quasisymplectic algorithm with great accuracy in the position evaluation exhibits better overall accuracy and stability than the other ones. On the other hand, the well-known BrunGer-Brooks-Karplus [Chem. Phys. Lett. 105, 495 (1982)] algorithm is found to be faster with limited accuracy loss but less stable. It is also found that using higher-order algorithms does not necessarily improve the accuracy. Moreover, they usually require more force evaluations per single step, thus leading to poorer performances.     

Electrical conductivity studies in sulphate glasses and the mixed alkali effect

Electrical conductivities of alkali sulphate-zinc sulphate glasses have been measured. The variation of conductivity with compositon confirms the presence of the mixed alkali effect. The origin of mixed alkali effect has been explained on the basis of structural considerations reported earlier by us. Communication No. 160 from the Solid State and Structural Chemistry Unit.     

Coarse-grained molecular-dynamics simulations of the self-assembly of pentablock copolymers into micelles

Multiblock polymers in aqueous solution, where one or several blocks are hydrophobic, exhibit a rich variety of phases and states of aggregation. In this paper, we investigate a pentablock system ABCBA, where the B block is always hydrophilic and the A and C blocks have varying degrees of hydrophobicity depending on external conditions. We report coarse-grained molecular-dynamics simulations where the solvent is included explicitly and monomers interact via a 6-9 Lennard Jones potential function. The hydrophobic interaction is modeled by tuning the parameter controlling the strength of the interaction between the hydrophobic monomers and the solvent. We investigate the structure and morphology of the micelles for two concrete situations representing changes in temperature and the pH level. The simulated system is directly relevant to a recently synthesized pentablock system consisting of a triblock Pluronic with an added pH-sensitive end group [B. C. Anderson et al., Macromolecules 36, 1670 (2003)].     

Polymer dynamics under nanoscopic constraints: the "corset effect" as revealed by NMR relaxometry and diffusometry.

A spinodal demixing technique was employed for the preparation of linear poly(ethylene oxide) (PEO) confined in nanoscopic strands, which in turn are embedded in a quasi-solid methacrylate matrix impenetrable to PEO. Both the molecular weight of the PEO and the mean diameter of the strands are variable to a certain degree. Chain dynamics of the PEO in the molten state were examined with the aid of field-gradient NMR diffusometry and field-cycling NMR relaxometry. The dominating mechanism for translational displacements in the nanoscopic strands is shown to be reptation. A formalism for the evaluation of NMR diffusometry is presented, which permits the estimation of the mean PEO strand diameter. Samples of different composition revealed diameters in the range 9-58 nm, in reasonable agreement with electron micrographs. The time scale of the diffusion measurements was 10-300 ms. On the much shorter time scale of field-cycling NMR relaxometry, 10(-9)-10(-4)s, a frequency dispersion of the spin-lattice relaxation time characteristic for reptation clearly showed up in all samples. An effective tube diameter of only 0.6 nm was found even when the strand diameter was larger than the radius of gyration of the PEO chain random coils. The finding that the tube diameter effective on the short time scale of field-cycling NMR relaxometry is much smaller than the diameter of the confining structure is termed the "corset effect", and is traced back to the lack of local free-volume fluctuation capacity under nanoscale confinements. The order of magnitude of the 'pore' diameter, at which the cross-over from confined to bulk chain dynamics is expected, is estimated.     

Speed up of dynamic observables in coarse-grained molecular-dynamics simulations of unentangled polymers

Coarse-grained models that preserve atomistic detail display faster dynamics than atomistic systems alone. We show that this " indirect speed up" is robust: coarse-grained dynamic observables computed with time scaled by a constant factor are in excellent agreement with their underlying atomistic counterparts. Borrowing from accelerated dynamics methods used in the field of rare events, we predict the scaling factor within 7%, based on reduced intermolecular attraction yielding faster neighbor cage escapes.     

Polarizability response in polar solvents: molecular-dynamics simulations of acetonitrile and chloroform

The relaxation of the many-body polarizability in liquid acetonitrile and chloroform at room temperature was studied by molecular-dynamics simulations. The collective polarizability induced by intermolecular interactions was included using first- and all-orders dipole-induced-dipole models and calculated considering both molecule-centered and distributed site polarizabilities. The anisotropic response was analyzed using a separation scheme that allows a decomposition of the total response in terms of orientational and collision-induced effects. We found the method effective in approximately separating the contributions of these relaxation mechanisms, although the orientational-collision-induced interference makes a non-negligible contribution to the total response. In both liquids the main contribution to the anisotropic response is due to orientational dynamics, but intermolecular collision-induced (or translational) effects are important, especially at short times. We found that higher-order interaction-induced effects were essentially negligible for both liquids. Larger differences were found between the center-center and site-site models, with the latter showing faster polarizability relaxation and better agreement with experiment. Isotropic and anisotropic spectra were computed from the corresponding time correlation functions. The lowest-frequency contributions are largely suppressed in the isotropic spectra and their overall shape is similar to the purely collision-induced contribution to the anisotropic spectra, but with an amplitude which is smaller by a factor of approximately 5 in acetonitrile and approximately 3 in chloroform.     

The mixed alkali effect in ionically conducting glasses revisited: a study by molecular dynamics simulation

When more than two kinds of mobile ions are mixed in ionic conducting glasses and crystals, there is a non-linear decrease of the transport coefficients of either type of ion. This phenomenon is known as the mixed mobile ion effect or Mixed Alkali Effect (MAE), and remains an unsolved problem. We use molecular dynamics simulation to study the complex ion dynamics in ionically conducting glasses including the MAE. In the mixed alkali lithium-potassium silicate glasses and related systems, a distinct part of the van Hove functions reveals that jumps from one kind of site to another are suppressed. Although, consensus for the existence of preferential jump paths for each kind of mobile ions seems to have been reached amongst researchers, the role of network formers and the number of unoccupied ion sites remain controversial in explaining the MAE. In principle, these factors when incorporated into a theory can generate the MAE, but in reality they are not essential for a viable explanation of the ion dynamics and the MAE. Instead, dynamical heterogeneity and "cooperativity blockage" originating from ion-ion interaction and correlation are fundamental for the observed ion dynamics and the MAE. Suppression of long range motion with increased back-correlated motions is shown to be a cause of the large decrease of the diffusivity especially in dilute foreign alkali regions. Support for our conclusion also comes from the fact that these features of ion dynamics are common to other ionic conductors, which have no glassy networks, and yet they all exhibit the MAE.     

Atomistic molecular-dynamics simulations of the size and shape of polyethylene in hexane at infinite dilution

Parameters characteristic of size and shape of single polyethylene chains consisting of 15-60 monomer units dissolved in hexane are calculated by use of molecular-dynamics simulations based on a fully atomistic representation of the system. Results are compared with corresponding calculations in vacuum as well as Monte Carlo simulations of coarse-grained chains. The major concern of the study is a careful check of actual limits and possibilities of atomistic simulations of global properties of polymers. As expected such simulations are still restricted to rather small chain lengths but are already large enough to obey the characteristics of polymer coils.