The connection between Mueller and Jones matrices of polarization optics

A necessary and sufficient condition for the 4 × 4 Mueller matrix to be derivable from the 2 × 2 Jones matrix is obtained. This condition allows one to determine if a given Mueller matrix describes a totally polarized system or a partially polarized (depolarizing) system. The result of Barakat is analysed in the light of this condition. A recently reported experimentally measured Mueller matrix is examined using this condition and is shown to represent a partially polarized system.

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.

Universality in the fast orientational relaxation near isotropic-nematic transition

Detailed molecular dynamics simulations of Lennard-Jones ellipsoids have been carried out to investigate the emergence of criticality in the single-particle orientational relaxation near the isotropic-nematic (IN) phase transition. The simulations show a sudden appearance of a power-law behavior in the decay of the second-rank orientational relaxation as the IN transition is approached. The simulated value of the power-law exponent is 0.56, which is larger than the mean-field value (0.5) but less than the observed value (0.63) and may be due to the finite size of the simulated system. The decay of the first-rank orientational time correlation function, on the other hand, is nearly exponential but its decay becomes very slow near the isotropic-nematic transition, The zero-frequency rotational friction, calculated from the simulated angular Velocity correlation function, shows a marked increase near the IN transition.

The sputtering temperature of a cooling cylindrical rod with an insulated core

By a standard application of Jones's method associated with the Wiener-Hopf technique an explicit solution is obtained for the temperature distribution inside a cylindrical rod with an insulated inner core when the rod is allowed to enter into a fluid of large extent with a uniform speed, and a simple integral expression is derived for the value of the sputtering temperature of the rod at the points of entry. Numerical results under certain special circumstances are also obtained and presented in the form of a table.

The Milky Way's Dark Matter Halo Appears To Be Lopsided

The atomic hydrogen gas (H I) disk in the outer region (beyond similar to 10 kpc from the center) of Milky Way can provide valuable information about the structure of the dark matter halo. The recent three-dimensional thickness map of the outer H I disk from the all sky 21 cm line Leiden/Argentine/Bonn survey, gives us a unique opportunity to investigate the structure of the dark matter halo of Milky Way in great detail. A striking feature of this new survey is the north-south (N-S) asymmetry in the thickness map of the atomic hydrogen gas. Assuming vertical hydrostatic equilibrium under the total potential of the Galaxy, we derive the model thickness map of the H I gas. We show that simple axisymmetric halo models, such as softened isothermal halo (producing a flat rotation curve with V-c similar to 220 km s(-1)) or any halo with density falling faster than the isothermal one, are not able to explain the observed radial variation of the gas thickness. We also show that such axisymmetric halos along with different H I velocity dispersion in the two halves, cannot explain the observed asymmetry in the thickness map. Amongst the nonaxisymmetric models, it is shown that a purely lopsided (m = 1, first harmonic) dark matter halo with reasonable H I velocity dispersion fails to explain the N-S asymmetry satisfactorily. However, we show that by superposing a second harmonic (m = 2) out of phase onto a purely lopsided halo, e. g., our best fit and more acceptable model A (with parameters epsilon(1)(h) = 0.2, epsilon(2)(h) = 0.18, and sigma(H I) = 8.5 km s(-1)) can provide an excellent fit to the observation and reproduce the N-S asymmetry naturally. The emerging picture of the asymmetric dark matter halo is supported by the. cold dark matter halos formed in the cosmological N-body simulation.

Asymptotic behaviour of point processes with Markov-dependent intervals

This splitting techniques for MARKOV chains developed by NUMMELIN (1978a) and ATHREYA and NEY (1978b) are used to derive an imbedded renewal process in WOLD's point process with MARKOV-correlated intervals. This leads to a simple proof of renewal theorems for such processes. In particular, a key renewal theorem is proved, from which analogues to both BLACKWELL's and BREIMAN's forms of the renewal theorem can be deduced.

The sputtering temperature of a cooling cylindrical rod with an insulated core

By a standard application of Jones's method associated with the Wiener-Hopf technique an explicit solution is obtained for the temperature distribution inside a cylindrical rod with an insulated inner core when the rod is allowed to enter into a fluid of large extent with a uniform speed, and a simple integral expression is derived for the value of the sputtering temperature of the rod at the points of entry. Numerical results under certain special circumstances are also obtained and presented in the form of a table.

Crystal and molecular structure of the quinoxaline antibiotic analog TANDEM (des-N-tetramethyltriostin A)

The crystal structure of TANDEM (des-N-tetramethyltriostin A), a synthetic analogue of the quinoxaline antibiotic triostin A, has been determined independently at -135 and 7 'C and refined to R values of 0.088 and 0.147, respectively. The molecule has approximate 2-fold symmetry, with the quinoxaline chromophores and the disulfide cross-bridge projecting from opposite sides of the peptide ring. The quinoxaline groups are nearly parallel to each other and separated by about 6.5 A. The peptide backbone resembles a distorted antiparallel 13 ribbon joined by intramolecular hydrogen bonds N-H(LVal)--O(L-Ala). At low temperatures, the TANDEM molecule is surrounded by a regular first- and second-order hydration sphere containing 14 independent water molecules. At room temperature, only the first-order hydration shell is maintained. Calculations of the interplanar separation of the quinoxaline groups as a function of their orientation with respect to the peptide ring support the viability of TANDEM to intercalate bifunctionally into DNA.

A simulation study of the mixing, atomic flow and velocity profiles of crystalline materials during sliding

Self-contained Non-Equilibrium Molecular Dynamics (NEMD) simulations using Lennard-Jones potentials were performed to identify the origin and mechanisms of atomic scale interfacial behavior between sliding metals. The mixing sequence and velocity profiles were compared via MD simulations for three cases, viz.: sell-mated, similar and hard-softvcrystal pairs. The results showed shear instability, atomic scale mixing, and generation of eddies at the sliding interface. Vorticity at the interface suggests that atomic flow during sliding is similar to fluid flow under Kelvin-Helmholtz instability and this is supported by velocity profiles from the simulations. The initial step-function velocity profile spreads during sliding. However the velocity profile does not change much at later stages of the simulation and it eventually stops spreading. The steady state friction coefficient during simulation was monitored as a function of sliding velocity. Frictional behavior can be explained on the basis of plastic deformation and adiabatic effects. The mixing layer growth kinetics was also investigated.

A Monte Carlo study of crystal structure transformations

The Metropolis algorithm has been generalized to allow for the variation of shape and size of the MC cell. A calculation using different potentials illustrates how the generalized method can be used for the study of crystal structure transformations. A restricted MC integration in the nine dimensional space of the cell components also leads to the stable structure for the Lennard-Jones potential.