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.

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).

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.

A preliminary genetic map of Zhikong scallop (Chlamys farreri Jones et Preston 1904)

Zhikong scallop (Chlamys farreri Jones et Preston 1904) is one of the most important aquaculture species in China. The development of a genetic linkage map would provide a powerful tool for the genetic improvement of this species. Amplified fragment length polymorphism (AFLP) is a PCR-based technique that has proven to be powerful in genome fingerprinting and mapping, and population analysis. Genetic maps of C. farreri were constructed using AFLP markers and a full-sib family with 60 progeny. A total of 503 segregating AFLP markers were obtained, with 472 following the Mendelian segregation ratio of 1:1 and 31 markers showing significant (P< 0.05) segregation distortion. The male map contained 166 informative AFLP markers in 23 linkage groups covering 2468 cM. The average distance between markers was 14.9 cM. The female genetic map consisted of 198 markers in 25 linkage groups spanning 3130 cM with an average inter-marker spacing of 15.8 cM. DNA polymorphisms that segregated in a 3:1 ratio as well as the AFLP markers that were heterozygous in both parents were included to construct combined linkage genetic map. Five shared linkage groups, ranging from 61.1 to 162.5 cM, were identified between the male and female maps, covering 431 cM. Amplified fragment length polymorphism markers appeared to be evenly distributed within the linkage groups. Although preliminary, these maps provide a starting point for the mapping of the functional genes and quantitative trait loci in C. farreri.

Construction of AFLP-based genetic linkage map for Zhikong scallop, Chlamys farreri Jones et Preston and mapping of sex-linked markers

Zhikong scallop (Chlamys farreri) is an economically important aquaculture species in China; however, frequent mass mortality seriously affects the development of its industry. Genetic linkage map is useful for genetic improvement and selective breeding of C. farreri. Linkage maps were constructed using an intraspecific F-1 cross and amplified fragment length polymorphism (AFLP) markers. Thirty-two selected AFLP primer combinations produced 545 AFLP markers that were polymorphic in either of the parents and segregated in the progeny. Of these segregating markers, 166 were mapped to 19 linkage groups of the female framework map, covering a total of 1503.9 cM, with an average marker spacing of 10.2 cM; and 197 markers were assigned to 20 linkage groups of the male map, covering a total of 1630.7 cM, with 9.2 cM per marker. A sex-linked marker was mapped on the female map with zero recombination and a LOD of 27.3. The genetic length of C farreri genome was estimated as 1889.0 cM for the female and 1995.9 cM for the male. The coverage of the framework map was calculated as 79.6% for the female and 81.7% for the male. When the triplets and doublets were considered, the observed length of the map was calculated as 1610.2 cM with coverage of 85.2% for the female, and 1880.5 cM with coverage of 94.2% for the male. The genetic maps presented here will serve as a basis for the construction of a high-resolution genetic map and mapping of economically important genes. (C) 2004 Published by Elsevier B.V.

Triploid and tetraploid Zhikong scallop, Chlamys farreri Jones et Preston, produced by inhibiting polar body I

Triploid scallops are valuable for aquaculture because of their enlarged adductor muscle, and tetraploids are important for the commercial production of triploids. We tested tetraploid induction in the zhikong scallop by inhibiting polar body I in newly fertilized eggs. The ploidy of resultant embryos was determined by chromosome counting at 2- to 4-cell stage and by flow cytometry thereafter. Embryos from the control groups were mostly diploids (79%), along with some aneuploids. Embryos from the treated groups were 13% diploids, 18% triploids, 26% tetraploids, 13% pentaploids, and 36% aneuploids. Tetraploids, pentaploids, and most aneuploids suffered heavy mortality during the first week and became undetectable among the larvae at day 14. Five tetraploids (2%) were found among a sample of 267 spat from one of the replicates, and none was detected at day 450. The adductor muscle of triploid scallops was 44% heavier (P < .01) than that of diploids, confirming the value of the triploid technology in this species.

Construction and characterization of two bacterial artificial chromosome libraries of Zhikong scallop, Chlamys farreri Jones et Preston, and identification of BAC clones containing the genes involved in its innate immune system

Large-insert bacterial artificial chromosome (BAC) libraries are necessary for advanced genetics and genomics research. To facilitate gene cloning and characterization, genome analysis, and physical mapping of scallop, two BAC libraries were constructed from nuclear DNA of Zhikong scallop, Chlamys farreri Jones et Preston. The libraries were constructed in the BamHI and MboI sites of the vector pECBAC1, respectively. The BamHI library consists of 73,728 clones, and approximately 99% of the clones contain scallop nuclear DNA inserts with an average size of 110 kb, covering 8.0x haploid genome equivalents. Similarly, the MboI library consists of 7680 clones, with an average insert of 145 kb and no insert-empty clones, thus providing a genome coverage of 1.1x. The combined libraries collectively contain a total of 81,408 BAC clones arrayed in 212 384-well microtiter plates, representing 9.1x haploid genome equivalents and having a probability of greater than 99% of discovering at least one positive clone with a single-copy sequence. High-density clone filters prepared from a subset of the two libraries were screened with nine pairs of Overgos designed from the cDNA or DNA sequences of six genes involved in the innate immune system of mollusks. Positive clones were identified for every gene, with an average of 5.3 BAC clones per gene probe. These results suggest that the two scallop BAC libraries provide useful tools for gene cloning, genome physical mapping, and large-scale sequencing in the species.

Growth, survival and immune activity of scallops, Chlamys farreri Jones et Preston, compared between suspended and bottom culture in Haizhou Bay, China

We examined the growth, survival and immune response of the scallop, Chlamys farreri, during a 1-year period in deep water of Haizhou Bay. Scallops were cultured using two methods: (1) in lantern nets at a 5 m depth and (2) in a bottom culture system (sleeves) on the seabed at about a 25 m depth. Shell heights, meat dry weight and immune activities in the haemolymph (superoxide dismutase and myeloperoxidase) were measured bimonthly or quarterly from July 2007 to June 2008. Survival was measured at the end of the study and environmental parameters in the experimental layers were monitored during the experiment. The growth and immune activities of scallops were lower when the water temperature was high, which was consistent with the main mortality occurring in summer. The growth and immunity of scallops were higher in the suspended culture than in the bottom culture during the experiment, with the exception of shell growth during the last study period. Survival of scallops in the suspended culture (54.6 +/- 12.3%) was significantly lower than that in the bottom culture (86.8 +/- 3.5%) at the end of this study. We conclude from our results that the high mortality of C. farreri can be prevented by culturing them in a bottom culture system before November of the first year, and then transferring them to a suspended culture to improve scallop production.

Studies on mass summer mortality of cultured zhikong scallops (Chlamys farreri Jones et Preston) in China

Mass mortalities of cultured zhikong scallops (Chlamys farreri) have occurred each summer in most culture areas of northern China since 1996. Among the hypothesized causes are high culture density, infectious disease and genetic inbreeding. To investigate these potential agents, C. farreri were deployed at three densities (low, medium and high) at three sites (Jiaonan, Penglai and Yantai) in the summer of 2000. Scallops were sampled for survival, growth and histopathology before, during and after a mortality episode. Most of the mortality occurred in July and August, during and toward the later part of the spawning season, when water temperature reached 23-26 degrees C. Final cumulative mortalities reached 85% to 90% at all three sites. Scallops in the medium and high densities had higher initial death rates than did those at the low density. High densities also inhibited growth. Ciliates from the genus Trichodina, larvae of various organisms and anomalous secretions were observed in sections of the gill cavity, with highest prevalence during and at the end of the mortality period. Prokaryotic inclusion bodies were found in the soft tissues, but their prevalence was low and apparently without correlation with mortalities. Genetic analysis with random amplified polymorphic DNA markers showed slightly lower heterozygosity in the cultured stocks (0.301) than in the wild stocks (0.331). It is possible that the mortalities are caused by a combination of several factors such as stress associated with reproduction, high temperature, overcrowding and poor circulation in the growout cages, opportunistic invaders or pathogens, and possibly inbreeding. (c) 2005 Elsevier B.V. All rights reserved.