Irreversible deformation of isotactic polypropylene in the pre-yield regime

In the modeling of the mechanical response of a polymer over a large strain range, the nonlinear viscoelastic and viscoplastic behavior must be considered. For many polymers, nonlinear behavior is observed at low loads, e.g. by a stress-dependence of the creep compliance for stresses above 2 MPa in case of the polypropylene used in this study. Additionally, plastic deformation has been observed at strains below the yield point for several polymers. In this study, the irreversible deformation by cavitation and shear yielding of polypropylene are characterized in the pre-yield regime in uniaxial tensile tests using digital image correlation. The recovery of strain after unloading at a prescribed strain level is measured and used to identify the evolution of the plastic strain during uniaxial tension. An experimental technique for simultaneous determination of the true stress–true strain curve and the degree of stress whitening, which relates to the amount of cavitation, is introduced and the initiation of cavitation is compared to the plastic deformation detected in strain recovery at various temperatures.     

Global optimization of binary Lennard-Jones clusters using three perturbation operators

Global optimization of binary Lennard-Jones clusters is a challenging problem in computational chemistry. The difficulty lies in not only that there are enormous local minima on the potential energy surface but also that we must determine both the coordinate position and the atom type for each atom and thus have to deal with both continuous and combinatorial optimization. This paper presents a heuristic algorithm (denoted by 3OP) which makes extensive use of three perturbation operators. With these operators, the proposed 3OP algorithm can efficiently move from a poor local minimum to another better local minimum and detect the global minimum through a sequence of local minima with decreasing energy. The proposed 3OP algorithm has been evaluated on a set of 96 × 6 instances with up to 100 atoms. We have found most putative global minima listed in the Cambridge Cluster Database as well as discovering 12 new global minima missed in previous research.     

Potential energy and free energy landscapes

Familiar concepts for small molecules may require reinterpretation for larger systems. For example, rearrangements between geometrical isomers are usually considered in terms of transitions between the corresponding local minima on the underlying potential energy surface, V. However, transitions between bulk phases such as solid and liquid, or between the denatured and native states of a protein, are normally addressed in terms of free energy minima. To reestablish a connection with the potential energy surface we must think in terms of representative samples of local minima of V, from which a free energy surface is projected by averaging over most of the coordinates. The present contribution outlines how this connection can be developed into a tool for quantitative calculations. In particular, stepping between the local minima of V provides powerful methods for locating the global potential energy minimum, and for calculating global thermodynamic properties. When the transition states that link local minima are also sampled we can exploit statistical rate theory to obtain insight into global dynamics and rare events. Visualizing the potential energy landscape helps to explain how the network of local minima and transition states determines properties such as heat capacity features, which signify transitions between free energy minima. The organization of the landscape also reveals how certain systems can reliably locate particular structures on the experimental time scale from among an exponentially large number of local minima. Such directed searches not only enable proteins to overcome Levinthal's paradox but may also underlie the formation of "magic numbers" in molecular beams, the self-assembly of macromolecular structures, and crystallization.     

Simulated annealing to locate various stationary points in semiempirical methods

In this work, an algorithm was developed to study the potential energy surfaces in the coordinate spaces of molecules by a nonlocal way, in contrast to classic energy minimizers as the BFGS or the DFP method. This algorithm, based on the specificities of semiempirical methods, mixes simulated annealing and local searches to reduce computation costs. By this technique, the global energy minimum can be localized. Moreover, local minima that are close in energy to the global minimum are also obtained. If the search is not only for minima but for all stationary points (minima, saddle points…), then the energy is replaced by the gradient norm, which reaches its minimum values at stationary points. The annealing process is stopped before having accurately reached the global minimum and generates a list of geometries whose energies (respectively, whose gradients) are optimized by local minimizers. This list of geometries is shortened from the nearly equivalent geometries by a dynamic single-clustering analysis. The energy/gradient local minimizers act on the clustered list to produce a set of minima/stationary points. A targeted search of these points and reduction of the costs are reached by the way of several penalty functions. They eliminate—without energy calculation—most of the points generated by random walks on the potential energy surface. These penalty functions (on the total moment of inertia or on interatomic distances) are specific to the class of problem studied. They account for the nonrupture of either specific chemical bonds or rings in cyclic molecules, they assure that molecular systems are kept bonded, and they avoid the collapsing of atoms. © 1992 John Wiley & Sons, Inc.     

A strategy to find minimal energy nanocluster structures

An unbiased strategy to search for the global and local minimal energy structures of free standing nanoclusters is presented. Our objectives are twofold: to find a diverse set of low lying local minima, as well as the global minimum. To do so, we use massively the fast inertial relaxation engine algorithm as an efficient local minimizer. This procedure turns out to be quite efficient to reach the global minimum, and also most of the local minima. We test the method with the Lennard–Jones (LJ) potential, for which an abundant literature does exist, and obtain novel results, which include a new local minimum for LJ₁₃, 10 new local minima for LJ₁₄, and thousands of new local minima for . Insights on how to choose the initial configurations, analyzing the effectiveness of the method in reaching low‐energy structures, including the global minimum, are developed as a function of the number of atoms of the cluster. Also, a novel characterization of the potential energy surface, analyzing properties of the local minima basins, is provided. The procedure constitutes a promising tool to generate a diverse set of cluster conformations, both two‐ and three‐dimensional, that can be used as an input for refinement by means of ab initio methods. © 2013 Wiley Periodicals, Inc.     

Inherent structures of crystalline tetracene

We have systematically sampled the potential energy surface of crystalline tetracene to identify its local minima. These minima represent all possible stable configurations and constitute the "inherent structures" of the system. The crystal is described in terms of rigid molecules with Coulombic and atom-atom interactions. Hundreds of distinct minima are identified, mostly belonging to the space groups P (triclinic) and P2(1)/c (monoclinic), with a variety of structural arrangements. The deepest minimum corresponds to the high temperature-low pressure polymorph. This is the only polymorph with a completely described X-ray structure, which is satisfactorily described by the calculations. The next deep minimum is likely to correspond to the low temperature-high pressure polymorph, which has been experimentally identified but not yet fully described.     

Theoretical model based on the memory effect for the strange photoisomerization kinetics of diarylethene derivatives dispersed on polymer films

In the present paper the authors present a theoretical model to explain the kinetics involving the induction period observed by Irie et al. [Nature (London) 420, 759 (2002)] for photoisomerization of diarylethene derivatives dispersed on polymer films at a single molecular level. In the model we assume that both ground state and excited state free energy landscapes which result from the interaction between the photochromic molecule and the surrounding polymer are rugged and have several local minima along the pathway to the critical point at which isomerization actually occurs. We assume that after one photoexcitation a fraction of the photochromic molecule moves to a new local minimum and stays there, although the other fraction returns to the original local minimum. The former effect is referred to as the memory effect. After repeated photoexcitations the photochromic molecule moves gradually from one local minimum to another in the pathway to the isomerization point. It finally reaches the isomerization point, where isomerization occurs. Their model successfully reproduces the kinetics of photoisomerization of diarylethene derivatives dispersed on polymer films observed experimentally.     

Global energy minimization by rotational energy embedding

Given a sufficiently good empirical potential function for the internal energy of molecules, prediction of the preferred conformations is nearly impossible for large molecules because of the enormous number of local energy minima. Energy embedding has been a promising method for locating extremely good local minima, if not always the global minimum. The algorithm starts by locating a very good local minimum when the molecule is in a high-dimensional Euclidean space, and then it gradually projects down to three dimensions while allowing the molecule to relax its energy throughout the process. Now we present a variation on the method, called rotational energy embedding, where the descent into three dimensions is carried out by a sequence of internal rotations that are the multidimensional generalization of varying torsion angles in three dimensions. The new method avoids certain kinds of difficulties experienced by ordinary energy embedding and enables us to locate conformations very near the native for avian pancreatic polypeptide and apamin, given only their amino acid sequences and a suitable potential function.     

Atomic mobility in strained glassy polymers: The role of fold catastrophes on the potential energy surface

Deformation is known to enhance the atomic mobility in disordered systems such as polymer materials. To elucidate the origin of this phenomenon, we carry out two types of simulations: molecular dynamics (MD) simulations, which determine the atomic trajectories at finite temperature, and quasi-static simulations, which determine the atomic trajectories in the limit of zero temperature (and in the limit of zero shear rate). The quasi-static simulations show discontinuous changes in properties, such as system energy and atomic mobility. We use a normal mode analysis to show that these discontinuous changes arise from fold catastrophes of the potential energy landscape, in which energy minima flatten out and the heights of energy barriers decrease to zero; this was demonstrated by normal mode frequencies following a power law with an exponent of 0.5 as the discontinuous change is approached. After the fold catastrophe, the system relaxes to a different energy minimum, giving rise to atomic displacements. These fold catastrophes are the only mechanism for diffusive atomic displacements in the quasi-static simulations, where there is no thermal energy. We compared the mean-squared displacements as a function of strain from the quasi-static simulations to those from MD simulations (which do include thermal effects)—the similarity of the values of the mean-squared displacements in these two types of simulations demonstrates that the fold catastrophes underlie the enhanced dynamics in strained polymer systems even at finite temperature. © 2012 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys, 2012     

Potential-hypersurface local minima and temperature

Summary The presence of several local energy minima on a potential hypersurface is treated in terms of geometry, energy, and harmonic vibrations. Partition functions of the minima are employed in order to treat temperature excitations of rotational-vibrational motions. Proportions of relative stabilities of the individual structures change with temperature (including interchanges of the relative stabilities so that the global energy minimum can even be less populated than a higher local energy minimum). Illustrative examples are given on B₂H₄ and Ga₂H₄ systems. The treatment is suggested as a standard complement of the local-minimum hypersurface representation (before whole potential hypersurfaces are constructed and employed in molecular dynamics treatments).Dedicated to Prof. Klaus Ruedenberg on the occasion of his 70th birthday     

Simulation of the plastic behavior of amorphous glassy bis-phenol-A-polycarbonate

A protocol for studying the plastic deformation of amorphous glassy polymers is presented. The protocol is based on a viable computational procedure which combines constant-stress molecular dynamics simulations and fixed-cell energy minimizations, followed by kinetic, configurational, and energy analyses. It is shown that the computational results can be accounted for within a "potential energy landscape" theoretical framework, in which the plastic transition is interpreted as a crossing between and a collapse onto each other of "ideal (thermodynamic) structures." The procedure is applied to bis-phenol-A-polycarbonate (BPA-PC), but is equally valid for a wide variety of polymeric species. Allowing for the limited size of the simulation cell, the high strain rate, and the fact that the simulation are conducted at low temperature, the values of the density, Young's modulus, yield strain, yield stress, activation energy, and activation volume are in fair agreement with the experimental data on BPA-PC. The analysis of the results shows that the plastic relaxation for this polymer has both a collective and cooperative character (as in classical percolation theories), involves a significant fraction of the simulation cell, and can be viewed as a "nanoscopic shear band."     

Surveying a potential energy surface by eigenvector-following

We have developed a method to search potential energy surfaces which avoids some of the difficulties associated with trapping in local minima. Steps are directly taken between minima using eigenvector-following. Exploration of this space by low temperature Metropolis Monte Carlo is a useful global optimisation tool. This method successfully finds the lowest energy icosahedral minima of Lennard- Jones clusters from random starting configurations, but cannot find the global minimum in a reasonable time for difficult cases such as the 38-atom Lennard-Jones cluster where the face-centred-cubic truncated octahedron is lowest in energy. However, by performing searches at higher temperatures, we have found a pathway between the truncated octahedron and the lowest energy icosahedral minima. Such a pathway may be illustrative of some of the structural transformations that are observed for supported metal clusters by electron microscopy.     

Atomic origin of hysteresis during cyclic loading of Si due to bond rearrangements at the crack surfaces

The atomistic origin of fatigue failure in micron-sized silicon devices is not fully understood. Two series of density-functional theory calculations on cubic diamond Si explore the effect of surface bond formation on crack healing in systems which exhibit strong surface reconstruction. Both series introduce a separation between Si(100) layers (i.e., the crack) and allow the ions to relax to their minimum-energy configuration. The initial surface ionic positions are either bulk terminated or 2 x 1 reconstructed. A plot of the energy versus the introduced separation reveals that once the surfaces reconstruct, the crack is no longer able to return to the equilibrium configuration. Rather, the healed crack interface contains defects which places the flawed energy minimum at a finite strain of 3% and an increased energy of 1.13 Jm2 relative to the equilibrium configuration. The irreversible plastic deformation supports the mechanism proposed by Kahn et al. [Science 298 1215 (2002)] that invokes mechanically induced subcritical cracking to explain the delayed onset of failure.     

A Comparison of radiationless transitions involving single-minimum and double minima potential energy surfaces

The quantum-statistical-mechanical (QSM) approach to molecular relaxation phenomena is employed to compare radiationless transitions originating from an electronic state characterized by a single minimum and double minima potential surface for a vibronically active, non-totally symmetric mode. The vibronic level dependence of the decay rates of these two cases has been investigated for both small and large energy gap transitions. It is shown that the behavior of a molecular system is quite different for an initial state possesing a double minima potential surface as compared to the case in which the initial state possesses a single minimum.     

On the investigation of the low-energy conformational space of molecules

A new perspective on traditional energy minimization problems is provided by a connection between statistical thermodynamics and combinatorial optimization (finding the minimum of a function depending on many variables). The joint use of a new method for uncovering the global minimum of intramolecular potential energy functions, based on following the asymptotic behavior of a system of stochastic differential equations, and an iterative-improvement technique, whereby a search for relative minima is made by carrying out local quasi-Newton minimizations starting from many distinct points of the energy hypersurface, proved most effective for investigating the low-energy conformational space of molecules.     

Force field development for poly(dimethylsilylenemethylene) with the aid of ab initio calculations

The molecular geometries, conformational energies, and zero-point energies of di(trimethylsilylene)methylene have been determined from high-level quantum chemistry calculations. The results are further used in the parametrization of a classical potential energy function suitable for performing simulations of the corresponding polymer, namely, poly(dimethylsilylenemethylene). Di(trimethylsilylene)methylene geometrical parameter optimizations for a proper location of the global minimum and other local minima, constrained at certain dihedral and bond angles, were performed at both the B3LYP/6-311G and MP2(full)/6-311G levels of theory. The global minimum configuration is slightly displaced from a perfectly staggered geometry, approximately by 16.0 degrees, at both levels of theory. Molecular mechanics and Monte Carlo calculations for isolated polymer chains together with molecular dynamics runs for the modeled dimer provide very good results in terms of conformational and thermodynamic properties.     

Global potential energy minima of C60(H2O)n clusters

Likely candidates for the global potential energy minima of C60(H2O)n clusters with n < or = 21 are found using basin-hopping global optimization. The potential energy surfaces are constructed using the TIP4P intermolecular potential for the water molecules, a Lennard-Jones water-fullerene potential, and a water-fullerene polarization potential, which depends on the first few nonvanishing C60 multipole polarizabilities. This combination produces a rather hydrophobic water-fullerene interaction. As a consequence, the water component of the lowest C60(H2O)n minima is quite closely related to low-lying minima of the corresponding TIP4P (H2O)n clusters. In most cases, the geometrical substructure of the water molecules in the C60(H2O)n global minimum coincides with that of the corresponding free water cluster. Exceptions occur when the interaction with C60 induces a change in geometry. This qualitative picture does not change significantly if we use the TIP3P model for the water-water interaction. Structures such as C60@(H2O)60, in which the water molecules surround the C60 fullerene, correspond to local minima with much higher potential energies. For such a structure to become the global minimum, the magnitude of the water-fullerene interaction must be increased to an unphysical value.     

Plasticity in polymeric honeycombs made by photo-polymerization and nozzle based 3D-printing

Study of the plastic deformation in polymeric honeycombs can pave the way for understanding the deformation localization in more complex cellular structures, which have received progressive attention in the past few years. This study compares the strain localization in deforming honeycombs made by two cost-effective 3D-printing technologies. Hexagonal honeycombs and their unit cell models were 3D-printed by both PolyJet™, using a photo-crosslinkable polymer, and fused deposition modeling (FDM) using a thermoplastic material. The state of the art digital image correlation (DIC) technique was employed as the experimental route in order to calculate the strain field during the deformation of manufactured parts. It was found that DIC is an effective tool to study the localization in 3D-printed honeycomb struts. Moreover, in comparison with FDM, PolyJet technology provides more homogeneous strain distributions in struts. In addition, FDM decreases the maximal strains generated on the side layers of the honeycomb struts. Accordingly, the ligament damage under plastic deformation can be postponed and the energy absorption capability of the product can be improved when PolyJet technology is utilized.     

Athermal simulation of plastic deformation in amorphous solids at constant pressure

An algorithm is introduced for the molecular simulation of constant-pressure plastic deformation in amorphous solids at zero temperature. This allows to directly study the volume changes associated with plastic deformation (dilatancy) in glassy solids. In particular, the dilatancy of polymer glasses is an important aspect of their mechanical behavior. The new method is closely related to Berendsen's barostat, which is widely used for molecular dynamics simulations at constant pressure. The new algorithm is applied to plane strain compression of a binary Lennard-Jones glass. Conditions of constant volume lead to an increase of pressure with strain, and to a concommitant increase in shear stress. At constant (zero) pressure, by contrast, the shear stress remains constant up to the largest strains investigated (ε = 1), while the system density decreases linearly with strain. The linearity of this decrease suggests that each elementary shear relaxation event brings about an increase in volume which is proportional to the amount of shear. In contrast to the stress–strain behavior, the strain-induced structural relaxation, as measured by the self-part of the intermediate structure factor, was found to be the same in both cases. This suggests that the energy barriers that must be overcome for their nucleation continually grow in the case of constant-volume deformation, but remain the same if the deformation is carried out at constant pressure. © 2004 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 42: 2057–2065, 2004     

Misorientation mapping for visualization of plastic deformation via electron back-scattered diffraction.

The ability to map plastic deformation around high strain gradient microstructural features is central in studying phenomena such as fatigue and stress corrosion cracking. A method for the visualization of plastic deformation in electron back-scattered diffraction (EBSD) data has been developed and is described in this article. This technique is based on mapping the intragrain misorientation in polycrystalline metals. The algorithm maps the scalar misorientation between a local minimum misorientation reference pixel and every other pixel within an individual grain. A map around the corner of a Vickers indentation in 304 stainless steel was used as a test case. Several algorithms for EBSD mapping were then applied to the deformation distributions around air fatigue and stress corrosion cracks in 304 stainless steel. Using this technique, clear visualization of a deformation zone around high strain gradient microstructural features (crack tips, indentations, etc.) is possible with standard EBSD data.