A comment on the consistency of truncated nonlinear integral equation based theories of freezing

We report the results of two studies of aspects of the consistency of truncated nonlinear integral equation based theories of freezing: (i) We show that the self-consistent solutions to these nonlinear equations are unfortunately sensitive to the level of truncation. For the hard sphere system, if the Wertheim–Thiele representation of the pair direct correlation function is used, the inclusion of part but not all of the triplet direct correlation function contribution, as has been common, worsens the predictions considerably. We also show that the convergence of the solutions found, with respect to number of reciprocal lattice vectors kept in the Fourier expansion of the crystal singlet density, is slow. These conclusions imply great sensitivity to the quality of the pair direct correlation function employed in the theory. (ii) We show the direct correlation function based and the pair correlation function based theories of freezing can be cast into a form which requires solution of isomorphous nonlinear integral equations. However, in the pair correlation function theory the usual neglect of the influence of inhomogeneity of the density distribution on the pair correlation function is shown to be inconsistent to the lowest order in the change of density on freezing, and to lead to erroneous predictions. The Journal of Chemical Physics is copyrighted by The American Institute of Physics.

A comment on the consistency of truncated nonlinear integral equation based theories of freezing

We report the results of two studies of aspects of the consistency of truncated nonlinear integral equation based theories of freezing: (i) We show that the self-consistent solutions to these nonlinear equations are unfortunately sensitive to the level of truncation. For the hard sphere system, if the Wertheim–Thiele representation of the pair direct correlation function is used, the inclusion of part but not all of the triplet direct correlation function contribution, as has been common, worsens the predictions considerably. We also show that the convergence of the solutions found, with respect to number of reciprocal lattice vectors kept in the Fourier expansion of the crystal singlet density, is slow. These conclusions imply great sensitivity to the quality of the pair direct correlation function employed in the theory. (ii) We show the direct correlation function based and the pair correlation function based theories of freezing can be cast into a form which requires solution of isomorphous nonlinear integral equations. However, in the pair correlation function theory the usual neglect of the influence of inhomogeneity of the density distribution on the pair correlation function is shown to be inconsistent to the lowest order in the change of density on freezing, and to lead to erroneous predictions. The Journal of Chemical Physics is copyrighted by The American Institute of Physics.

Solvation dynamics in a Brownian dipolar lattice. Comparison between computer simulation and various molecular theories of solvation dynamics

Several recent theoretical and computer simulation studies have considered solvation dynamics in a Brownian dipolar lattice which provides a simple model solvent for which detailed calculations can be carried out. In this article a fully microscopic calculation of the solvation dynamics of an ion in a Brownian dipolar lattice is presented. The calculation is based on the non‐Markovian molecular hydrodynamic theory developed recently. The main assumption of the present calculation is that the two‐particle orientational correlation functions of the solid can be replaced by those of the liquid state. It is shown that such a calculation provides an excellent agreement with the computer simulation results. More importantly, the present calculations clearly demonstrate that the frequency‐dependent dielectric friction plays an important role in the long time decay of the solvation time correlation function. We also find that the present calculation provides somewhat better agreement than either the dynamic mean spherical approximation (DMSA) or the Fried–Mukamel theory which use the simulated frequency‐dependent dielectric function. It is found that the dissipative kernels used in the molecular hydrodynamic approach and in the Fried–Mukamel theory are vastly different, especially at short times. However, in spite of this disagreement, the two theories still lead to comparable results in good agreement with computer simulation, which suggests that even a semiquantitatively accurate dissipative kernel may be sufficient to obtain a reliable solvation time correlation function. A new wave vector and frequency‐dependent dissipative kernel (or memory function) is proposed which correctly goes over to the appropriate expressions in both the single particle and the collective limits. This form is expected to lead to better results than all the existing descriptions.

Ultrafast underdamped solvation: Agreement between computer simulation and various theories of solvation dynamics

A theoretical analysis of the three currently popular microscopic theories of solvation dynamics, namely, the dynamic mean spherical approximation (DMSA), the molecular hydrodynamic theory (MHT), and the memory function theory (MFT) is carried out. It is shown that in the underdamped limit of momentum relaxation, all three theories lead to nearly identical results when the translational motions of both the solute ion and the solvent molecules are neglected. In this limit, the theoretical prediction is in almost perfect agreement with the computer simulation results of solvation dynamics in the model Stockmayer liquid. However, the situation changes significantly in the presence of the translational motion of the solvent molecules. In this case, DMSA breaks down but the other two theories correctly predict the acceleration of solvation in agreement with the simulation results. We find that the translational motion of a light solute ion can play an important role in its own solvation. None of the existing theories describe this aspect. A generalization of the extended hydrodynamic theory is presented which, for the first time, includes the contribution of solute motion towards its own solvation dynamics. The extended theory gives excellent agreement with the simulations where solute motion is allowed. It is further shown that in the absence of translation, the memory function theory of Fried and Mukamel can be recovered from the hydrodynamic equations if the wave vector dependent dissipative kernel in the hydrodynamic description is replaced by its long wavelength value. We suggest a convenient memory kernel which is superior to the limiting forms used in earlier descriptions. We also present an alternate, quite general, statistical mechanical expression for the time dependent solvation energy of an ion. This expression has remarkable similarity with that for the translational dielectric friction on a moving ion.

Structural and electrical transport properties of Al-Cu-Cr quasicrystals

Transport properties of quasicrystals in rapidly solidified as well as heat-treated Al65CU20Cr15 alloys were studied over a wide temperature range as a function of structure and microstructure. The characterization was done using x-ray diffraction, transmission electron microscopy and differential scanning calorimetry. Particular attention was paid to primitive to face-centered quasicrystalline transformation which occurs on annealing and the effect of microstructures on the transport behavior. The temperature dependence of resistivity is found to depend crucially on the microstructure of the alloy. Further, ordering enhances the negative temperature coefficient of resistivity. The low-temperature (T less than or equal to 25 K) resistivity of Al65Cu20Cr15 has been compared with that of Al63.5Cu24.5Fe12 alloy. In this region p(T) can be well described by a root T contribution arising from electron-electron interaction. We discuss our results in view of current theories.

Spontaneous generation of optical activity in urea inclusion compounds vis-a-vis current theories

A novel test of recent theories of the origin of optical activity has been designed based on the inclusion of certain alkyl 2-methylhexanoates into urea channels.

Structural and lithium ion transport studies in sol-gel-prepared lithium silicophosphate glasses

Lithium silicophosphate glasses have been prepared by a sol-gel route over a wide range of compositions. Their structural and electrical properties have been investigated. Infrared spectroscopic studies show the presence of hydroxyl groups attached to Si and P. MAS NMR investigations provide evidence for the presence of different phosphatic units in the structure. The variations of de conductivities at 423 K and activation energies have been studied as a function of composition, and both exhibit an increasing trend with the ratio of nonbridging oxygen to bridging oxygen in the structure. Ac conductivity behavior shows that the power law exponent, s, is temperature dependent and exhibits a minimum. Relaxation behavior has been examined in detail using an electrical modulus formalism, and modulus data were fitted to Kohlraush-William-Watts stretched exponential function. A structural model has been proposed and the unusual properties exhibited by this unique system of glasses have been rationalized using this model. Ion transport in these glasses appears to be confined to unidimensional conduits defined by modified phosphate chains and interspersed with unmodified silica units.

Structural correlations in GexSe1-x glasses - a neutron diffraction study

Neutron diffraction measurement is carried out on GexSe1-x glasses, where 0.1 less than or equal to x less than or equal to 0.4, in a Q interval of 0.55-13.8 Angstrom(-1). The first sharp diffraction peak (FSDP) in the structure factor, S(Q), shows a systematic increase in the intensity and shifts to a lower Q with increasing Ge concentration. The coherence length of FSDP increases with x and becomes maximum for 0.33 less than or equal to x less than or equal to 0.4. The Monte-Carlo method, due to Soper, is used to generate S(Q) and also the pair correlation function, g(r). The generated S(Q) is in agreement with the experimental data for all x. Analysis of the first four peaks in the total correlation function, T(r), shows that the short range order in GeSe2 glass is due to Ge(Se-1/2)(4) tetrahedra, in agreement with earlier reports. Se-rich glasses contain Se-chains which are cross-linked with Ge(Se-1/2)(4) tetrahedra. Ge-2(Se-1/2)(6) molecular units are the basic structural units in Ge-rich, x = 0.4, glass. For x = 0.2, 0.33 and 0.4 there is evidence for some of the tetrahedra being in an edge-shared configuration. The number of edge-shared tetrahedra in these glasses increase with increasing Ge content.

Orientational relaxation in a random dipolar lattice: Wave-number and frequency dependence

In order to understand the role of translational modes in the orientational relaxation in dense dipolar liquids, we have carried out a computer ''experiment'' where a random dipolar lattice was generated by quenching only the translational motion of the molecules of an equilibrated dipolar liquid. The lattice so generated was orientationally disordered and positionally random. The detailed study of orientational relaxation in this random dipolar lattice revealed interesting differences from those of the corresponding dipolar liquid. In particular, we found that the relaxation of the collective orientational correlation functions at the intermediate wave numbers was markedly slower at the long times for the random lattice than that of the liquid. This verified the important role of the translational modes in this regime, as predicted recently by the molecular theories. The single-particle orientational correlation functions of the random lattice also decayed significantly slowly at long times, compared to those of the dipolar liquid.

Analysis on conservation of disulphide bonds and their structural features in homologous protein domain families

Background: Disulphide bridges are well known to play key roles in stability, folding and functions of proteins. Introduction or deletion of disulphides by site-directed mutagenesis have produced varying effects on stability and folding depending upon the protein and location of disulphide in the 3-D structure. Given the lack of complete understanding it is worthwhile to learn from an analysis of extent of conservation of disulphides in homologous proteins. We have also addressed the question of what structural interactions replaces a disulphide in a homologue in another homologue. Results: Using a dataset involving 34,752 pairwise comparisons of homologous protein domains corresponding to 300 protein domain families of known 3-D structures, we provide a comprehensive analysis of extent of conservation of disulphide bridges and their structural features. We report that only 54% of all the disulphide bonds compared between the homologous pairs are conserved, even if, a small fraction of the non-conserved disulphides do include cytoplasmic proteins. Also, only about one fourth of the distinct disulphides are conserved in all the members in protein families. We note that while conservation of disulphide is common in many families, disulphide bond mutations are quite prevalent. Interestingly, we note that there is no clear relationship between sequence identity between two homologous proteins and disulphide bond conservation. Our analysis on structural features at the sites where cysteines forming disulphide in one homologue are replaced by non-Cys residues show that the elimination of a disulphide in a homologue need not always result in stabilizing interactions between equivalent residues. Conclusion: We observe that in the homologous proteins, disulphide bonds are conserved only to a modest extent. Very interestingly, we note that extent of conservation of disulphide in homologous proteins is unrelated to the overall sequence identity between homologues. The non-conserved disulphides are often associated with variable structural features that were recruited to be associated with differentiation or specialisation of protein function.

Synthesis, reactivity and structural characterisation of bicyclic tetraphosphapentazanes

The diphenoxy bicyclic tetraphosphapentazane derivatives (EtN)(5)P-4(OPh)(2) 2 and its monoxide (EtN)(5)P-4(O)(OPh)(2) 3 have been prepared. Both 2 and 3 exist as a mixture of two isomers. One isomer of (EtN)(5)P-4(O)(OPh)(2) 3a has been isolated and its reaction with tetrachloro-1,2-benzoquinone yielded (EtN)(5)P-4(O)(OPh)(2)(O2C6Cl4) 5 in which the junction phosphorus atom becomes five-co-ordinated. Treatment of 2 or 3a with [Mo(CO)(4)(nbd)] (nbd = norbornadiene, bicyclo[2.2.1]hepta-2,5-diene), on the other hand, yielded the chelate complex [Mo(CO)(4){(EtN)(5)P-4(O)(n)(OPh)(2)}] (n = 0 or 1; 6 or 7) in which the peripheral phosphorus atoms are bonded to the metal. The structures of 3a and 5-7 have been confirmed by single-crystal X-ray diffraction studies. The two P3N3 rings in 3a and 5 adopt twist/twist and irregular/twist conformations respectively; the phenoxy substituents occupy the 'pseudo axial' positions. However, an ideal chair conformation is observed for the P3N3 rings in 6 and 7 with the phenoxy substituents taking up the 'pseudo equatorial' positions. The NMR spectroscopic data for the compounds are discussed.

Micropillar compression studies on a bulk metallic glass in different structural states

Uniaxial compression experiments on 0.3, 1 and 3 mu m diameter micropillars of a Zr-based bulk metallic glass in as-cast, shot-peened and structurally relaxed conditions were conducted. Shear band formation and stable propagation is observed to be the plastic deformation mode in all cases, with no detectable difference in yield strength according to either size or condition. The limitations of uniaxial compression tests in assessing the influence of various material conditions on plasticity, when it is inhomogeneous in nature, are illustrated.

Boundary element analysis of reinforced concrete structural elements

A simple and practical technique for the discrete representation of reinforcement in two-dimensional boundary element analysis of reinforced concrete structural elements is presented. The bond developed over the surface of contact between the reinforcing steel and concrete is represented using fictitious one-dimensional spring elements. Potentials of the model developed are demonstrated using a number of numerical examples. The results are seen to be in good agreement with the results obtained using standard finite element software.

Structural characteristics of a giant tropical liana and its mode of canopy spread in an alien environment

To circumvent the practical difficulties in research on tropical rainforest lianas in their natural habitat due to prevailing weather conditions, dense camouflaging vegetation and problems in transporting equipment for experimental investigations, Entada pursaetha DC (syn. Entada scandens Benth., Leguminosae) was grown inside a research campus in a dry subtropical environment. A solitary genet has attained a gigantic size in 17 years, infesting crowns of semi-evergreen trees growing in an area roughly equivalent to 1.6 ha. It has used aerially formed, cable-like stolons for navigating and spreading its canopy across tree gaps. Some of its parts which had remained unseen in its natural habitat due to dense vegetation are described. The attained size of this liana in a climatically different environment raises the question as to why it is restricted to evergreen rainforests. Some research problems for which this liana will be useful are pointed out.

Different phenomenological theories and their abilities to describe the interdiffusion process in a binary system during multiphase growth

There are essentially two different phenomenological models available to describe the interdiffusion process in binary systems in the olid state. The first of these, which is used more frequently, is based on the theory of flux partitioning. The second model, developed much more recently, uses the theory of dissociation and reaction. Although the theory of flux partitioning has been widely used, we found that this theory does not account for the mobility of both species and therefore is not suitable for use in most interdiffusion systems. We have first modified this theory to take into account the mobility of both species and then further extended it to develop relations or the integrated diffusion coefficient and the ratio of diffusivities of the species. The versatility of these two different models is examined in the Co-Si system with respect to different end-member compositions. From our analysis, we found that the applicability of the theory of flux partitioning is rather limited but the theory of dissociation and reaction can be used in any binary system.