Isospin- and momentum-dependent interaction in isospin-dependent quantum molecular dynamics

An isospin-dependent quantum molecular dynamical model (IQMD) is developed, with the isospin degree of freedom in the momentum-dependent interaction(MDI) included in IQMD, to obtain an isospin- and momentum-dependent interaction (IMDI) in IQMD. We investigate the effect of IMDI on the isospin fractionation ratio and its dynamical mechanism in the intermediate energy heavy ion collisions. It is found that the IMDI induces the significant reductions in the isospin fractionation ratio for all of beam energies, impact parameters, neutron-proton ratios and mass number of colliding systems. However, the strong dependence of isospin fractionation ratio on the symmetrical potential is preserved, with the isospin degree of freedom included in the MDI, i.e. the isospin fractionation ratio is still a good probe for extracting the information about the equation of state of isospin asymmetrical nuclear matter.

Accurate and efficient non-adiabatic quantum molecular dynamics approach for laser-matter interactions

A non-adiabatic quantum molecular dynamics approach for treating the interaction of matter with intense, short-duration laser pulses is developed. This approach, which is parallelized to run on massively-parallel supercomputers, is shown to be both accurate and efficient. Illustrative results are presented for harmonic generation occurring in diatomic molecules using linearly polarized laser pulses.

Quantum molecular dynamics and ab initio studies of the crystal structure of hydrogen and deuterium

Hydrogen isotopes play a critical role both in inertial and magnetic confinemen Nuclear Fusion. Since the preferent fuel needed for this technology is a mixture of deuterium and tritium. The study of these isotopes particularly at very low temperatures carries a technological interest in other applications. The present line promotes a deep study on the structural configuration that hydrogen and deuterium adopt at cryogenic temperatures and at high pressures. Typical conditions occurring in present Inertial Fusion target designs. Our approach is aims to determine the crystal structure characteristics, phase transitions and other parameters strongly correlated to variations of temperature and pressure.

A molecular dynamics model of the atomic structure of dysprosium alumino-phosphate glass

Molecular dynamics (MD) has been used to identify the relative distribution of dysprosium in the phosphate glass DyAl0.30P3.05O9.62. The MD model has been compared directly with experimental data obtained from neutron diffraction to enable a detailed comparison beyond the total structure factor level. The MD simulation gives Dy ... Dy correlations at 3.80(5) and 6.40(5) angstrom with relative coordination numbers of 0.8(1) and 7.3(5), thus providing evidence of minority rare-earth clustering within these glasses. The nearest neighbour Dy-O peak occurs at 2.30 angstrom with each Dy atom having on average 5.8 nearest neighbour oxygen atoms. The MD simulation is consistent with the phosphate network model based on interlinked PO4 tetrahedra where the addition of network modifiers Dy3+ depolymerizes the phosphate network through the breakage of P-(O)-P bonds whilst leaving the tetrahedral units intact. The role of aluminium within the network has been taken into explicit account, and A1 is found to be predominantly (78 tetrahedrally coordinated. In fact all four A1 bonds are found to be to P (via an oxygen atom) with negligible amounts of Al-O-Dy bonds present. This provides an important insight into the role of Al additives in improving the mechanical properties of these glasses.

MOLECULAR DYNAMICS STUDIES OF SIMPLE MODEL FLUIDS AND WATER CONFINED IN CARBON NANOTUBE

Molecular Dynamics (MD) simulation is one of the most important computational techniques with broad applications in physics, chemistry, chemical engineering, materials design and biological science. Traditional computational chemistry refers to quantum calculations based on solving Schrodinger equations. Later developed Density Functional Theory (DFT) based on solving Kohn-Sham equations became the more popular ab initio calculation technique which could deal with ~1000 atoms by explicitly considering electron interactions. In contrast, MD simulation based on solving classical mechanics equations of motion is a totally different technique in the field of computational chemistry. Electron interactions were implicitly included in the empirical atom-based potential functions and the system size to be investigated can be extended to ~106 atoms. The thermodynamic properties of model fluids are mainly determined by macroscopic quantities, like temperature, pressure, density. The quantum effects on thermodynamic properties like melting point, surface tension are not dominant. In this work, we mainly investigated the melting point, surface tension (liquid-vapor and liquid-solid) of model fluids including Lennard-Jones model, Stockmayer model and a couple of water models (TIP4P/Ew, TIP5P/Ew) by means of MD simulation. In addition, some new structures of water confined in carbon nanotube were discovered and transport behaviors of water and ions through nano-channels were also revealed.

Length-scale dependent mechanical properties of Al-Cu eutectic alloy: Molecular dynamics based model and its experimental verification

This paper attempts to gain an understanding of the effect of lamellar length scale on the mechanical properties of two-phase metal-intermetallic eutectic structure. We first develop a molecular dynamics model for the in-situ grown eutectic interface followed by a model of deformation of Al-Al2Cu lamellar eutectic. Leveraging the insights obtained from the simulation on the behaviour of dislocations at different length scales of the eutectic, we present and explain the experimental results on Al-Al2Cu eutectic with various different lamellar spacing. The physics behind the mechanism is further quantified with help of atomic level energy model for different length scale as well as different strain. An atomic level energy partitioning of the lamellae and the interface regions reveals that the energy of the lamellae core are accumulated more due to dislocations irrespective of the length-scale. Whereas the energy of the interface is accumulated more due to dislocations when the length-scale is smaller, but the trend is reversed when the length-scale is large beyond a critical size of about 80 nm. (C) 2014 Author(s).

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.

Shell effect and capture cross sections in the synthesis of superheavy nuclei

The shell effect is included in the improved isospin dependent quantum molecular dynamics model in which the shell correction energy of the system is calculated by using the deformed two-center shell model. A switch function is introduced to connect the shell correction energy of the projectile and the target with that of the compound nucleus during the dynamical fusion process. It is found that the calculated capture cross sections reproduce the experimental data quantitatively at the energy near the Coulomb barrier. The capture cross sections for reaction (35) (80) Br + (82) (208) Pb -> (117) (288) X are also calculated and discussed.

Average property of fragmentation reaction and momentum dissipation induced by halo-nuclei in heavy ion collisions

We study systematically the average property of fragmentation reaction and momentum dissipation induced by halo-nuclei in intermediate energy heavy ion collisions for different colliding systems and different beam energies within the isospin dependent quantum molecular dynamics model (IQMD). This study is based on the extended halo-nucleus density distributions, which indicates the average property of loosely inner halo nucleus structure, because the interaction potential and in-medium nucleon-nucleon cross section in IQMD model depend on the density distribution. In order to study the average properties of fragmentation reaction and momentum dissipation induced by halo-nuclei we also compare the results for the halo-nuclear colliding systems with those for corresponding stable colliding systems with same mass under the same incident channel condition. We find that the effect of extended halo density distribution on the fragment multiplicity and nuclear stopping (momentum dissipation) are important for the different beam energies and different colliding systems. For example the extended halo density distributions increase the fragment multiplicity but decrease the nuclear stopping for all of incident channel conditions in this paper.

Dependence of isoscaling parameters on nucleon-nucleon cross section and momentum-dependent interaction

Influences of the isospin dependence of the in-medium nucleon-nucleon cross section and the momentum-dependant interaction (MDI) on the isotope scaling are investigated by using the isospin-dependent quantum molecular dynamics model (IQMD). The results show that both the isospin dependence of the in-medium nucleon-nucleon cross section and the momentum-dependent interaction affect the isoscaling parameters appreciably and independently. The influence caused by the isospin dependence of two-body collision is relatively larger than that from the MDI in the mean field. Aiming at exploring the implication of isoscaling behaviour, which the statistical equilibrium in the reaction is reached, the statistical properties in the mass distribution and the kinetic energy distribution of the fragments simulated by IQMD are presented.

Influence of isospin dependence of nucleon-nucleon cross section and momentum dependent interaction on the isoscaling in intermediate energy collisions

The influences of the isospin dependent in-medium nucleon-nucleon cross section and the MomentumDependent Interaction(MDI) on the isotope scaling have been investigated within the Isospin dependent Quantum Molecular Dynamics Model(IQMD). The results show that both the isospin dependent in-medium nucleon-nucleon cross section and the momentum interaction reduce the isoscaling parameter a appreciably, which means they decrease the dependence of yield ratios of two systems on the isospin difference between two systems.