Molecular Dynamics Simulation Of Adhesion Forces In A Dipalmitoylphosphatidylcholine Membrane

Adhesion forces of Dipalmitoylphosphatidylcholine ( DPPC) membrane in the gel phase are investigated by molecular dynamics ( MD) simulation. In the simulations, individual DPPC molecules are pulled out of DPPC membranes with different rates and we get the maximum adhesion forces of DPPC membrane. We find that the maximum adhesion forces increase with pull rate, from about 400 to 700 pN when pull rates are from 0.001 to 0.03 nm/ps. We analyze the relationship between pull rate and adhesion forces of different origins using Brownian dynamics and notice that viscosity of solvent plays an important role in adhesion forces. Then we simulate the motion of a single DPPC molecule in solvent and it elucidates that the maximum drag force is almost linear with respect to the pull rate. We use Stokes' relation to describe the motion of a single DPPC molecule and deduce the effective length of a DPPC molecule. Conformational analyses indicate that the free energy variation of a DPPC molecule inside and outside of the DPPC membrane is an essential part of adhesion energy.

Strengthening Mechanisms And Dislocation Dynamics In Twinned Metal Nanowires

Microtwins are frequently observed in face-centered-cubic (fcc) metal nanowires with low stacking fault energy. The authors have previously reported that the tensile Yield strength of copper nanowires can be increased by, the presence of twin boundaries. lit this work, simulations are carried out under both uniaxial tension and compression loading, to demonstrate that the strengthening effects are inherent to these nanowires, independent of the loading condition (tensile/compressive). It appears that the strengthening mechanism of the twinned nanowires can be attributed to stress redistribution due to the change of crystallographic orientations across twin boundaries, which requires larger external stress to make them Yield as compared to the twin-free wire.

Study Of Mechanical Properties Of Amorphous Copper With Molecular Dynamics Simulation

The formation and mechanical properties of amorphous copper are studied using molecular dynamics simulation. The simulations of tension and shearing show that more pronounced plasticity is found under shearing, compared to tension. Apparent strain hardening and strain rate effect are observed. Interestingly, the variations of number density of atoms during deformation indicate free volume creation, especially under higher strain rate. In particular, it is found that shear induced dilatation does appear in the amorphous metal.

A Note On The Numerical Simulations Of Flow Past A Wavy Square-Section Cylinder

The flow past a square-section cylinder with a geometric disturbance is investigated by numerical simulations. The extra terms, due to the introduction of mapping transformation simulating the effect of disturbance into the transformed Navier-Stokes equations, are correctly derived, and the incorrect ones in the previous literature are pointed out and analyzed. Furthermore, the relationship between the vorticity, especially on the cylinder surface, and the disturbance is derived and explained theoretically. The computations are performed at two Reynolds numbers of 100 and 180 and three amplitudes of waviness of 0.006, 0.025 and 0.167 with another aim to explore the effects of different Reynolds numbers and disturbance on the vortex dynamics in the wake and forces on the body. Numerical results have shown that, at the mild waviness of 0.025, the Karman vortex shedding is suppressed completely for Re = 100, while the forced vortex dislocation is appeared in the near wake at the Reynolds number of 180. The drag reduction is up to 21.6% at Re = 100 and 25.7% at Re = 180 for the high waviness of 0.167 compared with the non-wavy cylinder. The lift and the Strouhal number varied with different Reynolds numbers and the wave steepness are also obtained.

Adsorption of His-tagged peptide to Ni, Cu and Au (100) surfaces: Molecular dynamics simulation

Molecular dynamics (MD) simulations are performed to study the interaction of His-tagged peptide with three different metal surfaces in explicit water. The equilibrium properties are analyzed by using pair correlation functions (PCF) to give an insight into the behavior of the peptide adsorption to metal surfaces in water solvent. The intermolecular interactions between peptide residues and the metal surfaces are evaluated. By pulling the peptide away from the peptide in the presence of solvent water, peeling forces are obtained and reveal the binding strength of peptide adsorption on nickel, copper and gold. From the analysis of the dynamics properties of the peptide interaction with the metal surfaces, it is shown that the affinity of peptide to Ni surface is the strongest, while on Cu and An the affinity is a little weaker. In MD simulations including metals, the His-tagged region interacts with the substrate to an extent greater than the other regions. The work presented here reveals various interactions between His-tagged peptide and Ni/Cu/Au surfaces. The interesting affinities and dynamical properties of the peptide are also derived. The results give predictions for the structure of His-tagged peptide adsorbing on three different metal surfaces and show the different affinities between them, which assist the understanding of how peptides behave on metal surfaces and of how designers select amino sequences in molecule devices design. (c) 2007 Elsevier Ltd. All rights reserved.

Molecular dynamics studies on the local disordering of σ3 and σ11 grain boundaries in aluminium bicrystal

Molecular dynamics (MD) simulations using Morse interaction potential are performed in studies of [110] symmetrical tilt grain boundary (GB) structures with mis-orientation angles 50.5 degrees(Sigma 11), 129.5 degrees(Sigma 11), 70.5 degrees(Sigma 3) and 109.5 degrees(Sigma 3) at various tempratures. The GB structures are found to start local disordering at about 0.5T(m)(T-m is the melting point of aluminium) for 50.5 degrees(Sigma 11), 0.32T(m) for 129.5 degrees(Sigma 11) and 0.38T(m) for 70.5 degrees(Sigma 3), respectively. These results agree with conclusions deduced from the anelastic measurements. But, for twin-boundary structure 109.5 degrees(Sigma 3), this disordering has not been found even when temperature increases up to 0.9T(m).

Molecular Dynamics Simulation of Barnacle Cement

Barnacle cement is an underwater adhesive that is used for permanent settlement. Its main components are insoluble protein complexes that have not been fully studied. In present article, we chose two proteins of barnacle cement for study, 36-KD protein and Mrcp-100K protein. In order to investigate the characteristic of above two proteins, we introduced the method of molecular modeling. And the simulation package GROMACS was used to simulate the behavior of these proteins. In this article, before the simulations, we introduce some theories to predict the time scale for polymer relaxation. During the simulation, we mainly focus on two properties of these two proteins: structural stability and adhesive force to substrate. First, we simulate the structural stability of two proteins in water, and then the stability of 36-KD protein in seawater environment is investigated. We find that the stability varies in the different environments. Next, to study adhesive ability of two proteins, we simulate the process of peeling the two proteins from the substrate (graphite). Then, we analyze the main reasons of these results. We find that hydrogen bonds in proteins play an important role in the protein stability. In the process of the peeling, we use Lennard-Jones 12-6 potential to calculate the van der Waals interactions between proteins and substrate.

Size-dependent elastic properties of Ni nanofilms by molecular dynamics simulation

Size-dependent elastic properties of Ni nanofilms are investigated by molecular dynamics ( MD) simulations with embedded atom method (EAM). The surface effects are considered by calculating the surface relaxation, surface energy, and surface stress. The Young's modulus and yield stress are obtained as functions of thickness and crystallographic orientation. It is shown that the surface relaxation has important effects on the the elastic properties at nanoscale. When the surface relaxation is outward, the Young's modulus decreases with the film thickness decreasing, and vice versa. The results also show that the yield stresses of the films increase with the films becoming thinner. With the thickness of the nanofilms decreasing, the surface effects on the elastic properties become dominant.

Molecular dynamics simulation of adhesion forces in a dipalmitoylphosphatidylcholine membrane

Adhesion forces of Dipalmitoylphosphatidylcholine ( DPPC) membrane in the gel phase are investigated by molecular dynamics ( MD) simulation. In the simulations, individual DPPC molecules are pulled out of DPPC membranes with different rates and we get the maximum adhesion forces of DPPC membrane. We find that the maximum adhesion forces increase with pull rate, from about 400 to 700 pN when pull rates are from 0.001 to 0.03 nm/ps. We analyze the relationship between pull rate and adhesion forces of different origins using Brownian dynamics and notice that viscosity of solvent plays an important role in adhesion forces. Then we simulate the motion of a single DPPC molecule in solvent and it elucidates that the maximum drag force is almost linear with respect to the pull rate. We use Stokes' relation to describe the motion of a single DPPC molecule and deduce the effective length of a DPPC molecule. Conformational analyses indicate that the free energy variation of a DPPC molecule inside and outside of the DPPC membrane is an essential part of adhesion energy.

A comparative study of Young's modulus of single-walled carbon nanotube by CPMD, MD, and first principle simulations

Carbon nanotubes (CNTs), due to their exceptional magnetic, electrical and mechanical properties, are promising candidates for several technical applications ranging from nanoelectronic devices to composites. Young's modulus holds the special status in material properties and micro/nano-electromechanical systems (MEMS/NEMS) design. The excellently regular structures of CNTs facilitate accurate simulation of CNTs' behavior by applying a variety of theoretical methods. Here, three representative numerical methods, i.e., Car-Parrinello molecular dynamics (CPMD), density functional theory (DFT) and molecular dynamics (MD), were applied to calculate Young's modulus of single-walled carbon nanotube (SWCNT) with chirality (3,3). The comparative studies showed that the most accurate result is offered by time consuming DFT simulation. MID simulation produced a less accurate result due to neglecting electronic motions. Compared to the two preceding methods the best performance, with a balance between efficiency and precision, was deduced by CPMD.

Electron dynamics and harmonics emission spectra due to electron oscillation driven by intense laser pulses

The dynamics and harmonics emission spectra due to electron oscillation driven by intense laser pulses have been investigated considering a single electron model. The spectral and angular distributions of the harmonics radiation are numerically analyzed and demonstrate significantly different characteristics from those of the low-intensity field case. Higher-order harmonic radiation is possible for a sufficiently intense driving laser pulse. A complex shifting and broadening structure of the spectrum is observed and analyzed for different polarization. For a realistic pulsed photon beam, the spectrum of the radiation is redshifted for backward radiation and blueshifted for forward radiation, and spectral broadening is noticed. This is due to the changes in the longitudinal velocity of the electron during the laser pulse. These effects are much more pronounced at higher laser intensities giving rise to even higher-order harmonics that eventually leads to a continuous spectrum. Numerical simulations have further shown that broadening of the high harmonic radiation can be limited by increasing the laser pulse width. The complex shifting and broadening of the spectra can be employed to characterize the ultrashort and ultraintense laser pulses and to study the ultrafast dynamics of the electrons. (c) 2006 American Institute of Physics.

Ultrafast carrier dynamics in undoped and p-doped InAs/GaAs quantum dots characterized by pump-probe reflection measurements

We investigate the dependence of the differential reflection on the structure parameters of quantum dot (QD) heterostructures in pump-probe reflection measurements by both numerical simulations based on the finite-difference time-domain technique and theoretical calculations based on the theory of dielectric films. It is revealed that the value and sign of the differential reflection strongly depend on the thickness of the cap layer and the QD layer. In addition, a comparison between the carrier dynamics in undoped and p-doped InAs/GaAs QDs is carried out by pump-probe reflection measurements. The carrier capture time from the GaAs barrier into the InAs wetting layer and that from the InAs wetting layer into the InAs QDs are extracted by appropriately fitting differential reflection spectra. Moreover, the dependence of the carrier dynamics on the injected carrier density is identified. A detailed analysis of the carrier dynamics in the undoped and p-doped QDs based on the differential reflection spectra is presented, and its difference with that derived from the time-resolved photoluminescence is discussed. (C) 2008 American Institute of Physics.

Non-exponential photoluminescence decay dynamics of localized carriers in disordered InGaN/GaN quantum wells: the role of localization length

In this article, we report a combined experimental and theoretical study on the luminescence dynamics of localized carriers in disordered InGaN/GaN quantum wells. The luminescence intensity of localized carriers is found to exhibit an unusual non-exponential decay. Adopting a new model taking the radiative recombination and phonon-assisted hopping transition between different localized states into account, which was recently developed by Rubel et al., the non-exponential decay behavior of the carriers can be quantitatively interpreted. Combining with precise structure characterization, the theoretical simulations show that the localization length of localized carriers is a key parameter governing their luminescence decay dynamics. (c) 2006 Optical Society of America.

Dynamics of Incoherent Photovoltaic Spatial Solitons

Propagation properties of bright and dark incoherent beams are numerically studied in photovoltaic-photorefractive crystal by using coherent density approach for the first time. Numerical simulations not only exhibit that bright incoherent photovoltaic quasi-soliton, grey-like incoherent photovoltaic soliton, incoherent soliton doublet and triplet can be established under proper conditions, but also display that the spatial coherence properties of these incoherent beams can be significantly affected during propagation by the photovoltaic field.

Brownian dynamics simulation of the crystallization dynamics of charged colloidal particles

Crystal formation process of charged colloidal particles is investigated using Brownian dynamics (BD) simulations. The particles are assumed to interact with the pair-additive repulsive Yukawa potential. The time evolution of crystallization process and the crystal structure during the simulation are characterized by means of the radial distribution functions (RDF) and mean square displacement (MSD). The simulations show that when the interaction is featured with long-range, particles can spontaneously assemble into body-centered-cubic (BCC) arrays at relatively low particle number density. When the interaction is short-ranged, with increasing the number density particles become trapped into a stagnant disordered configuration before the crystallization could be actualized. The simulations further show that as long as the trapped configurations are bypassed, the face-centered-cubic (FCC) structures can be achieved and are actually more stable than BCC structures. (C) 2010 Elsevier Inc. All rights reserved.