Biomedical subjects

Kurt Binder

Publications and source records attributed to Kurt Binder.

8 recordsLinked to original sources

PubMed2004

Polymer nanodroplets forming liquid bridges in chemically structured slit pores: a computer simulation.

Using a coarse-grained bead-spring model of flexible polymer chains, the structure of a polymeric nanodroplet adsorbed on a chemically decorated flat wall is investigated by means of molecular dynamics simulation. We consider sessile drops on a lyophilic (attractive for the monomers) region of circular shape with radius R(D) while the remaining part of the substrate is lyophobic. The variation of the droplet shape, including its contact angle, with R(D) is studied, and the density profiles across these droplets also are obtained. In addition, the interaction of droplets adsorbed on two walls forming a slit pore with two lyophilic circular regions just opposite of one another is investigated, paying attention to the formation of a liquid bridge between both walls. A central result of our study is the measurement of the force between the two substrate walls at varying wall separation as well as the kinetics of droplet merging. Our results are compared to various phenomenological theories developed for liquid droplets of mesoscopic rather than nanoscopic size.

Journal Article↗

PubMed2004

Polymer translocation through a nanopore induced by adsorption: Monte Carlo simulation of a coarse-grained model.

Dynamic Monte Carlo simulation of a bead-spring model of flexible macromolecules threading through a very narrow pore in a very thin rigid membrane are presented, assuming at the cis side of the membrane a purely repulsive monomer-wall interaction, while the trans side is attractive. Two choices of monomer-wall attraction epsilon are considered, one choice is slightly below and the other slightly above the "mushroom to pancake" adsorption threshold epsilon(c) for an infinitely long chain. Studying chain lengths N=32, 64, 128, and 256 and varying the number of monomers N(trans) (time t=0) that have already passed the pore when the simulation started, over a wide range, we find for epsilon epsilon(c) a finite number N(trans)(t=0) suffices that the translocation probability is close to unity. In the case epsilon epsilon(c), we find that the translocation time scales as tau proportional, variant N(1.65+/-0.08). We suggest a tentative scaling explanation for this result. Also the distribution of translocation times is obtained and discussed.

Journal Article↗

PubMed2004

Stability of thin polymer films: influence of solvents.

The interface and surface properties and the wetting behavior of polymer-solvent mixtures are investigated using Monte Carlo simulations and self-consistent field calculations. We carry out Monte Carlo simulations in the framework of a coarse-grained bead-spring model using short chains (oligomers) of N(P)=5 beads and a monomeric solvent, N(S)=1. The self-consistent field calculations are based on a simple phenomenological equation of state for compressible binary mixtures and we employ Gaussian chain model. The bulk behavior of the polymer-solvent mixture belongs to type III in the classification of van Konynenburg and Scott [Phil. Trans. R. Soc. London, Ser. A 298, 495 (1980)]. It is characterized by a triple line on which the polymer-liquid coexists with solvent-vapor and a solvent-rich liquid. The solvent is not homogeneously distributed across the dense polymer film but tends to accumulate at the surface and the polymer-vapor interface. This solvent enrichment at the interface and surface becomes more pronounced upon increasing the vapor pressure and alters the surface and interface tensions. This effect gives rise to a nonmonotonic dependence of the contact angle on the vapor pressure and one might observe reentrant wetting. The results of the Monte Carlo simulations and the self-consistent field calculations qualitatively agree. The profiles of drops are investigated by Monte Carlo simulations and a pronounced solvent enrichment is observed at the wedge formed by the substrate and the liquid-vapor interface at the three-phase contact line.

Journal Article↗

PubMed2004

The evaporation/condensation transition of liquid droplets.

The condensation of a supersaturated vapor enclosed in a finite system is considered. A phenomenological analysis reveals that the vapor is found to be stable at densities well above coexistence. The system size at which the supersaturated vapor condenses into a droplet is found to be governed by a typical length scale which depends on the coexistence densities, temperature and surface tension. When fluctuations are neglected, the chemical potential is seen to show a discontinuity at an effective spinodal point, where the inhomogeneous state becomes more stable than the homogeneous state. If fluctuations are taken into account, the transition is rounded, but the slope of the chemical potential versus density isotherm develops a discontinuity in the thermodynamic limit. In order to test the theoretical predictions, we perform a simulation study of droplet condensation for a Lennard-Jones fluid and obtain loops in the chemical potential versus density and pressure. By computing probability distributions for the cluster size, chemical potential, and internal energy, we confirm that the effective spinodal point may be identified with the occurrence of a first order phase transition, resulting in the condensation of a droplet. An accurate equation of state is employed in order to estimate the droplet size and the coexisting vapor density and good quantitative agreement with the simulation data is obtained. The results highlight the need of an accurate equation of state data for the Laplace equation to have predictive power.

Journal Article↗

PubMed2004

Analysis of the phase transitions in alkyl-mica by density and pressure profiles.

In a previous work [Heinz, Castelijns, and Suter, J. Am. Chem. Soc. 115, 9500 (2003)], we developed an accurate force field and simulated the phase transitions in C18-mica (octadecyltrimethylammonium-mica) as well as the absence of such transitions in 2C18-mica (dioctadecyldimethylammonium-mica) between room temperature and 100 degrees C. Here we analyze (i) average z coordinates of the carbon atoms and interdigitation of the hydrocarbon bilayers, (ii) density profiles, and (iii) pressure profiles of the structures along all Cartesian axes. In C18-mica, the standard deviation in the z coordinate for the chain atoms is high and more than doubles in the disordered phase. The order-disorder transition is accompanied by a change in the orientation of the ammonium head group, as well as decreasing tensile and shear stress in the disordered phase. In 2C18-mica, the standard deviation in the z coordinate for the chain atoms is low and does not increase remarkably on heating. The backbones display a highly regular structure, which is slightly obscured by rotations in the C18 backbones and minor head group displacements at 100 degrees C. Close contacts between the bulky head groups with sidearms cause significant local pressure which is in part not relieved at 100 degrees C. An increase of the basal-plane spacing at higher temperature is found in both systems due to larger separation between the two hydrocarbon layers and an increased z spacing between adjacent chain atoms (=decreased tilt of the chains relative to the surface normal), and, in C18-mica only, a stronger upward orientation of the C18 chain at the ammonium head group. The likelihood for chain interdigitation between the two hydrocarbon layers is 24%-30% for C18-mica, and 65%-26% for 2C18-mica (for 20-100 degrees C).

Journal Article↗

PubMed2004

Structure and diffusion in amorphous aluminum silicate: a molecular dynamics computer simulation.

The amorphous aluminum silicate (Al2O3)2(SiO2) [AS2] is investigated by means of large scale molecular dynamics computer simulations. We consider fully equilibrated melts in the temperature range 6100 K> or =T> or =2300 K as well as glass configurations that were obtained from cooling runs from T=2300 to 300 K with a cooling rate of about 10(12) K/s. Already at temperatures as high as 4000 K, most of the Al and Si atoms are fourfold coordinated by oxygen atoms. Thus, the structure of AS2 is that of a disordered tetrahedral network. The packing of AlO4 tetrahedra is very different from that of SiO4 tetrahedra in that Al is involved with a relatively high probability in small-membered rings and in triclusters in which an O atom is surrounded by four cations. We find as typical configurations two-membered rings with two Al atoms in which the shared O atoms form a tricluster. On larger length scales, the system shows a microphase separation in which the Al-rich network structure percolates through the SiO2 network. The latter structure gives rise to a prepeak in the static structure factor at a wave number q=0.5 A(-1). A comparison of experimental x-ray data with the results from the simulation shows good agreement for the structure function. The diffusion dynamics in AS2 is found to be much faster than in SiO2. We show that the self-diffusion constants for O and Al are very similar and that they are by a factor of 2-3 larger than the one for Si.

Journal Article↗

PubMed2002

Surface-directed phase separation with off-critical composition: analytical and numerical results.

We study the interplay of wetting and phase separation in an unstable binary mixture (AB) with off-critical composition, placed in contact with a surface which prefers the component A. We consider surface potentials V(z) approximately z(-n), where z is the distance from the surface, and present analytical arguments and detailed numerical results to elucidate wetting-layer kinetics for arbitrary mixture compositions. If the preferred component is the minority phase, the wetting-layer thickness exhibits a potential-specific behavior at early times tau, R1 approximately tau(1/(n+2)), before crossing over to the universal growth law, R1 approximately tau(1/3). On the other hand, if the preferred component is the majority phase, there is a crossover from potential-specific growth (as before) to a slower growth regime.

Journal Article↗

PubMed2002

Dynamics of sodium in sodium disilicate: channel relaxation and sodium diffusion.

We use molecular-dynamics computer simulations to study the dynamics of amorphous (Na2O)2(SiO2). We find that the Na ions move in channels embedded in a SiO2 matrix. The characteristic distance between these channels gives rise to a prepeak in the structure factor at q approximately equal to 0.95 A(-1). The dynamics of sodium is given by a fast process which can be seen in the incoherent scattering function and a slow process which is seen in the coherent function. The relaxation time of the latter coincides with the alpha-relaxation time of the matrix. The Kohlrausch exponent of the fast process for q>1.6 A(-1) is the same as the von Schweidler exponent for the slow one. Thus the two processes are closely related.

Journal Article↗