Temperature and concentration dependence of adsorption properties of methane in NaY: a molecular dynamics study

Molecular dynamics calculations on methane sorbed in NaY (Si/Al = 3.0) employing realistic methane-methane and methane-zeolite intermolecular potential functions at different temperatures (50, 150, 220, and 300 K) and concentrations (2, 4, 6, and 8 molecules/cage) are reported. The thermodynamic results are in agreement with the available experimental data. Guest-guest and guest-host radial distribution functions (rdfs), energy distribution functions, distribution of cage occupancy, center-of-cage-center-of-mass (coc-com) rdfs, velocity autocorrelation functions for com and angular motion and the Fourier transformed power spectra, and diffusion coefficients are presented as a function of temperature and concentration. At 50 K, methane is localized near the adsorption site. Site-site migration and essentially free rotational motion are observed at 150 K. Molecules preferentially occupy the region near the inner surface of the alpha-cage. The vibrational frequencies for the com of methane shift toward higher values with decreasing temperature and increasing adsorbate concentration. The observed frequencies for com motion are 36, 53, and 85 cm-1 and for rotational motion at 50 K, 95 and 150 cm-1 in agreement with neutron scattering data. The diffusion coefficients show a type I behavior as a function of loading in agreement with NMR measurements. Cage-to-cage diffusion is found to be always mediated by the surface.

Importance of orientational rearrangement during vitrification of hydrocarbons: dependence on molecular shape

On the basis of Monte Carlo calculations of 2,2-dimethylpropane (neopentane), n-pentane, and 2,2-dimethylbutane (neohexane) at several temperatures, thermodynamic properties and radial distribution functions as well as dimerization and bonding energy distribution functions are reported for both liquid and glassy states. Changes in the radial distribution functions on cooling depend on whether the groups are accessible (peripheral) or inaccessible. Peaks in the radial distribution functions corresponding to peripheral groups do not shift to lower distances on cooling and at times display a large increase in the intensity of the first peak. The peaks due to inaccessible groups, on the other hand, shift to lower distances on cooling. The magnitude of the reorientational contribution in determining the resulting structure of the glass is estimated for the different hydrocarbon molecules investigated. The reorientational contribution is highest for neopentane (26%) followed by isopentane (24%), neohexane (22%), and n-pentane (0%). It appears that molecular geometry has an important role in determining the magnitude of the reorientational contribution to the structure of the glass.

Dimethyl sulfoxide induced structural transformations and non-monotonic concentration dependence of conformational fluctuation around active site of lysozyme

Experimental studies have observed significant changes in both structure and function of lysozyme (and other proteins) on addition of a small amount of dimethyl sulfoxide (DMSO) in aqueous solution. Our atomistic molecular dynamic simulations of lysozyme in water-DMSO reveal the following sequence of changes on increasing DMSO concentration. (i) At the initial stage (around 5% DMSO concentration) protein's conformational flexibility gets markedly suppressed. From study of radial distribution functions, we attribute this to the preferential solvation of exposed protein hydrophobic residues by the methyl groups of DMSO. (ii) In the next stage (10-15% DMSO concentration range), lysozome partially unfolds accompanied by an increase both in fluctuation and in exposed protein surface area. (iii) Between 15-20% concentration ranges, both conformational fluctuation and solvent accessible protein surface area suddenly decrease again indicating the formation of an intermediate collapse state. These results are in good agreement with near-UV circular dichroism (CD) and fluorescence studies. We explain this apparently surprising behavior in terms of a structural transformation which involves clustering among the methyl groups of DMSO. (iv) Beyond 20% concentration of DMSO, the protein starts its final sojourn towards the unfolding state with further increase in conformational fluctuation and loss in native contacts. Most importantly, analysis of contact map and fluctuation near the active site reveal that both partial unfolding and conformational fluctuations are centered mostly on the hydrophobic core of active site of lysozyme. Our results could offer a general explanation and universal picture of the anomalous behavior of protein structure-function observed in the presence of cosolvents (DMSO, ethanol, tertiary butyl alcohol, dioxane) at their low concentrations. (C) 2012 American Institute of Physics. [http://dx.doi.org/10.1063/1.3694268]

Molecular Origin of the Intrinsic Bending Force for Helical Morphology Observed in Chiral Amphiphilic Assemblies: Concentration and Size Dependence

The formation of the helical morphology in monolayers and bilayers of chiral amphiphilic assemblies is believed to be driven at least partly by the interactions at the chiral centers of the amphiphiles. However, a detailed microscopic understanding of these interactions and their relation with the helix formation is still not clear. In this article a study of the molecular origin of the chirality-driven helix formation is presented by calculating, for the first time, the effective pair potential between a pair of chiral molecules. This effective potential depends on the relative sizes of the groups attached to the two chiral centers, on the orientation of the amphiphile molecules, and also on the distance between them. We find that for the mirror-image isomers (in the racemic modification) the minimum energy conformation is a nearly parallel alignment of the molecules. On the other hand, the same for a pair of molecules of one kind of enantiomer favors a tilt angle between them, thus leading to the formation of a helical morphology of the aggregate. The tilt angle is determined by the size of the groups attached to the chiral centers of the pair of molecules considered and in many cases predicted it to be close to 45 degrees. The present study, therefore, provides a molecular origin of the intrinsic bending force, suggested by Helfrich (J. Chem. Phys. 1986, 85, 1085-1087), to be responsible for the formation of helical structure. This effective potential may explain many of the existing experimental results, such as the size and the concentration dependence of the formation of helical morphology. It is further found that the elastic forces can significantly modify the pitch predicted by the chiral interactions alone and that the modified real pitch is close to the experimentally observed value. The present study is expected to provide a starting point for future microscopic studies.

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.

Free Energy Gap Dependence of the Electron-Transfer Rate from the Inverted to the Normal Region

numerical study of the free energy gap (FEG) dependence of the electron-transfer rate in polar solvents is presented. This study is based on the generalized multidimensional hybrid model, which not only includes the solvent polarization and the molecular vibration modes, but also the biphasic polar response of the solvent. The free energy gap dependence is found to be sensitive to several factors, including the solvent relaxation rate, the electronic coupling between the surfaces, the frequency of the high-frequency quantum vibrational mode, and the magnitude of the solvent reorganization energy. It is shown that in some cases solvent relaxation can play an important role even in the Marcus normal regime. The minimal hybrid model involves a large number of parameters, giving rise to a diverse non-Marcus FEG behavior which is often determined collectively by these parameters. The model gives the linear free energy gap dependence of the logarithmic rate over a substantial range of FEG, spanning from the normal to the inverted regime. However, even for favorable values of the relevant parameters, a linear free energy gap dependence of the rate could be obtained only over a range of 5000-6000 cm(-1) (compared to the experimentally observed range of 10000 cm(-1) reported by Benniston et al.). The present work suggests several extensions/generalizations of the hybrid model which might be necessary to fully understand the observed free energy gap dependence.

Long range resonance energy transfer from a dye molecule to graphene has (distance)-4 dependence

In our previous report on resonance energy transfer from a dye molecule to graphene [J. Chem. Phys.129, 054703 (2008)], we had derived an expression for the rate of energy transfer from a dye to graphene. An integral in the expression for the rate was evaluated approximately. We found a Yuwaka-type dependence of the rate on the distance. We now present an exact evaluation of the integral involved, leading to very interesting results. For short distances (z < 20 A), the present rate and the previous rate are in good agreement. For larger distances, the rate is found to have a z(-4) dependence on the distance, exactly. Thus we predict that for the case of pyrene on graphene, it is possible to observe fluorescence quenching up to a distance of 300 A. This is in sharp contrast to the traditional fluorescence resonance energy transfer where the quenching is observable only up to 100 A.

How does contrasting dependence of impurity-atom diffusivity on the density of host disordered medium arise?

We report results of molecular dynamics investigations into neutral impurity diffusing within an amorphous solid as a function of the size of the diffusant and density of the host amorphous matrix. We find that self diffusivity exhibits an anomalous maximum as a function of the size of the impurity species. An analysis of properties of the impurity atom with maximum diffusivity shows that it is associated with lower mean square force, reduced backscattering of velocity autocorrelation function, near-exponential decay of the intermediate scattering function (as compared to stretched-exponential decay for other sizes of the impurity species) and lower activation energy. These results demonstrate the existence of size-dependent diffusivity maximum in disordered solids. Further, we show that the diffusivity maximum is observed at lower impurity diameters with increase in density. This is explained in terms of the Levitation parameter and the void structure of the amorphous solid. We demonstrate that these results imply contrasting dependence of self diffusivity (D) on the density of the amorphous matrix, p. D increases with p for small sizes of the impurity but shows an increase followed by a decrease for intermediate sizes of the impurity atom. For large sizes of the impurity atom, D decreases with increase in p. These contrasting dependence arises naturally from the existence of Levitation Effect.

Dependence of metal Auger intensity ratios on the effective charge on metal atoms

Metal Auger line intensity ratios were shown by Rao and others to be directly related to the occupancy of valence states. It is now shown that these intensity ratios are more generally related to the effective charge on the metal atom. The Auger intensity ratios are also directly proportional to valence band intensities of metals. Correlations of the intensity ratios with Auger parameters have also been examined.

Dependence of low-lying energy eigenvalues on the short- and long-range parts of confining potentials

The potential description of a quark-antiquark system seems to work very well in describing a number of hadronic properties. However, the precise form of the potential is unknown. The changes in the low-lying eigenvalues as a result of changes in the long-range part of the potential are investigated in a non-perturbative manner. It is shown by considering a variety of examples that the low-lying eigenvalues are insensitive to the long-range part of the potential.

Buckling in polymer monolayers: Molecular-weight dependence

We present systematic investigations of buckling in Langmuir monolayers of polyvinyl acetate formed at the air-water interface. On compression the polymer monolayers are converted to a continuous membrane with a thickness of ~2–3 nm of well-defined periodicity, lambdab. Above a certain surface concentration the membrane undergoes a morphological transition buckling, leading to the formation of striped patterns. The periodicity seems to depend on molecular weight as per the predictions of the gravity-bending buckling formalism of Milner et al. for fluidlike films on water. However anomalously low values of bending rigidity and Young's modulus are obtained using this formalism. Hence we have considered an alternative model of buckling-based solidlike films on viscoelastic substrates. The values of bending rigidity and Young's modulus obtained by this method, although lower than expected, are closer to the bulk values. Remarkably, no buckling is found to occur above a certain molecular weight. We have tried to explain the observed molecular-weight dependence in terms of the variation in isothermal compressive modulus of the monolayers with surface concentration as well as provided possible explanations for the obtained low values of mechanical properties similar to that observed for ultrathin polymer films.

Temperature dependence of the dielectric constant of diamond

The temperature dependence of the dielectric constant of diamond has been measured over the temperature range 50-2OO"c. The value of E-ldc dT over this range is + 1 x 10-j. Details of the method of measuring the temperature coefficient of dielectric constant are also given. The magnitude and sign of c-ldc, dT for diamond has been theoretically calculated using Maxwell's relationship and Kramers-Heisenberg theory. The agreement between theoretical and experimental values is extremely good.

Unusual stress and grain size dependence for creep in nanocrystalline materials

Analyses of diffusion and dislocation creep in nanocrystals needs to take into account the generally utilized low temperatures, high stresses and very fine grain sizes. In nanocrystals, diffusion creep may be associated with a nonlinear stress dependence and dislocation creep may involve a grain size dependence.

Temperature dependence of 35Cl NQR in 2,6-dichloropyridine

Chlorine NQR in 2,6-dichloropyridine has been investigated in the temperature range 77 K to room temperature and a single resonance line has been observed throughout. Using this data, torsional frequencies of the molecule have been evaluated on the basis of both the Bayer theory and the modified Bayer theory incorporating Tatsuzaki correction.