Two-Phase Thermodynamic Model for Efficient and Accurate Absolute Entropy of Water from Molecular Dynamics Simulations

Presented here is the two-phase thermodynamic (2PT) model for the calculation of energy and entropy of molecular fluids from the trajectory of molecular dynamics (MD) simulations. In this method, the density of state (DoS) functions (including the normal modes of translation, rotation, and intramolecular vibration motions) are determined from the Fourier transform of the corresponding velocity autocorrelation functions. A fluidicity parameter (f), extracted from the thermodynamic state of the system derived from the same MD, is used to partition the translation and rotation modes into a diffusive, gas-like component (with 3Nf degrees of freedom) and a nondiffusive, solid-like component. The thermodynamic properties, including the absolute value of entropy, are then obtained by applying quantum statistics to the solid component and applying hard sphere/rigid rotor thermodynamics to the gas component. The 2PT method produces exact thermodynamic properties of the system in two limiting states: the nondiffusive solid state (where the fluidicity is zero) and the ideal gas state (where the fluidicity becomes unity). We examine the 2PT entropy for various water models (F3C, SPC, SPC/E, TIP3P, and TIP4P-Ew) at ambient conditions and find good agreement with literature results obtained based on other simulation techniques. We also validate the entropy of water in the liquid and vapor phases along the vapor-liquid equilibrium curve from the triple point to the critical point. We show that this method produces converged liquid phase entropy in tens of picoseconds, making it an efficient means for extracting thermodynamic properties from MD simulations.

Rotation sensitivity of Lau fringes: an analysis based on coherence theory

An analysis based on coherence theory is presented, which explains the experimentally observed rotation sensitivity of the contrast of Lau fringes obtained under spatially incoherent illumination.

The intramolecular external heavy atom effect in bromo-, benzo-, and naphthonorbornenes

A theoretical analysis of the external heavy atom effect of a halogen atom on the radiative rate constant of phosphorescence is examined as a function of position of a bromine atom or atoms relative to a naphthalene or a benzene chromophore for a series of mono- and dibromo-, naphtho-, and benzonorbornenes. The theoretical results are then compared to experimentaldata and lead to the conclusion that the enhancement of the phosphorescence process takes place through the second-ordermixing of the triplet states of the chromophore with the singlet charge transfer states arising primarily from an electron transferfrom the orbitals of the heavy atom perturber to the unfilled x* orbitals of the chromophore.

Role of spatial coherence on the rotation sensitivity of Lau fringes: an experimental study

A simple technique involving the use of a rotating and a stationary diffuser has been developed to vary the spatial coherence of light from a He-Ne laser. Using this technique an experimental investigation of the dependence of rotation sensitivity of Lau fringes on the spatial coherence of the illuminating wavefield has been carried out. It is observed that (i) the rotation sensitivity of Lau fringes varies in a well-defined manner as a function of the spatial coherence of the light used; (ii) the extremely good rotation sensitivity of Lau fringes can be used to great advantage (compared to the conventional double slit method) in the measurement of the spatial coherence of a wavefield; (iii) Lau fringes are formed at various levels of spatial coherence and as such it appears that the Lau effect need not be associated with an incoherent optical field

An accurate numerical integration scheme for finite rotations using rotation vector parametrization

A numerical integration procedure for rotational motion using a rotation vector parametrization is explored from an engineering perspective by using rudimentary vector analysis. The incremental rotation vector, angular velocity and acceleration correspond to different tangent spaces of the rotation manifold at different times and have a non-vectorial character. We rewrite the equation of motion in terms of vectors lying in the same tangent space, facilitating vector space operations consistent with the underlying geometric structure. While any integration algorithm (that works within a vector space setting) may be used, we presently employ a family of explicit Runge-Kutta algorithms to solve this equation. While this work is primarily motivated out of a need for highly accurate numerical solutions of dissipative rotational systems of engineering interest, we also compare the numerical performance of the present scheme with some of the invariant preserving schemes, namely ALGO-C1, STW, LIEMIDEA] and SUBCYC-M. Numerical results show better local accuracy via the present approach vis-a-vis the preserving algorithms. It is also noted that the preserving algorithms do not simultaneously preserve all constants of motion. We incorporate adaptive time-stepping within the present scheme and this in turn enables still higher accuracy and a `near preservation' of constants of motion over significantly longer intervals. (C) 2010 The Franklin Institute. Published by Elsevier Ltd. All rights reserved.

An infrared spectroscopic study of C7 intramolecular hydrogen bonds in peptides

The ir-spectra in the N-H stretching region of Piv-Pro-NHMe and Boc-Pro-NHMe have been studied in carbon tetrachloride and chloroform solutions over a wide range of concentrations. Based on the concentration dependence of the N-H stretching bands, it has been shown that the characteristic N-H stretching band due to the C7 intramolecular hydrogen bond is around 3335 cm-'. Intermolecular hydrogen bonding also occurs to a small extent in these peptides, giving rise to a slight concentration dependence of the N-H stretching bands. The band around 3335 cm-* need not necessarily be due to C7 hydrogen bonds alone as proposed by Tsuboi et al. or to intermolecular hydrogen bonding alone as proposed by Maxfield et al.; this conclusion is supported by studies on Boc-Leu-NHMe, which undergoes only intermolecular hydrogen bonding We have shown that 2-Aib-Aib-OMe and Z-Aib- Ala-OMe form C7 intramolecular hydrogen bonds in addition to C5 intramolecular hydrogen bonds. The present studies also show that all the peptides studied exist in more than one conformation in solution.

On the relation between rotation increments in different tangent spaces

In computational mechanics, finite rotations are often represented by rotation vectors. Rotation vector increments corresponding to different tangent: spaces are generally related by a linear operator, known as the tangential transformation T. In this note, we derive the higher order terms that are usually left out in linear relation. The exact nonlinear relation is also presented. Errors via the linearized T are numerically estimated. While the concept of T arises out of the nonlinear characteristics of the rotation manifold, it has been derived via tensor analysis in the context of computational mechanics (Cardona and Geradin, 1988). We investigate the operator T from a Lie group perspective, which provides a better insight and a 1-1 correspondence between approaches based on tensor analysis and the standard matrix Lie group theory. (C) 2010 Elsevier Ltd. All rights reserved.

Stereochemical studies on cyclic peptides. Part XI. Conformation of cyclic pentapeptides having intramolecular 3 leads to 1 hydrogen bonds.

Conformational analysis of cyclic pentapeptides having two intra-ring 3 leads to 1 hydrogen bonds has been carried out. It is found that the structure can easily be formed with trans planar peptide units without causing significant angular strain at the alpha-carbon atoms. Four different types of conformations designated Types I--IV are possible for the backbone structure. Details of these four types of conformations and also the accommodating possibility of these types for allglycyl and all-alanyl residues are presented. Three of the four types have relatively low energies for glycyl residues whereas the other one has a slightly higher energy. When alanyl residues are introduced at the five alpha-carbon atoms, the types that are energetically favourable depend upon the sequence of isomers. Energy calculations have also been carried out for the combinations of glycyl, L- and D-alanyl residues. The theoretical results are compared with available experimental observations both from solution and solid state studies.

Time-Dependent Dynamics of a Planar Galaxy Model

Time-dependent models of collisionless stellar systems with harmonic potentials allowing for an essentially exact analytic description have recently been described. These include oscillating spheres and spheroids. This paper extends the analysis to time-dependent elliptic discs. Although restricted to two space dimensions, the systems are richer in that their parameters form a 10-dimensional phase space (in contrast to six for the earlier models). Apart from total energy and angular momentum, two additional conserved quantities emerge naturally. These can be chosen as the areas of extremal sections of the ellipsoidal region of phase space occupied by the system (their product gives the conserved volume). The present paper describes the construction of these models. An application to a tidal encounter is given which allows one to go beyond the impulse approximation and demonstrates the effects of rotation of the perturbed system on energy and angular-momentum transfer. The angular-momentum transfer is shown to scale inversely as the cube of the encounter velocity for an initial configuration of the perturbed galaxy with zero quadrupole moment.