The dynamics of solvation of an electron in the image potential state by a layer of polar adsorbates

Recently, ultrafast two-photon photoemission has been used to study electron solvation at a two-dimensional metal/polar adsorbate interfaces [A. Miller , Science 297, 1163 (2002)]. The electron is bound to the surface by the image interaction. Earlier we have suggested a theoretical description of the states of the electron interacting with a two-dimensional layer of the polar adsorbate [K. L. Sebastian , J. Chem. Phys. 119, 10350 (2003)]. In this paper we have analyzed the dynamics of electron solvation, assuming a trial wave function for the electron and the solvent polarization and then using the Dirac-Frenkel variational method to determine it. The electron is initially photoexcited to a delocalized state, which has a finite but large size, and causes the polar molecules to reorient. This reorientation acts back on the electron and causes its wave function to shrink, which will cause further reorientation of the polar molecules, and the process continues until the electron gets self-trapped. For reasonable values for the parameters, we are able to obtain fair agreement with the experimental observations. (c) 2005 American Institute of Physics.

Main chain and segmental dynamics of semi interpenetrating based on polyisoprene and poly(methyl methacrylate)

Main chain and segmental dynamics of polyisoprene (PI) and poly(methyl methacrylate)(PMMA) chains in semi IPNs were systematically studied over a wide range of temperatures (above and below T-g of both polymers) as a function of composition, crosslink density, and molecular weight. The immiscible polymers retained most of its characteristic molecular motion; however, the semi IPN synthesis resulted in dramatic changes in the motional behavior of both polymers due to the molecular level interpenetration between two polymer chains. ESR spin probe method was found to be sensitive to the concentration changes of PMMA in semi IPNs. Low temperature spectra showed the characteristics of rigid limit spectra, and in the range of 293-373 K.complex spectra were obtained with the slow component mostly arisingout of the PMMA rich regions and fast component from the PI phase. We found that the rigid PMMA chains closely interpenetrated into thehighly mobile PI network imparts motional restriction in nearby PI chains, and the highly mobile PI chains induce some degree of flexibility in highly rigid PMMA chains. Molecular level interchain mixing was found to be more efficient at a PMMA concentration of 35 wt.%. Moreover, the strong interphase formed in the above mentionedsemi IPN contributed to the large slow component in the ESR spectra at higher temperature. The shape of the spectra along with the data obtained from the simulations of spectra was correlated to the morphology of the semi IPNs. The correlation time measurement detected the motional region associated with the glass transition of PI and PMMA, and these regions were found to follow the same pattern of shifts in a-relaxation of PI and PMMA observed in DMA analysis. Activation energies associated with the T-g regions were also calculated. T-50G was found to correlate with the T-g of PMMA, and the volume of polymer segments undergoing glass transitional motion was calculated to be 1.7 nm(3).C-13 T-1 rho measurements of PMMA carbons indicate that the molecular level interactions were strong in semi IPN irrespective of the immiscible nature of polymers. The motional characteristics of H atoms attached to carbon atoms in both polymers were analyzed using 2D WISE NMR. Main relaxations of both components shifted inward, and both SEM and TEM analysis showed the development of a nanometer sized morphology in the case of highly crosslinked semi IPN. (C) 2010 Elsevier Ltd. All rights reserved.

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.

Effect of √-irradiation on the structure and dynamics of domains in triglycine selenate

A comparative study of the switching properties of pure and √-irradiated TGSe crystals has been carried out to see the effect of irradiation on the structure and dynamics of domains. The switching behaviour of √-irradiated TGSe has been found to be qualitatively similar to that of unirradiated crystal and this has been interpreted in terms of structural inhibition caused by the formation of radiolysis products as well as the difference between the domain structures of the unirradiated and irradiated samples. Confirmation of this has been obtained by studying the domain patterns using the etch method.

Dynamical cluster approximation: Nonlocal dynamics of correlated electron systems

We recently introduced the dynamical cluster approximation (DCA), a technique that includes short-ranged dynamical correlations in addition to the local dynamics of the dynamical mean-field approximation while preserving causality. The technique is based on an iterative self-consistency scheme on a finite-size periodic cluster. The dynamical mean-field approximation (exact result) is obtained by taking the cluster to a single site (the thermodynamic limit). Here, we provide details of our method, explicitly show that it is causal, systematic, Phi derivable, and that it becomes conserving as the cluster size increases. We demonstrate the DCA by applying it to a quantum Monte Carlo and exact enumeration study of the two-dimensional Falicov-Kimball model. The resulting spectral functions preserve causality, and the spectra and the charge-density-wave transition temperature converge quickly and systematically to the thermodynamic limit as the cluster size increases.

Monte carlo simulation for molecular gas dynamics

The dynamics of low-density flows is governed by the Boltzmann equation of the kinetic theory of gases. This is a nonlinear integro-differential equation and, in general, numerical methods must be used to obtain its solution. The present paper, after a brief review of Direct Simulation Monte Carlo (DSMC) methods due to Bird, and Belotserkovskii and Yanitskii, studies the details of theDSMC method of Deshpande for mono as well as multicomponent gases. The present method is a statistical particle-in-cell method and is based upon the Kac-Prigogine master equation which reduces to the Boltzmann equation under the hypothesis of molecular chaos. The proposed Markoff model simulating the collisions uses a Poisson distribution for the number of collisions allowed in cells into which the physical space is divided. The model is then extended to a binary mixture of gases and it is shown that it is necessary to perform the collisions in a certain sequence to obtain unbiased simulation.

Effect of √-irradiation on the structure and dynamics of domains in triglycine selenate

A comparative study of the switching properties of pure and √-irradiated TGSe crystals has been carried out to see the effect of irradiation on the structure and dynamics of domains. The switching behaviour of √-irradiated TGSe has been found to be qualitatively similar to that of unirradiated crystal and this has been interpreted in terms of structural inhibition caused by the formation of radiolysis products as well as the difference between the domain structures of the unirradiated and irradiated samples. Confirmation of this has been obtained by studying the domain patterns using the etch method.

Entropy and dynamics of water in hydration layers of a bilayer

We compute the entropy and transport properties of water in the hydration layer of dipalmitoylphosphatidylcholine bilayer by using a recently developed theoretical scheme two-phase thermodynamic model, termed as 2PT method; S.-T. Lin et al., J. Chem. Phys. 119, 11792 (2003)] based on the translational and rotational velocity autocorrelation functions and their power spectra. The weights of translational and rotational power spectra shift from higher to lower frequency as one goes from the bilayer interface to the bulk. Water molecules near the bilayer head groups have substantially lower entropy (48.36 J/mol/K) than water molecules in the intermediate region (51.36 J/mol/K), which have again lower entropy than the molecules (60.52 J/mol/K) in bulk. Thus, the entropic contribution to the free energy change (T Delta S) of transferring an interface water molecule to the bulk is 3.65 kJ/mol and of transferring intermediate water to the bulk is 2.75 kJ/mol at 300 K, which is to be compared with 6.03 kJ/mol for melting of ice at 273 K. The translational diffusion of water in the vicinity of the head groups is found to be in a subdiffusive regime and the rotational diffusion constant increases going away from the interface. This behavior is supported by the slower reorientational relaxation of the dipole vector and OH bond vector of interfacial water. The ratio of reorientational relaxation time for Legendre polynomials of order 1 and 2 is approximately 2 for interface, intermediate, and bulk water, indicating the presence of jump dynamics in these water molecules. (C) 2010 American Institute of Physics. doi:10.1063/1.3494115]

Communication pathways from the anti-codon region to the aminoacylation site in Methionyl tRNA Synthetase: Molecular dynamics simulations and the structure network analysis

The enzymes of the family of tRNA synthetases perform their functions with high precision by synchronously recognizing the anticodon region and the aminoacylation region, which are separated by ?70 in space. This precision in function is brought about by establishing good communication paths between the two regions. We have modeled the structure of the complex consisting of Escherichia coli methionyl-tRNA synthetase (MetRS), tRNA, and the activated methionine. Molecular dynamics simulations have been performed on the modeled structure to obtain the equilibrated structure of the complex and the cross-correlations between the residues in MetRS have been evaluated. Furthermore, the network analysis on these simulated structures has been carried out to elucidate the paths of communication between the activation site and the anticodon recognition site. This study has provided the detailed paths of communication, which are consistent with experimental results. Similar studies also have been carried out on the complexes (MetRS + activated methonine) and (MetRS + tRNA) along with ligand-free native enzyme. A comparison of the paths derived from the four simulations clearly has shown that the communication path is strongly correlated and unique to the enzyme complex, which is bound to both the tRNA and the activated methionine. The details of the method of our investigation and the biological implications of the results are presented in this article. The method developed here also could be used to investigate any protein system where the function takes place through long-distance communication.

Analysis of Hydrogen Bonding and Stability of Protein Secondary Structures in Molecular Dynamics Simulation

Molecular Dynamics (MD) simulations provide an atomic level account of the molecular motions and have proven to be immensely useful in the investigation of the dynamical structure of proteins. Once an MD trajectory is obtained, specific interactions at the molecular level can be directly studied by setting up appropriate combinations of distance and angle monitors. However, if a study of the dynamical behavior of secondary structures in proteins becomes important, this approach can become unwieldy. We present herein a method to study the dynamical stability of secondary structures in proteins, based on a relatively simple analysis of backbone hydrogen bonds. The method was developed for studying the thermal unfolding of beta-lactamases, but can be extended to other systems and adapted to study relevant properties.

Matrix method for eigenstructure assignment - The multi-input case with application

The eigenvalue and eigenstructure assignment procedure has found application in a wide variety of control problems. In this paper a method for assigning eigenstructure to a linear time invariant multi-input system is proposed. The algorithm determines a matrix that has eigenvalues and eigenvectors at the desired locations. It is obtained from the knowledge of the open-loop system and the desired eigenstructure. Solution of the matrix equation, involving unknown controller gams, open-loop system matrices, and desired eigenvalues and eigenvectors, results hi the state feedback controller. The proposed algorithm requires the closed-loop eigenvalues to be different from those of the open-loop case. This apparent constraint can easily be overcome by a negligible shift in the values. Application of the procedure is illustrated through the offset control of a satellite supported, from an orbiting platform, by a flexible tether.

Matrix method for eigenvalue assignment: The single input case

The eigenvalue assignment/pole placement procedure has found application in a wide variety of control problems. The associated literature is rather extensive with a number of techniques discussed to that end. In this paper a method for assigning eigenvalues to a Linear Time Invariant (LTI) single input system is proposed. The algorithm determines a matrix, which has eigenvalues at the desired locations. It is obtained from the knowledge of the open-loop system and the desired eigenvalues. Solution of the matrix equation, involving unknown controller gains, open-loop system matrices and desired eigenvalues, results in the state feedback controller. The proposed algorithm requires the closed-loop eigenvalues to be different from those of the open-loop case. This apparent constraint is easily overcome by a negligible shift in the values. Two examples are considered to verify the proposed algorithm. The first one pertains to the in-plane libration of a Tethered Satellite System (TSS) while the second is concerned with control of the short period dynamics of a flexible airplane. Finally, the method is extended to determine the Controllability Grammian, corresponding to the specified closed-loop eigenvalues, without computing the controller gains.

Matrix method for eigenstructure assignment - The multi-input case with application

The eigenvalue and eigenstructure assignment procedure has found application in a wide variety of control problems. In this paper a method for assigning eigenstructure to a Linear time invariant multi-input system is proposed. The algorithm determines a matrix that has eigenvalues and eigenvectors at the desired locations. It is obtained from the knowledge of the open-loop system and the desired eigenstructure. solution of the matrix equation, involving unknown controller gains, open-loop system matrices, and desired eigenvalues and eigenvectors, results in the state feedback controller. The proposed algorithm requires the closed-loop eigenvalues to be different from those of the open-loop case. This apparent constraint can easily be overcome by a negligible shift in the values. Application of the procedure is illustrated through the offset control of a satellite supported, from an orbiting platform, by a flexible tether,

Matrix method for eigenvalue assignment: The single input case

The eigenvalue assignment/pole placement procedure has found application in a wide variety of control problems. The associated literature is rather extensive with a number of techniques discussed to that end. In this paper a method for assigning eigenvalues to a Linear Time Invariant (LTI) single input system is proposed. The algorithm determines a matrix, which has eigenvalues at the desired locations. It is obtained from the knowledge of the open-loop system and the desired eigenvalues. Solution of the matrix equation, involving unknown controller gains, open-loop system matrices and desired eigenvalues, results in the state feedback controller. The proposed algorithm requires the closed-loop eigenvalues to be different from those of the open-loop case. This apparent constraint is easily overcome by a negligible shift in the values. Two examples are considered to verify the proposed algorithm. The first one pertains to the in-plane libration of a Tethered Satellite System (TSS) while the second is concerned with control of the short period dynamics of a flexible airplane. Finally, the method is extended to determine the Controllability Grammian, corresponding to the specified closed-loop eigenvalues, without computing the controller gains.

PVU and wave-particle splitting schemes for Euler equations of gas dynamics

A new way of flux-splitting, termed as the wave-particle splitting is presented for developing upwind methods for solving Euler equations of gas dynamics. Based on this splitting, two new upwind methods termed as Acoustic Flux Vector Splitting (AFVS) and Acoustic Flux Difference Splitting (AFDS) methods are developed. A new Boltzmann scheme, which closely resembles the wave-particle splitting, is developed using the kinetic theory of gases. This method, termed as Peculiar Velocity based Upwind (PVU) method, uses the concept of peculiar velocity for upwinding. A special feature of all these methods that the unidirectional and multidirectional parts of the flux vector are treated separately. Extensive computations done using these schemes demonstrate the soundness of the ideas.