Effect of Sorbate-Zeolite Interaction on Cluster Lifetime and Size of Sorbates: Xe in NaY

Molecular dynamics calculations are reported for Xe in sodium Y zeolite with varying strengths of sorbate-zeolite dispersion interaction. In the absence of any dispersion interaction between the sorbate and the zeolite, the presence of the zeolite has a purely geometrical role. Increase in the strength of the sorbate-zeolite interaction increases the monomer population and decreases the population of dimers and higher sized clusters. The lifetime of the monomers as well as dimers increases with the strength of the dispersion interaction. The observed variations in the lifetime and the population of the different sized clusters is explained in terms of the changes in the potential energy surface caused by the increase in the strength of the dispersion interaction.

Sorbate properties and cage-to-cage diffusion of argon in NaCaA: a molecular dynamics study

Detailed molecular dynamics simulations of argon in zeolite NaCaA are reported. Thermodynamic, structural, and dynamical properties of the sorbate as a function of temperature have been obtained. The properties calculated include various site-site radial distribution functions, different energy distribution functions, selfdiffusion coefficients, the power spectra, and properties relating to cage-to-cage diffusion. The results suggest that sorbate is delocalized above 300 K. Both modes of cage-to-cage diffusion-the surface-mediated and centralized diffusion-are associated with negative barrier heights. Surprisingly, rate of cage-to-cage diffusion is associated with negative and positive activation energies below and above 500 K. The observed differences in the behavior of the rate of cage-to-cage diffusion between Xe-NaY and Ar-NaCaA systems and the nature of the potential energy surface are discussed. Presence of sorbatezeolite interactions results in significant enhancement in the rate of cage-to-cage diffusion and rate of cage visits. It is shown that properties dependent on the long-time behavior such as the diffusion coefficient and the rate of cages visited exhibit the expected Arrhenius dependence on temperature.

Inherent structures of phase-separating binary mixtures: Nucleation, spinodal decomposition, and pattern formation

An energy landscape view of phase separation and nonideality in binary mixtures is developed by exploring their potential energy landscape (PEL) as functions of temperature and composition. We employ molecular dynamics simulations to study a model that promotes structure breaking in the solute-solvent parent binary liquid, at low temperatures. The PEL of the system captures the potential energy distribution of the inherent structures (IS) of the system and is obtained by removing the kinetic energy (including that of intermolecular vibrations). The broader distribution of the inherent structure energy for structure breaking liquid than that of the structure making liquid demonstrates the larger role of entropy in stabilizing the parent liquid of the structure breaking type of binary mixtures. At high temperature, although the parent structure of the structure breaking binary mixture is homogenous, the corresponding inherent structure is found to be always phase separated, with a density pattern that exhibits marked correlation with the energy of its inherent structure. Over a broad range of intermediate inherent structure energy, bicontinuous phase separation prevails with interpenetrating stripes as signatures of spinodal decomposition. At low inherent structure energy, the structure is largely phase separated with one interface where as at high inherent structure energy we find nucleation type growth. Interestingly, at low temperature, the average inherent structure energy (< EIS >) exhibits a drop with temperature which signals the onset of crystallization in one of the phases while the other remains in the liquid state. The nonideal composition dependence of viscosity is anticorrelated with average inherent structure energy.

Reactions of cyanogen radical with alkanes and an explanation for negative temperature dependence of rate constants

Reactions of cyanide radicals with alkanes have been investigated by ab initio methods. It is found that the potential energy surface for reaction of CN with a primary C-H bond in methane has a small positive barrier while reactions of CN with a secondary and a tertiary C-H bond in alkanes are barrierless at the correlated level. A simple explanation for the obtained negative temperature dependence of rate constants for reactions of CN with a secondary and a tertiary C-H bond in alkanes are given in terms of the collision theory of bimolecular reactions. It is shown that for barrierless reactions the negative temperature dependence of the rate constants is attributed to the variation of the pre-exponential factor with temperature.

A profile of adenosine triphosphate

An attempt is made to draw a profile of adenosine triphosphate (ATP) and to project its many actions. The amazing versatility of its participation in a number of synthetic reactions lies in the oligophosphate structure. Many proteins that use ATP have conserved binding 'P-loop' but this gives no clue what makes it so special. The energy transducing reactions leading to synthesis of the terminal phosphodiester had at least three strategies. Of these, direct dehydration and transfer of inorganic phosphate using respiratory energy operate through mechano-coupling in a multisubunit protein. This tripartite, knob-stalk-base structure provides a novel mechanism of rotational catalysis and the tiniest molecular motor, All the reactions occur in concert with no sign of energized chemical intermediate. With the new knowledge on the crystal structure of F-1-ATPase, proton translocation needs a relook. An alternative perspective is emerging on energy being received and stored in polypeptide structure by breaking hydrogen bonds. Membrane serves the purpose of mobilizing the constituent proteins and also as a potential energy carrier of proteins with little loss of energy.

Ultrasonic wave characteristics of a monolayer graphene on silicon substrate

The present work deals with an ultrasonic type of wave propagation characteristics of monolayer graphene on silicon (Si) substrate. An atomistic model of a hybrid lattice involving a hexagonal lattice of graphene and surface atoms of diamond lattice of Si is developed to identify the carbon-silicon bond stiffness. Properties of this hybrid lattice model is then mapped into a nonlocal continuum framework. Equivalent force constant due to Si substrate is obtained by minimizing the total potential energy of the system. For this equilibrium configuration, the nonlocal governing equations are derived to analyze the ultrasonic wave dispersion based on spectral analysis. From the present analysis we show that the silicon substrate affects only the flexural wave mode. The frequency band gap of flexural mode is also significantly affected by this substrate. The results also show that, the silicon substrate adds cushioning effect to the graphene and it makes the graphene more stable. The analysis also show that the frequency bang gap relations of in-plane (longitudinal and lateral) and out-of-plane (flexural) wave modes depends not only on the y-direction wavenumber but also on nonlocal scaling parameter. In the nonlocal analysis, at higher values of the y-directional wavenumber, a decrease in the frequency band gap is observed for all the three fundamental wave modes in the graphene-silicon system. The atoms movement in the graphene due to the wave propagation are also captured for all the tree fundamental wave modes. The results presented in this work are qualitatively different from those obtained based on the local analysis and thus, are important for the development of graphene based nanodevices such as strain sensor, mass and pressure sensors, atomic dust detectors and enhancer of surface image resolution that make use of the ultrasonic wave dispersion properties of graphene. (C) 2011 Elsevier Ltd. All rights reserved.

Monte Carlo and molecular dynamics simulation of argon clusters andn-alkanes in the confined regions of zeolites

Geometry and energy of argon clusters confined in zeolite NaCaA are compared with those of free clusters. Results indicate the possible existence of magic numbers among the confined clusters. Spectra obtained from instantaneous normal mode analysis of free and confined clusters give a larger percentage of imaginary frequencies for the latter indicating that the confined cluster atoms populate the saddle points of the potential energy surface significantly. The variation of the percentage of imaginary frequencies with temperature during melting is akin to the variation of other properties. It is shown that confined clusters might exhibit inverse surface melting, unlike medium-to-large-sized free clusters that exhibit surface melting. Configurational-bias Monte Carte (CBMC) simulations of n-alkanes in zeolites Y and A are reported. CBMC method gives reliable estimates of the properties relating to the conformation of molecules. Changes in the conformational properties of n-butane and other longer n-alkanes such as n-hexane and n-heptane when they are confined in different zeolites are presented. The changes in the conformational properties of n-butane and n-hexane with temperature and concentration is discussed. In general, in zeolite Y as well as A, there is significant enhancement of the gauche population as compared to the pure unconfined fluid.

Temperature dependence of the levitation effect - Implications for separation of multicomponent mixtures

Sufficiently long molecular dynamics simulations have been carried out on spherical monatomic sorbates in NaY zeolite, interacting via simple Lennard-Jones potentials, to investigate the dependence of the levitation effect on the temperature. Simulations carried out in the range 100-300 K suggest that the anomalous peak in the diffusion coefficient (observed when the levitation parameter, gamma, is near unity) decreases in intensity with increase in temperature. The rate of cage-to-cage migrations also exhibits a similar trend. The activation energy obtained from Arrhenius plots is found to exhibit a minimum when the diffusion coefficient is a maximum, corresponding to the gamma approximate to 1 sorbate diameter. In the linear or normal regime, the activation energy increases with increase in sorbate diameter until it shows a sharp decrease in the anomalous regime. Locations and energies of the adsorption sites and their dependence on the sorbate size gives interesting insight into the nature of the underlying potential-energy surface and further explain the observed trend in the activation energy with sorbate size. Cage residence times, tau(c), show little or no change with temperature for the sorbate with diameter corresponding to gamma approximate to 1, whereas there is a significant decrease in tau(c) with increase in temperature for sorbates in the linear regime. The implications of the present study for the separation of mixtures of sorbates are discussed.

Irreversibility-to-reversibility crossover in transient response of an optically trapped particle

We study the transient response of a colloidal bead which is released from different heights and allowed to relax in the potential well of an optical trap. Depending on the initial potential energy, the system's time evolution shows dramatically different behaviors. Starting from the short-time reversible to long-time irreversible transition, a stationary reversible state with zero net dissipation can be achieved as the release point energy is decreased. If the system starts with even lower energy, it progressively extracts useful work from thermal noise and exhibits an anomalous irreversibility. In addition, we have verified the Transient Fluctuation Theorem and the Integrated Transient Fluctuation Theorem even for the non-ergodic descriptions of our system. Copyright (C) EPLA, 2011

Nonlinear dynamics and chaotic motions in feedback-controlled two- and three-degree-of-freedom robots

The dynamics of a feedback-controlled rigid robot is most commonly described by a set of nonlinear ordinary differential equations. In this paper we analyze these equations, representing the feedback-controlled motion of two- and three-degrees-of-freedom rigid robots with revolute (R) and prismatic (P) joints in the absence of compliance, friction, and potential energy, for the possibility of chaotic motions. We first study the unforced or inertial motions of the robots, and show that when the Gaussian or Riemannian curvature of the configuration space of a robot is negative, the robot equations can exhibit chaos. If the curvature is zero or positive, then the robot equations cannot exhibit chaos. We show that among the two-degrees-of-freedom robots, the PP and the PR robot have zero Gaussian curvature while the RP and RR robots have negative Gaussian curvatures. For the three-degrees-of-freedom robots, we analyze the two well-known RRP and RRR configurations of the Stanford arm and the PUMA manipulator respectively, and derive the conditions for negative curvature and possible chaotic motions. The criteria of negative curvature cannot be used for the forced or feedback-controlled motions. For the forced motion, we resort to the well-known numerical techniques and compute chaos maps, Poincare maps, and bifurcation diagrams. Numerical results are presented for the two-degrees-of-freedom RP and RR robots, and we show that these robot equations can exhibit chaos for low controller gains and for large underestimated models. From the bifurcation diagrams, the route to chaos appears to be through period doubling.

[n]peristylanes and [n]oxa[n]peristylanes (n = 3-6): A theoretical study

Theoretical studies at the HF and Becke3LYP levels using 6-31G* basis sets were carried out on a series of [n]peristylanes and [n]oxa[n]peristylanes (n = 3-6) to understand their structure and energetics. The structures of the [3]- and [4]peristylanes (1, 2) and their era-derivatives (5, 6) were calculated to have the anticipated high symmetry, C-nv. In contrast, a C-s structure (9) at HF/6-31G* and another (25) at the Becke3LYP/6-31G* level were calculated for the [5]oxa[5]peristylane. The energy difference between them is extremely small even though there are major differences in the structures indicating every soft potential energy surface: On the other hand, the potential energy surface of [6]oxa[6]peristylane is not as soft. Similar structures were also calculated for the top rings. Calculations on the seco-compounds 11-14 and 15-19 (Table 4) indicate that there is no unusual strain involved in the formation of 27 from 19. The Li+ interaction energies of the [n]oxa[n]peristylanes are 61.7 (n = 3), 72.8 (n = 4), 84.2 (n = 5) and 91.7 (n = 6) kcal mol(-1) at the Becke3LYP/6-3IG* level. Dramatic differences between the C-C bond lengths obtained from the solid state X-ray diffraction studies and those from the calculations for the [n]oxa[n]peristylanes were also observed.

Investigation of short-time isomerization dynamics in p-nitroazobenzene from resonance Raman intensity analysis

Resonance Raman (RR) spectra are presented for p-nitroazobenzene dissolved in chloroform using 18 excitation Wavelengths, covering the region of (1)(n --> pi*) electronic transition. Raman intensities are observed for various totally symmetric fundamentals, namely, C-C, C-N, N=N, and N-O stretching vibrations, indicating that upon photoexcitation the excited-state evolution occurs along all of these vibrational coordinates. For a few fundamentals, interestingly, in p-nitroazobenzene, it is observed that the RR intensities decrease near the maxima of the resonant electronic (1)(n --> pi*) transition. This is attributed to the interference from preresonant scattering due to the strongly allowed (1)(pi --> pi*) electronic transition. The electronic absorption spectrum and the absolute Raman cross section for the nine Franck-Condon active fundamentals of p-nitroazobenzene have been successfully modeled using Heller's time-dependent formalism for Raman scattering. This employs harmonic description of the lowest energy (1)(n --> pi*) potential energy surface. The short-time isomerization dynamics is then examined from a priori knowledge of the ground-state normal mode descriptions of p-nitroazobenzene to convert the wave packet motion in dimensionless normal coordinates to internal coordinates. It is observed that within 20 fs after photoexcitation in p-nitroazobenzene, the N=N and C-N stretching vibrations undergo significant changes and the unsubstituted phenyl ring and the nitro stretching vibrations are also distorted considerably.

Through ring umbrella inversion of a molecular rattle

We report theoretical investigations on some [Ring]Li--(+) compounds, which can exhibit a through ring umbrella like inversion. Our studies predict cyclononatetraenyllithium to be molecular rattle, in which such inversions can occur. The potential energy for the motion is a double well, with an activation barrier of 11.50 kcal/mol. We find that the lithium should go through the ring easily by an excitation to nu = 17 vibrational level. (C) 2002 Elsevier Science B.V. All rights reserved.

Towards the concept of hydrodynamic cavitation control

A careful study of the existing literature available in the field of cavitation reveals the potential of ultrasonics as a tool for controlling and, if possible, eliminating certain types of hydrodynamic cavitation through the manipulation of nuclei size present in a flow. A glass venturi is taken to be an ideal device to study the cavitation phenomenon at its throat and its potential control. A piezoelectric transducer, driven at the crystal resonant frequency, is used to generate an acoustic pressure field and is termed an �ultrasonic nuclei manipulator (UNM)�. Electrolysis bubbles serve as artificial nuclei to produce travelling bubble cavitation at the venturi throat in the absence of a UNM but this cavitation is completely eliminated when a UNM is operative. This is made possible because the nuclei, which pass through the acoustic field first, cavitate, collapse violently and perhaps fragment and go into dissolution before reaching the venturi throat. Thus, the potential nuclei for travelling bubble cavitation at the venturi throat seem to be systematically destroyed through acoustic cavitation near the UNM. From the solution to the bubble dynamics equation, it has been shown that the potential energy of a bubble at its maximum radius due to an acoustic field is negligible compared to that for the hydrodynamic field. Hence, even though the control of hydrodynamic macro cavitation achieved in this way is at the expense of acoustic micro cavitation, it can still be considered to be a significant gain. These are some of the first results in this direction.

Identity, energetics, dynamics and environment of interfacial water molecules in a micellar solution

The structure and energetics of interfacial water molecules in the aqueous micelle of cesium perfluorooctanoate have been investigated, using large-scale atomistic molecular dynamics simulations, with the primary objective of classifying them. The simulations show that the water molecules at the interface fall into two broad classes: bound and free, present in a ratio of 9:1. The bound water molecules can be further categorized on the basis of the number of hydrogen bonds (one or two) that they form with the surfactant headgroups. The hydrogen bonds of the doubly hydrogen-bonded species are found to be, on the average, slightly weaker than those in the singly bonded species. The environment around interfacial water molecules is more ordered than that in the bulk. The surface water molecules have substantially lower potential energy, because of interaction with the micelle. In particular, both forms of bound water have energies that are lower by �2.5-4.0 kcal/ mol. Entropy is found to play an important role in determining the relative concentration of the species.