On the critical behavior of the magnetization in high-dimensional Ising models

We derive rigorously general results on the critical behavior of the magnetization in Ising models, as a function of the temperature and the external field. For the nearest-neighbor models it is shown that ind4 dimensions the magnetization is continuous atT _c and its critical exponents take the classical values=3 and=1/2, with possible logarithmic corrections atd=4. The continuity, and other explicit bounds, formally extend tod>3 1/2. Other systems to which the results apply include long-range models ind=1 dimension, with 1/|x–y| couplings, for which 2/(–1) replacesd in the above summary. The results are obtained by means of differential inequalities derived here using the random current representation, which is discussed in detail for the case of a nonvanishing magnetic field.Research supported in part by NSF grant PHY-8301493 A02, and by a John S. Guggenheim Foundation fellowship (M.A.).     

Heat conductivity of three periodic Hard Disks via nonequilibrium molecular dynamics

We use a driving field, of the type first suggested by Evans, to generate a steady heat current in the simplest possible system, a two-dimensional periodic fluid of three hard disks. Hard-disk motion equations can be conveniently derived from repulsive constant-force or linear-force potentials by considering the infinitely repulsive limit of these potentials. We show that the isoenergetic and isokinetic forms of the nonequilibrium equations of motion generate steady-state heat conductivities differing by terms of order 1/N, whereN is the number of particles. The resulting conductivities appear to vary as the logarithm of the driving field strength. Even at low fields, the three-body periodic-system results lie well below Enskog's infinite-system prediction.     

Anion substitution effects on structure and magnetism in the chromium chalcogenide Cr5Te8—Part I: Cluster glass behavior in trigonal Cr(1+x)Q2 with basic cell (Q=Te, Se; Te:Se=7:1)

The effect of substitution of the anion Te by Se in non-stoichiometric Cr₅Te₈ has been investigated with respect to its crystal structure, magnetic properties, and electronic structure. The compounds Cr_(1+x)Q₂ (Q=Te, Se; Te:Se=7:1; (1+x)=1.234(6), 1.264(6), 1.300(7)) were synthesized at elevated temperatures followed by quenching the samples to room temperature. The crystal structures have been refined with X-ray powder diffraction data with the Rietveld method in the trigonal space group with lattice parameters a=3.8651(1)-3.8831(1) Å and c=5.9917(2)-6.0528(2) Å. The structure is related to the NiAs structure with full and deficient metal layers stacking alternatively along the c-axis. The irreversibility in the field-cooled/zero-field-cooled magnetization suggests that the substitution effects of one Te by one Se is strong enough to cause cluster-glass behavior, from ferromagnetic Cr₅Te₈ to cluster-glass Cr_(1+x)Q₂. Non-saturation magnetizations at 5.5 T and the magnetic relaxation results further support the existence of cluster-glass behavior. Accompanying SPR-KKR (spin-polarized relativistic Korringa-Kohn-Rostoker) band structure calculations strongly support the observation that the Cr(1) sites are preferentially occupied by Cr atoms and predict that these compounds are metallic. Results for the spin-resolved DOS and magnetic moments on each crystallographic sites are presented.     

Parameterization of peptide 13C carbonyl chemical shielding anisotropy in molecular dynamics simulations.

NMR chemical shielding anisotropy (CSA) relaxation is an important tool in the study of dynamical processes in proteins and nucleic acids in solution. Herein, we investigate how dynamical variations in local geometry affect the chemical shielding anisotropy relaxation of the carbonyl carbon nucleus, using the following protocol: 1) Using density functional theory, the carbonyl (13)C' CSA is computed for 103 conformations of the model peptide group N-methylacetamide (NMA). 2) The variations in computed (13)C' CSA parameters are fitted against quadratic hypersurfaces containing cross terms between the variables. 3) The predictive quality of the CSA hypersurfaces is validated by comparing the predicted and de novo calculated (13)C' CSAs for 20 molecular dynamics snapshots. 4) The CSA fluctuations and their autocorrelation and cross correlation functions due to bond-length and bond-angle distortions are predicted for a chemistry Harvard molecular mechanics (CHARMM) molecular dynamics trajectory of Ca(2+)-saturated calmodulin and GB3 from the hypersurfaces, as well as for a molecular dynamics (MD) simulation of an NMA trimer using a quantum mechanically correct forcefield. We find that the fluctuations can be represented by a 0.93 scaling factor of the CSA tensor for both R(1) and R(2) relaxations for residues in helix, coil, and sheet alike. This result is important, as it establishes that (13)C' relaxation is a valid tool for measurement of interesting dynamical events in proteins.     

室温熔盐二(三氟甲基磺酸酰)亚胺锂-尿素体系的分子动力学模拟

应用分子动力学模拟了25 ℃和50 ℃时新型室温熔盐二(三氟甲基磺酸酰)亚胺锂[LiN(SO2CF3)2, LiTFSI]与尿素(摩尔比为1:3.6)体系的结构与动力学性质. 在两个温度下体系的微观结构基本相同, Li+的配位数约为5, 且都是与溶剂和阴离子中的氧原子发生配位. 对TFSI-的研究表明, 每个TFSI-只提供四个氧中的一个与Li+配位; 而且在Li+的配位层中, TFSI-具有顺、反和gauche 等不同的构象, 并且不同构象出现的几率会随着温度的改变而改变.     

Incorporating the effect of ionic strength in free energy calculations using explicit ions

The incorporation of explicit ions to mimic the effect of ionic strength or to neutralize the overall charge on a system in free energy calculations using molecular dynamics simulations is investigated. The difference in the free energy of hydration between two triosephosphate isomerase inhibitors calculated at five different ion concentrations is used as an example. We show that the free energy difference can be highly sensitive to the presence of explicit ions even in cases where the mutation itself does not involve a change in the overall charge. The effect is most significant if the molecule carries a net charge close to the site mutated. Furthermore, it is shown that the introduction of a small number of ions can lead to very severe sampling problems suggesting that in practical calculations convergence can best be achieved by incorporating either no counterions or by simulating at high ionic strength to ensure sufficient sampling of the ion distribution.     

Binding Ensembles of p53-MDM2 Peptide Inhibitors by Combining Bayesian Inference and Atomistic Simulations

Designing peptide inhibitors of the p53-MDM2 interaction against cancer is of wide interest. Computational modeling and virtual screening are a well established step in the rational design of small molecules. But they face challenges for binding flexible peptide molecules that fold upon binding. We look at the ability of five different peptides, three of which are intrinsically disordered, to bind to MDM2 with a new Bayesian inference approach (MELD × MD). The method is able to capture the folding upon binding mechanism and differentiate binding preferences between the five peptides. Processing the ensembles with statistical mechanics tools depicts the most likely bound conformations and hints at differences in the binding mechanism. Finally, the study shows the importance of capturing two driving forces to binding in this system: the ability of peptides to adopt bound conformations (ΔGconformation) and the interaction between interface residues (ΔGinteraction).     

A computational model of the nicotinic acetylcholine binding site

We have derived a model of the nicotinic acetylcholine binding site. This was accomplished by using three known agonists (acetylcholine, nicotine and epibatidine) as templates around which polypeptide side chains, found to be part of the receptor cavity from published molecular biology studies, are allowed to flow freely in molecular dynamics simulations and mold themselves around these templates. The resulting supramolecular complex should thus be a complement, both in terms of steric effects as well as electronic effects, to the agonists and it should be a good estimation of the true receptor cavity structure. The shapes of those minireceptor cavities equilibrated rapidly on the simulation time scale and their structural congruence is very high, implying that a satisfactory model of the nicotinic acetylcholine binding site has been achieved. The computational methodology was internally tested against two rigid and specific antagonists (dihydro--erytroidine and erysoidine), that are expected to give rise to a somewhat differently shaped binding site compared to that derived from the agonists. Using these antagonists as templates there were structural reorganizations of the initial receptor cavities leading to distinctly different cavities compared to agonists. This indicates that adequate times and temperatures were used in our computational protocols to achieve equilibrium structures for the agonists. Overall, both minireceptor geometries for agonists and antagonists are similar with the exception of one amino acid (ARG209).     

受限于纳米碳管中的乙醇分子的结构和扩散的分子模拟研究

利用分子模拟研究了常温常压下受限于(8,8) (管径1.081 nm)和(15,15) (管径 2.035 nm)单壁纳米碳管中的乙醇分子. 对受限分子的径向密度分布和氢键等静态性质以及扩散性质进行了分析. 结果显示在管内乙醇分子的平均氢键数目和主体相一致. 乙醇分子在(8,8)碳管内具有高度有序的结构, 而在(15,15)碳管内由于空间的增大导致结构有序度的降低, 其中分子取向已呈随机分布. 进一步对扩散系数的分析发现, 在管内乙醇分子的轴向扩散系数低于主体相, 特别在(8,8)碳管内乙醇分子几乎丧失了轴向扩散能力.     

Solvation of uranyl(II), europium(III) and europium(II) cations in "basic" room-temperature ionic liquids: a theoretical study

We report a molecular dynamics study of the solvation of UO2(2+), Eu3+ and Eu2+ ions in two "basic" (Lewis acidity) room-temperature ionic liquids (IL) composed of the 1-ethyl-3-methylimidazolium cation (EMI+) and a mixture of AlCl4- and Cl- anions, in which the Cl-/AlCl4- ratio is about 1 and 3, respectively. The study reveals the importance of the [UO2Cl4]2- species, which spontaneously form during most simulations, and that the first solvation shell of europium is filled with Cl- and AlCl4- ions embedded in a cationic EMI+ shell. The stability of the [UO2Cl4]2- and [Eu(III)Cl6]3- complexes is supported by quantum mechanical calculations, according to which the uranyl and europium cations intrinsically prefer Cl- to the AlCl4- ion. In the gas phase, however, [Eu(III)Cl6]3- and [Eu(II)Cl6]4- complexes are predicted to be metastable and to lose two to three Cl- ions. This contrasts with the results of simulations of complexes in ILs, in which the "solvation" of the europium complexes increases with the number of coordinated chlorides, leading to an equilibrium between different chloro species. The behavior of the hydrated [Eu(OH2)8]3+ complex is considered in the basic liquids; the complex exchanges H2O molecules with Cl- ions to form mixed [EuCl3(OH2)4] and [EuCl4(OH2)3]- complexes. The results of the simulations allow us to better understand the microscopic nature and solvation of lanthanide and actinide complexes in "basic" ionic liquids.