Fusion dynamics of symmetric systems near barrier energies

Within the framework of the improved isospin-dependent quantum molecular dynamics (ImIQMD) model,he fusion dynamics of symmetric reaction systems are investigated systematically. Calculations show that the number of nucleon transfer in the neck region is appreciably dependent on the incident energies, but strongly on he reaction systems. A comparison of the neck dynamics is performed for the symmetric reactions 58Ni+58Niand 64Ni+64Ni at energies in the vicinity of the Coulomb barrier. An increase of the ratios of the neutron to proton in the neck region at initial collision stage is observed and obvious for the latter system, which reduces the fusion barrier of two colliding nuclei. The distribution of the dynamical fusion barriers and the fusion excitation functions are calculated and compared with the available experimental data.

Nuclear reaction dynamics induced by halo-nuclei at intermediate energy heavy ion collisions

We studied systematically the reaction dynamics induced by neutron-halo nuclei and proton-halo nuclei within the isospin dependent quantum molecular dynamics, such as the effects of loose bound halo-nuclei on the fragmentation reaction and momentum dissipation for different colliding systems with different beam energies and different impact parameters. In order to emphasize the roles of neutron-halo nucleus B-19 and proton-halo nucleus Al-23 on the reaction dynamics we also calculated the the reaction dynamics induced by the stable nuclei F-19 and Na-23 with equal mass under identical incident channel conditions. Based on the comparison of results of reaction dynamics induced by halo-nucleus colliding systems and stable nucleus collidinmg systems we found that the roles of loose bound halo-nucleus structure on the fragmentation multiplicity and nuclear stopping (momentum dissipation) are important for all of colliding systems with different beam energies and minor impact parameters, such as, the loose bound halo-nuclei structure increases the fragmentation multiplicity, but reduces the nuclear stopping.

Dynamics of strangeness production in heavy-ion collisions near threshold energies

Within the framework of the improved isospin-dependent quantum molecular dynamics (ImIQMD) model, the dynamics of strangeness (K-0,K-+, Lambda, and Sigma(-,0,+)) production in heavy-ion collisions near threshold energies is investigated systematically, with the strange particles considered to be produced mainly by inelastic collisions of baryon-baryon and pion-baryon. Collisions in the region of suprasaturation densities of the dense baryonic matter formed in heavy-ion collisions dominate the yields of strangeness production. Total multiplicities as functions of incident energies and collision centralities are calculated with the Skyrme parameter SLy6. The excitation function of strangeness production is analyzed and also compared with the KaoS data for K+ production in the reactions C-12 + C-12 and Au-197 + Au-197.

Dynamics of pion production in heavy-ion collisions around 1A GeV energies

Within the framework of the improved isospin-dependent quantum molecular dynamics model, the dynamics of pion emission in heavy-ion collisions in the region of 1A GeV energies as a probe of nuclear symmetry energy at suprasaturation densities is investigated systematically. The total pion multiplicities and the pi(-)/pi(+) yields are calculated for selected Skyrme parameters SkP, SLy6, Ska, and SIII and also for the cases of different stiffness of symmetry energy with the parameter SLy6. The influence of Coulomb potential, symmetry energy, and in-medium pion potential on the pion production is investigated and compared to each other by analyzing the distributions of transverse momentum and longitudinal rapidity and also the excitation functions of the total pion and the pi(-)/pi(+) ratio. The directed flow, elliptic flow, and polar-angle distributions are calculated for the cases of different collision centralities and also the various stiffnesses of the symmetry energies. A comparison of the calculations with the available experimental data is performed.

Influence of molecular topology on the static and dynamic properties of single polymer chain in solution

The influence of molecular topology on the structural and dynamic properties of polymer chain in solution with ring structure, three-arm branched structure, and linear structure are studied by molecular dynamics simulation. At the same degree of polymerization (N), the ring-shaped chain possesses the smallest size and largest diffusion coefficient. With increasing N, the difference of the radii of gyration between the three types of polymer chains increases, whereas the difference of the diffusion coefficients among them decreases. However, the influence of the molecular topology on the static and the dynamic scaling exponents is small. The static scaling exponents decrease slightly, and the dynamic scaling exponents increase slightly, when the topology of the polymer chain is changed from linear to ring-shaped or three-arm branched architecture. The dynamics of these three types of polymer chain in solution is Zimm-like according to the dynamic scaling exponents and the dynamic structure factors.

Exploring the mechanism of flexible biomolecular recognition with single molecule dynamics

Combining a single-molecule study of protein binding with a coarse grained molecular dynamics model including solvent (water molecules) effects, we find that biomolecular recognition is determined by flexibilities in addition to structures. Our single-molecule study shows that binding of CBD (a fragment of Wiskott-Aldrich syndrome protein) to Cdc42 involves bound and loosely bound states, which can be quantitatively explained in our model as a result of binding with large conformational changes. Our model identified certain key residues for binding consistent with mutational experiments. Our study reveals the role of flexibility and a new scenario of dimeric binding between the monomers: first bind and then fold.

Energetic basis of molecular recognition in a DNA aptamer

The thermal stability and ligand binding properties of the L-argininamide-binding DNA aptamer (5'-GATCGAAACGTAGCGCCTTCGATC3') were studied by spectroscopic and calorimetric methods. Differential calorimetric studies showed that the uncomplexed aptamer melted in a two-state reaction with a melting temperature T-m = 50.2 +/- 0.2 degrees C and a folding enthalpy Delta H degrees(fold) = -49.0 +/- 2.1 kcal mol(-1). These values agree with values of T-m = 49.6 degrees C and Delta H degrees(fold) = -51.2 kcal mol(-1) predicted for a simple hairpin structure. Melting of the uncomplexed aptamer was dependent upon salt concentration, but independent of strand concentration. The T of aptamer melting was found to increase as L-argininamide concentrations increased. Analysis of circular dichroism titration data using a single-site binding model resulted in the determination of a binding free energy Delta G degrees(bind) = -5.1 kcal mol(-1). Isothermal titration calorimetry studies revealed an exothermic binding reaction with Delta H degrees(bind) = -8.7 kcal mol(-1). Combination of enthalpy and free energy produce ail unfavorable entropy of -T Delta S degrees = +3.6 kcal mol(-1). A molar heat capacity change of -116 cal mol(-1) K-1 was determined from calorimetric measurements at four temperatures over the range of 15-40 degrees C. Molecular dynamics simulations were used to explore the structures of the unligated and ligated aptamer structures.

Single-molecule dynamics reveals cooperative binding-folding in protein recognition

The study of associations between two biomolecules is the key to understanding molecular function and recognition. Molecular function is often thought to be determined by underlying structures. Here, combining a single-molecule study of protein binding with an energy-landscape-inspired microscopic model, we found strong evidence that biomolecular recognition is determined by flexibilities in addition to structures. Our model is based on coarse-grained molecular dynamics on the residue level with the energy function biased toward the native binding structure ( the Go model). With our model, the underlying free-energy landscape of the binding can be explored. There are two distinct conformational states at the free-energy minimum, one with partial folding of CBD itself and significant interface binding of CBD to Cdc42, and the other with native folding of CBD itself and native interface binding of CBD to Cdc42. This shows that the binding process proceeds with a significant interface binding of CBD with Cdc42 first, without a complete folding of CBD itself, and that binding and folding are then coupled to reach the native binding state.

高温高压下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.

Dynamics of Delta/Notch signaling on endomesoderm segregation in the sea urchin embryo.

Endomesoderm is the common progenitor of endoderm and mesoderm early in the development of many animals. In the sea urchin embryo, the Delta/Notch pathway is necessary for the diversification of this tissue, as are two early transcription factors, Gcm and FoxA, which are expressed in mesoderm and endoderm, respectively. Here, we provide a detailed lineage analysis of the cleavages leading to endomesoderm segregation, and examine the expression patterns and the regulatory relationships of three known regulators of this cell fate dichotomy in the context of the lineages. We observed that endomesoderm segregation first occurs at hatched blastula stage. Prior to this stage, Gcm and FoxA are co-expressed in the same cells, whereas at hatching these genes are detected in two distinct cell populations. Gcm remains expressed in the most vegetal endomesoderm descendant cells, while FoxA is downregulated in those cells and activated in the above neighboring cells. Initially, Delta is expressed exclusively in the micromeres, where it is necessary for the most vegetal endomesoderm cell descendants to express Gcm and become mesoderm. Our experiments show a requirement for a continuous Delta input for more than two cleavages (or about 2.5 hours) before Gcm expression continues in those cells independently of further Delta input. Thus, this study provides new insights into the timing mechanisms and the molecular dynamics of endomesoderm segregation during sea urchin embryogenesis and into the mode of action of the Delta/Notch pathway in mediating mesoderm fate.

On the structure and dynamics of ionic liquids

The structure and dynamics of the ionic liquid 1-ethyl-3-methylimidazolium nitrate is studied by molecular dynamics simulations. We find long-range spatial correlations between the ions and a three-dimensional local structure that reflects the asymmetry of the cations. The main contribution to the configurational energy comes from the electrostatic interactions which leads to charge-ordering effects. Radial screening and threedimensional distribution of charge are also analyzed. The motion of a single ion is studied via velocity and reorientational correlation functions. It is found that ions "rattle" in a long-lived cage, while the orientational structure relaxes on a time scale longer than 200 ps. As in a supercooled liquid, the mean square displacements reveal a subdiffusive dynamics. In addition, the presence of dynamic heterogeneities can be detected by analyzing the non-Gaussian behavior of the van Hove correlation function and the spatial arrangement of the most mobile ions. The short-time collective dynamics is also studied through the electric current time correlation function.

Correlated electron-ion dynamics: the excitation of atomic motion by energetic electrons

Correlated electron-ion dynamics (CEID) is an extension of molecular dynamics that allows us to introduce in a correct manner the exchange of energy between electrons and ions. The formalism is based on a systematic approximation: small amplitude moment expansion. This formalism is extended here to include the explicit quantum spread of the ions and a generalization of the Hartree-Fock approximation for incoherent sums of Slater determinants. We demonstrate that the resultant dynamical equations reproduce analytically the selection rules for inelastic electron-phonon scattering from perturbation theory, which control the mutually driven excitations of the two interacting subsystems. We then use CEID to make direct numerical simulations of inelastic current-voltage spectroscopy in atomic wires, and to exhibit the crossover from ionic cooling to heating as a function of the relative degree of excitation of the electronic and ionic subsystems.

Salt-bridge dynamics control substrate-induced conformational change in the membrane transporter GlpT.

Active transport of substrates across cytoplasmic membranes is of great physiological, medical and pharmaceutical importance. The glycerol-3-phosphate (G3P) transporter (GlpT) of the E. coli inner membrane is a secondary active antiporter from the ubiquitous major facilitator superfamily that couples the import of G3P to the efflux of inorganic phosphate (Pi) down its concentration gradient. Integrating information from a novel combination of structural, molecular dynamics simulations and biochemical studies, we identify the residues involved directly in binding of substrate to the inward-facing conformation of GlpT, thus defining the structural basis for the substrate-specificity of this transporter. The substrate binding mechanism involves protonation of a histidine residue at the binding site. Furthermore, our data suggest that the formation and breaking of inter- and intradomain salt bridges control the conformational change of the transporter that accompanies substrate translocation across the membrane. The mechanism we propose may be a paradigm for organophosphate:phosphate antiporters.