Translational and rotational motion in molecular liquids: A computer simulation study of Lennard-Jones ellipsoids

In order to understand the translational and rotational motion in dense molecular liquids, detailed molecular dynamics simulations of Lennard-Jones ellipsoids have been carried out for three different values of the aspect ratio kappa. For ellipsoids with an aspect ratio equal to 2, the product of the translational diffusion coefficient (D-T) and the average orientational correlation time of the l-th rank harmonics (tau(lR)), converges to a nearly constant value at high density. Surprisingly, this density independent value of D-T tau(lR) is within 5% of the hydrodynamic prediction with the slip boundary condition. This is despite the fact that both D-T and tau(lR) themselves change nearly by an order of magnitude in the density range considered, and the rotational correlation function itself is strongly nonexponential. For small aspect ratios (kappa less than or equal to 1.5), the rotational correlation function remains largely Gaussian even at a very large density, while for a large aspect ratio (kappa greater than or equal to 3), the transition to the nematic liquid-crystalline phase precludes the hydrodynamic regime. Thus, the rotational dynamics of ellipsoids show great sensitivity to the aspect ratio. At low density, tau(lR) goes through a minimum value, indicating the role of interactions in enhancing the rate of orientational relaxation. (C) 1997 American Institute of Physics. [S0021-9606(97)50142-5].

Dependence of diffusivity on density and solute diameter in liquid phase: A molecular dynamics study of Lennard-Jones system

Investigations into the variation of self-diffusivity with solute radius, density, and degree of disorder of the host medium is explored. The system consists of a binary mixture of a relatively smaller sized solute, whose size is varied and a larger sized solvent interacting via Lennard-Jones potential. Calculations have been performed at three different reduced densities of 0.7, 0.8, and 0.933. These simulations show that diffusivity exhibits a maximum for some intermediate size of the solute when the solute diameter is varied. The maximum is found at the same size of the solute at all densities which is at variance with the prediction of the levitation effect. In order to understand this anomaly, additional simulations were carried out in which the degree of disorder has been varied while keeping the density constant. The results show that the diffusivity maximum gradually disappears with increase in disorder. Disorder has been characterized by means of the minimal spanning tree. Simulations have also been carried out in which the degree of disorder is constant and only the density is altered. The results from these simulations show that the maximum in diffusivity now shifts to larger distances with decrease in density. This is in agreement with the changes in void and neck distribution with density of the host medium. These results are in excellent agreement with the predictions of the levitation effect. They suggest that the effect of disorder is to shift the maximum in diffusivity towards smaller solute radius while that of the decrease in density is to shift it towards larger solute radius. Thus, in real systems where the degree of disorder is lower at higher density and vice versa, the effect due to density and disorder have opposing influences. These are confirmed by the changes seen in the velocity autocorrelation function, self part of the intermediate scattering function and activation energy. (C) 2012 American Institute of Physics. http://dx.doi.org/10.1063/1.3701619]

Solid-liquid transition in polydisperse Lennard-Jones systems

We study melting of a face-centered crystalline solid consisting of polydisperse Lennard-Jones spheres with Gaussian polydispersity in size. The phase diagram reproduces the existence of a nearly temperature invariant terminal polydispersity (delta(t) similar or equal to 0.11), with no signature of reentrant melting. The absence of reentrant melting can be attributed to the influence of the attractive part of the potential upon melting. We find that at terminal polydispersity the fractional density change approaches zero, which seems to arise from vanishingly small compressibility of the disordered phase. At constant temperature and volume fraction the system undergoes a sharp transition from crystalline solid to the disordered amorphous or fluid state with increasing polydispersity. This has been quantified by second- and third-order rotational invariant bond orientational order, as well as by the average inherent structure energy. The translational order parameter also indicates a similar sharp structural change at delta similar or equal to 0.09 in case of T* = 1.0, phi = 0.58. The free energy calculation further supports the sharp nature of the transition. The third-order rotationally invariant bond order shows that with increasing polydispersity, the local cluster favors a more icosahedral arrangement and the system loses its local crystalline symmetry. Interestingly, the value of structure factor S(k) of the amorphous phase at delta similar or equal to 0.10 (just beyond the solid-liquid transition density at T* = 1) becomes 2.75, which is below the value of 2.85 required for freezing given by the empirical Hansen-Verlet rule of crystallization, well known in the theory of freezing.

Use of polydispersity index as control parameter to study melting/freezing of Lennard-Jones system: Comparison among predictions of bifurcation theory with Lindemann criterion, inherent structure analysis and Hansen-Verlet rule

Using polydispersity index as an additional order parameter we investigate freezing/melting transition of Lennard-Jones polydisperse systems (with Gaussian polydispersity in size), especially to gain insight into the origin of the terminal polydispersity. The average inherent structure (IS) energy and root mean square displacement (RMSD) of the solid before melting both exhibit quite similar polydispersity dependence including a discontinuity at solid-liquid transition point. Lindemann ratio, obtained from RMSD, is found to be dependent on temperature. At a given number density, there exists a value of polydispersity index (delta (P)) above which no crystalline solid is stable. This transition value of polydispersity(termed as transition polydispersity, delta (P) ) is found to depend strongly on temperature, a feature missed in hard sphere model systems. Additionally, for a particular temperature when number density is increased, delta (P) shifts to higher values. This temperature and number density dependent value of delta (P) saturates surprisingly to a value which is found to be nearly the same for all temperatures, known as terminal polydispersity (delta (TP)). This value (delta (TP) similar to 0.11) is in excellent agreement with the experimental value of 0.12, but differs from hard sphere transition where this limiting value is only 0.048. Terminal polydispersity (delta (TP)) thus has a quasiuniversal character. Interestingly, the bifurcation diagram obtained from non-linear integral equation theories of freezing seems to provide an explanation of the existence of unique terminal polydispersity in polydisperse systems. Global bond orientational order parameter is calculated to obtain further insights into mechanism for melting.

A Lennard-Jones Embedded-Atom Potential and Its Application to the Study of Melting

The nearest-neighbour Lennard-Jones potential from the embedded-atom method is extended to a form that includes more than nearest neighbours. The model has been applied to study melting with molecular dynamics. The calculated melting point, fractional volume change on melting, heat of fusion and linear coefficients of thermal expansion are in good agreement with experimental data. We have found that the second and third neighbours influence the melting point distinctly.

1. Ultrasonic studies of binary liquid structure in the critical region. Theory and experiment for the 2,6-lutidine/water system. 2. Hartree-Fock calculations of electric polarizabilities of some simple atoms and molecules, and their practicality. 3. Calculation of vibrational transition probabilities in collinear atom-diatom and diatom-diatom collisions with Lennard-Jones interaction

Part 1. Many interesting visual and mechanical phenomena occur in the critical region of fluids, both for the gas-liquid and liquid-liquid transitions. The precise thermodynamic and transport behavior here has some broad consequences for the molecular theory of liquids. Previous studies in this laboratory on a liquid-liquid critical mixture via ultrasonics supported a basically classical analysis of fluid behavior by M. Fixman (e. g., the free energy is assumed analytic in intensive variables in the thermodynamics)--at least when the fluid is not too close to critical. A breakdown in classical concepts is evidenced close to critical, in some well-defined ways. We have studied herein a liquid-liquid critical system of complementary nature (possessing a lower critical mixing or consolute temperature) to all previous mixtures, to look for new qualitative critical behavior. We did not find such new behavior in the ultrasonic absorption ascribable to the critical fluctuations, but we did find extra absorption due to chemical processes (yet these are related to the mixing behavior generating the lower consolute point). We rederived, corrected, and extended Fixman's analysis to interpret our experimental results in these more complex circumstances. The entire account of theory and experiment is prefaced by an extensive introduction recounting the general status of liquid state theory. The introduction provides a context for our present work, and also points out problems deserving attention. Interest in these problems was stimulated by this work but also by work in Part 3.

Part 2. Among variational theories of electronic structure, the Hartree-Fock theory has proved particularly valuable for a practical understanding of such properties as chemical binding, electric multipole moments, and X-ray scattering intensity. It also provides the most tractable method of calculating first-order properties under external or internal one-electron perturbations, either developed explicitly in orders of perturbation theory or in the fully self-consistent method. The accuracy and consistency of first-order properties are poorer than those of zero-order properties, but this is most often due to the use of explicit approximations in solving the perturbed equations, or to inadequacy of the variational basis in size or composition. We have calculated the electric polarizabilities of H₂, He, Li, Be, LiH, and N₂ by Hartree-Fock theory, using exact perturbation theory or the fully self-consistent method, as dictated by convenience. By careful studies on total basis set composition, we obtained good approximations to limiting Hartree-Fock values of polarizabilities with bases of reasonable size. The values for all species, and for each direction in the molecular cases, are within 8% of experiment, or of best theoretical values in the absence of the former. Our results support the use of unadorned Hartree-Pock theory for static polarizabilities needed in interpreting electron-molecule scattering data, collision-induced light scattering experiments, and other phenomena involving experimentally inaccessible polarizabilities.

Part 3. Numerical integration of the close-coupled scattering equations has been carried out to obtain vibrational transition probabilities for some models of the electronically adiabatic H₂-H₂ collision. All the models use a Lennard-Jones interaction potential between nearest atoms in the collision partners. We have analyzed the results for some insight into the vibrational excitation process in its dependence on the energy of collision, the nature of the vibrational binding potential, and other factors. We conclude also that replacement of earlier, simpler models of the interaction potential by the Lennard-Jones form adds very little realism for all the complication it introduces. A brief introduction precedes the presentation of our work and places it in the context of attempts to understand the collisional activation process in chemical reactions as well as some other chemical dynamics.

Learning Multirobot Hose Transportation and Deployment by Distributed Round-Robin Q-Learning

Multi-Agent Reinforcement Learning (MARL) algorithms face two main difficulties: the curse of dimensionality, and environment non-stationarity due to the independent learning processes carried out by the agents concurrently. In this paper we formalize and prove the convergence of a Distributed Round Robin Q-learning (D-RR-QL) algorithm for cooperative systems. The computational complexity of this algorithm increases linearly with the number of agents. Moreover, it eliminates environment non sta tionarity by carrying a round-robin scheduling of the action selection and execution. That this learning scheme allows the implementation of Modular State-Action Vetoes (MSAV) in cooperative multi-agent systems, which speeds up learning convergence in over-constrained systems by vetoing state-action pairs which lead to undesired termination states (UTS) in the relevant state-action subspace. Each agent's local state-action value function learning is an independent process, including the MSAV policies. Coordination of locally optimal policies to obtain the global optimal joint policy is achieved by a greedy selection procedure using message passing. We show that D-RR-QL improves over state-of-the-art approaches, such as Distributed Q-Learning, Team Q-Learning and Coordinated Reinforcement Learning in a paradigmatic Linked Multi-Component Robotic System (L-MCRS) control problem: the hose transportation task. L-MCRS are over-constrained systems with many UTS induced by the interaction of the passive linking element and the active mobile robots.

A new species and a new record of Sphaeromopsis Holdich and Jones, 1973 from Pakistan waters (Isopoda, Sphaeromatidae)

Sphaeromopsis minutus, sp. nov., the sixth species of Splitteromopsis Holdich and Jones is described and figured in detail. The genus Sphaeromopsis and species S. serriguberna Holdich and Harrison, 1981 is recorded for the first time from northern Arabian Sea, Karachi (Pakistan).

Room-temperature dual-wavelength erbium-doped fibre laser based on a sampled fibre Bragg grating and a photonic Robin Hood

With the assistance of a kind of photonic Robin Hood that is originated from four-wave mixing in a dispersion-flattened high-nonlinearity photonic-crystal fibre, a novel dual-wavelength erbium-doped fibre (EDF) laser is proposed and demonstrated by using a sampled fibre Bragg grating. The experiments show that, due to the contribution of the photonic Robin Hood, the proposed fibre laser has the advantage of excellent uniformity, high stability and stable operation at room temperature. Our dual-wavelength EDF laser has the unique merit that the wavelength spacing remains unchanged when tuning the two wavelengths of laser, and this laser is simpler and more stable than the laser reported by Liu et al. [Opt. Express, 13 142 (2005)].

Chromosomal localization and molecular marker development of the lipopolysaccharide and beta-1,3-glucan binding protein gene in the Zhikong scallop Chlamys farreri (Jones et Preston) (Pectinoida, Pectinidae)

Zhikong scallop Chlamys farreri(Jones et Preston) is an economically important species in China. Understanding its immune system would be of great help in controlling diseases. In the present study, an important immunity-related gene, the Lipopolysaccharide and Beta-1,3-glucan Binding Protein (LGBP) gene, was located on C. farreri chromosomes by mapping several lgbp-containing BAC clones through fluorescence in situ hybridization (FISH). Through the localization of various BAC clones, it was shown that only one locus of this gene existed in the genome of C. farreri, and that this was located on the long arm of a pair of homologous chromosomes. Molecular markers, consisting of eight single nucleotide polymorphism (SNPs) markers and one insertion-deletion (indel), were developed from the LGBP gene. Indel marker testing in an F1 family revealed slightly distorted segregation (p = 0.0472). These markers can be used to map the LGBP gene to the linkage map and assign the linkage group to the corresponding chromosome. Segregation distortion of the indel marker indicated genes with deleterious alleles might exist in the surrounding region of the LGBP gene.