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.

Magnetotransport of Dirac fermions on the surface of a topological insulator

We study the properties of Dirac fermions on the surface of a topological insulator in the presence of crossed electric and magnetic fields. We provide an exact solution to this problem and demonstrate that, in contrast to their counterparts in graphene, these Dirac fermions allow relative tuning of the orbital and Zeeman effects of an applied magnetic field by a crossed electric field along the surface. We also elaborate and extend our earlier results on normal-metal-magnetic film-normal metal (NMN) and normal-metal-barrier-magnetic film (NBM) junctions of topological insulators [S. Mondal, D. Sen, K. Sengupta, and R. Shankar, Phys. Rev. Lett. 104, 046403 (2010)]. For NMN junctions, we show that for Dirac fermions with Fermi velocity vF, the transport can be controlled using the exchange field J of a ferromagnetic film over a region of width d. The conductance of such a junction changes from oscillatory to a monotonically decreasing function of d beyond a critical J which leads to the possible realization of magnetic switches using these junctions. For NBM junctions with a potential barrier of width d and potential V-0, we find that beyond a critical J, the criteria of conductance maxima changes from chi=eV(0)d/h upsilon(F)=n pi to chi=(n+1/2)pi for integer n. Finally, we compute the subgap tunneling conductance of a normal-metal-magnetic film-superconductor junctions on the surface of a topological insulator and show that the position of the peaks of the zero-bias tunneling conductance can be tuned using the magnetization of the ferromagnetic film. We point out that these phenomena have no analogs in either conventional two-dimensional materials or Dirac electrons in graphene and suggest experiments to test our theory.

Study of magnetic field effect on the specific heat and entropy in DyBa2Cu3O7-x

The change in the specific heat by the application of magnetic field up to 161 for high temperature superconductor system for DyBa2Cu3O7-x by Revaz et al. [23] is examined through the phenomenological Ginzburg-Landau(G-L) theory of anisotropic Type-II superconductors. The observed specific heat anomaly near T-c with magnetic field is explained qualitatively through the expression <Delta C > = (B-a/T-c) t/(1 - t)(alpha Theta(gamma)lambda(2)(m)(0)), which is the anisotropic formulation of the G-L theory in the London limit developed by Kogan and coworkers; relating to the change in specific heat Delta C for the variation of applied magnetic field for different orientations with c-axis. The analysis of this equation explains satisfactorily the specific heat anomaly near T-c and determines the anisotropic ratio gamma as 5.608, which is close to the experimental value 5.3 +/- 0.5given in the paper of Revaz et al. for this system. (C) 2010 Elsevier B.V. All rights reserved.

Topological analysis off charge density distribution in concomitant polymorphs of 3-acetylcoumarin, a case off packing polymorphism

Detailed investigation of the charge density distribution in concomitant polymorphs of 3-acetylcoumarin in terms of experimental and theoretical densities shows significant differences in the intermolecular features when analyzed based on the topological properties via the quantum theory of atoms in molecules. The two forms, triclinic and monoclinic (Form A and Form B), pack in the crystal lattice via weak C-H---O and C-H---pi interactions. Form A results in a head-to-head molecular stack, while Form B generates a head-to-tail stack. Form A crystallizes in PI (Z' = 2) and Form B crystallizes in P2(1)/n (Z = 1). The electron density maps of the polymorphs demonstrate the differences in the nature of the charge density distribution in general. The charges derived from experimental and theoretical analysis show significant differences with respect to the polymorphic forms. The molecular dipole moments differ significantly for the two forms. The lattice energies evaluated at the HF and DFT (B3LYP) methods with 6-31G** basis set for the two forms clearly suggest that Form A is the thermodynamically stable form as compared to Form B. Mapping of electrostatic potential over the molecular surface shows dominant variations in the electronegative region, which bring out the differences between the two forms.

Topological defects in crystals: A density-wave theory

A new approach for describing dislocations and other topological defects in crystals, based on the density wave theory of Ramakrishnan and Yussouff is presented. Quantitative calculations are discussed in brief for the order parameter profiles, the atomic configuration and the free energy of a screw dislocation with Burgers vector b = (a/2, a/2,a/2 ) in a bcc solid. Our results for the free energy of the dislocation in a crystal of sizeR, when expressed as (λb 2/4π) ln (αR/|b|) whereλ is the shear elastic constant, yield, for example, the valueα ⋍ 1·85 for sodium at its freezing temperature (371°K). The density distribution in the presence of the dislocation shows that the dislocation core has a columnar character. To our knowledge, this study represents the first calculation of dislocation structure, including the core, within the framework of an order parameter theory incorporating thermal effects.

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]

Path integral quantisation of topological field theories

Conditions for quantum topological invariance of classically topological field theories in the path integral formulation are discussed. Both the three-dimensional Chern-Simons system and a Witten-type topological field theory are shown to satisfy these conditions.

Euclideanization, topological theories, higher dimensions and all that

It is shown that the euclideanized Yukawa theory, with the Dirac fermion belonging to an irreducible representation of the Lorentz group, is not bounded from below. A one parameter family of supersymmetric actions is presented which continuously interpolates between the N = 2 SSYM and the N = 2 supersymmetric topological theory. In order to obtain a theory which is bounded from below and satisfies Osterwalder-Schrader positivity, the Dirac fermion should belong to a reducible representation of the Lorentz group and the scalar fields have to be reinterpreted as the extra components of a higher dimensional vector field.

Maximum-entropy calculation of the electron density at 4 A resolution of Pf1 filamentous bacteriophage

A 4 A electron-density map of Pf1 filamentous bacterial virus has been calculated from x-ray fiber diffraction data by using the maximum-entropy method. This method produces a map that is free of features due to noise in the data and enables incomplete isomorphous-derivative phase information to be supplemented by information about the nature of the solution. The map shows gently curved (banana-shaped) rods of density about 70 A long, oriented roughly parallel to the virion axis but slewing by about 1/6th turn while running from a radius of 28 A to one of 13 A. Within these rods, there is a helical periodicity with a pitch of 5 to 6 A. We interpret these rods to be the helical subunits of the virion. The position of strongly diffracted intensity on the x-ray fiber pattern shows that the basic helix of the virion is right handed and that neighboring nearly parallel protein helices cross one another in an unusual negative sense.

Evidence of gate-tunable topological excitations in two-dimensional electron systems

We report experimental observations of a new mechanism of charge transport in two-dimensional electron systems (2DESs) in the presence of strong Coulomb interaction and disorder. We show that at low enough temperature the conductivity tends to zero at a nonzero carrier density, which represents the point of essential singularity in a Berezinskii-Kosterlitz-Thouless-like transition. Our experiments with many 2DESs in GaAs/AlGaAs heterostructures suggest that the charge transport at low carrier densities is due to the melting of an underlying ordered ground state through proliferation of topological defects. Independent measurement of low-frequency conductivity noise supports this scenario.

Quenching across quantum critical points: Role of topological patterns

We introduce a one-dimensional version of the Kitaev model consisting of spins on a two-legged ladder and characterized by Z(2) invariants on the plaquettes of the ladder. We map the model to a fermionic system and identify the topological sectors associated with different Z2 patterns in terms of fermion occupation numbers. Within these different sectors, we investigate the effect of a linear quench across a quantum critical point. We study the dominant behavior of the system by employing a Landau-Zener-type analysis of the effective Hamiltonian in the low-energy subspace for which the effective quenching can sometimes be non-linear. We show that the quenching leads to a residual energy which scales as a power of the quenching rate, and that the power depends on the topological sectors and their symmetry properties in a non-trivial way. This behavior is consistent with the general theory of quantum quenching, but with the correlation length exponent nu being different in different sectors. Copyright (C) EPLA, 2010

Interdiffusion study on Co(W) solid solution and topological close-packed mu phase in Co-W system

An interdiffusion study is conducted on the Co-W system by a diffusion couple technique. The interdiffusion coefficient of the Co(W) solid solution and the Co7W6 mu phase is determined. The activation energy is found to increase with the W content of the Co(W) solid solution. (C) 2010 Elsevier Ltd. All rights reserved.

Fluorescence and topological gap of conjugated phenylene polymers

We propose that strong fluorescence in conjugated polymers requires a dipole-allowed state to be the lowest singlet. Hückel theory for para-conjugated phenyl rings yields an extended, topologically one-dimensional ?-system with increased alternation, states localized on each ring, and charge-transfer excitations between them. Exact Pariser�Parr�Pople results and molecular spectra for oligomers support a topological contribution and a lowest dipole-allowed singlet in phenylene polymers.

Molecular interpretation of the linear relationship between the entropy and the enthalpy of activation of charge transfer reactions in polar liquids

The well-known linear relationship (T?S# =??H# +?, where 1 >? > 0,? > 0) between the entropy (?S#) and the enthalpy (?H#) of activation for reactions in polar liquids is investigated by using a molecular theory. An explicit derivation of this linear relation from first principles is presented for an outersphere charge transfer reaction. The derivation offers microscopic interpretation for the quantities? and?. It has also been possible to make connection with and justify the arguments of Bell put forward many years ago.

Studies in crystal engineering: Effect of fluorine substitution in crystal packing and topological photodimerization of styryl coumarins in the solid state

Styryl coumarins generally yield centrosymmetric (alpha-mode, anti-HT) photodimers when subjected to irradiation in the solid state, However, the substitution of fluorine dramatically alters the packing mode and steers the molecules 4-(4-fluorostyryl)coumarin 1 and 4-(2-fluorostyryl)coumarin 2 to form a stereospecific photodimer, beta-mode, syn-HH across the styrenic double bond (yield 78-85%). The stereochemistry of the photodimer 2a has been established by X-ray crystallography. There is no evidence for the presence of C-H ... F interactions. The true nature of the weak atom-atom interactions called into play when fluorine is substituted is not clear, It is observed that the fluoro substituted compounds have greater crystal density than the corresponding unsubstituted ones.