Effect of quantum nuclear motion on hydrogen bonding

This work considers how the properties of hydrogen bonded complexes, X-H center dot center dot center dot Y, are modified by the quantum motion of the shared proton. Using a simple two-diabatic state model Hamiltonian, the analysis of the symmetric case, where the donor (X) and acceptor (Y) have the same proton affinity, is carried out. For quantitative comparisons, a parametrization specific to the O-H center dot center dot center dot O complexes is used. The vibrational energy levels of the one-dimensional ground state adiabatic potential of the model are used to make quantitative comparisons with a vast body of condensed phase data, spanning a donor-acceptor separation (R) range of about 2.4-3.0 angstrom, i.e., from strong to weak hydrogen bonds. The position of the proton (which determines the X-H bond length) and its longitudinal vibrational frequency, along with the isotope effects in both are described quantitatively. An analysis of the secondary geometric isotope effect, using a simple extension of the two-state model, yields an improved agreement of the predicted variation with R of frequency isotope effects. The role of bending modes is also considered: their quantum effects compete with those of the stretching mode for weak to moderate H-bond strengths. In spite of the economy in the parametrization of the model used, it offers key insights into the defining features of H-bonds, and semi-quantitatively captures several trends. (C) 2014 AIP Publishing LLC.

Structural crossover from nonmodulated to long-period modulated tetragonal phase and anomalous change in ferroelectric properties in the lead-free piezoelectric Na1/2Bi1/2TiO3-BaTiO3

The highly complex structure-property interrelationship in the lead-free piezoelectric (x) Na1/2Bi1/2TiO3 - (1 - x) BaTiO3 is a subject of considerable contemporary debate. Using comprehensive x-ray, neutron diffraction, dielectric, and ferroelectric studies, we have shown the existence of a new criticality in this system at x = 0.80, i.e., well within the conventional tetragonal phase field. This criticality manifests as a nonmonotonic variation of the tetragonality and coercivity and is shown to be associated with a crossover from a nonmodulated tetragonal phase (for x < 0.8) to a long-period modulated tetragonal phase (for x > 0.80). It is shown that the stabilization of long-period modulation introduces a characteristic depolarization temperature in the system. While differing qualitatively from the two-phase model often suggested for the critical compositions of this system, our results support the view with regard to the tendency in perovskites to stabilize long-period modulated structures as a result of complex interplay of antiferrodistortive modes Bellaiche and Iniguez, Phys. Rev. B 88, 014104 ( 2013); Prosandeev, Wang, Ren, Iniguez, ands Bellaiche, Adv. Funct. Mater. 23, 234 (2013)].

Theory of the non-Fermi-liquid transition point in the two-impurity Kondo model

We present an explicit solution of the problem of two coupled spin-1/2 impurities, interacting with a band of conduction electrons. We obtain an exact effective bosonized Hamiltonian, which is then treated by two different methods (low-energy theory and mean-field approach). Scale invariance is explicitly shown at the quantum critical point. The staggered susceptibility behaves like ln(T(K)/T) at low T, whereas the magnetic susceptibility and [S1.S2] are well behaved at the transition. The divergence of C(T)/T when approaching the transition point is also studied. The non-Fermi-liquid (actually marginal-Fermi-liquid) critical point is shown to arise because of the existence of anomalous correlations, which lead to degeneracies between bosonic and fermionic states of the system. The methods developed in this paper are of interest for studying more physically relevant models, for instance, for high-T(c) cuprates.

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.

Absolute Entropy and Energy of Carbon Dioxide Using the Two-Phase Thermodynamic Model

The two-phase thermodynamic (2PT) model is used to determine the absolute entropy and energy of carbon dioxide over a wide range of conditions from molecular dynamics trajectories. The 2PT method determines the thermodynamic properties by applying the proper statistical mechanical partition function to the normal modes of a fluid. The vibrational density of state (DoS), obtained from the Fourier transform of the velocity autocorrelation function, converges quickly, allowing the free energy, entropy, and other thermodynamic properties to be determined from short 20-ps MD trajectories. The anharmonic effects in the vibrations are accounted for by the broadening of the normal modes into bands from sampling the velocities over the trajectory. The low frequency diffusive modes, which lead to finite DoS at zero frequency, are accounted for by considering the DoS as a superposition of gas-phase and solid-phase components (two phases). The analytical decomposition of the DoS allows for an evaluation of properties contributed by different types of molecular motions. We show that this 2PT analysis leads to accurate predictions of entropy and energy of CO2 over a wide range of conditions (from the triple point to the critical point of both the vapor and the liquid phases along the saturation line). This allows the equation of state of CO2 to be determined, which is limited only by the accuracy of the force field. We also validated that the 2PT entropy agrees with that determined from thermodynamic integration, but 2PT requires only a fraction of the time. A complication for CO2 is that its equilibrium configuration is linear, which would have only two rotational modes, but during the dynamics it is never exactly linear, so that there is a third mode from rotational about the axis. In this work, we show how to treat such linear molecules in the 2PT framework.

A numerical model of a liquid-feed solid polymer electrolyte DMFC and its experimental validation

A one-dimensional, biphasic, multicomponent steady-state model based on phenomenological transport equations for the catalyst layer, diffusion layer, and polymeric electrolyte membrane has been developed for a liquid-feed solid polymer electrolyte direct methanol fuel cell (SPE- DMFC). The model employs three important requisites: (i) implementation of analytical treatment of nonlinear terms to obtain a faster numerical solution as also to render the iterative scheme easier to converge, (ii) an appropriate description of two-phase transport phenomena in the diffusive region of the cell to account for flooding and water condensation/evaporation effects, and (iii) treatment of polarization effects due to methanol crossover. An improved numerical solution has been achieved by coupling analytical integration of kinetics and transport equations in the reaction layer, which explicitly include the effect of concentration and pressure gradient on cell polarization within the bulk catalyst layer. In particular, the integrated kinetic treatment explicitly accounts for the nonhomogeneous porous structure of the catalyst layer and the diffusion of reactants within and between the pores in the cathode. At the anode, the analytical integration of electrode kinetics has been obtained within the assumption of macrohomogeneous electrode porous structure, because methanol transport in a liquid-feed SPE- DMFC is essentially a single-phase process because of the high miscibility of methanol with water and its higher concentration in relation to gaseous reactants. A simple empirical model accounts for the effect of capillary forces on liquid-phase saturation in the diffusion layer. Consequently, diffusive and convective flow equations, comprising Nernst-Plank relation for solutes, Darcy law for liquid water, and Stefan-Maxwell equation for gaseous species, have been modified to include the capillary flow contribution to transport. To understand fully the role of model parameters in simulating the performance of the DMCF, we have carried out its parametric study. An experimental validation of model has also been carried out. (C) 2003 The Electrochemical Society.

Model of resistances in systems of Tomonaga-Luttinger liquid wires

In a recent paper, we combined the technique of bosonization with the concept of a Rayleigh dissipation function to develop a model for resistances in one-dimensional systems of interacting spinless electrons Europhys. Lett. 93, 57007 (2011)]. We also studied the conductance of a system of three wires by using a current splitting matrix M at the junction. In this paper, we extend our earlier work in several ways. The power dissipated in a three-wire system is calculated as a function of M and the voltages applied in the leads. By combining two junctions of three wires, we examine a system consisting of two parallel resistances. We study the conductance of this system as a function of the M matrices and the two resistances; we find that the total resistance is generally quite different from what one expects for a classical system of parallel resistances. We do a sum over paths to compute the conductance of this system when one of the two resistances is taken to be infinitely large. We study the conductance of a three-wire system of interacting spin-1/2 electrons, and show that the charge and spin conductances can generally be different from each other. Finally, we consider a system of two wires that are coupled by a dissipation function, and we show that this leads to a current in one wire when a voltage bias is applied across the other wire.

Interference of Two Closely Spaced Strip Footings on Sand Using Model Tests

By using small scale model tests, the interference effect on the ultimate bearing capacity of two closely spaced strip footings, placed on the surface of dry sand, was investigated. At any time, the footings were assumed to (1) carry exactly the same magnitude of load; and (2) settle to the same extent. No tilt of the footing was allowed. The effect of clear spacing (s) between two footings was explicitly studied. An interference of footings leads to a significant increase in their bearing capacity; the interference effect becomes even more substantial with an increase in the relative density of sand. The bearing capacity attains a peak magnitude at a certain (critical) spacing between two footings. The experimental observations presented in this technical note were similar to those given by different available theories. However, in a quantitative sense, the difference between the experiments and theories was seen to be still significant and it emphasizes the need of doing a further rigorous analysis in which the effect of stress level on the shear strength parameters of soil mass can be incorporated properly.

Testing approximate theories of first-order phase transitions on the two-dimensional Potts model

The two-dimensional,q-state (q>4) Potts model is used as a testing ground for approximate theories of first-order phase transitions. In particular, the predictions of a theory analogous to the Ramakrishnan-Yussouff theory of freezing are compared with those of ordinary mean-field (Curie-Wiess) theory. It is found that the Curie-Weiss theory is a better approximation than the Ramakrishnan-Yussouff theory, even though the former neglects all fluctuations. It is shown that the Ramakrishnan-Yussouff theory overestimates the effects of fluctuations in this system. The reasons behind the failure of the Ramakrishnan-Yussouff approximation and the suitability of using the two-dimensional Potts model as a testing ground for these theories are discussed.

Coalescence of drops in stirred dispersion. A white noise model for coalescence

Coalescence between two droplets in a turbulent liquid-liquid dispersion is generally viewed as a consequence of forces exerted on the drop-pair squeezing out the intervening continuous phase to a critical thickness. A new synthesis is proposed herein which models the film drainage as a stochastic process driven by a suitably idealized random process for the fluctuating force. While the true test of the model lies in detailed parameter estimations with measurement of drop-size distributions in coalescing dispersions, experimental measurements on average coalescence frequencies lend preliminary support to the model.

The hole spin model for oxides

A model of mobile 0-holes hybrized with Cu-spins on a square lattice is examined. A variational groundstate wavefunction which interpolates smoothly between n.n. RVB and Néel limits gives a Néellike minimum. A hole in an AF lattice polarizes it locally and becomes quite mobile. Two n.n. holes attract. Finally we speculate how holes can stabilize a spin liquid state.

Thermodynamic scaling theory for the two-impurity Anderson model

We present results of a study of the two-impurity Anderson model using a thermodynamic scaling theory developed recently. The model is characterized by the Coulomb energy U, the orbital energy epsilond, the d-level width Gamma, and the separation between impurities R. If Gamma<<−epsilond<~Gamma. Here we find that the single-impurity physics dominates the low-temperature behavior, and impurity-impurity interactions are perturbative. The qualitative features of the temperature-dependent susceptibility are discussed. Journal of Applied Physics is copyrighted by The American Institute of Physics.

A new kinetic model for bulk polymerization of vinyl chloride based on two-phase hypothesis

A kinetic model has been developed for the bulk polymerization of vinyl chloride using Talamini's hypothesis of two-phase polymerization and a new concept of kinetic solubility which assumes that rapidly growing polymer chains have considerably greater solubility than the thermodynamic solubility of preformed polymer molecules of the same size and so can remain in solution even under thermodynamically unfavourable conditions. It is further assumed that this kinetic solubility is a function of chain length. The model yields a rate expression consistent with the experimental data for vinyl chloride bulk polymerization and moreover is able to explain several characteristic kinetic features of this system. Application of the model rate expression to the available rate data has yielded 2.36 × 108l mol−1 sec−1 for the termination rate constant in the polymer-rich phase; as expected, this value is smaller than that reported for homogenous polymerization by a factor of 10–30.

Line tension of a two dimensional gas-liquid interface

In two dimensional (2D) gas-liquid systems, the reported simulation values of line tension are known to disagree with the existing theoretical estimates. We find that while the simulation erred in truncating the range of the interaction potential, and as a result grossly underestimated the actual value, the earlier theoretical calculation was also limited by several approximations. When both the simulation and the theory are improved, we find that the estimate of line tension is in better agreement with each other. The small value of surface tension suggests increased influence of noncircular clusters in 2D gas-liquid nucleation, as indeed observed in a recent simulation.

Rate Model and Mechanism of Liquid-Phase Oxidation of Propionaldehyde

A rate equation is developed for the liquid-phase oxidation of propionaldehyde with oxygen in the presence of manganese propionate catalyst in a sparged reactor. The equation takes into account diffusional limitations based on Brian's solution for mass transfer accompanied by a pseudo m-. nth-order reaction. Sauter-mean bubble diameter, gas holdup, interfacial area, and bubble rise velocity are measured, and rates of mass transfer within the gas phase and across the gas-liquid interface are computed. Statistically designed experiments show the adequacy of the equation. The oxidation reaction is zero order with respect to oxygen concentration, 3/2 order with respect to aldehyde concentration, and order with respect to catalyst concentration. The activation energy is 12.1 kcal/g mole.