Bimodality and re-entrant behaviour in the hierarchical self-assembly of polymeric nanoparticles

We show that a film of a suspension of polymer grafted nanoparticles on a liquid substrate can be employed to create two-dimensional nanostructures with a remarkable variation in the pattern length scales. The presented experiments also reveal the emergence of concentration-dependent bimodal patterns as well as re-entrant behaviour that involves length scales due to dewetting and compositional instabilities. The experimental observations are explained through a gradient dynamics model consisting of coupled evolution equations for the height of the suspension film and the concentration of polymer. Using a Flory-Huggins free energy functional for the polymer solution, we show in a linear stability analysis that the thin film undergoes dewetting and/or compositional instabilities depending on the concentration of the polymer in the solution. We argue that the formation via `hierarchical self-assembly' of various functional nanostructures observed in different systems can be explained as resulting from such an interplay of instabilities.

Ab Initio MD Simulations of the Bronsted Acidity of Glutathione in Aqueous Solutions: Predicting pK(a) Shifts of the Cysteine Residue

The tripeptide glutathione (GSH) is one of the most abundant peptides and the major repository for nonprotein sulfur in both animal and plant cells. It plays a critical role in intracellular oxidative stress management by the reversible formation of glutathione disulfide with the thiol-disulfide pair acting as a redox buffer. The state of charge of the ionizable groups of GSH can influence the redox couple, and hence the pK(a) value of the cysteine residue of GSH is critical to its functioning. Here we report ab initio Car-Parrinello molecular dynamics simulations of glutathione solvated by 200 water molecules, all of which are considered in the simulation. We show that the free-energy landscape for the protonation-deprotonation reaction of the cysteine residue of GSH computed using metadynamics sampling provides shift in the dissociation constant values as compared with the isolated accurate estimates of the pK(a) and correctly predicts the cysteine amino acid.

Doping dependence of fluctuation diamagnetism in high T-c superconductors

Using a recently proposed Ginzburg-Landau-like lattice free energy functional due to Banerjee et al. (2011) we calculate the fluctuation diamagnetism of high -T-c superconductors as a function of doping, magnetic field and temperature. We analyse the pairing fluctuations above the superconducting transition temperature in the cuprates, ranging from the strong phase fluctuation dominated underdoped limit to the more conventional amplitude fluctuation dominated overdoped regime. We show that a model where the pairing scale increases and the superfluid density decreases with underdoping produces features of the observed magnetization in the pseudogap region, in good qualitative and reasonable quantitative agreement with the experimental data. In particular, we explicitly show that even when the pseudogap has a pairing origin the magnetization actually tracks the superconducting dome instead of the pseudogap temperature, as seen in experiment. We discuss the doping dependence of the `onset' temperature for fluctuation diamagnetism and comment on the role of vortex core -energy jn our model. (C) 2015 Elsevier Inc. All rights reserved.

Structure-rheology relationship in a sheared lamellar fluid

The structure-rheology relationship in the shear alignment of a lamellar fluid is studied using a mesoscale model which provides access to the lamellar configurations and the rheology. Based on the equations and free energy functional, the complete set of dimensionless groups that characterize the system are the Reynolds number (rho gamma L-2/mu), the Schmidt number (mu/rho D), the Ericksen number (mu(gamma)/B), the interface sharpness parameter r, the ratio of the viscosities of the hydrophilic and hydrophobic parts mu(r), and the ratio of the system size and layer spacing (L/lambda). Here, rho and mu are the fluid density and average viscosity, (gamma) over dot is the applied strain rate, D is the coefficient of diffusion, B is the compression modulus, mu(r) is the maximum difference in the viscosity of the hydrophilic and hydrophobic parts divided by the average viscosity, and L is the system size in the cross-stream direction. The lattice Boltzmann method is used to solve the concentration and momentum equations for a two dimensional system of moderate size (L/lambda = 32) and for a low Reynolds number, and the other parameters are systematically varied to examine the qualitative features of the structure and viscosity evolution in different regimes. At low Schmidt numbers where mass diffusion is faster than momentum diffusion, there is fast local formation of randomly aligned domains with ``grain boundaries,'' which are rotated by the shear flow to align along the extensional axis as time increases. This configuration offers a high resistance to flow, and the layers do not align in the flow direction even after 1000 strain units, resulting in a viscosity higher than that for an aligned lamellar phase. At high Schmidt numbers where momentum diffusion is fast, the shear flow disrupts layers before they are fully formed by diffusion, and alignment takes place by the breakage and reformation of layers by shear, resulting in defects (edge dislocations) embedded in a background of nearly aligned layers. At high Ericksen number where the viscous forces are large compared to the restoring forces due to layer compression and bending, shear tends to homogenize the concentration field, and the viscosity decreases significantly. At very high Ericksen number, shear even disrupts the layering of the lamellar phase. At low Ericksen number, shear results in the formation of well aligned layers with edge dislocations. However, these edge dislocations take a long time to anneal; the relatively small misalignment due to the defects results in a large increase in viscosity due to high layer stiffness and due to shear localization, because the layers between defects get pinned and move as a plug with no shear. An increase in the viscosity contrast between the hydrophilic and hydrophobic parts does not alter the structural characteristics during alignment. However, there is a significant increase in the viscosity, due to pinning of the layers between defects, which results in a plug flow between defects and a localization of the shear to a part of the domain.

Wave Propagation Analysis in a Pressure-Wave-Refrigerator

Pressure wave refrigerators (PWR) refrigerate the gas through periodical expansion waves. Due to its simple structure and robustness, PWR may have many potential applications if the efficiency becomes competitive with existing alternative devices. In order to improve the efficiency, the characteristics of wave propagation in a PWR are studied by experiment, numerical simulation and theoretical analysis. Based on the experimental results and numerical simulation, a simplified model is suggested, which includes the assumptions of flux-equilibrium and conservation of the free energy. This allows the independent analysis of the operation parameters and design specifics. Furthermore, the optimum operation condition can be deduced. Some considerations to improve the PWR efficiency are also given.

Surface attachment of protein fibrils via covalent modification strategies.

Chemical control of surface functionality and topography is an essential requirement for many technological purposes. In particular, the covalent attachment of monomeric proteins to surfaces has been the object of intense studies in recent years, for applications as varied as electrochemistry, immuno-sensing, and the production of biocompatible coatings. Little is known, however, about the characteristics and requirements underlying surface attachment of supramolecular protein nanostructures. Amyloid fibrils formed by the self-assembly of peptide and protein molecules represent one important class of such structures. These highly organized beta-sheet-rich assemblies are a hallmark of a range of neurodegenerative disorders, including Alzheimer's disease and type II diabetes, but recent findings suggest that they have much broader significance, potentially representing the global free energy minima of the energy landscapes of proteins and having potential applications in material science. In this paper, we describe strategies for attaching amyloid fibrils formed from different proteins to gold surfaces under different solution conditions. Our methods involve the reaction of sulfur containing small molecules (cystamine and 2-iminothiolane) with the amyloid fibrils, enabling their covalent linkage to gold surfaces. We demonstrate that irreversible attachment using these approaches makes possible quantitative analysis of experiments using biosensor techniques, such as quartz crystal microbalance (QCM) assays that are revolutionizing our understanding of the mechanisms of amyloid growth and the factors that determine its kinetic behavior. Moreover, our results shed light on the nature and relative importance of covalent versus noncovalent forces acting on protein superstructures at metal surfaces.

Experimental Investigation of the Velocity Effect on Adhesion Forces with an Atomic Force Microscope Atomic Force Microscope

Capillary forces are significantly dominant in adhesive forces measured with an atomic force microscope (AFM) in ambient air, which are always thought to be dependent on water film thickness, relative humidity, and the free energy of water film. We study the nature of the pull-off force on a variety of surfaces as a function of tip velocity. It is found that the capillary forces are of relatively strong dependence on tip velocity. The present experiment is expected to provide a better understanding of the work mechanism of AFM in ambient air.

Population of nonnative states of lysozyme variants drives amyloid fibril formation.

The propensity of protein molecules to self-assemble into highly ordered, fibrillar aggregates lies at the heart of understanding many disorders ranging from Alzheimer's disease to systemic lysozyme amyloidosis. In this paper we use highly accurate kinetic measurements of amyloid fibril growth in combination with spectroscopic tools to quantify the effect of modifications in solution conditions and in the amino acid sequence of human lysozyme on its propensity to form amyloid fibrils under acidic conditions. We elucidate and quantify the correlation between the rate of amyloid growth and the population of nonnative states, and we show that changes in amyloidogenicity are almost entirely due to alterations in the stability of the native state, while other regions of the global free-energy surface remain largely unmodified. These results provide insight into the complex dynamics of a macromolecule on a multidimensional energy landscape and point the way for a better understanding of amyloid diseases.

Tolman Effect On Fluid Dynamics In Carbon Nanotubes

The numerical solutions of or(R) given by two different methods (Samsonov et al., 2003; and Lu et al., 2005) are compared with the result that they are coincident closely (the difference is within 4%). We conclude that it is necessary to consider the Tolman correction in the calculation of fluid dynamics in carbon nanotubes. Although our conclusion is the same as that of Prylutskyy et al. (2005), the sign of our Tolman correction is opposite to theirs, and the difference can be attributed to the errors appeared in the paper of Prylutskyy et al.

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.

Theory Of Phase Transformation And Reorientation In Single Crystalline Shape Memory Alloys

A constitutive model, based on an (n + 1)-phase mixture of the Mori-Tanaka average theory, has been developed for stress-induced martensitic transformation and reorientation in single crystalline shape memory alloys. Volume fractions of different martensite lattice correspondence variants are chosen as internal variables to describe microstructural evolution. Macroscopic Gibbs free energy for the phase transformation is derived with thermodynamics principles and the ensemble average method of micro-mechanics. The critical condition and the evolution equation are proposed for both the phase transition and reorientation. This model can also simulate interior hysteresis loops during loading/unloading by switching the critical driving forces when an opposite transition takes place.

Surface-induced melting of metal nanoclusters

We investigate the size effect on melting of metal nanoclusters by molecular dynamics simulation and thermo dynamic theory based on Kofman's melt model. By the minimization of the free energy of metal nanoclusters with respect to the thickness of the surface liquid layer, it has been found that the nanoclusters of the same metal have the same premelting temperature T-pre = T-0 - T-0(gamma(su) - gamma(lv) - gamma(sl))/(rhoLxi) (T-0 is the melting point of bulk metal, gamma(sv) the solid-vapour interfacial free energy, gamma(sl) the liquid-vapour interfacial free energy, gamma(sl),l the solid-liquid interfacial free energy, p the density of metal, L the latent heat of bulk metal, and xi the characteristic length of surface-interface interaction) to be independent of the size of nanoclusters, so that the characteristic length of a metal can be obtained easily by T-pre, which can be obtained by experiments or molecular dynamics (MD) simulations. The premelting temperature T-pre of Cu is obtained by AID simulations, then xi is obtained. The melting point T-cm is further predicted by free energy analysis and is in good agreement with the result of our MD simulations. We also predict the maximum premelting-liquid width of Cu nanoclusters with various sizes and the critical size, below which there is no premelting.

Growth of liquid bridge in AFM

Capillary forces are dominant in adhesive forces measured with an atomic force microscope (AFM) in ambient air, which are thought to be dependent on water film thickness, relative humidity and the free energy of the water film. In this paper, besides these factors, we study the nature of the 'pull-off' force on a variety of atmospheres as a function of the contact time. It is found that capillary forces strongly depend on the contact time. In lower relative humidity atmosphere, the adhesion force is almost independent of the contact time. However, in higher relative humidity, the adhesion force increases with the contact time. Based on the experiment and a model that we present in this paper, the growth of the liquid bridge can be seen as undergoing two processes: one is water vapour condensation; the other is the motion of the thin liquid film that is absorbed on the substrate. The experiment and the growth model presented in this paper have direct relevance to the working mechanism of AFM in ambient air.

A new dynamic model for study of dislocation pattern formation

Based on the principle given in nonlinear diffusion-reaction dynamics, a new dynamic model for dislocation patterning is proposed by introducing a relaxation time to the relation between dislocation density and dislocation flux. The so-called chemical potential like quantities, which appear in the model can be derived from variation principle for free energy functional of dislocated media, where the free energy density function is expressed in terms of not only the dislocation density itself but also their spatial gradients. The Linear stability analysis on the governing equations of a simple dislocation density shows that there exists an intrinsic wave number leading to bifurcation of space structure of dislocation density. At the same time, the numerical results also demonstrate the coexistence and transition between different dislocation patterns.

Wave Propagation Analysis in a Pressure-Wave-Refrigerator

Pressure wave refrigerators (PWR) refrigerate the gas through periodical expansion waves. Due to its simple structure and robustness, PWR may have many potential applications if the efficiency becomes competitive with existing alternative devices. In order to improve the efficiency, the characteristics of wave propagation in a PWR are studied by experiment, numerical simulation and theoretical analysis. Based on the experimental results and numerical simulation, a simplified model is suggested, which includes the assumptions of flux-equilibrium and conservation of the free energy. This allows the independent analysis of the operation parameters and design specifics. Furthermore, the optimum operation condition can be deduced. Some considerations to improve the PWR efficiency are also given.