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.

Entanglement and nonclassicality for multimode radiation-field states

Nonclassicality in the sense of quantum optics is a prerequisite for entanglement in multimode radiation states. In this work we bring out the possibilities of passing from the former to the latter, via action of classicality preserving systems like beam splitters, in a transparent manner. For single-mode states, a complete description of nonclassicality is available via the classical theory of moments, as a set of necessary and sufficient conditions on the photon number distribution. We show that when the mode is coupled to an ancilla in any coherent state, and the system is then acted upon by a beam splitter, these conditions turn exactly into signatures of negativity under partial transpose (NPT) entanglement of the output state. Since the classical moment problem does not generalize to two or more modes, we turn in these cases to other familiar sufficient but not necessary conditions for nonclassicality, namely the Mandel parameter criterion and its extensions. We generalize the Mandel matrix from one-mode states to the two-mode situation, leading to a natural classification of states with varying levels of nonclassicality. For two-mode states we present a single test that can, if successful, simultaneously show nonclassicality as well as NPT entanglement. We also develop a test for NPT entanglement after beam-splitter action on a nonclassical state, tracing carefully the way in which it goes beyond the Mandel nonclassicality test. The result of three-mode beam-splitter action after coupling to an ancilla in the ground state is treated in the same spirit. The concept of genuine tripartite entanglement, and scalar measures of nonclassicality at the Mandel level for two-mode systems, are discussed. Numerous examples illustrating all these concepts are presented.

Inertial mass of a vortex in cuprate superconductors

We present here a calculation of the inertial mass of a moving vortex in cuprate superconductors. This is a poorly known basic quantity of obvious interest in vortex dynamics. The motion of a vortex causes a dipolar density distortion and an associated electric field which is screened. The energy cost of the density distortion as well as the related screened electric field contributes to the vortex mass, which is small because of efficient screening. As a preliminary, we present a discussion and calculation of the vortex mass using a microscopically derivable phase-only action functional for the far region which shows that the contribution from the far region is negligible and that most of it arises from the (small) core region of the vortex. A calculation based on a phenomenological Ginzburg-Landau functional is performed in the core region. Unfortunately such a calculation is unreliable; the reasons for it are discussed. A credible calculation of the vortex mass thus requires a fully microscopic non-coarse-grained theory. This is developed, and results are presented for an s-wave BCS-like gap, with parameters appropriate to the cuprates. The mass, about 0.5m(e) per layer, for a magnetic field along the c axis arises from deformation of quasiparticle states bound in the core and screening effects mentioned above. We discuss earlier results, possible extensions to d-wave symmetry, and observability of effects dependent on the inertial mass. [S0163-1829(97)05534-3].

Self-Organization of n-Alkane Chains in Water: Length Dependent Crossover from Helix and Toroid to Molten Globule

We demonstrate a chain length dependent crossover in the structural properties of linear hydrocarbon (n-alkane) chains using detailed atomistic simulations in explicit water. We identify a number of exotic structures of the polymer chain through energy minimization of representative snapshots collected from molecular dynamics trajectory. While the collapsed state is ring-like (circular) for small chains (CnH2n+2; n <= 20) and spherical for very long ones (n = 100), we find the emergence of ordered helical structures at intermediate lengths (n similar to 40). We find different types of disordered helices and toroid-like structures at n = 60. We also report a sharp transition in the stability of the collapsed state as a function of the chain length through relevant free energy calculations. While the collapsed state is only marginally metastable for C20H42, a clear bistable free energy surface emerges only when the chain is about 30 monomers long. For n = 30, the polymer exhibits an intermittent oscillation between the collapsed and the coil structures, characteristic of two stable states separated by a small barrier.

Jump Reorientation of Water Molecules Confined in Narrow Carbon Nanotubes

We used molecular dynamics (MD) simulations to study the reorientational dynamics of water molecules confined inside narrow carbon nanotubes immersed in a bath of water. Our simulations show that the confined water molecules exhibit bistability in their reorientational relaxation, which proceeds by angular jumps between the two stable states. The angular jump of a water molecule in the bulk involves the breaking of a hydrogen bond with one of its neighbors and the formation of a hydrogen bond with a different neighbor. In contrast, the angular jump of a confined water molecule corresponds to an interchange of the two hydrogen atoms that can form a hydrogen bond with the same neighbor. The free energy barrier between these two states is a few k(B)T. The analytic solution of a simplified two-state jump model that qualitatively explains the reorientational behavior observed in simulations is also presented.

The Mycobacterial MsDps2 Protein Is a Nucleoid-Forming DNA Binding Protein Regulated by Sigma Factors sigma(A) and sigma(B)

The Dps (DNA-binding protein from starved cells) proteins from Mycobacterium smegmatis MsDps1 and MsDps2 are both DNA-binding proteins with some differences. While MsDps1 has two oligomeric states, with one of them responsible for DNA binding, MsDps2 has only one DNA-binding oligomeric state. Both the proteins however, show iron-binding activity. The MsDps1 protein has been shown previously to be induced under conditions of starvation and osmotic stress and is regulated by the extra cellular sigma factors sigma(H) and sigma(F). We show here, that the second Dps homologue in M. smegmatis, namely MsDps2, is purified in a DNA-bound form and exhibits nucleoid-like structures under the atomic force microscope. It appears that the N-terminal sequence of Dps2 plays a role in nucleoid formation. MsDps2, unlike MsDps1, does not show elevated expression in nutritionally starved or stationary phase conditions; rather its promoter is recognized by RNA polymerase containing sigma(A) or sigma(B), under in vitro conditions. We propose that due to the nucleoid-condensing ability, the expression of MsDps2 is tightly regulated inside the cells.

Optimum design of cantilever sheet pile walls in sandy soils using inverse reliability approach

In this paper an approach for obtaining depth and section modulus of the cantilever sheet pile wall using inverse reliability method is described. The proposed procedure employs inverse first order reliability method to obtain the design penetration depth and section modulus of the steel sheet pile wall in order that the reliability of the wall against failure modes must meet a desired level of safety. Sensitivity analysis is conducted to assess the effect of uncertainties in design parameters on the reliability of cantilever sheet pile walls. The analysis is performed by treating back fill soil properties, depth of the water table from the top of the sheet pile wall, yield strength of steel and section modulus of steel pile as random variables. Two limit states, viz., rotational and flexural failure of sheet pile wall are considered. The results using this approach are used to develop a set of reliability based design charts for different coefficients of variation of friction angle of the backfill (5%, 10% and 15%). System reliability considerations in terms of series and parallel systems are also studied.

Manifestation of conical intersections in resonance Raman intensities

This paper deals with the manifestations of conical intersections (CIs), unequivocal spectroscopic signatures of which are still elusive, in the resonance Raman intensities. In particular, the results of our calculations on the `two state-two vibrational mode' and the `two state-three vibrational mode' models are presented. The models comprise two excited states of different spatial symmetry, one bright and one dark, which are coupled by a nontotally symmetric mode while the energy gap between them is tuned by one/two totally symmetric modes. Time dependent theory for vibronically coupled states is employed for the calculation and analysis of Raman excitation profiles (REPs). The manifestation of intersections in REPs is studied by extensive modelm calculations and the results of two specific models are presented. Themfeasibility of using REPs to probe the role of CIs in polyatomic systems is ascertained by multimode calculations on two polyatomic systems viz., pyrazine and trans-azobenzene. The study also notes the importance of the pump excitation wavelength dependence in a femtosecond time-resolved experiment probing the intersection-induced nonadiabatic dynamics. Copyright (C) 2009 John Wiley & Sons, Ltd.

Solution structures of conformationally equilibrium forms of holo-acyl carrier protein (PfACP) from Plasmodium falciparum provides insight into the mechanism of activation of ACPs

Acyl Carrier Protein (ACP) from the malaria parasite, Plasmodium falciparum (PfACP) in its holo form is found to exist in two conformational states in solution. Unique 3D solution structures of holo-PfACP have been determined for both equilibrium conformations, using high-resolution NMR methods. Twenty high-resolution solution structures for each of the two forms of holo-PfACP have been determined on the basis of 1226 and 1218 unambiguously assigned NOEs (including NOEs between 4 '-phosphopantetheine prosthetic group (4 '-PP) and protein), 55 backbone dihedral angles and 26 hydrogen bonds. The atomic rmsd values of the determined structures of two equilibrium forms, about the mean coordinates of the backbone and heavy atoms, are 0.48 +/- 0.09 and 0.92 +/- 0.10 and 0.49 +/- 0.08 and 0.97 +/- 0.11 angstrom, respectively. The interaction of 4 '-PP with the polypeptide backbone is reported here for the first time for any of the ACPs. The structures of holo-PfACP consist of three well-defined helices that are tightly packed. The structured regions of the molecule are stabilized by extensive hydrophobic interactions. The difference between the two forms arises from a reorientation of the 4 '-PP group. The enthalpy difference between the two forms, although small, implies that a conformational switch is essential for the activation of holo-ACP. Sequence and structures of holo-PfACP have been compared with those of the ACPs from type I and type II fatty acid biosynthesis pathways (FAS), in particular with the ACP from rat and the butyryl-ACP from E. coli. The PfACP structure, thus determined has several novel features hitherto not seen in other ACPs.

A class of one-dimensional steady-state flows in radiation-gas-dynamics

The study of steady-state flows in radiation-gas-dynamics, when radiation pressure is negligible in comparison with gas pressure, can be reduced to the study of a single first-order ordinary differential equation in particle velocity and radiation pressure. The class of steady flows, determined by the fact that the velocities in two uniform states are real, i.e. the Rankine-Hugoniot points are real, has been discussed in detail in a previous paper by one of us, when the Mach number M of the flow in one of the uniform states (at x=+∞) is greater than one and the flow direction is in the negative direction of the x-axis. In this paper we have discussed the case when M is less than or equal to one and the flow direction is still in the negative direction of the x-axis. We have drawn the various phase planes and the integral curves in each phase plane give various steady flows. We have also discussed the appearance of discontinuities in these flows.

A density matrix renormalization group study of low-lying excitations of polythiophene within a Pariser-Parr-Pople model

Symmetrized density-matrix-renormalization-group calculations have been carried out, within Pariser-Parr-Pople Hamiltonian, to explore the nature of the ground and low-lying excited states of long polythiophene oligomers. We have exploited C-2 symmetry and spin parity of the system to obtain excited states of experimental interest, and studied the lowest dipole allowed excited state and lowest dipole forbidden two photon state, for different oligomer sizes. In the long system limit, the dipole allowed excited state always lies below the lowest dipole forbidden two-photon state which implies, by Kasha rule, that polythiophene fluoresces strongly. The lowest triplet state lies below two-photon state as usual in conjugated polymers. We have doped the system with a hole and an electron and obtained the charge excitation gap and the binding energy of the 1(1)B(u)(-) exciton. We have calculated the charge density of the ground, one-photon and two-photon states for the longer system size of 10 thiophene rings to characterize these states. We have studied bond order in these states to get an idea about the equilibrium excited state geometry of the system. We have also studied the charge density distribution of the singly and doubly doped polarons for longer system size, and observe that polythiophenes do not support bipolarons.

Transport properties of phosphomolybdate and phosphotungstate glasses

The d.c. conductivity of phosphomolybdate and phosphotungstate glasses is discussed. The conductivity of these glasses is due to the hopping of electrons between two valence states (Mo5+ to Mo6+ or W5+ W6+). In some of the glasses, the activation energy itself is found to be a function of temperature. This appears to be due to thermally activated and variable-range hopping mechanisms operating in different temperature regimes. The relation between conductivity and the [M5+]/[Mtotal](M ≡ Mo, W) ratio does not show any systematic variation. This anomaly can be understood using the structural models of these glasses. In contrast, Mott's theory and the Triberis and Friedman model have been used to obtain conductivity parameters such as the percolation distance Rij and 2agrRij (agr is the tunnelling probability). The conductivity parameter 2agrRij is quite useful to resolve the controversy regarding the tunnelling term exp(2agrRij) existing in the literature. For low values of 2agrRij, it is shown that the exp (2agrRij) term is very significant.

X-ray studies on crystalline complexes involving amino acids and peptides XIX. Effects of change in chirality in the complexes of succinic acid with dl- and l-arginine

Crystals of dl-arginine hemisuccinate dihydrate (I)(monoclinic; P21/c; a = 5.292, b = 16.296, c = 15.203 Å; α= 92.89°; Z = 4) and l-arginine hemisuccinate hemisuccinic acid monohydrate (II) (triclinic; P1; a = 5.099; b = 10.222, c = 14.626 Å; α= 77.31, β= 89.46, γ= 78.42°; Z = 2) were grown under identical conditions from aqueous solutions of the components in molar proportions. The structures were solved by direct methods and refined to R = 0.068 for 2585 observed reflections in the case of (I) and R = 0.036 for 2154 observed reflections in the case of (11). Two of the three crystallographically independent arginine molecules in the complexes have conformations different from those observed so far in the crystal structures containing arginine. The succinic acid molecules and the succinate ions in the structures are centrosymmetric and planar. The crystal structure of (II) is highly pseudosymmetric. Arginine-succinate interactions in both the complexes involve specific guanidyl-carboxylate interactions. The basic elements of aggregation in both the structures are ribbons made up of alternating arginine dimers and succinate ions. However, the ribbons pack in different ways in the two structures. (II) presents an interesting case in which two ionisation states of the same molecule coexist in a crystal. The two complexes provide a good example of the effect of change in chirality on stoichiometry, conformation, aggregation, and ionisation state in the solid state.

Negative differential capacitance in n-GaN/p-Si heterojunctions

Negative differential capacitance (NDC) has been observed in n-GaN/p-Si heterojunctions grown by plasma assisted molecular beam epitaxy (PAMBE). The NDC is observed at low frequencies 1 and 10 kilohertz (kHz) and disappeared at a higher testing frequency of 100 kHz. The NDC is also studied with temperature and found that it has disappeared above 323 degrees C. Current-Voltage (I-V) characteristics of n-GaN /p-Si heterojunction were measured at different temperatures and are attributed to the space-charge-limited current (SCLC). A simple model involving two quantum states is proposed to explain the observed NDC behavior. (C) 2010 Elsevier Ltd. All rights reserved.

Photodissociation dynamics of CH2ICl at 222, 236, 266, 280, and-304 nm

Dynamics of I*(P-2(1/2)) formation from CH2ICl dissociation has-been investigated at five different ultraviolet excitation wavelengths, e.g., 222, 236, 266, 280, and similar to304 nm. The quantum yield of I*((2)p(1/2)) production, phi*, has been measured by monitoring nascent I(P-2(3/2)) and I* concentrations using a resonance enhanced multiphoton ionization detection scheme. The measured quantum yield as a function of excitation energy follows the same trend as that of methyl iodide except at 236 run. The photodissociation dynamics of CH2ICl also involves three upper states similar to methyl iodide, and a qualitative correlation diagram has been constructed to account for the observed quantum yield. From the difference in behavior at 236 nm, it appears that the crossing region between the two excited states ((3)Q(0) and (1)Q(1)) is located near the exit valley away from the Franck Condon excitation region. The B- and C-band transitions do not participate in the dynamics, and the perturbation of the methyl iodide states due to Cl-I interaction is relatively weak at the photolysis wavelengths employed in this investigation.