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.

Calculating the equation of state parameters and predicting the spinodal curve of isotactic polypropylene/poly(ethylene-co-octene) blend by molecular dynamics simulations combined with Sanchez-Lacombe lattice fluid theory

Molecular dynamics simulations are adopted to calculate the equation of state characteristic parameters P*, rho*, and T* of isotactic polypropylene (iPP) and poly(ethylene-co-octene) (PEOC), which can be further used in the Sanchez-Lacombe lattice fluid theory (SLLFT) to describe the respective physical properties. The calculated T* is a function of the temperature, which was also found in the literature. To solve this problem, we propose a Boltzmann fitting of the data and obtain T* at the high-temperature limit. With these characteristic parameters, the pressure-volume-temperature (PVT) data of iPP and PEOC are predicted by the SLLFT equation of state. To justify the correctness of our results, we also obtain the PVT data for iPP and PEOC by experiments. Good agreement is found between the two sets of data. By integrating the Euler-Lagrange equation and the Cahn-Hilliard relation, we predict the density profiles and the surface tensions for iPP and PEOC, respectively. Furthermore, a recursive method is proposed to obtain the characteristic interaction energy parameter between iPP and PEOC. This method, which does not require fitting to the experimental phase equilibrium data, suggests an alternative way to predict the phase diagrams that are not easily obtained in experiments.

Study on the two class glass transition and ion transport of a comb polymer electrolyte

The thermal behaviour and ion-transport properties of a comb polymer electrolyte CP350/LiSCN based on methyl vinyl ether/maleic anhydride copolymer with oligo-oxyethylene side chains were studied by means of DSC and ac impedance method. The two glass transition temperatures which can be attributed to side chains and main chains respectively were found to increase with increasing salt concentration. Conductivities which displayed non-Arrhenius behaviour were analyzed by using Vogel-Tammann-Fulcher equation and interpreted on the basis of the configurational entropy model derived by Gibbs and coworkers. The optimum ionic conductivity at 25 degrees C achieved was 2.19x10(-5)S/cm.

Symmetry-breaking and chaos in electron transport in semiconductor superlattices

We study the motion of electrons in a single miniband of a semiconductor superlattice driven by THz electric field polarized along the growth direction. We work in the semiclassical balance-equation model, including different elastic and inelastic scattering rates, and incorporating the self-consistent electric field generated by electron motion. We explore regions of complex dynamics, which can include chaotic behaviour and symmetry-breaking. We estimate the magnitudes of dc current and dc voltage that spontaneously appear in regions of broken-symmetry for parameters characteristic of modern semiconductor superlattices. This work complements PRL 80(1998)2669 [ cond-mat/9709026 ].

Water transport in resin-modified glass-ionomer dental cement

Water uptake and water loss have been studied in a commercial resin-modified glass-ionomer cement, Fuji II LC, under a variety of conditions. Uptake was generally non-Fickian, but affected by temperature. At room temperature, the equilibrium water uptake values varied from 2.47 to 2.78% whereas at low temperature (12 degrees C), it varied from 0.85 to 1.18%. Cure time affected uptake values significantly. Water uptake was much lower than in conventional glass-ionomer restorative cements exposed to water vapor. Loss of water under desiccating conditions was found to be Fickian for the first 5 h loss at both 22 and 12 degrees C. Diffusion coefficients were between 0.45 and 0.76 x 10( -7) cm(2)/s, with low temperature diffusion coefficients slightly greater than those at room temperature. Plotting water loss as percentage versus s(-(1/2)) allowed activation energies to be determined from the Arrhenius equation and these were found to be 65.6, 79.8, and 7.7 kJ/mol respectively for 30, 20, and 10 s cure times. The overall conclusion is that the main advantage of incorporating HEMA into resin-modified-glass-ionomers is to alter water loss behavior. Rate of water loss and total amount lost are both reduced. Hence, resin-modified glass-ionomers are less sensitive to water loss than conventional glass-ionomers.

Solving the advection-dispersion equation using the discrete puff particle method

A flexible, mass-conservative numerical technique for solving the advection-dispersion equation for miscible contaminant transport is presented. The method combines features of puff transport models from air pollution studies with features from the random walk particle method used in water resources studies, providing a deterministic time-marching algorithm which is independent of the grid Peclet number and scales from one to higher dimensions simply. The concentration field is discretised into a number of particles, each of which is treated as a point release which advects and disperses over the time interval. The dispersed puff is itself discretised into a spatial distribution of particles whose masses can be pre-calculated. Concentration within the simulation domain is then calculated from the mass distribution as an average over some small volume. Comparison with analytical solutions for a one-dimensional fixed-duration concentration pulse and for two-dimensional transport in an axisymmetric flow field indicate that the algorithm performs well. For a given level of accuracy the new method has lower computation times than the random walk particle method.

A Van der Pol-Mathieu equation for the dynamics of dust grain charge in dusty plasmas

The chaotic profile of dust grain dynamics associated with dust-acoustic oscillations in a dusty plasma is considered. The collective behaviour of the dust plasma component is described via a multi-fluid model, comprising Boltzmann distributed electrons and ions, as well as an equation of continuity possessing a source term for the dust grains, the dust momentum and Poisson's equations. A Van der Pol–Mathieu-type nonlinear ordinary differential equation for the dust grain density dynamics is derived. The dynamical system is cast into an autonomous form by employing an averaging method. Critical stability boundaries for a particular trivial solution of the governing equation with varying parameters are specified. The equation is analysed to determine the resonance region, and finally numerically solved by using a fourth-order Runge–Kutta method. The presence of chaotic limit cycles is pointed out.

Transport properties of protic ionic liquids, pure and in aqueous solutions: Effects of the anion and cation structure

Ionic conductivities of twelve protic ionic liquids (PILs) and their mixtures with water over the whole composition range are reported at 298.15 K and atmospheric pressure. The selected PILs are the pyrrolidinium-based PILs containing nitrate, acetate or formate anions; the formate-based PILs containing diisopropylethylammonium, amilaminium, quinolinium, lutidinium or collidinium cations; and the pyrrolidinium alkylcarboxylates, [Pyrr][CnH2n+1COO] with n = 5–8. This study was performed in order to investigate the influence of molecular structures of the ions on the ionic conductivities in aqueous solutions. The ionic conductivities of the aqueous solutions are 2–30 times higher than the conductivities of pure PILs. The maximum in conductivity varies from ww=0.41???to???0.74 and is related to the nature of cations and anions. The molar conductance and the molar conductance at infinite dilution for (PIL + water) solutions are then determined. Self-diffusion coefficients of the twelve protic ionic liquids in water at infinite dilution and at 298.15 K are calculated by using the Nernst–Haskell, the original and the modified Wilke–Chang equations. These calculations show that similar values are obtained using the modified Wilke–Chang and the Nernst–Haskell equations. Finally, the effective hydrodynamic (or Stokes) radius of the PILs was determined by using the Stokes–Einstein equation. A linear relationship was established in order to predict this radius as a function of the anion alkyl chain length in the case of the pyrrolidinium alkylcarboxylates PILs.

Transport Properties Investigation of Aqueous Protic Ionic Liquid Solutions through Conductivity, Viscosity, and NMR Self-Diffusion Measurements

We present a study on the transport properties through conductivity (s), viscosity (?), and self-diffusion coefficient (D) measurements of two pure protic ionic liquids—pyrrolidinium hydrogen sulfate, [Pyrr][HSO4], and pyrrolidinium trifluoroacetate, [Pyrr][CF3COO]—and their mixtures with water over the whole composition range at 298.15 K and atmospheric pressure. Based on these experimental results, transport mobilities of ions have been then investigated in each case through the Stokes–Einstein equation. From this, the proton conduction in these PILs follows a combination of Grotthuss and vehicle-type mechanisms, which depends also on the water composition in solution. In each case, the displacement of the NMR peak attributed to the labile proton on the pyrrolidinium cation with the PILs concentration in aqueous solution indicates that this proton is located between the cation and the anion for a water weight fraction lower than 8%. In other words, for such compositions, it appears that this labile proton is not solvated by water molecules. However, for higher water content, the labile protons are in solution as H3O+. This water weight fraction appears to be the solvation limit of the H+ ions by water molecules in these two PILs solutions. However, [Pyrr][HSO4] and [Pyrr][CF3COO] PILs present opposed comportment in aqueous solution. In the case of [Pyrr][CF3COO], ?, s, D, and the attractive potential, Epot, between ions indicate clearly that the diffusion of each ion is similar. In other words, these ions are tightly bound together as ion pairs, reflecting in fact the importance of the hydrophobicity of the trifluoroacetate anion, whereas, in the case of the [Pyrr][HSO4], the strong H-bond between the HSO4– anion and water promotes a drastic change in the viscosity of the aqueous solution, as well as on the conductivity which is up to 187 mS·cm–1 for water weight fraction close to 60% at 298 K.

Modelling solute and microorganism transport through a peri-alpine gravel aquifer by considering the aquifer as a dual porosity fractured medium.

Approach:
In-situ passive gradient comparative artificial tracer testing, undertaken using solutes (Uranine and Iodide), Bacteria (E.coli and P.putida) and bacteriophage (H40/1), permitted comparison of the mobility of different sized microorganisms relative to solutes in the sand and gravel aquifer underlying Dornach, Germany.
Tracer breakthrough curves reveal that even though uranine initially arrived at observation wells at the same time as microbiological tracers, maximum relative concentrations were sometimes less than those of microbiological tracers, while solute breakthrough curves proved more disperse.
Monitoring uranine breakthrough with depth suggested tracers arrived in observation wells in discrete 0.5m-1m thick intervals, over the aquifer’s 12m saturated thickness. Nearby exposures of aquifer material suggested that the aquifer consisted of sandy gravels enveloping sequences of open framework (OW) gravel up to 1m thick. Detailed examination of OW units revealed that they contained lenses of silty sand up to 1m long x 30cm thick., while granulometric data suggested that the gravel was two to three orders of magnitude more permeable than the enveloping sandy gravel.
Solute and microorganism tracer responses could not be simulated using conventional advective-dispersive equation solutions employing the same velocity and dispersion terms. By contrast solute tracer responses, modelled using a dual porosity approach for fractured media (DP-1D) corresponded well to observed field data. Simulating microorganism responses using the same transport terms, but no dual porosity term, generated good model fits and explained the higher relative concentration of the bacteria, compared to the non-reactive solute, even with first order removal to account for lower RR. Geologically, model results indicate that the silty units within open framework gravels are accessible to solute tracers, but not to microorganisms.
Importance:
Results highlight the benefits of geological observations developing appropriate conceptual models of solute and micro organism transport and in developing suitable numerical approaches to quantifying microorganism mobility at scales appropriate for the development of groundwater supply (wellhead) protection zones.

Comparative study on transport properties for LiFAP and LiPF6 in alkyl-carbonates as electrolytes through conductivity, viscosity and NMR self-diffusion measurements

We present in this work a comparative study on density and transport properties, such as the conductivity (sigma), viscosity (eta) and self-diffusion coefficients (D), for electrolytes based on the lithium hexafluorophosphate, LiPF6; or on the lithium tris(pentafluoroethane)-trifluorophosphate, LiFAP dissolved in a binary mixture of ethylene carbonate (EC) and dimethylcarbonate (DMC) (50:50 wt%). For each electrolyte, the temperature dependence on transport properties over a temperature range from 10 to 80 degrees C and 20 to 70 degrees C for viscosity and conductivity, respectively, exhibits a non-Arrhenius behavior. However, this dependence is correctly correlated by using the Vogel-Tamman-Fulcher (VTF) type fitting equation. In each case, the best-fit parameters, such as the pseudo activation energy and ideal glass transition temperature were then extracted. The self-diffusion coefficients (D) of the Li+ cation and PF6- or FAP(-) anions species, in each studied electrolyte, were then independently determined by observing Li-3, F-19 and P-31 nuclei with the pulsed-gradient spin-echo (PGSE) NMR technique over the same temperature range from 20 to 80 degrees C. Results show that even if the diffusion of the lithium cation is quite similar in both electrolytes, the anions diffusion differs notably. In the case of the LiPF6-based electrolyte, for example at T approximate to 75 degrees C (high temperature), the self-diffusion coefficients of Li+ cations in solution (D (Li+)approximate to 5 x 10(-19) m(2) s(-1)) is 1.6 times smaller than that of PF6- anions (D (PF6-) = 8.5 x 10(-19) m(2) s(-1)), whereas in the case of the LiFAP-based electrolyte, FAP(-) anions diffuse at same rate as the Li+ cations (D (FAP(-)) = 5 x 10(-1) m(2) s(-1)). Based on these experimental results, the transport mobility of ions were then investigated through Stokes-Einstein and Nernst-Einstein equations to determine the transport number of lithium t(Li)(+), effective radius of solvated Li+ and of PF6- and FAP(-) anions, and the degree of dissociation of these lithium salts in the selected EC/DMC (50:50 wt%) mixture over a the temperature range from 20 to 80 degrees C. This study demonstrates the conflicting nature of the requirements and the advantage of the well-balanced properties as ionic mobility and dissociation constant of the selected electrolytes. (C) 2013 Elsevier Ltd. All rights reserved.

Generalized Langevin equation: An efficient approach to nonequilibrium molecular dynamics of open systems

The generalized Langevin equation (GLE) has been recently suggested to simulate the time evolution of classical solid and molecular systems when considering general nonequilibrium processes. In this approach, a part of the whole system (an open system), which interacts and exchanges energy with its dissipative environment, is studied. Because the GLE is derived by projecting out exactly the harmonic environment, the coupling to it is realistic, while the equations of motion are non-Markovian. Although the GLE formalism has already found promising applications, e. g., in nanotribology and as a powerful thermostat for equilibration in classical molecular dynamics simulations, efficient algorithms to solve the GLE for realistic memory kernels are highly nontrivial, especially if the memory kernels decay nonexponentially. This is due to the fact that one has to generate a colored noise and take account of the memory effects in a consistent manner. In this paper, we present a simple, yet efficient, algorithm for solving the GLE for practical memory kernels and we demonstrate its capability for the exactly solvable case of a harmonic oscillator coupled to a Debye bath.

Applications of the generalized Langevin equation: Towards a realistic description of the baths

The generalized Langevin equation (GLE) method, as developed previously [L. Stella et al., Phys. Rev. B 89, 134303 (2014)], is used to calculate the dissipative dynamics of systems described at the atomic level. The GLE scheme goes beyond the commonly used bilinear coupling between the central system and the bath, and permits us to have a realistic description of both the dissipative central system and its surrounding bath. We show how to obtain the vibrational properties of a realistic bath and how to convey such properties into an extended Langevin dynamics by the use of the mapping of the bath vibrational properties onto a set of auxiliary variables. Our calculations for a model of a Lennard-Jones solid show that our GLE scheme provides a stable dynamics, with the dissipative/relaxation processes properly described. The total kinetic energy of the central system always thermalizes toward the expected bath temperature, with appropriate fluctuation around the mean value. More importantly, we obtain a velocity distribution for the individual atoms in the central system which follows the expected canonical distribution at the corresponding temperature. This confirms that both our GLE scheme and our mapping procedure onto an extended Langevin dynamics provide the correct thermostat. We also examined the velocity autocorrelation functions and compare our results with more conventional Langevin dynamics.

Abrasion, transport and distribution of sediment in selected streams of Southern Ontario and Western New York /

The rate of decrease in mean sediment size and weight per square metre along a 54 km reach of the Credit River was found to depend on variations in the channel geometry. The distribution of a specific sediment size consist of: (1) a transport zone; (2) an accumulation zone; and (3) a depletion zone. These zones shift downstream in response to downcurrent decreases in stream competence. Along a .285 km man-made pond, within the Credit River study area, the sediment is also characterized by downstream shifting accumulation zones for each finer clast size. The discharge required to initiate movement of 8 cm and 6 cm blocks in Cazenovia Creek is closely approximated by Baker and Ritter's equation. Incipient motion of blocks in Twenty Mile Creek is best predicted by Yalin's relation which is more efficient in deeper flows. The transport distance of blocks in both streams depends on channel roughness and geometry. Natural abrasion and distribution of clasts may depend on the size of the surrounding sediment and variations in flow competence. The cumulative percent weight loss with distance of laboratory abraded dolostone is defined by a power function. The decrease in weight of dolostone follows a negative exponential. In the abrasion mill, chipping causes the high initial weight loss of dolostone; crushing and grinding produce most of the subsequent weight loss. Clast size was found to have little effect on the abrasion of dolostone within the diameter range considered. Increasing the speed of the mill increased the initial amount of weight loss but decreased the rate of abrasion. The abrasion mill was found to produce more weight loss than stream action. The maximum percent weight loss determined from laboratory and field abrasion data is approximately 40 percent of the weight loss observed along the Credit River. Selective sorting of sediment explains the remaining percentage, not accounted for by abrasion.

A semi-phenomenological approach to the structure and transport properties of macromolecules in solution

Motiviert durch die Lebenswissenschaften (Life sciences) haben sich Untersuchungen zur Dynamik von Makromolekülen in Lösungen in den vergangenen Jahren zu einem zukunftsweisenden Forschungsgebiet etabliert, dessen Anwendungen von der Biophysik über die physikalische Chemie bis hin zu den Materialwissenschaften reichen. Neben zahlreichen experimentellen Forschungsprogrammen zur räumlichen Struktur und den Transporteigenschaften grosser MolekÄule, wie sie heute praktisch an allen (Synchrotron-) Strahlungsquellen und den Laboren der Biophysik anzutreffen sind, werden gegenwärtig daher auch umfangreiche theoretische Anstrengungen unternommen, um das Diffusionsverhalten von Makromolekülen besser zu erklären. Um neue Wege für eine quantitative Vorhersagen des Translations- und Rotationsverhaltens grosser Moleküle zu erkunden, wurde in dieser Arbeit ein semiphänomenologischer Ansatz verfolgt. Dieser Ansatz erlaubte es, ausgehend von der Hamiltonschen Mechanik des Gesamtsystems 'Molekül + Lösung', eine Mastergleichung für die Phasenraumdichte der Makromoleküle herzuleiten, die den Einfluss der Lösung mittels effektiver Reibungstensoren erfasst. Im Rahmen dieses Ansatzes gelingt es z.B. (i) sowohl den Einfluss der Wechselwirkung zwischen den makromolekularen Gruppen (den sogenannten molekularen beads) und den Lösungsteilchen zu analysieren als auch (ii) die Diffusionseigen schaften für veschiedene thermodynamische Umgebungen zu untersuchen. Ferner gelang es auf der Basis dieser Näherung, die Rotationsbewegung von grossen Molekülen zu beschreiben, die einseitig auf einer Oberfläche festgeheftet sind. Im Vergleich zu den aufwendigen molekulardynamischen (MD) Simulationen grosser Moleküle zeichnet sich die hier dargestellte Methode vor allem durch ihren hohen `Effizienzgewinn' aus, der für komplexe Systeme leicht mehr als fünf Grössenordnungen betragen kann. Dieser Gewinn an Rechenzeit erlaubt bspw. Anwendungen, wie sie mit MD Simulationen wohl auch zukünftig nicht oder nur sehr zögerlich aufgegriffen werden können. Denkbare Anwendungsgebiete dieser Näherung betreffen dabei nicht nur dichte Lösungen, in denen auch die Wechselwirkungen der molekularen beads zu benachbarten Makromolekülen eine Rolle spielt, sondern auch Untersuchungen zu ionischen Flüssigkeiten oder zur Topologie grosser Moleküle.