First-principle molecular dynamics calculation of selenium clusters

The equiliblium structures of small selenium clusters are obtained via first-principle molecular dynamics calculations based on the linearized-augmented-plane-wave (LAPW) method. Resulting equiliblium structures show a good agreement with experimental data and other firstprinciple calculations.     

The dynamics of phase changes in molecular clusters

The significant advantages offered by systems of molecular clusters in the study of homogeneous nucleation are discussed. Determinations of nucleation rates in clusters can be followed experimentally in supersonic jets or computationally in molecular dynamics simulations. Extraordinarily high rates may be encountered, both in freezing and in solid-state transitions. From such information can be inferred the interfacial free energies, σ_sl and σ_SS, mechanisms of solid-state transitions. and an explanation of why certain crystalline phases not found in bulk systems can be seen in large molecular clusters.     

Ligand binding inside the cavities of lacunar and saddle-shaped cyclidene complexes: molecular mechanics and molecular dynamics studies

Cobalt(II) complexes with tetradentate macrocyclic cyclidene ligands are known to coordinate one additional axial base molecule, leaving the sixth vacant coordination site at the metal ion available for small ligand (e.g., O2) binding. Molecular mechanics and molecular dynamics simulations provide a microscopic view of 1-methylimidazole (MeIm) binding within the cavities of several lacunar (bridged) and saddle-shaped (unbridged) cyclidenes and uncover the roles of the bridges and the walls of the clefts in steric protection of the cobalt(II) coordination site. Short bridges (C3 and C6) prevent inside-the-cavity MeIm binding because of severe ligand distortions leading to high-energy penalties (58 and 25 kcal/mol, respectively), while long bridges (C8 and C12) flip away from the MeIm binding site, allowing for penalty-free MeIm inclusion. In the unbridged saddle-shaped complex, there is no energy difference between inside- and outside-the-cavity MeIm binding. The preferential existence of the coordinatively unsaturated, five-coordinate species Co(unbrCyc)(MeIm)2+ should therefore be explained by electronic, rather than steric, factors. Molecular dynamics and free energy simulations reveal the presence of a weak (ca. 4 kcal/mol in the gas phase and ca. 2 kcal/mol in methanol solution) noncovalent MeIm binding site at the entrance of the cleft of cobalt(II) unbridged cyclidene, at a distance of about 4 A from the metal ion. The macrocycle geometry remains undistorted at such large Co-N(MeIm) separations, while the cavity opens up by 0.9 A upon covalent MeIm binding (Co-N(MeIm) distance of 2 A). An increase in macrocycle strain energy upon MeIm inclusion is compensated by favorable nonbonded interactions between the incoming base and the walls of the unbridged cyclidene.     

Far-infrared absorption of water clusters by first-principles molecular dynamics

Based on first-principle molecular dynamic simulations, we calculate the far-infrared spectra of small water clusters (H(2)O)(n) (n = 2, 4, 6) at frequencies below 1000 cm(-1) and at 80 K and at atmospheric temperature (T>200 K). We find that cluster size and temperature affect the spectra significantly. The effect of the cluster size is similar to the one reported for confined water. Temperature changes not only the shape of the spectra but also the total strength of the absorption, a consequence of the complete anharmonic nature of the classical dynamics at high temperature. In particular, we find that in the frequency region up to 320 cm(-1), the absorption strength per molecule of the water dimer at 220 K is significantly larger than that of bulk liquid water, while tetramer and hexamer show bulklike strengths. However, the absorption strength of the dimer throughout the far-infrared region is too small to explain the measured vapor absorption continuum, which must therefore be dominated by other mechanisms.     

Polar solvation dynamics of lysozyme from molecular dynamics studies

The solvation dynamics of a protein are believed to be sensitive to its secondary structures. We have explored such sensitivity in this article by performing room temperature molecular dynamics simulation of an aqueous solution of lysozyme. Nonuniform long-time relaxation patterns of the solvation time correlation function for different segments of the protein have been observed. It is found that relatively slower long-time solvation components of the α-helices and β-sheets of the protein are correlated with lower exposure of their polar probe residues to bulk solvent and hence stronger interactions with the dynamically restricted surface water molecules. These findings can be verified by appropriate experimental studies.     

Improving the performance of molecular dynamics simulations on parallel clusters

In this article a procedure is derived to obtain a performance gain for molecular dynamics (MD) simulations on existing parallel clusters. Parallel clusters use a wide array of interconnection technologies to connect multiple processors together, often at different speeds, such as multiple processor computers and networking. It is demonstrated how to configure existing programs for MD simulations to efficiently handle collective communication on parallel clusters with processor interconnections of different speeds.     

Tight‐binding molecular dynamics simulation of Si?H bond dissociation in silicon clusters

New model of Si? H bond dissociation is proposed and tested in the cluster Si₁₀H₁₆ by the simulation approach that combines classic molecular dynamics method and the self‐consistent tight‐binding electronic and total energy calculation one. It is shown that the monohydride Si? H bond is unstable with respect to silicon dangling bond and bend‐bridge Si? H? Si bond formation when this cluster traps the single positive charge and that hydrogen migrates through a path involving rather rotation around the Si? Si bond than the center of this bond (the bond‐centered position). These results can be useful for understanding hydrogen‐related phenomena at surfaces, interfaces, and internal voids of various hydrogenated silicon systems: electronic devices, silicon solar cells, and nanocrystalline and porous silicon. © 2003 Wiley Periodicals, Inc. Int J Quantum Chem 93: 351–359, 2003     

Ab initio model potential and molecular dynamics simulation of Ag6 clusters

We present an approach to generate a model potential with parameters fitted to ab initio energetic surfaces. The potential includes two-, three- and four-body terms. Each of them consists of an exponential exchange and dispersion terms. The analytical form of the latter was taken from perturbation theory up to fourth order. We illustrate the present approach by constructing an ab initio model potential for the Ag₆ cluster. A molecular dynamics simulation using this potential reveals interesting features in the isomerization of the C _2v structure. A two step isomerization transition is observed: First, at temperatures around 350 K, the cluster structures fluctuate between two-dimensional isomers. At higher temperatures (450 K), fast transitions occur between two- and three-dimensional cluster configurations.     

Ab initio molecular dynamics study of methanol adsorption on copper clusters

The preferential structures of small copper clusters Cun (n=2-9) and the adsorption of methanol molecules on these clusters are examined with first principles, molecular dynamics simulations. The results show that the copper clusters undergo systematic changes in bond length and bond order associated with altering their preferential structures from one-dimensional structures, to two-dimensional and three-dimensional structures. The results also indicate that low coordination number sites on the copper clusters are both the most favorable for methanol adsorption and have the greatest localization of electronic charge. The simulations predict that charge transfer between the neutral copper clusters and the incident methanol molecules is a key process by which adsorption is stabilized. Importantly, the changes in the dimensionality of the copper clusters do not significantly influence methanol adsorption.     

Simulation of silicon clusters from “quantum” Langevin molecular dynamics

We present a simulation method to determine from first principles the structure of low symmetry atomic systems. Our method is based on Langevin molecular dynamics and quantum mechanical interactions derived fromab initio pseudopotential calculations. The molecular dynamics time step with this approach can be one to two orders of magnitude larger than in the Car-Parrinello method, compensating for the time required for self-consistency at each step. Moreover, because the simulation is constrained to reside on the Born-Oppenheimer surface, this method can be used for insulating as well as metallic and charged systems. Application will be made to small silicon clusters.     

Pt-Au合金团簇热力学性质分子动力学模拟研究

采用分子动力学方法和原子嵌入法模型势模拟了Pt原子和Au原子合金纳米团簇的熔化过程,研究了这些金属原子纳米团簇熔点与团簇组分的关系,发现不同组分纳米团簇的熔点不是单调变化的,同时均出现了负热容现象.通过对各种团簇溶化前后结构的比较研究,分析了导致这种现象的原因.     

Electron hydration dynamics in water clusters: A direct ab initio molecular dynamics approach

Electron attachment dynamics of excess electron in water cluster (H2O)n (n = 2 and 3) have been investigated by means of full-dimensional direct ab initio molecular dynamics (MD) method at the MP26-311++G(d,p) level. It was found that the hydrogen bond breaking due to the excess electron is an important process in the first stage of electron capture in water trimer. Time scale of electron localization and hydrogen bond breaking were determined by the direct ab initio MD simulation. The initial process of hydration in water cluster is clearly visualized in the present study. In n = 3, an excess electron is first trapped around the cyclic water trimer with a triangular form, where the excess electron is equivalently distributed on the three water molecules at time zero. After 50 fs, the excess electron is concentrated into two water molecules, while the potential energy of the system decreases by -1.5 kcal/mol from the vertical point. After 100 fs, the excess electron is localized in one of the water molecules and the potential energy decreases by -5.3 kcal/mol, but the triangular form still remained. After that, one of the hydrogen bonds in the triangular form is gradually broken by the excess electron, while the structure becomes linear at 100-300 fs after electron capture. The time scale of hydrogen bond breaking due to the excess electron is calculated to be about 300 fs. Finally, a dipole bound state is formed by the linear form of three water molecules. In the case of n = 2, the dipole bound anion is formed directly. The mechanism of electron hydration dynamics was discussed on the basis of theoretical results.     

Reactivity of molecular oxygen with aluminum clusters: Density functional and Ab Initio molecular dynamics simulation study

Dissociative adsorption of molecular oxygen (O₂) on aluminum (Al) clusters has attracted much interest in the field of surface science and catalysis, but theoretical predictions of the reactivity of this reaction in terms of barrier height is still challenging. In this regard, we systematically investigate the reactivity of O₂ with Al clusters using density functional theory (DFT) and atom‐centered density matrix propagation (ADMP) simulations. We also calculate potential energy surfaces (PESs) of the reaction between O₂ and Al clusters to estimate the barrier energy of this reaction. The M06‐2X functional gives the barrier energy in agreement with the one calculated by coupled cluster singles and doubles with perturbed triples (CCSD(T)) while the TPSSh functional significantly underestimates the barrier height. The ADMP simulation using the M06‐2X functional predicts the reactivity of O₂ with the Al cluster in agreement with the experimental findings, that is, singlet O₂ readily reacts with Al clusters but triplet O₂ is less reactive. We found that the ability of a DFT functional to describe the charge transfer appropriately is critical for calculating the barrier energy and the reactivity of the reaction of O₂ with Al clusters. The M06‐2X functional is relevant for investigating chemical reactions involving Al and O₂. © 2016 Wiley Periodicals, Inc.     

Electron and nuclear dynamics of molecular clusters in ultraintense laser fields. III. Coulomb explosion of deuterium clusters

In this paper we present a theoretical and computational study of the energetics and temporal dynamics of Coulomb explosion of molecular clusters of deuterium (D2)n/2 (n = 480 - 7.6 x 10(4), cluster radius R0 = 13.1 - 70 A) in ultraintense laser fields (laser peak intensity I = 10(15) - 10(20)W cm(-2)). The energetics of Coulomb explosion was inferred from the dependence of the maximal energy EM and the average energy Eav of the product D+ ions on the laser intensity, the laser pulse shape, the cluster radius, and the laser frequency. Electron dynamics of outer cluster ionization and nuclear dynamics of Coulomb explosion were investigated by molecular dynamics simulations. Several distinct laser pulse shape envelopes, involving a rectangular field, a Gaussian field, and a truncated Gaussian field, were employed to determine the validity range of the cluster vertical ionization (CVI) approximation. The CVI predicts that Eav, EM proportional to R0(2) and that the energy distribution is P(E) proportional to E1/2. For a rectangular laser pulse the CVI conditions are satisfied when complete outer ionization is obtained, with the outer ionization time toi being shorter than both the pulse width and the cluster radius doubling time tau2. By increasing toi, due to the increase of R0 or the decrease of I, we have shown that the deviation of Eav from the corresponding CVI value (Eav(CVI)) is (Eav(CVI) - Eav)/Eav(CVI) approximately (toi/2.91tau2)2. The Gaussian pulses trigger outer ionization induced by adiabatic following of the laser field and of the cluster size, providing a pseudo-CVI behavior at sufficiently large laser fields. The energetics manifest the existence of a finite range of CVI size dependence, with the validity range for the applicability of the CVI being R0 < or = (R0)I, with (R0)I representing an intensity dependent boundary radius. Relating electron dynamics of outer ionization to nuclear dynamics for Coulomb explosion induced by a Gaussian pulse, the boundary radius (R0)I and the corresponding ion average energy (Eav)I were inferred from simulations and described in terms of an electrostatic model. Two independent estimates of (R0)I, which involve the cluster size where the CVI relation breaks down and the cluster size for the attainment of complete outer ionization, are in good agreement with each other, as well as with the electrostatic model for cluster barrier suppression. The relation (Eav)I proportional to (R0)I(2) provides the validity range of the pseudo-CVI domain for the cluster sizes and laser intensities, where the energetics of D+ ions produced by Coulomb explosion of (D)n clusters is optimized. The currently available experimental data [Madison et al., Phys. Plasmas 11, 1 (2004)] for the energetics of Coulomb explosion of (D)n clusters (Eav = 5 - 7 keV at I = 2 x 10(18) W cm(-2)), together with our simulation data, lead to the estimates of R0 = 51 - 60 A, which exceed the experimental estimate of R0 = 45 A. The predicted anisotropy of the D+ ion energies in the Coulomb explosion at I = 10(18) W cm(-2) is in accord with experiment. We also explored the laser frequency dependence of the energetics of Coulomb explosion in the range nu = 0.1 - 2.1 fs(-1) (lambda = 3000 - 140 nm), which can be rationalized in terms of the electrostatic model.     

Characterization of supramolecular polyphenol-chromium(III) clusters by molecular dynamics simulations

The binding of Cr(III) with (2R,3S,4R)-(+)-3,3',4,4',7-flavanpentol in aqueous solution is investigated by atomistic molecular dynamics simulations concentrating the analysis of the sampled data on the polyphenol ability to chelate metal ions and to form large noncovalently bonded molecular and supramolecular architectures.     

Kinematic effects associated with molecular frames in structural isomerization dynamics of clusters

Kinematic effects associated with movements of molecular frames, which specify instantaneous orientation of molecules, is investigated in structural isomerization dynamics of a triatomic cluster whose total angular momentum is zero. The principal-axis frame is employed to introduce the so-called principal-axis hyperspherical coordinates, with which the mechanism of structural isomerization dynamics of the cluster is systematically analyzed. A force called "democratic centrifugal force" is extracted from the associated kinematics. This force arises from an intrinsic non-Euclidean metric in the internal space and has an effect of distorting the triatomic cluster to a collapsed shape and of trapping the system around collinear transition states. The latter effect is particularly important in that the kinematics effectively makes a basin at the saddle (transition state) on the potential surface. Based on this framework, we study the effect of the gauge field associated with the Eckart frame in internal space, which has not been carefully examined in the conventional reaction rate theories. Numerical comparison between the dynamics with and without the gauge field has revealed that this field has an effect to suppress the rate of isomerization reaction to a considerable amount. Thus a theory neglecting this effect will significantly overestimate the rate of isomerization. We show the physical origin of this suppressing effect.     

Energy minimization and molecular dynamics studies of Asn-102 elastase

Summary Four isomeric forms of the Asn-102 PPE (D 102 N mutant according to the emerging protocol, [Knowles, Science, 236 (1987) 1252–1258]) have been investigated using energy minimization (EM) and molecular dynamics (MD) techniques. MD simulation data for 175 ps are reported for each form (in total 700 ps for about 2500 atoms). The His-57 N-protonated forms are calculated to be more stable than the N-protonated ones. The active site region of the most stable form is very similar to that found in the D102N rat trypsin enzyme [Craik et al., Science, 237 (1987) 909–913]. Conformations of the active sites and their hydrogen bond patterns are presented for each of these forms and are compared with the structure of the native enzyme active site. The pH dependent activity of the D102N derivative is discussed.