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.

Temperature- and pressure-dependent study of Cl-35 NQR frequency and spin lattice relaxation time in 2,3-dichloroanisole

The temperature and pressure dependence of Cl-35 NQR frequency and spin lattice relaxation time (T-1) were investigated in 2,3-dichloroanisole. Two NQR signals were observed throughout the temperature and pressure range studied. T-1 were measured in the temperature range from 77 to 300 K and from atmospheric pressure to 5 kbar. Relaxation was found to be due to the torsional motion of the molecule and also reorientation f motion of the CH3 group. T-1 versus temperature data were analyzed on the basis of Woessner and Gutowsky model, and the activation energy for the reorientation of the CH3 group was estimated. The temperature dependence of the average torsional lifetimes of the molecules and the transition probabilities were also obtained. NQR frequency shows a nonlinear behavior with pressure, indicating both dynamic and static effects of pressure. The pressure coefficients were observed to be positive for both the lines. A thermodynamic analysis of the data was carried out to determine the constant volume temperature coefficients of the NQR frequency. The variation of spin lattice time with pressure was very small, showing that the relaxation is mainly due to the torsional motions of the molecules. Copyright (C) 2010 John Wiley & Sons, Ltd.

Critical micelle concentration from a lattice gas model

The temperature dependence of the critical micelle concentration (CMC) and a closed-loop coexistence curve are obtained, via Monte Carlo simulations, in the water surfactant limit of a two-dimensional version of a statistical mechanical model for micro-emulsions, The CMC and the coexistence curve reproduce various experimental trends as functions of the couplings. In the oil-surfactant limit, there is a conventional coexistence cure with an upper consolute point that allows for a region of three-phase coexistence between oil-rich, water-rich and microemulsion phases.

A kinetic Ising model study of dynamical correlations in confined fluids: Emergence of both fast and slow time scales

Experiments and computer simulation studies have revealed existence of rich dynamics in the orientational relaxation of molecules in confined systems such as water in reverse micelles, cyclodextrin cavities, and nanotubes. Here we introduce a novel finite length one dimensional Ising model to investigate the propagation and the annihilation of dynamical correlations in finite systems and to understand the intriguing shortening of the orientational relaxation time that has been reported for small sized reverse micelles. In our finite sized model, the two spins at the two end cells are oriented in the opposite directions to mimic the effects of surface that in real system fixes water orientation in the opposite directions. This produces opposite polarizations to propagate inside from the surface and to produce bulklike condition at the center. This model can be solved analytically for short chains. For long chains, we solve the model numerically with Glauber spin flip dynamics (and also with Metropolis single-spin flip Monte Carlo algorithm). We show that model nicely reproduces many of the features observed in experiments. Due to the destructive interference among correlations that propagate from the surface to the core, one of the rotational relaxation time components decays faster than the bulk. In general, the relaxation of spins is nonexponential due to the interplay between various interactions. In the limit of strong coupling between the spins or in the limit of low temperature, the nature of relaxation of the spins undergoes a qualitative change with the emergence of a homogeneous dynamics where decay is predominantly exponential, again in agreement with experiments. (C) 2010 American Institute of Physics. doi: 10.1063/1.3474948]

High-temperature dielectric response in pulsed laser deposited Bi1.5Zn1.0Nb1.5O7 thin films

Cubic pyrochlore Bi1.5Zn1.0Nb1.5O7 thin films were deposited by pulsed laser ablation on Pt(200)/SiO2/Si at 500, 550, 600, and 650 degrees C. The thin films with (222) preferred orientation were found to grow at 650 degrees C with better crystallinity which was established by the lowest full-width half maxima of similar to 0.38. The dielectric response of the thin films grown at 650 degrees C have been characterized within a temperature range of 270-650 K and a frequency window of 0.1-100 kHz. The dielectric dispersion in the thin films shows a Maxwell-Wagner type relaxation with two different kinds of response confirmed by temperature dependent Nyquist plots. The ac conduction of the films showed a varied behavior in two different frequency regions. The power law exponent values of more than 1 at high frequency are explained by a jump-relaxation-model. The possibility of grain boundary related large polaronic hopping, due to two different power law exponents and transformation of double to single response in Nyquist plots at high temperature, has been excluded. The ``attempt jump frequency'' obtained from temperature dependent tangent loss and real part of dielectric constants, has been found to lie in the range of their lattice vibronic frequencies (10(12)-10(13) Hz). The activation energy arising from a large polaronic hopping due to trapped charge at low frequency region has been calculated from the ac conduction behavior. The range of activation energies (0.26-0.59. eV) suggests that the polaronic hopping at low frequency is mostly due to oxygen vacancies. (C) 2010 American Institute of Physics. doi:10.106311.3457335]

The Blackett Memorial Lecture, 1991: Chemical Insights into High-Temperature Superconductors

The high-temperature superconductors are complex oxides, generally containing two-dimensional CuO2 sheets. Various families of the cuprate superconductors are described, paying special attention to aspects related to oxygen stoichiometry, phase stability, synthesis and chemical manipulation of charge carriers. Other aspects discussed are chemical applications of cuprates, possibly as gas sensors and copper-free oxide superconductors. All but the substituted Nd and Pr cuprates are hole-superconductors. Several families of cuprates show a nearly constant n(h) at maximum T(c). Besides this universality, the cuprates exhibit a number of striking common features. Based on Cu(2p) photoemission studies, it is found that the Cu-O charge-transfer energy, DELTA, and the Cu(3d)-O(2p) hybridization strength, t(pd), are key factors in the superconductivity of cuprates. The relative intensity of the satellite in the Cu(2p) core-level spectra, the polarizability of the CuO2 sheets as well as the hole concentration are related to DELTA/t(pd). These chemical bonding factors have to be explicitly taken into account in any model for superconductivity of the cuprates.

Sensitive high-resolution ion microprobe U-Pb dating of prograde and retrograde ultrahigh-temperature metamorphism as exemplified by Sri Lankan granulites

Ultrahigh-temperature (UHT) granulites of the central Highland Complex, Sri Lanka, underwent some of the highest known peak temperatures of crustal metamorphism. Zircon and monazite U-Pb systems in granulites near Kandy, the highest grade region (similar to 1050 degrees C; 0.9 GPa), preserve both a record of the timing of prograde and retrograde phases of UHT metamorphism and evidence for the ages of older protolith components. Zircon grains from a quartz-saturated granulite containing relics of the peak UHT assemblage have remnant detrital cores with dates of ca. 2.5-0.83 Ga. Date clusters of ca. 1.7 and 1.04-0.83 Ga record episodes of zircon growth in the source region of the protolith sediment. Two generations of overgrowths with contrasting Th/U record metamorphic zircon growth at 569 +/- 5 and 551 +/- 7 Ma, probably in the absence and presence of monazite, respectively. The age of coexisting metamorphic monazite (547 +/- 7 Ma) is indistinguishable from that of the younger, low-Th/U zircon overgrowths. Zircon from a quartz-undersaturated monazite-absent UHT granulite with a mainly retrograde assemblage is mostly metamorphic (551 +/- 5 Ma). The ca. 570 Ma zircon overgrowths in the quartz-saturated granulite probably record partial melting just before or at the metamorphic peak. The ca. 550 Ma zircon in both rocks, and the ca. 550 Ma monazite in the quartz-saturated sample, record post-peak isothermal decompression. A possible model for this pressure-temperature-time evolution is ultrahot collisional orogeny during the assembly of Gondwana, locally superheated by basaltic underplating, followed by fast extensional exhumation.

Temperature dependent transport behavior of n-InN nanodot/p-Si heterojunction structures

The present work explores the temperature dependent transport behavior of n-InN nanodot/p-Si(100) heterojunction diodes. InN nanodot (ND) structures were grown on a 20 nm InN buffer layer on p-Si(100) substrates. These dots were found to be single crystalline and grown along 001] direction. The junction between these two materials exhibits a strong rectifying behavior at low temperatures. The average barrier height (BH) was determined to be 0.7 eV from current-voltage-temperature, capacitance-voltage, and flat band considerations. The band offsets derived from built-in potential were found to be Delta E-C=1.8 eV and Delta E-V=1.3 eV and are in close agreement with Anderson's model. (C) 2010 American Institute of Physics. doi:10.1063/1.3517489]

A spiral growth model for quasicrystals in two dimensions

Spiral space filling geometrical constructions using rhombuses in two dimensions are considered as plausible mechanisms for quasicrystal growth. These models will show staircase-like features which may be observed experimentally.

Two-Fermi liquid model for high Tc superconductivity involving suppression of tunnelling by decoherence

We consider a model system of two interacting Fermi-liquids, one of which is light and the other much heavier. In the normal state the lighter component provides a quantum mechanical bath coupled 'ohmically' to the heavier component in the sense of Caldeira and Leggett, suppressing thereby the band (tunnelling) matrix elements of the heavier component. Thus we lose the energy of delocalization. On the other hand, a superconducting ordering stiffens the bath spectral function at low energies and so restores the tunnelling. The resulting regain of the delocalization energy bootstraps so as to stabilize the superconducting order that caused it. It is conceivable that the motions parallel to the easy ab-plane and along the hard c-axis may also effectively correspond to the light and the heavy Fermi-liquids, respectively.

A Theory of the Lifted Temperature Minimum on Calm Clear Nights

Numerous reports from several parts of the world have confirmed that on calm clear nights a minimum in air temperature can occur just above ground, at heights of the order of $\frac{1}{2}$ m or less. This phenomenon, first observed by Ramdas & Atmanathan (1932), carries the associated paradox of an apparently unstable layer that sustains itself for several hours, and has not so far been satisfactorily explained. We formulate here a theory that considers energy balance between radiation, conduction and free or forced convection in humid air, with surface temperature, humidity and wind incorporated into an appropriate mathematical model as parameters. A complete numerical solution of the coupled air-soil problem is used to validate an approach that specifies the surface temperature boundary condition through a cooling rate parameter. Utilizing a flux-emissivity scheme for computing radiative transfer, the model is numerically solved for various values of turbulent friction velocity. It is shown that a lifted minimum is predicted by the model for values of ground emissivity not too close to unity, and for sufficiently low surface cooling rates and eddy transport. Agreement with observation for reasonable values of the parameters is demonstrated. A heuristic argument is offered to show that radiation substantially increases the critical Rayleigh number for convection, thus circumventing or weakening Rayleigh-Benard instability. The model highlights the key role played by two parameters generally ignored in explanations of the phenomenon, namely surface emissivity and soil thermal conductivity, and shows that it is unnecessary to invoke the presence of such particulate constituents as haze to produce a lifted minimum.

Phase diagram of the two-dimensional disordered Hubbard model in the Hartree-Fock approximation

Ground-state properties of the two-dimensional Hubbard model with point-defect disorder are investigated numerically in the Hartree-Fock approximation. The phase diagram in the p(point defect concentration)-delta(deviation from half filling) plane exhibits antiferromagnetic, spin-density-wave, paramagnetic, and spin-glass-like phases. The disorder stabilizes the antiferromagnetic phase relative to the spin-density-wave phase. The presence of U strongly enhances the localization in the antiferromagnetic phase. The spin-density-wave and spin-glass-like phases are weakly localized.

Dense granular flow down an inclined plane A comparison between the hard particle model and soft particle simulations

The granular flow down an inclined plane is simulated using the discrete element (DE) technique to examine the extent to which the dynamics of an unconfined dense granular flow can be well described by a hard particle model First, we examine the average coordination number for the particles in the flow down an inclined plane using the DE technique using the linear contact model with and without friction, and the Hertzian contact model with friction The simulations show that the average coordination number decreases below 1 for values of the spring stiffness corresponding to real materials, such as sand and glass, even when the angle of inclination is only 10 larger than the angle of repose Additional measures of correlations in the system, such as the fraction of particles with multibody contact, the force ratio (average ratio of the magnitudes of the largest and the second largest force on a particle), and the angle between the two largest forces on the particle, show no evidence of force chains or other correlated motions in the system An analysis of the bond-orientational order parameter indicates that the flow is in the random state, as in event-driven (ED) simulations V Kumaran, J Fluid Mech 632, 107 (2009), J Fluid Mech 632, 145 (2009)] The results of the two simulation techniques for the Bagnold coefficients (ratio of stress and square of the strain rate) and the granular temperature (mean square of the fluctuating velocity) are compared with the theory V Kumaran, J Fluid Mech 632, 107 (2009), J Fluid Mech 632, 145 (2009)] and are found to be in quantitative agreement In addition, we also conduct a comparison of the collision frequency and the distribution of the precollisional relative velocities of particles in contact The strong correlation effects exhibited by these two quantities in event-driven simulations V Kumaran, J Fluid Mech 632, 145 (2009)] are also found in the DE simulations (C) 2010 American Institute of Physics doi 10 1063/1 3504660]

Resolution of conflicting views on thermodynamics of glass transition: A unified model

Classical description of thermodynamic properties during glass transition has been questioned by the entropy-loss model. The uncompensated loss of entropy at the glass transition temperature and zero residual entropy is at the heart of the controversy. Both the models are critically reviewed. A unified model is presented which incorporates features of both entropy loss and residual entropy. It implies two different types of contributions to the entropy of the supercooled liquid, one of which vanishes at the transition and the other which contributes to residual entropy. Entropy gain during spontaneous relaxation of glass, and the nature of heat capacity `hysteresis' during cooling and heating through the glass transition range support the proposed model. Experiments are outlined for differentiating between the models.