Hybrid Qm/Mm Simulation Of The Hydration Phenomena Of Dipalmitoylphosphatidylcholine Headgroup

The polar headgroup of dipalmitoylphosphatidylcholine (DPPC) molecule both in gas phase and aqueous Solution is investigated by the hybrid quantum mechanical/molecular mechanical (QM/MM) method, in which the polar head of DPPC molecule and the bound water molecules are treated with density functional theory (DFT), while the apolar hydrocarbon chain of DPPC molecule is treated with MM method. It is demonstrated that the hybrid QM/MM method is both accurate and efficient to describe the conformations of DPPC headgroup. Folded structures of headgroup are found in gas phase calculations. In this work, both monohydration and polyhydration phenomena are investigated. In monohydration, different water association sites are studied. Both the hydration energy and the quantum properties of DPPC and water molecules are calculated at the DFT level of theory after geometry optimization. The binding force of monohydration is estimated by using the scan method. In polyhydration, more extended conformations are found and hydration energies in different polyhydration styles are estimated. (C) 2008 Elsevier Inc. All rights reserved.

Spectroscopic studies of gas-phase molecular clusters

Spectroscopic investigations of hydrogen-bonding and van der Waals' interactions m molecular clusters were studied by the techniques of infrared predissociation and resonance-enhanced multiphoton ionization spectroscopies (REMPI). Ab initio calculations were applied in conjunction for data interpretation.

The infrared predissociation spectroscopy of CN^-•(H_2O)_n (n = 2 - 6) clusters was reported in the region of 2950 - 3850 cm^(-1). The hydrogen bondings for the C-site and N-site binding, and among the water molecules were identified for n = 2 to 4. A spectral transition was observed for n = 5 and 6, implying that the anion was surface-bound onto the water aggregates in larger clusters.

The infrared predissociation spectroscopy of Br^-•(NH_3) and I^-•(NH_3)_n (n =1-3) clusters was reported in the region of 3050-3450 cm^(-1). For the Br^-•(NH_3) complex, a dominating ionic NH stretch appeared at 3175 cm^(-1), and the weaker free NH stretch appeared at 3348 cm^(-1). The observed spectrum was consistent to the structure in which there was one nearly linear hydrogen bond between Br^- and the NH_3 moiety. For the I^- •(NH_3) complex, five distinct IR absorption bands were observed in the spectrum. The spectrum was not consistent with basic frequency patterns of three geometries considered in the ab initio calculations - complex with one, two and three hydrogen bondings between I^- and the NH_3 moiety. Substantial inhomogenous broadening were displayed in the spectra for I^-•(NH_3)_n (n =2-3), suggesting the presence of multiple isomers.

The REMPI spectroscopy of the bound 4p ^2П 1/2 and ^2П 3/2 states, and the dissociative 3d ^2Σ^+ 1/2 state in the Al•Ar complex was reported. The dissociative spectrum at Al^+ channel suggested the coupling of the 4p ^2П 1/2,3/2 states to the repulsive 3d ^2Σ^+ 1/2 state. The spin-electronic coupling was further manifested in the dissociative Al^+ spectrum of the 3d ^2Σ^+ 1/2 state. Using the potential energy curves obtained from ab initio calculations, a bound → continuum Franck-Condon-intensity simulation was performed and compared with the one-photon 3d ^2Σ^+ 1/2 profile. The agreement provided evidence for the petturbation above the Al(3d)Ar dissociation limit, and the repulsive character of the 3d ^2Σ^+ 1/2 state.

Tensile and compressive mechanical behavior of twinned silicon carbide nanowires

Molecular dynamics simulations with the Tersoff potential were used to study the response of twinned SiC nanowires under tensile and compressive strain. The critical strain of the twinned nanowires can be enhanced by twin stacking faults, and their critical strains are larger than those of perfect nanowires with the same diameters. Under axial tensile strain, the bonds of the nanowires are stretched just before failure. The failure behavior is found to depend on the twin segment thickness and the diameter of the nanowires. An atomic chain is observed for thin nanowires with small twin segment thickness under tension strain. Under axial compressive strain, the collapse of twinned SiC nanowires exhibits two different failure modes, depending on the length and diameter of the nanowires, i.e., shell buckling for short nanowires and columnar buckling for longer nanowires.

2-cyanoprop-2-yl dithiobenzoate mediated reversible addition-fragmentation chain transfer polymerization of acrylonitrile targeting a polymer with a higher molecular weight

The reversible addition-fragmentation chain transfer (RAFT) polymerization of acrylonitrile (AN) mediated by 2-cyanoprop-2-yl dithiobenzoate was first applied to synthesize polyacrylonitrile (PAN) with a high molecular weight up to 32,800 and a polydispersity index as low as 1.29. The key to success was ascribed to the optimization of the experimental conditions to increase the fragmentation reaction efficiency of the intermediate radical. In accordance with the atom transfer radical polymerization of AN, ethylene carbonate was also a better solvent candidate for providing higher controlled/living RAFT polymerization behaviors than dimethylformamide and dimethyl sulfoxide. The various experimental parameters, including the temperature, the molar ratio of dithiobenzoate to the initiator, the molar ratio of the monomer to dithiobenzoate, the monomer concentration, and the addition of the comonomer, were varied to improve the control of the molecular weight and polydispersity index. The molecular weights of PANS were validated by gel permeation chromatography along with a universal calibration procedure and intrinsic viscosity measurements. H-1 NMR analysis confirmed the high chain-end functionality of the resultant polymers.

高温高压下H2O-CH4体系PVTx性质的分子动力学模拟研究

Geological fluids are important components in the earth system. To study thephysical chemistry properties and the evolution of fluid system turns out to be one of the most challenging issues in geosciences. Besides the conventional experimental approaches and theoretical or semi-theoretical modeling, molecular level computer simulation(MLCS) emerges as an alternative tool to quantificationally study the physico-chemical properties of fluid under extreme conditions in order to find out the characteristics and interaction of geological fluids in and around earth. Based on our previous study of the intermolecular potential for pure H2O and thestrict evaluation of the competitive potential models for pure CH4 and the ab initio fitting potential surface across H2O-CH4 molecules in this study, we carried out more than two thousand molecular dynamics simulations for the PVTx properties of pure CH4 and the H2O-CH4 mixtures. Comparison of 1941 simulations with experimental PVT data for pure CH4 shows an average deviation of 0.96% and a maximum deviation of 2.82%. The comparison of the results of 519 simulations of the mixtures with the experimental measurements reveals that the PVTx properties of the H2O-CH4 mixtures generally agree with the extensive experimental data with an average deviation of 0.83% and 4% in maximum, which is equivalent to the experimental uncertainty. Moreover, the maximum deviation between the experimental data and the simulation results decreases to about 2% as temperature and pressure increase,indicating that the high accuracy of the simulation is well retained in the high temperature and pressure region. After the validation of the simulation method and the intermolecular potential models, we systematically simulated the PVTx properties of this binary system from 673 K and 0.05 GPa to 2573 K and 10 GPa. In order to integrate all the simulation results and the experimental data for the calculation of thermodynamic properties, an equation of state (EOS) is developed for the H2O-CH4 system covering 673 to 2573 K and 0.01 to 10 GPa. Isochores for compositions < 4 mol% CH4 up to 773 K and 600 MPa are also determined in this thesis.

分子水平上研究地质流体的物理化学性质

Geological fluids exist in every geosphere of the Earth and play important roles in many processes of material transformations, energetic interchanges and geochemical interactions. To study the physicochemical properties and geochemical behaviors of geological fluids turn Girt to be one of the challenging issues in geosciences. Compared with conventional approaches of experiments and semi-theoretical modeling, computer simulation on molecular level shows its advantages on quantitative predictions of the physicochemical properties of geological fluids under extreme conditions and emerges as a promising approach to find the characteristics of geological fluids and their interactions in different geospheres of the Earth interior.This dissertation systematically discusses the physicochemical properties of typical geological fluids with state-of-the-art computer simulation techniques. The main results can be summarized as follows: (1) The experimental phase behaviors of the systems CH4-C2H6 and. CO2 have been successfully reproduced with Monte Carlo simulations. (2) Through comprehensive isothermal-isobaric molecular dynamics simulations, the PVT data of water hia^e been extended beyond experimental range to about 2000 K and 20 GPa and an improved equation of state for water has been established. (3) Based on extensive computer simulations, am optimized molecular potential for carbon dioxide have been proposed, this model is expected to predict different properties of carbon dioxide (volumetric properties, phase equilibria, heat of vaporization, structural and dynamical properties) with improved accuracies. (4) On the basis of the above researches of the end-members, a set of parameters for unlike interactions has been proposed by non-linear fitting to the ab initio potential surface of CO2-H2O and is superior to the common used mixing rule and the results of prior workers vs/Ith remarkable accuracies, then a number of simulations of the mixture have been carried out to generate data under high temperatures and pressures as an important complement to the limited experiments. (5) With molecular dynamics simulations, various structural, dynamical and thermodynamical properties of ionic solvations and associations have been oomprehensively analyzed, these results not only agree well with experimental data and first principle calculation results, but also reveal some new insights into the microscopic ionic solvation and association processes.

Time-dependent transport through molecular junctions.

We investigate transport properties of molecular junctions under two types of bias--a short time pulse or an ac bias--by combining a solution for Green's functions in the time domain with electronic structure information coming from ab initio density functional calculations. We find that the short time response depends on lead structure, bias voltage, and barrier heights both at the molecule-lead contacts and within molecules. Under a low frequency ac bias, the electron flow either tracks or leads the bias signal (resistive or capacitive response) depending on whether the junction is perfectly conducting or not. For high frequency, the current lags the bias signal due to the kinetic inductance. The transition frequency is an intrinsic property of the junctions.

Molecular electrostatic properties of ions in an ionic liquid

We have analysed the electronic wave functions from an ab initio simulation of the ionic liquid (room temperature molten salt) dimethyl imidazolium chloride ([dmim][Cl] or [C1mim][Cl]) using localized Wannier orbitals. This allows us to assign electron density to individual ions. The probability distributions of the ionic dipole moments for an isolated ion and for ions in solution are compared. The liquid environment is found to polarize the cation by about 0.7 D and to increase the amplitude of the fluctuations in the dipole moments of both cation and anion. The relative changes in nuclear and electronic contributions are shown. The implications for classical force fields are discussed.

Suggested Binding Mechanism of the HIV-gp120 to its CD4 Receptor

The molecular recognition and attachment of the CD4 molecule and the HIV envelope glycoprotein (gp120) might be described as a consecutive three-step molecular recognition process. 1. (a) Long range interaction: electrostatic pre-orientation, 2. (b) short range interaction: electronic attachment followed by a ‘Locking-in’ (via aromatic ring orientation) and 3. (c) internal interaction (induced fit): conformational readjustment of the protein molecules. On the basis of the preliminary investigations (X-ray structures of CD4 and biological studies of CD4 and gp120 point mutants) we described a computational model. This approach consists of empirical calculations as well as ab initio level of quantum chemistry. The conformational analysis of the wild type and mutant CD4 molecules was supported by molecular mechanics and dynamics (Amber force field). The latter analysis involves the application of a novel method, the Amino Acid Conformation Assignment of Proteins (ACAP) software, developed for the notation of secondary protein structures. According to the cardinal role of the electrostatic factors during this interaction, several ab initio investigations were performed for better understanding of the recognition process on submolecular level. Using the above mentioned computational model, we could interpret the basic behaviours and predict some additional features of CD4-gp120 interaction, in spite of the missing gp120 X-ray structure.

Molecular Modeling on Inhibitor Complexes and Active-Site Dynamics of Cytochrome P450 C17, a Target for Prostate Cancer Therapy

A molecular model for the P450 enzyme cytochrome P450 C17 (CYP17) is presented based on sequence alignments of multiple template structures and homology modeling. This enzyme plays a central role in the biosynthesis of testosterone and is emerging as a major target in prostate cancer, with the recently developed inhibitor abiraterone currently in advanced clinical trials. The model is described in detail, together with its validation, by providing structural explanations to available site-directed mutagenesis data. The CYP17 molecule in this model is in the form of a triangular prism, with an edge of similar to 55 angstrom and a thickness of similar to 37 angstrom. It is predominantly helical, comprising 13 alpha helices interspersed by six 3(10) helices and 11 beta-sheets. Multinanosecond molecular dynamics simulations in explicit solvent have been carried out, and principal components analysis has been used to reveal the details of dynamics around the active site. Coarse-grained methods have also been used to verify low-frequency motions, which have been correlated with active-site gating. The work also describes the results of docking synthetic inhibitors, including the drug abiraterone and the natural substrate pregnenolone, in the CYP17 active site together with molecular dynamics simulations on the complexes. (C) 2010 Elsevier Ltd. All rights reserved.

Pinning down the solid-state polymorphism of the ionic liquid [bmim][PF6]

The solid-state polymorphism of the ionic liquid 1-butyl-3-methylimidazolium hexafluorophosphate, [bmim][PF6], has been investigated via low-temperature and high-pressure crystallisation experiments. The samples have been characterised by single-crystal X-ray diffraction, optical microscopy and Raman spectroscopy. The solid-state phase behaviour of the compound is confirmed and clarified with respect to previous phase diagrams. The structures of the previously reported gamma-form, which essentially exhibits a G'T cation conformation, as well as those of the elusive beta- and alpha-forms, are reported. Crystals of the beta-phase are twinned and the structure is heavily disordered; the cation conformation in this form is predominantly TT, though significant contributions from other less frequently encountered conformers are also observed at low temperature and high pressure. The cation conformation in the alpha-form is GT; the presence of the G'T conformer at 193 K in this phase can be eliminated on cooling to 100 K. Whilst X-ray structural data are overall in good agreement with previous interpretations based on Raman and NMR studies, they also reveal a more subtle interplay of intermolecular interactions, which give rise to a wider range of conformers than previously considered.

Positron scattering and annihilation on noble-gas atoms

Positron scattering and annihilation on noble-gas atoms is studied ab initio using many-body theory methods for positron energies below the positronium formation threshold. We show that in this energy range, the many-body theory yields accurate numerical results and provides a near-complete understanding of the positron–noble-gas atom system. It accounts for positron-atom and electron-positron correlations, including the polarization of the atom by the positron and the nonperturbative effect of virtual positronium formation. These correlations have a large influence on the scattering dynamics and result in a strong enhancement of the annihilation rates compared to the independent-particle mean-field description. Computed elastic scattering cross sections are found to be in good agreement with recent experimental results and Kohn variational and convergent close-coupling calculations. The calculated values of the annihilation rate parameter Zeff (effective number of electrons participating in annihilation) rise steeply along the sequence of noble-gas atoms due to the increasing strength of the correlation effects, and agree well with experimental data.

ORBITAL MIXING IN CO CHEMISORPTION ON TRANSITION-METAL SURFACES

The chemisorption of CO on metal surfaces is widely accepted as a model for understanding chemical bonding between molecules and solid surfaces, but is nevertheless still a controversial subject. Ab initio total energy calculations using density functional theory with gradient corrections for CO chemisorption on an extended Pd{110} slab yield good agreement with experimental adsorption energies. Examination of the spatial distribution of individual Bloch states demonstrates that the conventional model for CO chemisorption involving charge donation from CO 5 sigma states to metal states and back-donation from metal states into CO 2 pi states is too simplistic, but the computational results provide direct insight into the chemical bonding within the framework of orbital mixing (or hybridisation). The results provide a sound basis for understanding the bonding between molecules and metal surfaces.

GRADIENT CORRECTIONS IN DENSITY-FUNCTIONAL THEORY CALCULATIONS FOR SURFACES - CO ON PD(110)

Ab initio total energy calculations have been performed for CO chemisorption on Pd(110). Local density approximation (LDA) calculations yield chemisorption energies which are significantly higher than experimental values but inclusion of the generalised gradient approximation (GGA) gives better agreement. In general, sites with higher coordination of the adsorbate to surface atoms lead to a larger degree of overbinding with LDA, and give larger corrections with GGA. The reason is discussed using a first-order perturbation approximation. It is concluded that this may be a general failure of LDA for chemisorption energy calculations. This conclusion may be extended to many surface calculations, such as potential energy surfaces for diffusion.

Extraordinary mechanical properties of monatomic C3N2 chain

One-dimensional monatomic chains are promising candidates for technical applications in the field of nanoelectronics due to their unique mechanical, electrical and optical properties. In particular, we investigate the mechanical properties including Young's modulus, ultimate strength and ultimate strain, which are necessities for the stability of the materials by the Car-Parrinello molecular dynamics in this work. The comparative studies for the alternating carbon-nitrogen (C₃N₂) chain and carbon chains (carbyne) of different lengths show that the carbon-nitrogen (C-N) chain is obviously stronger and stiffer than carbynes. Thus the C-N chain, which has been found in decomposition products of the nitromethane explosive simulations, could be a superior nano-mechanical material than the carbyne chain. Furthermore, it is found that the bond order of weakest bond in monatomic chains is positively correlated with Young's modulus and ultimate strength of materials.