Gas-liquid nucleation at large metastability: unusual features and a new formalism

Nucleation at large metastability is still largely an unsolved problem, even though it is a problem of tremendous current interest, with wide-ranging practical value, from atmospheric research to materials science. It is now well accepted that the classical nucleation theory (CNT) fails to provide a qualitative picture and gives incorrect quantitative values for such quantities as activation-free energy barrier and supersaturation dependence of nucleation rate, especially at large metastability. In this paper, we present an alternative formalism to treat nucleation at large supersaturation by introducing an extended set of order parameters in terms of the kth largest liquid-like clusters, where k = 1 is the largest cluster in the system, k = 2 is the second largest cluster and so on. At low supersaturation, the size of the largest liquid-like cluster acts as a suitable order parameter. At large supersaturation, the free energy barrier for the largest liquid-like cluster disappears. We identify this supersaturation as the one at the onset of kinetic spinodal. The kinetic spinodal is system-size-dependent. Beyond kinetic spinodal many clusters grow simultaneously and competitively and hence the nucleation and growth become collective. In order to describe collective growth, we need to consider the full set of order parameters. We derive an analytic expression for the free energy of formation of the kth largest cluster. The expression predicts that, at large metastability (beyond kinetic spinodal), the barrier of growth for several largest liquid-like clusters disappears, and all these clusters grow simultaneously. The approach to the critical size occurs by barrierless diffusion in the cluster size space. The expression for the rate of barrier crossing predicts weaker supersaturation dependence than what is predicted by CNT at large metastability. Such a crossover behavior has indeed been observed in recent experiments (but eluded an explanation till now). In order to understand the large numerical discrepancy between simulation predictions and experimental results, we carried out a study of the dependence on the range of intermolecular interactions of both the surface tension of an equilibrium planar gas-liquid interface and the free energy barrier of nucleation. Both are found to depend significantly on the range of interaction for the Lennard-Jones potential, both in two and three dimensions. The value of surface tension and also the free energy difference between the gas and the liquid phase increase significantly and converge only when the range of interaction is extended beyond 6-7 molecular diameters. We find, with the full range of interaction potential, that the surface tension shows only a weak dependence on supersaturation, so the reason for the breakdown of CNT (with simulated values of surface tension and free energy gap) cannot be attributed to the supersaturation dependence of surface tension. This remains an unsettled issue at present because of the use of the value of surface tension obtained at coexistence.

Japanese encephalitis virus infection of mouse cell lines: ability to prime mice for generation of virus specific cytotoxic T lymphocytes and differences in CTL recognisable viral determinants

Ten different mouse cell lines were examined for Japanese encephalitis virus (JEV) infection in vitro and then tested for their ability to generate virus specific cytotoxic T lymphocytes (CTL). Among all cell lines examined, Neuro La (a neuroblastoma) was readily infected with JEV as examined by immunofluorescence and viral replication. Among other cells, P388D1, RAW 264.7 (Macrophage origin), Sp2/0 (B-cell Hybridoma), YAC-1 (T-cell lymphoma), and L929 (Fibroblast) were semipermissive to JEV infection. The cytopathic effects caused by progressive JEV infection varied from cell line to cell line. In the case of YAC-1 cells long-term viral antigen expression was observed without significant alterations in cell viability. Intermediate degrees of cytopathicity are seen in RAW 264.7 and L929 cells while infection of PS, Neuro 2a, P388D1 and Sp2/0 caused major viability losses. All infected cell lines were able to prime adult BALB/c (H-2(d)) mice for the generation of secondary JEV specific CTL. In contrast to YAC-1, the permissive neuroblastoma cell line Neuro 2a (H-2K(k)D(d)) was found to be least efficient in its ability to stimulate anti-viral CTL generation. Cold target competition studies demonstrated that both Neuro 2a and YAC-1 (H-2K(k)D(d)) cells expressed similar viral determinants that are recognised by CTL, suggesting that the reason for the lower ability of Neuro 2a to stimulate anti-viral CTL was not due to lack of viral CTL determinants. These findings demonstrate that a variety of mouse cell lines can be infected with Japanese encephalitis virus, and that these infected cells could be utilised to generate virus specific CTL in BALB/c mice.

Solvation dynamics in monohydroxy alcohols: Agreement between theory and different experiments

Recently three different experimental studies on ultrafast solvation dynamics in monohydroxy straight-chain alcohols (C-1-C-4) have been carried out, with an aim to quantify the time constant (and the amplitude) of the ultrafast component. The results reported are, however, rather different from different experiments. In order to understand the reason for these differences, we have carried out a detailed theoretical study to investigate the time dependent progress of solvation of both an ionic and a dipolar solute probe in these alcohols. For methanol, the agreement between the theoretical predictions and the experimental results [Bingemann and Ernsting J. Chem. Phys. 1995, 102, 2691 and Horng et al. J: Phys, Chern, 1995, 99, 17311] is excellent. For ethanol, propanol, and butanol, we find no ultrafast component of the time constant of 70 fs or so. For these three liquids, the theoretical results are in almost complete agreement with the experimental results of Horng et al. For ethanol and propanol, the theoretical prediction for ionic solvation is not significantly different from that of dipolar solvation. Thus, the theory suggests that the experiments of Bingemann and Ernsting and those of Horng et al. studied essentially the polar solvation dynamics. The theoretical studies also suggest that the experimental investigations of Joo et al. which report a much faster and larger ultrafast component in the same series of solvents (J. Chem. Phys. 1996, 104, 6089) might have been more sensitive to the nonpolar part of solvation dynamics than the polar part. In addition, a discussion on the validity of the present theoretical approach is presented. In this theory the ultrafast component arises from almost frictionless inertial motion of the individual solvent molecules in the force field of its neighbors.

Isomerization dynamics in viscous liquids: Microscopic investigation of the coupling and decoupling of the rate to and from solvent viscosity and dependence on the intermolecular potential

A detailed investigation of viscosity dependence of the isomerization rate is carried out for continuous potentials by using a fully microscopic, self-consistent mode-coupling theory calculation of both the friction on the reactant and the viscosity of the medium. In this calculation we avoid approximating the short time response by the Enskog limit, which overestimates the friction at high frequencies. The isomerization rate is obtained by using the Grote-Hynes formula. The viscosity dependence of the rate has been investigated for a large number of thermodynamic state points. Since the activated barrier crossing dynamics probes the high-frequency frictional response of the liquid, the barrier crossing rate is found to be sensitive to the nature of the reactant-solvent interaction potential. When the solute-solvent interaction is modeled by a 6-12 Lennard-Jones potential, we find that over a large variation of viscosity (eta), the rate (k) can indeed be fitted very well to a fractional viscosity dependence: (k similar to eta(-alpha)), with the exponent alpha in the range 1 greater than or equal to alpha >0. The calculated values of the exponent appear to be in very good agreement with many experimental results. In particular, the theory, for the first time, explains the experimentally observed high value of alpha even at the barrier frequency, omega(b). similar or equal to 9 X 10(12) s(-1) for the isomerization reaction of 2-(2'-propenyl)anthracene in liquid eta-alkanes. The present study can also explain the reason for the very low value of vb observed in another study for the isomerization reaction of trans-stilbene in liquid n-alkanes. For omega(b) greater than or equal to 2.0 X 10(13) s(-1), we obtain alpha similar or equal to 0, which implies that the barrier crossing rate becomes identical to the transition-state theory predictions. A careful analysis of isomerization reaction dynamics involving large amplitude motion suggests that the barrier crossing dynamics itself may become irrelevant in highly viscous liquids and the rate might again be coupled directly to the viscosity. This crossover is predicted to be strongly temperature dependent and could be studied by changing the solvent viscosity by the application of pressure. (C) 1999 American Institute of Physics. [S0021-9606(9950514-X].

Structural characterisation of a uracil containing hairpin DNA by NMR and molecular dynamics

Three-dimensional (3D) structure of a hairpin DNA d-CTAGAGGATCCTTTUGGATCCT (22mer; abbreviated as U4-hairpin), which has a uracil nucleotide unit at the fourth position from the 5' end of the tetra-loop has been solved by NMR spectroscopy. The H-1 resonances of this hairpin have been assigned almost completely. NMR restrained molecular dynamics and energy minimisation procedures have been used to describe the 3D structure of the U4 hairpin. This study establishes that the stem of the hairpin adopts a right handed B-DNA conformation while the T-12 and U-15 nucleotide stack upon 3' and 5' ends of the stem, respectively. Further, T-14 stacks upon both T-12 and U-15 while T-13 partially stacks upon T-14. Very weak stacking interaction is observed between T-13 and T-12. All the individual nucleotide bases adopt 'anti' conformation with respect to their sugar moiety. The turning phosphate in the loop is located between T-13 and T-14. The stereochemistry of U-15 mimics the situation where uracil would stack in a B-DNA conformation. This could be the reason as to why the U4-hairpin is found to be the best substrate for its interaction with uracil DNA glycosylase (UDG) compared to the other substrates in which the uracil is at the first, second and third positions of the tetra-loop from its 5' end, as reported previously.

Mutagenicity associated with O6-methylguanine-DNA damage and mechanism of nucleotide flipping by AGT during repair

Methylated guanine damage at O6 position (i.e. O6MG) is dangerous due to its mutagenic and carcinogenic character that often gives rise to G:C-A:T mutation. However, the reason for this mutagenicity is not known precisely and has been a matter of controversy. Further, although it is known that O6-alkylguanine-DNA alkyltransferase (AGT) repairs O6MG paired with cytosine in DNA, the complete mechanism of target recognition and repair is not known completely. All these aspects of DNA damage and repair have been addressed here by employing high level density functional theory in gas phase and aqueous medium. It is found that the actual cause of O6MG mediated mutation may arise due to the fact that DNA polymerases incorporate thymine opposite to O6MG, misreading the resulting O6MG:T complex as an A:T base pair due to their analogous binding energies and structural alignments. It is further revealed that AGT mediated nucleotide flipping occurs in two successive steps. The intercalation of the finger residue Arg 128 into the DNA double helix and its interaction with the O6MG: C base pair followed by rotation of the O6MG nucleotide are found to be crucial for the damage recognition and nucleotide flipping.

Theoretical Studies on the Structure and Bonding of Metallacyclocumulenes, -cyclopentynes, and -cycloallenes

In this paper, we compare the electronic structure of the hafnacycloallene complex Cp(2)HfC(4)Rr'(2)R `' (5Hf), which was previously described by Erker et al., with those of the titanium, zirconium, and hafnium complexes Cp(2)M(eta(4)-RHC(4)HR) (3M; i.e. metallacyclopent-2,3,4-trienes, metallacyclocumulenes) and Cp(2)M(eta(2)-R(2)C(4)R(2)) (4M; i.e. 1-metallacyclopent-3-ynes) using density functional theory (BP86/LANL2DZ) calculations. Moreover, the eta(3)-phenylallenyl zirconocene complex 7Zr, which was synthesized by Wojcicki et al., is included for the comparison. These calculations and extended Huckel calculations show that the bonding in complex 5Hf is remarkably similar to that of complexes 4M and 7Zr. An analysis of the structural parameters and bonding reveals that the unique interaction of the internal carbon atoms along with the terminal carbon atoms with the bent-metallocene moiety is the reason for the unusual stability of these metallacycles. The molecular orbital analysis further suggests that complex 5Hf can react with another metal fragment to give the bimetallic complexes 9 and 10. The electronic structures of complexes 3M, 4M, 5Hf, and 7Zr have been comparatively studied to get a general understanding of the bonding in these metallacycles.

Fluctuation quench and hydrophobic collapse of polymers and biopolymers in water-DMSO binary mixture at low DMSO concentration

Binary mixtures have strong influence on activities of polymers and biopolymers even at low cosolvent concentration. Among the several aqueous binary mixtures studied, water-DMSO especially stands out for its unusual behavior at certain specific concentrations of DMSO. In the present work, we study the effect of water-DMSO binary mixture on polymers and biopolymers by taking a simple linear hydrocarbon chain of intermediate length (n = 30) and the protein lysozyme, respectively. We find that at a mole fraction of 0.05 of DMSO (x(DMSO) = 0.05) in aqueous solution, the hydrocarbon chain adopts the collapsed conformation as the most stable and rigid state. In this case of 0.05 mole fraction of DMSO in bulk, the DMSO concentration in the first hydration layer around the polymer is found to be as large as 17%. Formation of such hydrophobic environment around the polymer is the reason for the collapsed state gaining so much stability. Interestingly, similar quench of conformational fluctuation is also observed for the protein investigated. It is observed that in the case of alkane polymer chains, long wavelength fluctuation gets easily quenched, the polymer being purely hydrophobic. However, in case of the protein, quench of fluctuation is prominent only at the hydrophobic surface, and quench of long wavelength fluctuation becomes insignificant for the full protein. As protein contains both hydrophobic and hydrophilic moieties, the extent of quench of conformational fluctuation with respect to that in pure water is almost half for the biopolymer complex (16.83%) than the same for pure hydrophobic polymer chain (32.43%).