Simulated annealing of silicon atom clusters in Langevin molecular dynamics

Extraction of kinetics from equilibrium distributions of states using the Metropolis Monte Carlo method

The Metropolis Monte Carlo (MMC) method is used to extract reaction kinetics from a given equilibrium distribution of states of a complex system. The approach is illustrated by the folding/unfolding reaction for two proteins: a model β-hairpin and α-helical protein α_{3}D. For the β-hairpin, the free energy surfaces (FESs) and free energy profiles (FEPs) are employed as the equilibrium distributions of states, playing a role of the potentials of mean force to determine the acceptance probabilities of new states in the MMC simulations. Based on the FESs and PESs for a set of temperatures that were simulated with the molecular dynamics (MD) method, the MMC simulations are performed to extract folding/unfolding rates. It has been found that the rate constants and first-passage time (FPT) distributions obtained in the MMC simulations change with temperature in good agreement with those from the MD simulations. For α_{3}D protein, whose equilibrium folding/unfolding was studied with the single-molecule FRET method [Chung et al., J. Phys. Chem. A 115, 3642 (2011)1089-563910.1021/jp1009669], the FRET-efficiency histograms at different denaturant concentrations were used as the equilibrium distributions of protein states. It has been found that the rate constants for folding and unfolding obtained in the MMC simulations change with denaturant concentration in reasonable agreement with the constants that were extracted from the photon trajectories on the basis of theoretical models. The simulated FPT distributions are single-exponential, which is consistent with the assumption of two-state kinetics that was made in the theoretical models. The promising feature of the present approach is that it is based solely on the equilibrium distributions of states, without introducing any additional parameters to perform simulations, which suggests its applicability to other complex systems.

First-passage times in protein folding: exploring the native-like states vs. overcoming the free energy barrier

Using a model β-hairpin protein as a representative example of simple two-state folders whose kinetics are uncomplicated by the presence of on- and off-pathway intermediates, it is studied how the search for the protein's native state among native-like states affects the folding kinetics. It is revealed that the first-passage time (FPT) distributions are essentially single-exponential not only for the times to overcome the free energy barrier between the unfolded and native-like states but also for the times to find the native state among the native-like ones. The FPT distributions of this type are observed through all studied two-state-like regimes of protein folding, varying from a regime close to two-state folding to a regime close to downhill folding. If the protein explores native-like states for a time much longer than the time to overcome the free energy barrier, which is characteristic of high temperatures, the resulting FPT distribution to reach the native state remains close to exponential but the mean FPT (MFPT) is determined not by the height of the free energy barrier but by the time to explore native-like states. In particular, the mean time to overcome the free energy barrier is in reasonable agreement with the Kramers rate formula and generally far shorter than the overall MFPT to reach the native state. The observed increase of the overall MFPT, as a result of longer exploration of native-like states, may lead to an overestimate of the height of the free energy barrier between the unfolded and folded states when it is calculated from the overall MFPT.

Synthesis of multiple longitudinal polarization vortex structures and its application in sorting chiral nanoparticles

As the essential properties of organisms, detection and characterization of chirality are of supreme importance in physiology and pharmacology. In this work, we propose an optical technique to sort chiral materials by use of longitudinal polarization vortex (LPV) structures, which is generated with tightly focusing Pancharatnam-Berry tailored Laguerre-Gaussian beam. The nonparaxial propagation of the focusing field leads to the creation of multiple pairs of dual LPV structures with arbitrary topological charge and location, which can be independently controlled by the spatial phase modulation applied on the illumination. More importantly, the opposite spin angular momentums carried by each pair of dual foci lead to different energy flow directions, making it suitable to sort nanoparticles by their handedness. In addition, the LPV structures would also bring different dynamic behaviors to the enantiomers, providing a feasible route toward all-optical enantiopure chemical syntheses and enantiomer separations in pharmaceuticals.

Efficient hierarchical models for reactivity of organic layers on semiconductor surfaces

Computational modeling of organic interface formation on semiconductors poses a challenge to a density functional theory-based description due to structural and chemical complexity. A hierarchical approach is presented, where parts of the interface are successively removed in order to increase computational efficiency while maintaining the necessary accuracy. First, a benchmark is performed to probe the validity of this approach for three model reactions and five dispersion corrected density functionals. Reaction energies are generally well reproduced by generalized gradient approximation-type functionals but accurate reaction barriers require the use of hybrid functionals. Best performance is found for the model system that does not explicitly consider the substrate but includes its templating effects. Finally, this efficient model is used to provide coverage dependent reaction energies and suggest synthetic principles for the prevention of unwanted growth termination reactions for organic layers on semiconductor surfaces.

Electrostatics-Assisted Building-Up Procedure for Capturing Energy Minima of Metal Clusters: Test Case of Agn Clusters

Global geometry optimization of metal clusters is an important problem in nanophysics. The starting geometries of the clusters generated with empirical or other model potentials are generally optimized further by density functional theory (DFT)-based energy minimization. For this purpose, several algorithms such as simulated annealing, genetic algorithms, basin hopping, etc. are used. Our building-up procedure generates putative lower-energy structures of metal (M) clusters, M_n+1, M_n+2, etc., by anchoring one or more metal atoms in the vicinity of the minima of the molecular electrostatic potential (MESP) of M_n. Here, we report an application of this method to Ag_n clusters, for 5 ≤ n ≤ 20, followed up by DFT-based geometry optimization, generating several lower-energy structures than those reported in the literature. New low-energy isomers are obtained by applying the same procedure to the test case of mixed-metal clusters, Ni_nAg_m, for n + m = 4 and 5. In conclusion, our MESP-based building-up procedure offers a new general methodology for generating lower-energy geometries of metal clusters.

How the dyes affect folding of small proteins in single-molecule FRET experiments: A simulation study

A key question in the application of the single-molecule Förster resonance energy transfer (smFRET) technique to study protein folding is how the dyes affect the protein behavior. Understanding of these effects is particularly important for small proteins, for which the dyes, along with their linkers, can be comparable in size (mass) with the protein. Using a coarse-grained model, we simulated folding of BBL protein and two of its FRET constructs. The obtained results suggest that even for small proteins, such as the 45-residue BBL, the appearance of the excluded volume in the protein conformation space due to the presence of dyes does not change the overall picture of folding. At the same time, some deviations from folding of the original protein are observed, in particular, the FRET constructs fold considerably slower than the original protein because the protein collapse in the initial state of folding is slowed down due to the protein loading with relatively massive dyes.

Alternation of phases of regular and irregular dynamics in protein folding

The regularity of the dynamics in different phases of protein folding is investigated for a set of proteins which undergo a cooperative, two-state folding transition. To determine the degree of regularity of the dynamics, the fractal dimension of probability fluxes is calculated on the basis of simulated folding trajectories. It has been found that the phases of regular and irregular dynamics alternate as follows. In the initial (collapse) phase of folding, the dynamics are essentially regular. Then, as the protein comes to the basin of semicompact states that precedes the transition state, the dynamics become irregular. At the transition state, the dynamics are regularized again but become less regular when the nativelike states are explored. Depending on the specific conditions at which the protein folding was considered, some phases of the dynamics could not be well resolved, but no significant deviation from this general picture has been observed. The regularization of the dynamics at the transition state is discussed in relation to the recent studies of the Hamiltonian dynamics of small clusters, where both regular and chaotic dynamics were observed depending on the flatness of the energy surface at the transition state.

Modeling of Multicolor Single-Molecule Förster Resonance Energy-Transfer Experiments on Protein Folding

Using a coarse-grained, C_α-model of BBL protein, a multicolor single-molecule Förster resonance energy transfer (FRET) experiment is modeled. Three fluorophores are introduced, which, for simplicity, are associated with C_α beads. Two fluorophores are placed at the ends of protein chain and the third one at the middle of the chain. The free-energy surfaces (FESs) depending on the interfluorophore distances and on the FRET efficiencies corresponding to these distances have been constructed and compared with the FESs depending on the conventional collective variables, such as the fraction of native contacts and radius of gyration. It has been found that multicolor experiments can successfully resolve all essential BBL states that are revealed by the conventional FESs. The resolution of these states with the FRET-efficiency histogram is found to be successful if the energy transfer is measured between the fluorophores at the BBL ends. We also show that, although the present model construct of BBL is very simple, it captures some characteristic features of the single-molecule FRET experiments, such as the pattern of the FRET-efficiency histograms and their evolution with the denaturant concentration.

Expedite random structure searching using objects from Wyckoff positions

Random structure searching has been proved to be a powerful approach to search and find the global minimum and the metastable structures. A true random sampling is in principle needed yet it would be highly time-consuming and/or practically impossible to find the global minimum for the complicated systems in their high-dimensional configuration space. Thus the implementations of reasonable constraints, such as adopting system symmetries to reduce the independent dimension in structural space and/or imposing chemical information to reach and relax into low-energy regions, are the most essential issues in the approach. In this paper, we propose the concept of "object" which is either an atom or composed of a set of atoms (such as molecules or carbonates) carrying a symmetry defined by one of the Wyckoff positions of space group and through this process it allows the searching of global minimum for a complicated system to be confined in a greatly reduced structural space and becomes accessible in practice. We examined several representative materials, including Cd₃As₂ crystal, solid methanol, high-pressure carbonates (FeCO₃), and Si(111)-7 × 7 reconstructed surface, to demonstrate the power and the advantages of using "object" concept in random structure searching.

An electrodynamics-Langevin dynamics (ED-LD) approach to simulate metal nanoparticle interactions and motion

Understanding the formation of electrodynamically interacting assemblies of metal nanoparticles requires accurate computational methods for determining the forces and propagating trajectories. However, since computation of electromagnetic forces occurs on attosecond to femtosecond timescales, simulating the motion of colloidal nanoparticles on milliseconds to seconds timescales is a challenging multi-scale computational problem. Here, we present a computational technique for performing accurate simulations of laser-illuminated metal nanoparticles. In the simulation, we self-consistently combine the finite-difference time-domain method for electrodynamics (ED) with Langevin dynamics (LD) for the particle motions. We demonstrate the ED-LD method by calculating the 3D trajectories of a single 100-nm-diameter Ag nanoparticle and optical trapping and optical binding of two and three 150-nm-diameter Ag nanoparticles in simulated optical tweezers. We show that surface charge on the colloidal metal nanoparticles plays an important role in their optically driven self-organization. In fact, these simulations provide a more complete understanding of the assembly of different structures of two and three Ag nanoparticles that have been observed experimentally, demonstrating that the ED-LD method will be a very useful tool for understanding the self-organization of optical matter.

Protein folding as a complex reaction: a two-component potential for the driving force of folding and its variation with folding scenario

The Helmholtz decomposition of the vector field of probability fluxes in a two-dimensional space of collective variables makes it possible to introduce a potential for the driving force of protein folding [Chekmarev, J. Chem. Phys. 139 (2013) 145103]. The potential has two components: one component (Φ) is responsible for the source and sink of the folding flow, which represent, respectively, the unfolded and native state of the protein, and the other (Ψ) accounts for the flow vorticity inherently generated at the periphery of the flow field and provides the canalization of the flow between the source and sink. Both components obey Poisson’s equations with the corresponding source/sink terms. In the present paper, we consider how the shape of the potential changes depending on the scenario of protein folding. To mimic protein folding dynamics projected onto a two-dimensional space of collective variables, the two-dimensional Müller and Brown potential is employed. Three characteristic scenarios are considered: a single pathway from the unfolded to the native state without intermediates, two parallel pathways without intermediates, and a single pathway with an off-pathway intermediate. To determine the probability fluxes, the hydrodynamic description of the folding reaction is used, in which the first-passage folding is viewed as a steady flow of the representative points of the protein from the unfolded to the native state. We show that despite the possible complexity of the folding process, the Φ-component is simple and universal in shape. The Ψ-component is more complex and reveals characteristic features of the process of folding. The present approach is potentially applicable to other complex reactions, for which the transition from the reactant to the product can be described in a space of two (collective) variables.

Protein folding: complex potential for the driving force in a two-dimensional space of collective variables

Using the Helmholtz decomposition of the vector field of folding fluxes in a two-dimensional space of collective variables, a potential of the driving force for protein folding is introduced. The potential has two components. One component is responsible for the source and sink of the folding flows, which represent respectively, the unfolded states and the native state of the protein, and the other, which accounts for the flow vorticity inherently generated at the periphery of the flow field, is responsible for the canalization of the flow between the source and sink. The theoretical consideration is illustrated by calculations for a model β-hairpin protein.

Characterizing plastic depinning dynamics with the fluctuation theorem

We demonstrate that the fluctuation theorem can be used to characterize plastic flow phases in collectively interacting particle assemblies driven over quenched disorder when strong fluctuations and crackling noise with 1/f(α) character occur. By measuring the frequency of entropy-destroying trajectories and the diffusivity near the threshold for motion, we map out the different dynamic phases and demonstrate that the fluctuation theorem holds in the strongly fluctuating plastic flow regime which was previously shown to be chaotic. For different driving rates and disorder strength, we find that it is possible to define an effective temperature which decreases with increasing drive, as expected for this type of system. When the size of the pinning sites is large, we identify specific regimes where the fluctuation theorem holds only at long times due to an excess of negative entropy events that occur when particles undergo circular motions within the traps. We discuss how the fluctuation theorem could be applied to plastic flow in other driven nonthermal systems with quenched disorder such as superconducting vortices, magnetic domain walls, Coulomb glasses, and earthquake models.

Atomic and Electronic Structure of Germanium Clusters at Finite Temperature Using Finite Difference Methods

Determining the electronic and structural properties of semiconductor clusters is one of the outstanding problems in materials science. The existence of numerous structures with nearly identical energies makes it very difficult to determine a realistic ground state structure. Moreover, even if an effective procedure can be devised to predict the ground state structure, questions can arise about the relevancy of the structure at finite temperatures. Kinetic effects and non-equilibrium structures may dominate the structural configurations present in clusters created under laboratory conditions. We illustrate theoretical techniques for predicting the structure and electronic properties of small germanium clusters. Spefically, we illustate that the detailed agreement between theoretical and experimental features can be exploited to identify the relevant isomers present under experimental conditions.

Ab Initio simulations of nonstoichiometric CdxTe1−x liquids

We present ab initio molecular-dynamics simulations for Cd(x)Te(1-x) liquids where the composition is nonstoichiometric. The simulations are performed following Born-Oppenheimer molecular dynamics. The required forces are obtained from a solution of the Kohn-Sham equation using ab initio pseudopotentials. We consider stoichiometries of the form: Cd(x)Te(1-x), where x=0.2, 0.4, 0.6, and 0.8. For each composition of the melt, we consider a range of temperatures near the experimentally determined liquid temperatures. We examine the microstructural properties of the melt, the viscosity, and self-diffusion properties of the liquid as a function of the stoichiometry and temperature. We also perform an analysis of the distribution of the electronic density of states in these liquids. We find that structural changes in the local order, experimentally predicted to occur when the concentration of Cd is increased, are closely related to changes in the electronic properties of the melt.

Taboo search by successive confinement: surveying a potential energy surface

A taboo search for minima on a potential energy surface (PES) is performed by means of confinement molecular dynamics: the molecular dynamics trajectory of the system is successively confined to various basins on the PES that have not been sampled yet. The approach is illustrated for a 13-atom Lennard-Jones cluster. It is shown that the taboo search radically accelerates the process of surveying the PES, with the probability of finding a new minimum defined by a propagating Fermi-like distribution.

Evidence of a reentrant peierls distortion in liquid GeTe

The local atomic order of semiconducting liquid GeTe is studied using first-principles molecular-dynamics simulations. Our work points out a high degree of alternating chemical order in the liquid and demonstrates the presence of a Peierls distortion close above the melting temperature. This distortion, absent in the high temperature crystalline structure of NaCl type, is a remnant of the atomic arrangement in the A7 low temperature crystalline phase. It disappears slowly with temperature, as the liquid evolves from a semiconducting to a metallic state.

COMPUTATIONAL STUDIES OF CLUSTERS:Methods and Results

▪ Abstract  The experimental and computational study of clusters has been an active field of research for over a decade. This review provides an overview of some of the methods that have been developed to study clusters and some of the results that have been obtained. Included are computational approaches to explore the potential energy surface for clusters, methods to extract thermodynamic properties from the potential surface information and approaches to insure simulation studies are performed in an ergodic fashion. The methods have proved to be useful in studying the structural transition from clusters to bulk phases, phase changes in small clusters and the importance of quantum effects. The review ends with a discussion of problems in cluster chemistry and physics that are of interest for future investigations.