Baohua Ji

Materials become insensitive to flaws at nanoscale: Lessons from nature

Natural materials such as bone, tooth, and nacre are nanocomposites of proteins and minerals with superior strength. Why is the nanometer scale so important to such materials? Can we learn from this to produce superior nanomaterials in the laboratory? These questions motivate the present study where we show that the nanocomposites in nature exhibit a generic mechanical structure in which the nanometer size of mineral particles is selected to ensure optimum strength and maximum tolerance of flaws (robustness). We further show that the widely used engineering concept of stress concentration at flaws is no longer valid for nanomaterial design.

Nano to Micro Structural Hierarchy Is Crucial for Stable Superhydrophobic and Water-Repellent Surfaces

Water-repellent biological systems such as lotus leaves and water striders legs exhibit two-level hierarchical surface structures with the smallest characteristic size on the order of a few hundreds nanometers. Here we show that such nano to micro structural hierarchy is crucial for a superhydrophobic and water-repellent surface. The first level structure at the scale of a few hundred nanometers allows the surface to sustain the highest pressure found in the natural environment of plants and insects in order to maintain a stable Cassie state. The second level structure leads to dramatic reduction in contact area, hence minimizing adhesion between water and the solid surface. The two level hierarchy further stabilizes the superhydrophobic state by enlarging the energy difference between the Cassie and the Wenzel states. The stability of Cassie state at the nanostructural scale also allows the higher level structures to restore superhydrophobicity easily after the impact of a rainfall.

Cracking and adhesion at small scales: atomistic and continuum studies of flaw tolerant nanostructures

Once the characteristic size of materials reaches nanoscale, the mechanical properties may change drastically and classical mechanisms of materials failure may cease to hold. In this paper, we focus on joint atomistic-continuum studies of failure and deformation of nanoscale materials. In the first part of the paper, we discuss the size dependence of brittle fracture. We illustrate that if the characteristic dimension of a material is below a critical length scale that can be on the order of several nanometres, the classical Griffith theory of fracture no longer holds. An important consequence of this finding is that materials with nano-substructures may become flaw-tolerant, as the stress concentration at crack tips disappears and failure always occurs at the theoretical strength of materials, regardless of defects. Our atomistic simulations complement recent continuum analysis (Gao et al 2003 Proc. Natl Acad. Sci. USA 100 5597–600) and reveal a smooth transition between Griffith modes of failure via crack propagation to uniform bond rupture at theoretical strength below a nanometre critical length. Our results may have consequences for understanding failure of many small-scale materials. In the second part of this paper, we focus on the size dependence of adhesion systems. We demonstrate that optimal adhesion can be achieved by either length scale reduction, or by optimization of the shape of the surface of the adhesion element. We find that whereas change in shape can lead to optimal adhesion strength, those systems are not robust against small deviations from the optimal shape. In contrast, reducing the dimensions of the adhesion system results in robust adhesion devices that fail at their theoretical strength, regardless of the presence of flaws. An important consequence of this finding is that even under the presence of surface roughness, optimal adhesion is possible provided the size of contact elements is sufficiently small. Our atomistic results corroborate earlier theoretical modelling at the continuum scale (Gao and Yao 2004 Proc. Natl Acad. Sci. USA 101 7851–6). We discuss the relevance of our studies with respect to natures design of bone nanostructures and nanoscale adhesion elements in geckos.

On the strength of β-sheet crystallites of Bombyx mori silk fibroin

Silk fibroin, a natural multi-domain protein, has attracted great attention due to its superior mechanical properties such as ultra-high strength and stretchability, biocompatibility, as well as its versatile biodegradability and processability. It is mainly composed of β-sheet crystallites and amorphous domains. Although its strength is well known to be controlled by the dissociation of protein chains from β-sheet crystallites, the way that water as the solvent affects its strength and the reason that its theoretically predicted strength is several times higher than experimental measurement remain unclear. We perform all-atom molecular dynamics simulations on a β-sheet crystallite of Bombyx mori silk. We find that water solvent reduces the number and strength of hydrogen bonds between β-chains, and thus greatly weakens the strength of silk fibroin. By dissociating protein chains at different locations from the crystallite, we also find that the pulling strength for the interior chains is several times higher than that for the surface/corner chains, with the former being consistent with the theoretically predicted value, while the latter on par with the experimental value. It is shown that the weakest rupture strength controls the failure strength of silk fibre. Hence, this work sheds light on the role of water in the strength of silk fibroin and also provides clues on the origin of the strength difference between theory and experiment.

PERSPECTIVES IN MECHANICS OF HETEROGENEOUS SOLIDS

The Micro- and Nano-mechanics Working Group of the Chinese Society of Theoretical and Applied Mechanics organized a forum to discuss the perspectives, trends, and directions in mechanics of heterogeneous materials in January 2010. The international journal, Acta Mechanica Solida Sinica, is devoted to all fields of solid mechanics and relevant disciplines in science, technology, and engineering, with a balanced coverage on analytical, experimental, numerical and applied investigations. On the occasion of the 30th anniversary of Acta Mechanica Solida Sinica, its editor-in-chief, Professor Q.S. Zheng invited some of the forum participants to review the state-of-the-art of mechanics of heterogeneous solids, with a particular emphasis on the recent research development results of Chinese scientists. Their reviews are organized into five research areas as reported in different sections of this paper. §I firstly brings in focus on micro- and nano-mechanics, with regards to several selective topics, including multiscale coupled models and computational methods, nanocrystal superlattices, surface effects, micromechanical damage mechanics, and microstructural evolution of metals and shape memory alloys. §II shows discussions on multifield coupled mechanical phenomena, e.g., multi-fields actuations of liquid crystal polymer networks, mechanical behavior of materials under radiations, and micromechanics of heterogeneous materials. In § III, we mainly address the multiscale mechanics of biological nanocomposites, biological adhesive surface mechanics, wetting and dewetting phenomena on microstructured solid surfaces. The phononic crystals and manipulation of elastic waves were elaborated in § IV. Finally, we conclude with a series of perspectives on solid mechanics. This review will set a primary goal of future science research and engineering application on solid mechanics with the effort of social and economic development.

Stability of Adhesion Clusters and Cell Reorientation under Lateral Cyclic Tension

This work is motivated by experimental observations that cells on stretched substrate exhibit different responses to static and dynamic loads. A model of focal adhesion that can consider the mechanics of stress fiber, adhesion bonds, and substrate was developed at the molecular level by treating the focal adhesion as an adhesion cluster. The stability of the cluster under dynamic load was studied by applying cyclic external strain on the substrate. We show that a threshold value of external strain amplitude exists beyond which the adhesion cluster disrupts quickly. In addition, our results show that the adhesion cluster is prone to losing stability under high-frequency loading, because the receptors and ligands cannot get enough contact time to form bonds due to the high-speed deformation of the substrate. At the same time, the viscoelastic stress fiber becomes rigid at high frequency, which leads to significant deformation of the bonds. Furthermore, we find that the stiffness and relaxation time of stress fibers play important roles in the stability of the adhesion cluster. The essence of this work is to connect the dynamics of the adhesion bonds (molecular level) with the cells behavior during reorientation (cell level) through the mechanics of stress fiber. The predictions of the cluster model are consistent with experimental observations.

Structure-based design of carbon nanotubes as HIV-1 protease inhibitors: atomistic and coarse-grained simulations.

Nanoparticles such as fullerenes and carbon nanotubes have been extensively studied for biomedical applications. In this paper, we report the design of carbon nanotubes as HIV-1 protease inhibitors. Docking and molecular dynamics calculations are performed using an atomistic model to explore the optimal interaction structure and free energy between the nanotube and HIV-1 protease. A coarse-grained model is then developed based on the atomistic model, allowing us to investigate the dynamic behaviors of the protease in the bound and unbound states. The dynamic process reveals that the carbon nanotube is able to bind to the active site of the protease and prevent the active flaps from opening up, thus blocking the function of the protease. This process is strongly influenced by the size of the nanotube. The binding of carbon nanotubes to an alternative binding site other than the active site is also explored. Therefore, carbon nanotube-based inhibitors have great potential for application as HIV-1 protease inhibitors.

Stabilizing to disruptive transition of focal adhesion response to mechanical forces

Strong mechanical forces can, obviously, disrupt cell-cell and cell-matrix adhesions, e.g., cyclic uniaxial stretch induces instability of cell adhesion, which then causes the reorientation of cells away from the stretching direction. However, recent experiments also demonstrated the existence of force dependent adhesion growth (rather than dissociation). To provide a quantitative explanation for the two seemingly contradictory phenomena, a microscopic model that includes both integrin-integrin interaction and integrin-ligand interaction is developed at molecular level by treating the focal adhesion as an adhesion cluster. The integrin clustering dynamics and integrin-ligand binding dynamics are then simulated within one unified theoretical frame with Monte Carlo simulation. We find that the focal adhesion will grow when the traction force is higher than a relative small threshold value, and the growth is dominated by the reduction of local chemical potential energy by the traction force. In contrast, the focal adhesion will rupture when the traction force exceeds a second threshold value, and the rupture is dominated by the breaking of integrin-ligand bonds. Consistent with the experiments, these results suggest a force map for various responses of cell adhesion to different scales of mechanical force.

Strength of Hydrogen Bond Network Takes Crucial Roles in the Dissociation Process of Inhibitors from the HIV-1 Protease Binding Pocket

To understand the underlying mechanisms of significant differences in dissociation rate constant among different inhibitors for HIV-1 protease, we performed steered molecular dynamics (SMD) simulations to analyze the entire dissociation processes of inhibitors from the binding pocket of protease at atomistic details. We found that the strength of hydrogen bond network between inhibitor and the protease takes crucial roles in the dissociation process. We showed that the hydrogen bond network in the cyclic urea inhibitors AHA001/XK263 is less stable than that of the approved inhibitor ABT538 because of their large differences in the structures of the networks. In the cyclic urea inhibitor bound complex, the hydrogen bonds often distribute at the flap tips and the active site. In contrast, there are additional accessorial hydrogen bonds formed at the lateral sides of the flaps and the active site in the ABT538 bound complex, which take crucial roles in stabilizing the hydrogen bond network. In addition, the water molecule W301 also plays important roles in stabilizing the hydrogen bond network through its flexible movement by acting as a collision buffer and helping the rebinding of hydrogen bonds at the flap tips. Because of its high stability, the hydrogen bond network of ABT538 complex can work together with the hydrophobic clusters to resist the dissociation, resulting in much lower dissociation rate constant than those of cyclic urea inhibitor complexes. This study may provide useful guidelines for design of novel potent inhibitors with optimized interactions.

Coarse-grained molecular dynamics of ligands binding into protein: The case of HIV-1 protease inhibitors

Binding dynamics and pathways of ligands or inhibitors to target proteins are challenging both experimental and theoretical biologists. A dynamics understanding of inhibitors interacting with protein is essential for the design of novel potent drugs. In this work we applied a coarse-grained molecular dynamics method for simulating inhibitors entering the binding cavity of human immunodeficiency virus type 1 protease (PR). It shows that the coarse-grained dynamics, consistent with the experimental results, can capture the essential molecular dynamics of various inhibitors binding into PR. The primary driving force for the binding processes is the nonbond interaction between inhibitors and PR. The size and topology of inhibitors and the interacting strength between inhibitors and PR have great influence on the binding mode and processes. The interaction strength between the PR and various inhibitors is also analyzed by atomistic molecular mechanics and Poisson-Boltzmann solvation area method.