Solvation dynamics, energy distribution and trapping of a light solute ion

In the theoretical treatments of the dynamics of solvation of a newly created ion in a dipolar solvent, the self-motion of the solute is usually ignored. Recently, it has been shown that for a light ion the translational motion of the ion can significantly enhance its own rate of solvation. Therefore, solvation itself may not be the rate determining step in the equilibration. Instead, the rate determining step is the search of the low energy configuration which serves to localize the light ion. In this article a microscopic calculation of the probability distribution of the interaction energy of the nascent charge with the dipolar solvent molecules is presented in order to address this problem of solute trapping. It is found that to a good approximation, this distribution is Gaussian and the second moment of this distribution is exactly equal to the half of its own solvation energy. It is shown that this is in excellent agreement with the simulation results that are available for the model Brownian dipolar lattice and for liquid acetonitrile. If the distortion of the solvent by the ion is negligible then the same relation gives the energy distribution for the solvated ion, with the average centered at the final equilibrium solvation energy. These results are expected to be useful in understanding various chemical processes in dipolar liquids. Another interesting outcome of the present study is a simple dynamic argument that supports Onsager's ''inverse snow-ball'' conjecture of solvation of a light ion. A simple derivation of the semi-phenomenological relation between the solvation time correlation function and the single particle orientation, reported recently by Maroncelli et al. (J. Phys. Chem. 97 (1993) 13), is also presented.

Gas-liquid nucleation at large metastability: unusual features and a new formalism

Nucleation at large metastability is still largely an unsolved problem, even though it is a problem of tremendous current interest, with wide-ranging practical value, from atmospheric research to materials science. It is now well accepted that the classical nucleation theory (CNT) fails to provide a qualitative picture and gives incorrect quantitative values for such quantities as activation-free energy barrier and supersaturation dependence of nucleation rate, especially at large metastability. In this paper, we present an alternative formalism to treat nucleation at large supersaturation by introducing an extended set of order parameters in terms of the kth largest liquid-like clusters, where k = 1 is the largest cluster in the system, k = 2 is the second largest cluster and so on. At low supersaturation, the size of the largest liquid-like cluster acts as a suitable order parameter. At large supersaturation, the free energy barrier for the largest liquid-like cluster disappears. We identify this supersaturation as the one at the onset of kinetic spinodal. The kinetic spinodal is system-size-dependent. Beyond kinetic spinodal many clusters grow simultaneously and competitively and hence the nucleation and growth become collective. In order to describe collective growth, we need to consider the full set of order parameters. We derive an analytic expression for the free energy of formation of the kth largest cluster. The expression predicts that, at large metastability (beyond kinetic spinodal), the barrier of growth for several largest liquid-like clusters disappears, and all these clusters grow simultaneously. The approach to the critical size occurs by barrierless diffusion in the cluster size space. The expression for the rate of barrier crossing predicts weaker supersaturation dependence than what is predicted by CNT at large metastability. Such a crossover behavior has indeed been observed in recent experiments (but eluded an explanation till now). In order to understand the large numerical discrepancy between simulation predictions and experimental results, we carried out a study of the dependence on the range of intermolecular interactions of both the surface tension of an equilibrium planar gas-liquid interface and the free energy barrier of nucleation. Both are found to depend significantly on the range of interaction for the Lennard-Jones potential, both in two and three dimensions. The value of surface tension and also the free energy difference between the gas and the liquid phase increase significantly and converge only when the range of interaction is extended beyond 6-7 molecular diameters. We find, with the full range of interaction potential, that the surface tension shows only a weak dependence on supersaturation, so the reason for the breakdown of CNT (with simulated values of surface tension and free energy gap) cannot be attributed to the supersaturation dependence of surface tension. This remains an unsettled issue at present because of the use of the value of surface tension obtained at coexistence.

Liquid alloy semiconductors revisited - Discussion

In 1990 Enderby and Barnes reviewed the electrical properties of liquid alloys which show features associated with semiconducting behaviour. They proposed an empirical classification scheme based on the notion that some liquid semiconductors are characterized by a finite gap in σ (E), the energy-dependent conductivity (narrow definition), whereas in others σ (E) is continuous (broad definition). Interesting behaviour occurs for systems at the narrow/broad boundary and further analysis of these liquid alloys will form the subject matter of this paper. Particular attention will be focused on liquid silver chalcogenides as these offer a severe test of current theories.

Measurement of Gibbs energy of formation of Ca2PbO4 using a solid-state cell with three electrodes

Phase relations in the system Ca-Pb-O at 1100 K have been determined by equilibrating 18 compositions in the ternary and identifying the phases present in quenched samples by X-ray diffraction and energy dispersive X-ray analysis (EDX). Only one ternary compound Ca2PbO4 was found to be present. The compound coexists with CaO and PbO. The intermetallic compounds Ca2Pb, Ca5Pb3 and CaPb and liquid alloys are in equilibrium with CaO. The standard Gibbs energies of formation of Ca2PbO4 (880 - 1100 K) and Pb3O4 (770 - 910 K) were determined using solid-state cells based on yttria-stabilized zirconia as the solid electrolyte. Pure oxygen gas at 0.1 MPa was used as the reference electrode. For measurements on Ca2PbO4, a novel cell design with three electrodes in series, separated by solid electrolyte membranes, was used to avoid polarization of the electrode containing three solid phases. Two three-phase electrodes were used. The first absorbs the electrochemical flux of oxygen from the reference electrode to the measuring electrode. The other three-phase electrode, which is unaffected by the oxygen flux through the solid electrolyte, is used for electromotive force (EMF) measurement. The results from EMF studies were cross-checked using thermogravimetry (TG) under controlled oxygen partial pressures. The stability of Pb3O4 was investigated using a conventional solid-state cell with RuO2 electrodes. The results can be summarized by the following equations: 2CaO + PbO +1/2O(2) --> Ca2PbO4 Delta(r)G degrees/J mol(-1) = (- 128340 + 93.21 T/K) +/- 200 3PbO + 1/2O(2) --> Pb3O4 Delta(r)G degrees/J mol(-1) = (- 70060 + 77.5 T/K) +/- 150

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.

Saturated liquid thermal conductivity correlations for some new refrigerants

This paper presents a set of linear equations describing the temperature dependence of the saturated liquid thermal conductivity covering the region of engineering importance for the new hydrofluorocarbons (HFC) 32, 125, 134a, 143a, 152a and hydrochlorofluorocarbons (HCFC) 123, 124, 141b and 142b. Available experimental data in the literature have been considered to arrive at a correlation of the form lambda = A - BT. It is observed that there exists an appreciable discrepancy between various sources of data in spite of the same purity of samples used and the same measurement technique being adopted. The correlations obtained here could be useful in engineering design applications. (C) 1998 John Wiley & Sons, Ltd.

Subquadratic quantum number dependence and other anomalies of vibrational dephasing in liquid nitrogen: Molecular dynamics simulation study from the triple point to the critical point and beyond

Vibrational phase relaxation near gas-liquid and liquid-solid phase coexistence has been studied by molecular dynamics simulations of N-N stretch in N-2. Experimentally observed pronounced insensitivity of phase relaxation from the triple point to beyond the boiling point is found to originate from a competition between density relaxation and resonant-energy transfer terms. The sharp rise in relaxation rate near the critical point (CP) can be attributed at least partly to the sharp, rise in vibration-rotation coupling contribution. Substantial subquadratic quantum number dependence of overtone dephasing rate is found near the CP and in supercritical fluids. [S0031-9007 (99)09318-7].

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.

Size-dependent cohesive energy and melting of non-spherical nanoparticles

All most all theoretical models assume spherical nanoparticles. However, thermodynamic properties of non-spherical nanoparticles are the subject of recent interests. In this article, we have discussed the size-dependent cohesive energy and melting of non-spherical nanoparticles based on liquid-drop model. The surface to volume ratio is different for different shapes of nanoparticles and as a consequence, the variation of cohesive energy and melting of non-spherical nanoparticles is different from that of spherical case. By analyzing the reported experimental results, it has been observed that liquid-drop model can be used to understand the size-dependent cohesive energy and melting of non-spherical nanoparticles.

Energy-insensitive Guidance of Solid Motor Propelled Long Range Flight Vehicles Using MPSP and Dynamic Inversion

A energy-insensitive explicit guidance design is proposed in this paper by appending newlydeveloped nonlinear model predictive static programming technique with dynamic inversion, which render a closed form solution of the necessary guidance command update. The closed form nature of the proposed optimal guidance scheme suppressed the computational difficulties, and facilitate realtime solution. The guidance law is successfully verified in a solid motor propelled long range flight vehicle, for which developing an effective guidance law is more difficult as compared to a liquid engine propelled vehicle, mainly because of the absence of thrust cutoff facility. The scheme guides the vehicle appropriately so that it completes the mission within a tight error bound assuming that the starting point of the second stage to be a deterministic point beyond the atmosphere. The simulation results demonstrate its ability to intercept the target, even with an uncertainty of greater than 10% in the burnout time

A Hybrid Energy-insensitive Explicit Guidance Scheme for Long Range Flight Vehicles with Solid Motors

Combining the newly developed nonlinear model predictive static programming technique with null range direction concept, a novel explicit energy-insensitive guidance design method is presented in this paper for long range flight vehicles, which leads to a closed form solution of the necessary guidance command update. Owing to the closed form nature, it does not lead to computational difficulties and the proposed optimal guidance algorithm can be implemented online. The guidance law is verified in a solid motor propelled long range flight vehicle, for which coming up with an effective guidance law is more difficult as compared to a liquid engine propelled vehicle (mainly because of the absence of thrust cutoff facility). Assuming the starting point of the second stage to be a deterministic point beyond the atmosphere, the scheme guides the vehicle properly so that it completes the mission within a tight error bound. The simulation results demonstrate its ability to intercept the target, even with an uncertainty of greater than 10% in burnout time.

Solubility and activity of oxygen in liquid gallium and gallium copper alloys

The solubility of oxygen in liquid gallium in the temperature range 775 –1125 °C and in liquid gallium-copper alloys at 1100 °C, in equilibrium with β-Ga2O3, has been measured by an isopiestic equilibrium technique. The solubility of oxygen in pure gallium is given by the equation log (at.% O) = −7380/T + 4.264 (±0.03) Using recently measured values for the standard free energy of formation of β-Ga2O3 and assuming that oxygen obeys Sievert's law up to the saturation limit, the standard free energy of solution of oxygen in liquid gallium may be calculated : View the MathML sourceΔ°298 = −52 680 + 6.53T (±200) cal where the standard state for dissolved oxygen is an infinitely dilute solution in which the activity is equal to atomic per cent. The effect of copper on the activity of oxygen dissolved in liquid gallium is found to be in good agreement with that predicted by a recent quasichemical model in which it was assumed that each oxygen is interstitially coordinated to four metal atoms and that the nearest neighbour metal atoms lose approximately half their metallic cohesive energies.

Solubility and activity of oxygen in liquid indium and copper indium alloys

The solubility of oxygen in liquid indium in the temperature range 650–820 °C and in liquid copper-indium alloys at 1100 °C in equilibrium with indium sesquioxide has been measured by a phase equilibration technique. The solubility of oxygen in pure indium is given by the relation log(at.% O) = −4726/T + 3.73 (±0.08) Using the recently measured values for the standard free energy of formation of In2O3 and assuming that oxygen obeys Sievert's law up to saturation, the standard free energy of solution of molecular oxygen in liquid indium is calculated as View the MathML sourceΔG°= −51 440 + 8.07 T (±500) cal where the standard state for dissolved oxygen is an infinitely dilute solution in which activity is equal to atomic per cent. The effect of indium additions on the activity coefficient of oxygen dissolved in liquid copper was measured by a solid oxide galvanic cell. The interaction parameter ϵ0In is given by View the MathML source The experimentally determined variation of the activity coefficient of oxygen in dilute solution in Cu-In alloys is in fair agreement with that predicted by a quasichemical model in which each oxygen atom is assumed to be interstitially coordinated to four metal atoms and the nearest neighbour metal atoms are assumed to lose approximately half their metallic cohesive energies.

Solubility and activity of oxygen in liquid germanium and germanium copper alloys

The solubility of oxygen in liquid germanium in the temperature range 1233 to 1397 K, and in liquid germanium-copper alloys at 1373 K, in equilibrium with GeO2 has been measured by the phase equilibration technique. The solubility of oxygen in pure germanium is given by the relation R470 log(at. pct 0)=-6470/T+4.24 (±0.07). The standard free energy of solution of oxygen in liquid germanium is calculated from the saturation solubility, and recently measured values for the free energy of formation of GeO2, assuming that oxygen obeys Sievert’s law up to the saturation limit. For the reaction, 1/2 O2(g)→ OGe ΔG° =-39,000 + 3.21T (±500) cal = -163,200 + 13.43T (±2100) J. where the standard state for dissolved oxygen is that which makes the value of activity equal to the concentration (in at. pct), in the limit, as concentration approaches zero. The effect of copper on the activity of oxygen dissolved in liquid germanium is found to be in good agreement with that predicted by a quasichemical model in which each oxygen was assumed to be bonded to four metal atoms and the nearest neighbor metal atoms to an oxygen atom are assumed to lose approximately half of their metallic bonds.

Vapor pressure and dissociation energy of (In2O)

The vapor pressure of pure liquid indium, and the sum of pressures of (In) and (In2O) species over the condensed phase mixture {In} + , contained in a silica vessel, have been measured by Knudsen effusion and Langmuir free vaporization methods in the temperatue range 600 to 950°C. Mass spectrometric studies reported in the literature show that (In) and (In2O) are the important species in the vapor phase over the {In} + ; mixture. The vapor pressure of (In2O) corresponding to the reaction, deduced from the present measurements is given by the equation, The “apparent evaporation coefficient” for the condensed phase mixture is approximately 0.8. The energy for the dissociation (In2O) molecule into atoms calculated from the above equation is D°0 = 180.0 (± 1.0) kcal mol−1.