Droplet dynamics for asymmetric Ising model

Nucleation from a metastable state is studied for an anisotropic Ising model at very low temperatures. It turns out that the critical nucleus as well as configurations on a typical path to it differ from the Wulff shape of an equilibrium droplet.     

Theoretical studies on the reaction path and dynamics of the reaction CH_3+H_2O→CH_4+OH

The reaction path, the dynamical properties along the reaction path and CVT rate constants are computed by the ab initio MO method, the reaction path Hamiltonian theory and the variational transition state theory. The results show that the effect of the electron correlation energy on activation barrier is large, the recrossing and tunneling effects exist in the reaction.     

Atomistic Simulation Study of Calcium,Phosphate and Fluoride Ion Association to the Teleopeptide‐Tails of Collagen – Initial Steps to Biomineral Formation

The attachment of single ions to putative adsorption sites in the tails of collagen fibers is investigated by means of molecular dynamics simulations and discussed with respect to the very early steps of apatite/collagen biomineral formation. Our studies clearly demonstrate an increased flexibility of the tails of the triple‐helical collagen protein. Apart from the termini of the backbone, several side chains were also observed to be freely accessible to ion attachment from aqueous solution. The teleopeptide was systematically scanned for suitable adsorption sites for calcium, phosphate and fluoride ions. Association of these ions was then explored from potential of mean force calculations. The resulting energy profiles reveal a variety of favorable protein‐ion bonds and hint at the suitability of the collagen tails to promote apatite aggregation.     

Using molecular dynamics simulations to identify the key factors responsible for chiral recognition by an amino acid-based molecular micelle

Molecular dynamics (MD) simulations were used to investigate the binding of six chiral compounds to the amino acid-based molecular micelle (MM) poly-(sodium undecyl-(L)-leucine-leucine) or poly(SULL). The MM investigated is used as a chiral selector in capillary electrophoresis. The project goal was to characterize the chiral recognition mechanism in these separations and to move toward predictive models to identify the best amino acid-based MM for a given separation. Poly(SULL) was found to contain six binding sites into which chiral compounds could insert. Four sites had similar sizes, shapes, and electrostatic properties. Enantiomers of alprenolol, propranolol, 1,1′-bi-2-naphthyl-2,2′-diyl hydrogen phosphate, 1,1′-bi-2-naphthol, chlorthalidone, or lorazepam were separately docked into each binding pocket and MD simulations with the resulting intermolecular complexes were performed. Solvent-accessible surface area calculations showed the compounds preferentially associated with binding sites where they penetrated into the MM core and shielded their non-polar atoms from solvent. Furthermore, with five of the six compounds the enantiomer with the most favorable free energy of MM association also experienced the most favorable intermolecular hydrogen bonding interactions with the MM. This result suggests that stereoselective intermolecular hydrogen bonds play an important role in chiral discrimination in separations using amino acid-based MMs.GRAPHICAL ABSTRACT     

Influence of static and dynamical structural changes on ultrafast processes mediated by conical intersections

The technological needs imposed by the exponential miniaturization trend of conventional electronic devices has drawn attention towards the development of smaller and faster devices like ultrafast molecular switches. In recent years molecular switches emerge again in the focus of active and innovative research with state-of-the-art optical tools recording their dynamics in real time. Still many questions about the underlying microscopic mechanism are left open, including potential factors that effect the switching process in either way, improve or worsen it. Due to the complexity of such molecules it is difficult to obtain a global answer from experiment alone. On the other side molecular switches are generally too large for a complete quantum chemical and quantum dynamical calculation. In our group we therefore developed an ab initio based modular model to handle the laser induced quantum dynamics in molecular switches like fulgides. It enables us to study the effect of internal molecular coupling and of the molecular response to external fields. We can investigate the related wave packet dynamics, the switching efficiency and the controllability. Our results focus on the laser induced ring opening in fulgides, which equals one direction of the switching process. Presented are the influence of a conical intersection seam and of time-dependent potentials, mimicking the mean interaction with the environment. Furthermore the relation of controllability and the wave packet's momentum is studied and the influence of potential barriers on the switching dynamics is shown.     

Considerations on the recognition of the D1 receptor by agonists

Summary A model of the mechanism for recognition of the D1 receptor has been developed. Conformational analysis for 10 agonists from diverse chemical families was carried out as a first step toward the characterization of the bioactive form. First, maximum structural overlap of the features common to all ligands allowed a simple identification of the candidate bioactive form for each ligand. At a second level of characterization, steric and electronic properties were computed for all accessible structures to analyze those properties that may modulate receptor recognition.     

ZnCl2熔盐的分子动力学模拟

近年来，若干作者根据几种粒子间函数，对ZnCl2熔盐结构做过分子动力学模拟［1－3］其出的Zn－Cl和Cl－Cl离子间的偏径向分布函数与中子衍射实测值符合较好，但Zn－Zn离子间距以及Zn－Zn间配位数计算值多偏高．鉴于Busing势函数在多价卤化物馆盐的分子动力学计算中应用效果较好[4],我们试用Busing势函数为基础对ZnCl2熔盐结构和能量做分子动力学计算．1研究方法计算所用粒子势为Busing势函数此处，Zi为离子的电行数（ZZn。＋＝2，Zcl－＝1），几；是离子有效半径，人为＊离子的“硬度”参数·据文献问，f二0．005071，尸zn。十二0．…     

Structure and Dynamics of Methyl-substituted Beryllocene: [Be(C5Me5)2]

In this paper we study theoretically the molecular structure of [Be(C₅Me₅)₂], with special interest in similarities and differences found in the computed geometric parameters, depending on the treatment of the electron correlation used. Given the low energy differences found between the different configurations studied (less than 4 kcal mol⁻¹), and the high fluxionality found in experimental studies for this compound, we analyzed the dynamics of the system by means of first principles molecular dynamics calculations. A complex dynamics is found and analyzed in terms of two molecular rearrangement processes: 1,2-sigmatropic intraring rearrangement and a ring inversion mechanism that interchanges the roles (with regard to their coordination to the central Be atom) of the two rings.Dedicated to Dr. Jeal Paul Malrieu on occasion of his 65th birthday.     

On a New Formulation of the Real-Time Propagator

The development of theoretical tools for the study of dynamical phenomena of many-particle systems on the quantum level is a fundamental challenge since many decades. A lot of efforts have been invested on Feynman's path integral approach, however, no computationally tractable method for investigating realistic systems could be developed up to now. In this paper we propose an alternative representation of the real-time many-body evolution operator formulated within the framework of the auxiliary field formalism. Our goal is to derive a new auxiliary field functional integral representation, in which the large oscillations of the functional integrand are reduced, in order to render the auxiliary field approach more attractive for real-time computation. This objective is attained using a generalized version of the method of Gaussian equivalent representation of Efimov and Ganbold [Phys. Stat. Sol. 168 (1991) 165], which eliminates the low-order fluctuations of the auxiliary field from the interaction functional.     

Effects of Poly(ethylene glycol) Doping on the Behavior of Pyrene, Rhodamine 6G, and Acrylodan-Labeled Bovine Serum Albumin Sequestered within Tetramethylorthosilane-Derived Sol-Gel-Processed Composites

We investigate the effects of controlled poly(ethylene glycol) (PEG) doping on the behavior of pyrene, rhodamine 6G (R6G), and acrylodan-labeled bovine serum albumin (BSA-Ac) sequestered within tetramethylorthosilicate (TMOS)-derived sol-gel-processed materials. To probe the dipolarity of the local environment within the composite we performed static fluorescence measurements on pyrene as the composites aged. We found that small levels of PEG loading effected significant enhancements in the local dipolarity surrounding the average pyrene molecule. Time-resolved fluorescence anisotropy measurements were used to follow the rotational reorientation dynamics of R6G as the composites aged. As the PEG loading increased, the R6G reorientational mobility increased. Nitrogen adsorption techniques were used to quantify the effects of PEG doping level on the surface area and final xerogel pore features. A large reduction in surface area was observed with PEG doping, but no detectable change in pore size was noted. The effects of PEG doping on a biomolecule were probed by following the time-resolved fluorescence anisotropy decay of BSA-Ac. These results showed that PEG doping resulted in increased biomolecule dynamics relative to that found for a neat, undoped TMOS-derived composites. Together these results show that PEG doping can be used to tune the sol-gel-processed composite dipolarity, alter the mobility of dopants sequestered within the composite, control analyte acessibility to the sensing chemistry, and modulate the internal dynamics within a biodopant.