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.

Structural stability of slender aerospace vehicles: Part I Mathematical modeling

A detailed mechanics based model is developed to analyze the problem of structural instability in slender aerospace vehicles. Coupling among the rigid-body modes, the longitudinal vibrational modes and the transverse vibrational modes due to asymmetric lifting-body cross-section are considered. The model also incorporates the effects of aerodynamic pressure and the propulsive thrust of the vehicle. The model is one-dimensional, and it can be employed to idealized slender vehicles with complex shapes. Condition under which a flexible body with internal stress waves behaves like a perfect rigid body is derived. Two methods are developed for finite element discretization of the system: (1) A time-frequency Fourier spectral finite element method and (2) h-p finite element method. Numerical results using the above methods are presented in Part II of this paper. (C) 2010 Elsevier Ltd. All rights reserved.

Extending lagrange interpolation to develop a 4-node quadrilateral element

Today finite element method is a well established tool in engineering analysis and design. Though there axe many two and three dimensional finite elements available, it is rare that a single element performs satisfactorily in majority of practical problems. The present work deals with the development of 4-node quadrilateral element using extended Lagrange interpolation functions. The classical univariate Lagrange interpolation is well developed for 1-D and is used for obtaining shape functions. We propose a new approach to extend the Lagrange interpolation to several variables. When variables axe more than one the method also gives the set of feasible bubble functions. We use the two to generate shape function for the 4-node arbitrary quadrilateral. It will require the incorporation of the condition of rigid body motion, constant strain and Navier equation by imposing necessary constraints. The procedure obviates the need for isoparametric transformation since interpolation functions are generated for arbitrary quadrilateral shapes. While generating the element stiffness matrix, integration can be carried out to the accuracy desired by dividing the quadrilateral into triangles. To validate the performance of the element which we call EXLQUAD4, we conduct several pathological tests available in the literature. EXLQUAD4 predicts both stresses and displacements accurately at every point in the element in all the constant stress fields. In tests involving higher order stress fields the element is assured to converge in the limit of discretisation. A method thus becomes available to generate shape functions directly for arbitrary quadrilateral. The method is applicable also for hexahedra. The approach should find use for development of finite elements for use with other field equations also.

Effect of coordinated ions on structure and flexibility of parallel G-quadruplexes: A molecular dynamics study

ingle tract guanine residues can associate to form stable parallel quadruplex structures in the presence of certain cations. Nanosecond scale molecular dynamics simulations have been performed on fully solvated fibre model of parallel d(G(7)) quadruplex structures with Na+ or K+ ions coordinated in the cavity formed by the O6 atoms of the guanine bases. The AMBER 4.1 force field and Particle Mesh Ewald technique for electrostatic interactions have been used in all simulations. There quadruplex structures are stable during the simulation, with the middle four base tetrads showing root mean square deviation values between 0.5 to 0.8 Angstrom from the initial structure as well the high resolution crystal structure. Even in the absence of any coordinated ion in the initial structure, the G-quadruplex structure remains intact throughout the simulation. During the 1.1 ns MD simulation, one Nai counter ion from the solvent as well as several water molecules enter the central cavity to occupy the empty coordination sites within the parallel quadruplex and help stabilize the structure. Hydrogen bonding pattern depends on the nature of the coordinated ion, with the G-tetrad undergoing local structural variation to accommodate cations of different sizes. in the absence of any coordinated ion. due to strong mutual repulsion, O6 atoms within G-tetrad are forced farther apart from each other, which leads to a considerably different hydrogen bonding scheme within the G-tetrads and very favourable interaction energy between the guanine bases constituting a G-tetrad. However, a coordinated ion between G-tetrads provides extra stacking energy for the G-tetrads and makes the quadruplex structure more rigid. Na+ ions, within the quadruplex cavity, are more mobile than coordinated K+ ions. A number of hydrogen bonded water molecules are observed within the grooves of all quadruplex structures.

Probing Fluorine Interactions in a Polyhydroxylated Environment: Conservation of a C-F center dot center dot center dot H-C Recognition Motif in Presence of O-H center dot center dot center dot O Hydrogen Bonds

Three conformationally locked fluorinated polycyclitols have been specially crafted on a rigid trans-decalin backbone, employing a surprisingly facile pyridine-poly(hydrogen fluoride)-mediated stereospecific epoxide ring opening as the key reaction. Molecula design of the three fluorinated probes under study focused on providing an efficient platform for (a) evaluating the ability of covalently bonded fluorine, vis-a-vis the isosteric hydroxy group, to act as a H-bond acceptor and (b) examining the possibility for an organic fluorine moiety, placed suitably in a spatially invariant position, to engage an 1,3-diaxial OH functionality in a purported intramolecular O-H center dot center dot center dot F hydrogen bond. The present endeavour reveals that C(sp(3))-F center dot center dot center dot H-C(sp(3)) hydrogen bonds, though weak and lesser investigated, can indeed be observed and supramolecular recognition motifs, involving such interactions, can be conserved even in crystal structures laden with stronger O-H center dot center dot center dot O hydrogen bonds.

On the L2 stability of linear multidimensional time varying systems

This paper analyzes the L2 stability of solutions of systems with time-varying coefficients of the form [A + C(t)]x′ = [B + D(t)]x + u, where A, B, C, D are matrices. Following proof of a lemma, the main result is derived, according to which the system is L2 stable if the eigenvalues of the coefficient matrices are related in a simple way. A corollary of the theorem dealing with small periodic perturbations of constant coefficient systems is then proved. The paper concludes with two illustrative examples, both of which deal with the attitude dynamics of a rigid, axisymmetric, spinning satellite in an eccentric orbit, subject to gravity gradient torques.

Adhesion-Induced Instability in Asperities

Adhesive forces between two approaching asperities will deform the asperities, and under certain conditions this will result in a sudden runaway deformations leading to a jump-to-contact instability. We present finite element-based numerical studies on adhesion-induced deformation and instability in asperities. We consider the adhesive force acting on an asperity, when it is brought near a rigid half-space, due to van der Waals interaction between the asperity and the half-space. The adhesive force is considered to be distributed over the volume of the asperity (body force), thus resulting in more realistic simulations for the length scales considered. Iteration scheme based on a ``residual stress update'' algorithm is used to capture the effect of deformation on the adhesion force, and thereby the equilibrium configuration and the corresponding force. The numerical results are compared with the previous approximate analytical solutions for adhesion force, deformation of the asperity and adhesion-induced mechanical instability (jump-to-contact). It is observed that the instability can occur at separations much higher,and could possibly explain the higher value of instability separation observed in experiments. The stresses in asperities, particularly in case of small ones, are found to be high enough to cause yielding before jump -to-contact. The effect of roughness is considered by modeling a spherical protrusion on the hemispherical asperity.This small-scale roughness at the tip of the asperities is found to control the deformation behavior at small separations, and hence are important in determining the friction and wear due to the jump-to-contact instability.

Helicopter blade flapping with and without small angle assumption in the presence of dynamic stall

The flapping equation for a rotating rigid helicopter blade is typically derived by considering (1)small flap angle, (2) small induced angle of attack and (3) linear aerodynamics. However, the use of nonlinear aerodynamics such as dynamic stall can make the assumptions of small angles suspect as shown in this paper. A general equation describing helicopter blade flap dynamics for large flap angle and large induced inflow angle of attack is derived. A semi-empirical dynamic stall aerodynamics model (ONERA model) is used. Numerical simulations are performed by solving the nonlinear flapping ordinary differential equation for steady state conditions and the validity of the small angle approximations are examined. It is shown that the small flapping assumption, and to a lesser extent, the small induced angle ofattack assumption, can lead to inaccurate predictions of the blade flap response in certain flight conditions for some rotors when nonlinear aerodynamics is considered. (C) 2010 Elsevier Inc. All rights reserved.

Carbon nanotube reinforced supramolecular gels with electrically conducting, viscoelastic and near-infrared sensitive properties

Pristine and long-chain functionalized single-walled carbon nanotubes (SWNTs) were incorporated successfully in supramolecular organogels formed by an all-trans tri(p-phenylenevinylene) bis-aldoxime to give rise to new nanocomposites with interesting mechanical, thermal and electrical properties. Variable-temperature UV-vis and fluorescence spectra reveal both pristine and functionalized SWNTs promote aggregation of the gelator molecules and result in quenching of the UV-vis and fluorescence intensity. Electron microscopy and confocal microscopy show the existence of a densely packed and directionally aligned fibrous network in the resulting nanocomposites. Differential scanning calorimetry (DSC) of the composites shows that incorporation of SWNTs increases the gel formation temperature. The DSC of the xerogels of 1-SWNT composites indicates formation of different thermotropic mesophases which is also evident from polarized optical microscopy. The reinforced aggregation of the gelators on SWNT doping was reflected in the mechanical properties of the composites. Rheology of the composites demonstrates the formation of a rigid and viscoelastic solid-like assembly on SWNT incorporation. The composites from gel-SWNTs were found to be semiconducting in nature and showed enhanced electrical conductivity compared to that of the native organogel. Upon irradiation with a near IR laser at 1064 nm for 5 min it was possible to selectively induce a gel-to-sol phase transition of the nanocomposites, while irradiation for even 30 min of the native organogel under identical conditions did not cause any gel-to-sol conversion.

Analysis of clearance fit pin joints

The analysis of clearance fit joints falls within the realm of mixed boundary problems with moving boundaries. In this paper, this problem is solved by a simple continuum method of analysis applying an inverse technique; the region of contact is specified and the corresponding causative load is evaluated. Illustrations are given for a rigid clearance fit pin in a large elastic plate with smooth zero-shear interface between pin and plate, under biaxial plate stress at infinity and due to load transfer through pin.

Theoretical studies on alpha-helix--DNA interactions

The interaction energies between (Ala)10 and alpha-helix fragment and different nucleotide sequences in right-handed B-form have been optimized using semi-empirical potential energy functions. The energies are calculated for two different orientations of the alpha-helix, viz., when the alpha-helix axis taken in the N----C direction is (i) parallel and (ii) antiparallel to the 5'-3' ascending strand of DNA, proximal to it. When both the DNA molecule as well as the alpha-helix are treated as rigid molecules it is found that a polyalanine alpha-helix has slightly more favourable contacts when it is in the proximity of a four nucleotide sequence of 5'-(N-A-T-N)-3' type, where N is either a purine or a pyrimidine. However, when the two interacting molecules are allowed to undergo local structural variations then the interaction energy appears to be independent of the base sequence confirming the non-specific nature of these interactions.

Edge reinforced circular elastic inclusions in cylindrical shells

The effect of having an edge reinforcement around a circular elastic inclusion in a cylindrical shell is studied. The influence of various parameters of the reinforcement such as area of cross section and moment of inertia on the stress concentrations around the inclusion is investigated. It is found that for certain inclusion parameters it is possible to get an optimum reinforcement, which gives minimum stress concentration around the inclusion. The effect of moment of inertia of the reinforcement of SCF is found to be negligible. The results are plotted in a non-dimensional form and a comparison with flat plate results is made which show the curvature effect. In the limiting case of a rigid reinforcement the results tend to those of a rigid circular inclusion. Results are also presented for different values of μe the ratio of extensional rigidity of shell to that of the inclusion.

X-ray and molecular-dynamics studies on Mycobacterium leprae single-stranded DNA-binding protein and comparison with other eubacterial SSB structures

The crystal structures of two forms of Mycobacterium leprae single-stranded DNA-binding protein (SSB) have been determined at 2.05 and 2.8 A resolution. Comparison of these structures with the structures of other eubacterial SSBs indicates considerable variation in their quaternary association, although the DNA-binding domains in all of them exhibit the same OB-fold. This variation has no linear correlation with sequence variation, but could be related to variation in protein stability. Molecular-dynamics simulations have been carried out on tetrameric molecules derived from the two forms and the prototype Escherichia coli SSB and the individual subunits of both proteins. Together, the X-ray studies and molecular-dynamics simulations yield information on the relatively rigid and flexible regions of the molecule and on the effect of oligomerization on flexibility. The simulations provide insight into the changes in subunit structure on oligomerization. They also provide insight into the stability and time evolution of the hydrogen bonds/water bridges that connect the two pairs of monomers in the tetramer.

Transition metal oxide perovskites by photoelectron and x-ray absorption spectroscopy

X-ray and ultraviolet photoelectron spectroscopy as well as x-ray absorption spectroscopy have been employed to investigate transition metal oxide perovskites of the general formula ABOs (A=La or rare-earth ion, B=trivalent transition metalion). Systematics in the core levels and in the valence bands in the series of LaBOa compounds have been discussed. Lanthanum chemical shifts in the x-ray absorption spectra in this series show interesting trends. Photoelectron spectra of the solid solutions, LaNil_x Coxes, LaNix_x FexO8 and LaFel_x Coxes show that the rigid band model is applicable to these systems. It is shown that x-ray photoelectron spectroscopy can be employed to identify multiple oxidation states of transition metal ions in oxide perovskites.

Studies on simultaneous nhibition of trypsin and chymotrypsin by horsegram Bowman-Birk inhibitor

Bowman-Birk inhibitors (BBI) isolated from plant seeds are small proteins active against trypsin and/or chymotrypsin. These inhibitors have been extensively studied in terms of their structure, interactions, function and evolution. Examination of the known three-dimensional structures of BBIs revealed similarities and subtle differences.The hydrophobic core, deduced from surface accessibility and hydrophobicity plots, corresponding to the two tandem structural domains of the double headed BBI are related by an almost exact two-fold, in contrast to the reactive site loops which depart appreciably from the two-fold symmetry. Also, the orientations of inhibitory loops in soybean and peanut inhibitors were different with respect to the rigid core. Based on the structure of Adzuki bean BBI-trypsin complex, models of trypsin and chymotryspin bound to the monomeric soybean BBI (SBI) were constructed. There were minor short contacts between the two enzymes bound to the inhibitor suggesting near independence of binding. Binding studies revealed that the inhibition of one enzyme in the presence of the other is associated with a minor negative cooperativity. In order to assess the functional significance of the reported oligomeric forms of BBI, binding of proteases to the crystallographic and non-crystallographic dimers as found in the crystal structure of peanut inhibitor were examined. It was found that all the active sites in these oligomers cannot simultaneously participate in inhibition.