Biomedical subjects

Zhong-Yuan Lu

Publications and source records attributed to Zhong-Yuan Lu.

10 recordsLinked to original sources

PubMed2007

Molecular dynamics simulation study on controlling the adsorption behavior of polyethylene by fine tuning the surface nanodecoration of graphite.

Molecular dynamics simulations are applied to study the adsorption of polyethylene with different chain lengths on patterned graphite surfaces that contain nanoscale protrusions. The influence of the nanostructure on the strong attractive interaction inherently in the hydrophobic polyethylene and hydrophobic graphite system is investigated by modifying the top surface area and the height and the shape of the protrusions. The results are analyzed in terms of the chain configuration, the adsorption energy, the global orientational order parameter, and the normalized surface-chain contacting pair number in the first adsorption layer. When the size of the protrusion increases, the adsorption energy, the order parameter, and the normalized surface-chain contacting pair number decrease at a fixed chain length. When the size of the protrusion is fixed, the average adsorption energy per monomer and the order parameter decrease with increasing chain length because of the stronger intramolecular interactions between the monomers. Changing the protrusion shape in a suitable way will effectively reduce the strong surface-chain interaction.

Journal Article↗

PubMed2006

An automatic coarse-graining and fine-graining simulation method: application on polyethylene.

Multiscale modeling of a polymeric system is a challenging task in polymer physics. Here we introduce a bottom-up and then top-down scheme for the simulation of polyethylene (PE). The coarse-grained numerical potential for PE is derived through an automatic updating program by mapping its radial distribution function (RDF) from the Lowe-Andersen temperature controlling (LA) simulation onto the one from detailed molecular dynamics (MD) simulation. This coarse-grained numerical potential can be applied in larger systems under the same thermodynamic conditions. We have tested the reliability of the derived potential in two ways. First, the blends of different linear low-density polyethylene (LLDPE) with high-density polyethylene (HDPE) have been simulated in LA with the coarse-grained numerical potentials and reasonable results are obtained. Moreover, Rouse scaling behavior is reproduced for monodispersed polymeric systems with different chain lengths. The atomistic details of the beads can be reintroduced into the coarse-grained HDPE and LLDPE/HDPE models, followed by a few MD runs to alleviate the local tension induced by this fine-graining procedure. The equilibrated large atomistic system can then be used for further studies.

Journal Article↗

PubMed2006

Dissipative particle dynamics study on the morphology changes of diblock copolymer lamellar microdomains due to steady shear.

The morphology changes of linear diblock copolymer lamellar microdomains under uniform simple shear are studied via the dissipative particle dynamics technique. The parallel and perpendicular reorientations of the lamellae are observed in the simulations, and two different reorientation mechanisms, under small and large shear rates respectively, are proposed. The parallel-to-perpendicular transition is also observed and the kinetics is discussed. Sinusoidal and chevron instabilities due to the shear are found. After relaxation the peculiar "bidirectionally undulating" instability is obtained.

Journal Article↗

PubMed2006

The dependence of nanostructures on the molecule rigidity of A2(B4)2-type miktoarm block copolymer.

Using the dissipative particle dynamics simulation technique, we have studied the influence of the molecule rigidity on the nanostructures of the A2(B4)2-type miktoarm block copolymers. A typical spherical micellar ordered structure is obtained for a coil-coil miktoarm block copolymer in melt. By introducing a bond angle potential in our model to enhance the molecule rigidity systematically, we find, respectively, a hexagonal cylindrical structure and a parallel ellipsoid in lamellae structure which is discovered for the first time.

Journal Article↗

PubMed2005

A molecular-dynamics simulation study of diffusion of a single model carbonic chain on a graphite (001) surface.

Molecular-dynamics simulations have been used to study the diffusion of a short single model carbonic chain on the graphite (001) surface. The calculated diffusion coefficient (D) first increases, then decreases with increasing chain length (N). This abnormal behavior is similar to polymer lateral diffusion at the solid-liquid interface. Furthermore, we have studied the relation between the mean-square gyration radius and N. [Figure: see text].

Carbon↗

PubMed2005

A molecular-dynamics simulation study on the dependence of Lennard-Jones gas-liquid phase diagram on the long-range part of the interactions.

The particle-transfer molecular-dynamics technique is adopted to construct the Lennard-Jones fluid gas-liquid phase diagram. Detailed study of the dependence of the simulation results on the system size and the cutoff distance is performed to test the validity of the simulation technique. Both the traditional cutoff plus long-range correction (CPC) and Ewald summation methods are used in the simulations to calculate the interactions. In the intermediate range of temperatures, the results with the Ewald summation method are almost the same as those with the CPC method. However, in the range close to the critical point, the results with the CPC method deviate from those with the Ewald summation. Compared with the results obtained via the Ewald summation in a smaller system, simply increasing the system size in the CPC scheme may not give better results.

Journal Article↗

PubMed2005

Molecular dynamics simulation study on adsorption and diffusion processes of a hydrophilic chain on a hydrophobic surface.

Molecular dynamics simulations are applied to investigate the adsorption and diffusion processes of a single hydrophilic poly(vinyl alcohol) (PVA) chain with different chain lengths on a hydrophobic graphite surface. It is expected that the chain and the surface "dislike" each other because one is hydrophilic and the other is hydrophobic. But surprisingly, a short PVA chain is well adsorbed on the surface, accompanied by large changes in the chain configuration. With increasing degree of polymerization (N), the chain turns gradually from two-dimensional adsorption to possessing certain height in the direction perpendicular to the surface. Moreover, the adsorption energy increases and the diffusion coefficient decreases with increasing N. In particular, for N = 20 in equilibrium, the hydroxyls of this short chain are close to the graphite surface in the stable adsorption configuration. In addition, we change the effective dielectric constant to 76.0 to mimic good solvent condition. The chain configurations and the diffusion coefficients both vary in contrast to the foregoing results.

Journal Article↗

PubMed2005

Dissipative particle dynamics study on the interfaces in incompatible A/B homopolymer blends and with their block copolymers.

Dissipative particle dynamics, a simulation technique appropriate at mesoscopic scales, has been applied to investigate the interfaces in immiscible binary A/B homopolymer blends and in the ternary systems with their block copolymers. For the binary blends, the interfacial tension increases and the interface thickness decreases with increasing Flory-Huggins interaction parameter chi while the homopolymer chain length is fixed. However, when the chi parameter and one of the homopolymer chain length is fixed, increasing another homopolymer chain length will induce only a small increase on interfacial tension and slight decrease on interface thickness. For the ternary blends, adding the A-b-B block copolymer will reduce the interfacial tension. When the mole number of the block copolymer is fixed, longer block chains have higher efficiency on reducing the interfacial tension than the shorter ones. But for the block copolymers with fixed volume fraction, shorter chains will be more efficient than the longer ones on reducing the interfacial tension. Increasing the block copolymer concentration reduces interfacial tension. This effect is more prominent for shorter block copolymer chains.

Journal Article↗

PubMed2005

Stability of two-dimensional tessellation ice on the hydroxylated beta-cristobalite (100) surface.

Monolayer adsorbed water on the beta-cristobalite (100) surface is studied via classical molecular dynamics simulations. The ordered two-dimensional (2D) tessellation ice structure (i.e., the four-membered and the eight-membered rings appear alternatively) is justified at low temperatures in the simulations. The stability of this possible new ice phase is further investigated by heating the system from 5 to 300 K. An order-disorder structural transition is observed between 100 and 200 K, featuring the melting process of the tessellation ice. This process is characterized by the water oxygen-oxygen radial distribution function, the coordination number, the distance vector between the center of mass of the oxygen and the hydrogen atoms in water, the mean square displacement of oxygen in water, and the vibrational density of state. The above techniques show consistency on that the order-disorder transition temperature of the 2D tessellation ice is far below 300 K. The 2D tessellation ice structure is also obtained via density functional calculations with different generalized gradient approximations. By comparing the calculated adsorption and the lateral energies between different methods, we find that the melting temperature of the specific 2D ice structure is strongly method dependent. Therefore, further experimental works are urged to justify this possible new ice phase and probe its stability.

Computer Simulation↗

PubMed2005

The effects of Lowe-Andersen temperature controlling method on the polymer properties in mesoscopic simulations.

Lowe-Andersen (LA) temperature controlling method [C. P. Lowe, Europhys. Lett. 47, 145 (1999)] is applied in a series of mesoscopic polymer simulations to test its validity and efficiency. The method is an alternative for dissipative particle dynamics simulation (DPD) technique which is also Galilean invariant. It shows excellent temperature control and gives correct radial distribution function as that from DPD simulation. The efficiency of LA method is compared with other typical DPD integration schemes and is proved to be moderately efficient. Moreover, we apply this approach to diblock copolymer microphase separation simulations. With LA method, we are able to reproduce all the results from the conventional DPD simulations. The calculated structure factors of the microphases are consistent with the experiments. We also study the microphase evolution dynamics with increasing chiN and find that the bath collision frequency Gamma does not affect the order of appearing phases. Although the thermostat does not affect the surface tension, the order-disorder transition (ODT) is somewhat sensitive to the values of Gamma, i.e., the ODT is nonmonotonic with increasing Gamma. The dynamic scaling law is also tested, showing that the relation obeys the Rouse theory with various Gamma.

Journal Article↗