Temperature and concentration dependence of adsorption properties of methane in NaY: a molecular dynamics study

Molecular dynamics calculations on methane sorbed in NaY (Si/Al = 3.0) employing realistic methane-methane and methane-zeolite intermolecular potential functions at different temperatures (50, 150, 220, and 300 K) and concentrations (2, 4, 6, and 8 molecules/cage) are reported. The thermodynamic results are in agreement with the available experimental data. Guest-guest and guest-host radial distribution functions (rdfs), energy distribution functions, distribution of cage occupancy, center-of-cage-center-of-mass (coc-com) rdfs, velocity autocorrelation functions for com and angular motion and the Fourier transformed power spectra, and diffusion coefficients are presented as a function of temperature and concentration. At 50 K, methane is localized near the adsorption site. Site-site migration and essentially free rotational motion are observed at 150 K. Molecules preferentially occupy the region near the inner surface of the alpha-cage. The vibrational frequencies for the com of methane shift toward higher values with decreasing temperature and increasing adsorbate concentration. The observed frequencies for com motion are 36, 53, and 85 cm-1 and for rotational motion at 50 K, 95 and 150 cm-1 in agreement with neutron scattering data. The diffusion coefficients show a type I behavior as a function of loading in agreement with NMR measurements. Cage-to-cage diffusion is found to be always mediated by the surface.

Sorbate properties and cage-to-cage diffusion of argon in NaCaA: a molecular dynamics study

Detailed molecular dynamics simulations of argon in zeolite NaCaA are reported. Thermodynamic, structural, and dynamical properties of the sorbate as a function of temperature have been obtained. The properties calculated include various site-site radial distribution functions, different energy distribution functions, selfdiffusion coefficients, the power spectra, and properties relating to cage-to-cage diffusion. The results suggest that sorbate is delocalized above 300 K. Both modes of cage-to-cage diffusion-the surface-mediated and centralized diffusion-are associated with negative barrier heights. Surprisingly, rate of cage-to-cage diffusion is associated with negative and positive activation energies below and above 500 K. The observed differences in the behavior of the rate of cage-to-cage diffusion between Xe-NaY and Ar-NaCaA systems and the nature of the potential energy surface are discussed. Presence of sorbatezeolite interactions results in significant enhancement in the rate of cage-to-cage diffusion and rate of cage visits. It is shown that properties dependent on the long-time behavior such as the diffusion coefficient and the rate of cages visited exhibit the expected Arrhenius dependence on temperature.

Thermodynamic and structural consequences of changing a sulfur atom to a methylene group in the M13Nle mutation in ribonuclease-S

Two fragments of pancreatic ribonuclease A, a truncated version of S-peptide (residues 1-15) and S-protein (residues 21-124), combine to give a catalytically active complex. We have substituted the wild-type residue at position 13, methionine (Met), with norleucine (Nle), where the only covalent change is the replacement of the sulfur atom with a methylene group. The thermodynamic parameters associated with the binding of this variant to S-protein, determined by titration calorimetry in the temperature range 10-40 degrees C, are reported and compared to values previously reported [Varadarajan, R., Connelly, P. R., Sturtevant, J. M., & Richards, F. M. (1992) Biochemistry 31, 1421-1426] for other position 13 analogs. The differences in the free energy and enthalpy of binding between the Met and Nle peptides are 0.6 and 7.9 kcal/mol at 25 degrees C, respectively. These differences are slightly larger than, but comparable to, the differences in the values for the Met/Ile and Met/Leu pairs. The structure of the mutant complex was determined to 1.85 Angstrom resolution and refined to an R-factor of 17.4% The structures of mutant and wild-type complexes are practically identical although the Nle side chain has a significantly higher average B-factor than the corresponding Met side chain. In contrast, the B-factors of the atoms of the cage of residues surrounding position 13 are all somewhat lower in the Nle variant than in the Met wild-type. Thus, the large differences in the binding enthalpy appear to reside entirely in the difference in chemical properties or dynamic behavior of the -S- and -CH2- groups and not in differences in the geometry of the side chains or the internal cavity surface. In addition, a novel method of obtaining protein stability data by means of isothermal titration calorimetry is introduced.

Dynamic and thermodynamic anomalies of water at low temperatures: from bulk water to reverse micelles and DNA hydration layer

Liquid water is known to exhibit remarkable thermodynamic and dynamic anomalies, ranging from solvation properties in supercritical state to an apparent divergence of the linear response functions at a low temperature. Anomalies in various dynamic properties of water have also been observed in the hydration layer of proteins, DNA grooves and inside the nanocavity, such as reverse micelles and nanotubes. Here we report studies on the molecular origin of these anomalies in supercooled water, in the grooves of DNA double helix and reverse micelles. The anomalies have been discussed in terms of growing correlation length and intermittent population fluctuation of 4- and 5-coordinated species. We establish correlation between thermodynamic response functions and mean squared species number fluctuation. Lifetime analysis of 4- and 5-coordinated species reveals interesting differences between the role of the two species in supercooled and constrained water. The nature and manifestations of the apparent and much discussed liquid-liquid transition under confinement are found to be markedly different from that in the bulk. We find an interesting `faster than bulk' relaxation in reverse micelles which we attribute to frustration effects created by competition between the correlations imposed by surface interactions and that imposed by hydrogen bond network of water.

A thermodynamic study of liquid Cu-Cr alloys and metastable liquid immiscibility

he thermodynamic acitivity of chromium in liquid Cu-Cr alloys is measured in the temperature range from 1473 to 1873 K using the solid state cell: Pt, W, Cr + Cr2O3 |(Y2O3) ThO2|Cu - Cr + Cr2O3, Pt The activity of copper and the Gibbs energy of mixing of the liquid alloy are derived. Activities exhibit large positive deviations from Raoult's law. The mixing properties can be represented by a pseudo-subregular solution model in which the excess entropy has the same type of functional dependence on composition as the enthalpy of mixing: ΔGE = XCr(1 - XCr)[60880 - 18750 XCr)-- T(16.25 - 7.55 XCr)]J mol-1 Pure liquid Cu and Cr are taken as the reference states. The results predict a liquid-liquid metastable miscibility gap, with TC = 1787 (±3) K and XCr = 0.436 (±0.02), lying below the liquidus. The results obtained in this study are in general agreement with experimental information reported in the literature, but provide further refinement of the thermodynamic parameters.

Thermodynamic, Raman and electrical switching studies on Si15Te85-xAgx (4 <= x <= 20) glasses

We have investigated thermal properties of bulk Si15Te85-xAgx (4 <= x <= 20) glasses in detail, through alternating differential scanning calorimetry experiments. The composition dependence of thermal parameters reveal the signatures of rigidity percolation and chemical threshold at compositions x = 12 and x = 19, respectively. The stability and glass forming ability of these glasses have also been determined using the data obtained from different thermodynamic quantities and it is found that the Si15Te85-xAgx glasses in the region 12 <= x <= 17 are more stable when compared to other glasses of the same series. Further, the blueshift observed in Raman spectroscopy investigations, in the composition range 12 <= x <= 13, support the occurrence of stiffness threshold in this composition range. All Si15Te85-xAgx (4 <= x <= 20) glasses are found to exhibit memory type switching (for sample thickness 0.25 mm) in the input current range 3-9 mA. The effect of rigidity percolation and chemical thresholds on switching voltages are observed at x = 12 and 19, respectively. (C) 2012 American Institute of Physics. [doi:10.1063/1.3682759]

Experiments on and thermodynamic analysis of a turbocharged engine with producer gas as fuel

This article presents the studies conducted on turbocharged producer gas engines designed originally for natural gas (NG) as the fuel. Producer gas, whose properties like stoichiometric ratio, calorific value, laminar flame speed, adiabatic flame temperature, and related parameters that differ from those of NG, is used as the fuel. Two engines having similar turbochargers are evaluated for performance. Detailed measurements on the mass flowrates of fuel and air, pressures and temperatures at various locations on the turbocharger were carried out. On both the engines, the pressure ratio across the compressor was measured to be 1.40 +/- 0.05 and the density ratio to be 1.35 +/- 0.05 across the turbocharger with after-cooler. Thermodynamic analysis of the data on both the engines suggests a compressor efficiency of 70 per cent. The specific energy consumption at the peak load is found to be 13.1 MJ/kWh with producer gas as the fuel. Compared with the naturally aspirated mode, the mass flow and the peak load in the turbocharged after-cooled condition increased by 35 per cent and 30 per cent, respectively. The pressure ratios obtained with the use of NG and producer gas are compared with corrected mass flow on the compressor map.

Role of electrolyte chemistry on electronic and in vitro electrochemical properties of micro-arc oxidized titania films on Cp Ti

The primary objective of the present work was to study the electronic and in vitro electrochemical properties of micro-arc oxidized titania films on Cp Ti, fabricated independently in various electrolyte solutions consisting of anions such as phosphate (PO43-), borate (B4O72-), citrate (C6H5O73-) and silicate (SiO32-). Further the role of anions on the structural, morphological and compositional properties of the fabricated films was studied. All the titania films were developed by micro-arc oxidation (MAO) technique for a fixed treatment time of 8 min under constant current mode. The surface morphology, elemental distribution, composition and structural characteristics of the films were assessed by scanning electron microscope (SEM) equipped with energy dispersive spectroscopy (EDS) and X-ray diffraction (XRD) techniques. The thermodynamic and kinetic corrosion properties of the films were studied under simulated body fluid (SBF) conditions (pH 7.4 and 37 degrees C) by conducting chronopotentiometric and potentiodynamic polarization tests. Electrochemical impedance spectroscopy (EIS) coupled with equivalent circuit modelling was carried out to analyse the frequency response and Mott-Schottky analysis was performed to study the semiconducting (electronic) properties of the films. Salt spray fog accelerated corrosion test was conducted for 168h as per ASTM B117 standard to corroborate the corrosion and semiconducting properties of the samples based on the visual examination. The XRD results showed that the transformation from the metastable anatase phase to the thermodynamically stable rutile phase and the crystalline growth of the respective phases were strongly influenced by the addition of anions. The SEM-EDS results demonstrated that the phosphorous (P) content in the films varied from 2.4 at% to 5.0 at% indicating that the amount of P in the films could be modified by adding an appropriate electrolyte additive. The electrochemical corrosion test results showed that the film fabricated in citrate (C6H5O73-) containing electrolyte is thermodynamically and kinetically more stable compared to that of all the others. The results of the Mott-Schottky analysis indicated that all the fabricated films showed an n-type semiconducting behaviour and the film developed in citrate (C6H5O73-) containing electrolyte exhibited the lowest donor concentration and the most negative flat band potential that contributed to its highest corrosion resistance in SBF solution. The results of the salt spray accelerated corrosion tests were in agreement with those obtained from the electrochemical and Mott-Schottky analysis.

Local structural disorder and its influence on the average global structure and polar properties in Na0.5Bi0.5TiO3

Na0.5Bi0.5TiO3 (NBT) and its derivatives have prompted a great surge in interest owing to their potential as lead-free piezoelectrics. In spite of five decades since its discovery, there is still a lack of clarity on crucial issues such as the origin of significant dielectric relaxation at room temperature, structural factors influencing its depoling, and the status of the recently proposed monoclinic (Cc) structure vis-a-vis the nanosized structural heterogeneities. In this work, these issues are resolved by comparative analysis of local and global structures on poled and unpoled NBT specimens using electron, x-ray, and neutron diffraction in conjunction with first-principles calculation, dielectric, ferroelectric, and piezoelectric measurements. The reported global monoclinic (Cc) distortion is shown not to correspond to the thermodynamic equilibrium state at room temperature. The global monocliniclike appearance rather owes its origin to the presence of local structural and strain heterogeneities. Poling removes the structural inhomogeneities and establishes a long-range rhombohedral distortion. In the process the system gets irreversibly transformed from a nonergodic relaxor to a normal ferroelectric state. The thermal depoling is shown to be associated with the onset of incompatible in-phase tilted octahedral regions in the field-stabilized long range rhombohedral distortion.

Origin of enhanced thermoelectric properties of doped CrSi2

Using first principles density functional theory, we report for CrSi2, a linear relationship between thermodynamic charge state transition levels of defects and maxima of thermopower T-m, thus proposing a unique way of tuning thermoelectric properties. We show for doped CrSi2 that the peak of thermopower occurs at the temperature which corresponds to the position of the defect transition level. Therefore, by modifying the defect transition level, a thermoelectric material with a given operational temperature can be designed.

Improving the Rigor and Consistency of the Thermodynamic Theory for Clathrate Hydrates through Incorporation of Movement of Water Molecules of Hydrate Lattice

Current applications of statistical thermodynamic theories for clathrate hydrates do not incorporate the translational and rotational movement of water molecules of the hydrate lattice,in a rigorous manner. Previous studies have shown that the movement of water molecules has a significant effect on the properties of clathrate hydrates. In this Article, a method is presented to incorporate the effect of water movement with as much rigor as possible. This method is then used to calculate the Langmuir constant of the guest species in a clathrate hydrate. Unlike previous studies on modeling of clathrate hydrate thermodynamics, the method presented in this paper does not regress either the intermolecular potentials or the properties of the empty hydrate from clathrate phase equilibria data. Also the properties of empty hydrate used in the theory do not depend on the nature and composition of the guest molecules. The predicted phase equilibria from the resulting theory are shown to be highly accurate and thermodynamically consistent by comparing them with the phase equilibria computed directly from molecular simulations.

Optical properties of excited silicon and germanium at low temperatures

Part I of this thesis deals with 3 topics concerning the luminescence from bound multi-exciton complexes in Si. Part II presents a model for the decay of electron-hole droplets in pure and doped Ge.

Part I.

We present high resolution photoluminescence data for Si doped With Al, Ga, and In. We observe emission lines due to recombination of electron-hole pairs in bound excitons and satellite lines which have been interpreted in terms of complexes of several excitons bound to an impurity. The bound exciton luminescence in Si:Ga and Si:Al consists of three emission lines due to transitions from the ground state and two low lying excited states. In Si:Ga, we observe a second triplet of emission lines which precisely mirror the triplet due to the bound exciton. This second triplet is interpreted as due to decay of a two exciton complex into the bound exciton. The observation of the second complete triplet in Si:Ga conclusively demonstrates that more than one exciton will bind to an impurity. Similar results are found for Si:Al. The energy of the lines show that the second exciton is less tightly bound than the first in Si:Ga. Other lines are observed at lower energies. The assumption of ground state to ground-state transitions for the lower energy lines is shown to produce a complicated dependence of binding energy of the last exciton on the number of excitons in a complex. No line attributable to the decay of a two exciton complex is observed in Si:In.

We present measurements of the bound exciton lifetimes for the four common acceptors in Si and for the first two bound multi-exciton complexes in Si:Ga and Si:Al. These results are shown to be in agreement with a calculation by Osbourn and Smith of Auger transition rates for acceptor bound excitons in Si. Kinetics determine the relative populations of complexes of various sizes and work functions, at temperatures which do not allow them to thermalize with respect to one another. It is shown that kinetic limitations may make it impossible to form two-exciton complexes in Si:In from a gas of free excitons.

We present direct thermodynamic measurements of the work functions of bound multi-exciton complexes in Al, B, P and Li doped Si. We find that in general the work functions are smaller than previously believed. These data remove one obstacle to the bound multi-exciton complex picture which has been the need to explain the very large apparent work functions for the larger complexes obtained by assuming that some of the observed lines are ground-state to ground-state transitions. None of the measured work functions exceed that of the electron-hole liquid.

Part II.

A new model for the decay of electron-hole-droplets in Ge is presented. The model is based on the existence of a cloud of droplets within the crystal and incorporates exciton flow among the drops in the cloud and the diffusion of excitons away from the cloud. It is able to fit the experimental luminescence decays for pure Ge at different temperatures and pump powers while retaining physically reasonable parameters for the drops. It predicts the shrinkage of the cloud at higher temperatures which has been verified by spatially and temporally resolved infrared absorption experiments. The model also accounts for the nearly exponential decay of electron-hole-droplets in lightly doped Ge at higher temperatures.

Effects of self energy of the ions on the double layer structure and properties at the dielectric interface

Although numerous theoretical efforts have been put forth, a systematic, unified and predictive theoretical framework that is able to capture all the essential physics of the interfacial behaviors of ions, such as the Hofmeister series effect, Jones-Ray effect and the salt effect on the bubble coalescence remain an outstanding challenge. The most common approach to treating electrostatic interactions in the presence of salt ions is the Poisson-Boltzmann (PB) theory. However, there are many systems for which the PB theory fails to offer even a qualitative explanation of the behavior, especially for ions distributed in the vicinity of an interface with dielectric contrast between the two media (like the water-vapor/oil interface). A key factor missing in the PB theory is the self energy of the ion.

In this thesis, we develop a self-consistent theory that treats the electrostatic self energy (including both the short-range Born solvation energy and the long-range image charge interactions), the nonelectrostatic contribution of the self energy, the ion-ion correlation and the screening effect systematically in a single framework. By assuming a finite charge spread of the ion instead of using the point-charge model, the self energy obtained by our theory is free of the divergence problems and gives a continuous self energy across the interface. This continuous feature allows ions on the water side and the vapor/oil side of the interface to be treated in a unified framework. The theory involves a minimum set of parameters of the ion, such as the valency, radius, polarizability of the ions, and the dielectric constants of the medium, that are both intrinsic and readily available. The general theory is first applied to study the thermodynamic property of the bulk electrolyte solution, which shows good agreement with the experiment result for predicting the activity coefficient and osmotic coefficient.

Next, we address the effect of local Born solvation energy on the bulk thermodynamics and interfacial properties of electrolyte solution mixtures. We show that difference in the solvation energy between the cations and anions naturally gives rise to local charge separation near the interface, and a finite Galvani potential between two coexisting solutions. The miscibility of the mixture can either increases or decreases depending on the competition between the solvation energy and translation entropy of the ions. The interfacial tension shows a non-monotonic dependence on the salt concentration: it increases linearly with the salt concentration at higher concentrations, and decreases approximately as the square root of the salt concentration for dilute solutions, which is in agreement with the Jones-Ray effect observed in experiment.

Next, we investigate the image effects on the double layer structure and interfacial properties near a single charged plate. We show that the image charge repulsion creates a depletion boundary layer that cannot be captured by a regular perturbation approach. The correct weak-coupling theory must include the self-energy of the ion due to the image charge interaction. The image force qualitatively alters the double layer structure and properties, and gives rise to many non-PB effects, such as nonmonotonic dependence of the surface energy on concentration and charge inversion. The image charge effect is then studied for electrolyte solutions between two plates. For two neutral plates, we show that depletion of the salt ions by the image charge repulsion results in short-range attractive and long-range repulsive forces. If cations and anions are of different valency, the asymmetric depletion leads to the formation of an induced electrical double layer. For two charged plates, the competition between the surface charge and the image charge effect can give rise to like- charge attraction.

Then, we study the inhomogeneous screening effect near the dielectric interface due to the anisotropic and nonuniform ion distribution. We show that the double layer structure and interfacial properties is drastically affected by the inhomogeneous screening if the bulk Debye screening length is comparable or smaller than the Bjerrum length. The width of the depletion layer is characterized by the Bjerrum length, independent of the salt concentration. We predict that the negative adsorption of ions at the interface increases linearly with the salt concentration, which cannot be captured by either the bulk screening approximation or the WKB approximation. For asymmetric salt, the inhomogeneous screening enhances the charge separation in the induced double layer and significantly increases the value of the surface potential.

Finally, to account for the ion specificity, we study the self energy of a single ion across the dielectric interface. The ion is considered to be polarizable: its charge distribution can be self-adjusted to the local dielectric environment to minimize the self energy. Using intrinsic parameters of the ions, such as the valency, radius, and polarizability, we predict the specific ion effect on the interfacial affinity of halogen anions at the water/air interface, and the strong adsorption of hydrophobic ions at the water/oil interface, in agreement with experiments and atomistic simulations.

The theory developed in this work represents the most systematic theoretical technique for weak-coupling electrolytes. We expect the theory to be more useful for studying a wide range of structural and dynamic properties in physicochemical, colloidal, soft-matter and biophysical systems.

A molecular dynamics study of the microscopic properties of simple dense fluids

The microscopic properties of a two-dimensional model dense fluid of Lennard-Jones disks have been studied using the so-called "molecular dynamics" method. Analyses of the computer-generated simulation data in terms of "conventional" thermodynamic and distribution functions verify the physical validity of the model and the simulation technique.

The radial distribution functions g(r) computed from the simulation data exhibit several subsidiary features rather similar to those appearing in some of the g(r) functions obtained by X-ray and thermal neutron diffraction measurements on real simple liquids. In the case of the model fluid, these "anomalous" features are thought to reflect the existence of two or more alternative configurations for local ordering.

Graphical display techniques have been used extensively to provide some intuitive insight into the various microscopic phenomena occurring in the model. For example, "snapshots" of the instantaneous system configurations for different times show that the "excess" area allotted to the fluid is collected into relatively large, irregular, and surprisingly persistent "holes". Plots of the particle trajectories over intervals of 2.0 to 6.0 x 10^-12 sec indicate that the mechanism for diffusion in the dense model fluid is "cooperative" in nature, and that extensive diffusive migration is generally restricted to groups of particles in the vicinity of a hole.

A quantitative analysis of diffusion in the model fluid shows that the cooperative mechanism is not inconsistent with the statistical predictions of existing theories of singlet, or self-diffusion in liquids. The relative diffusion of proximate particles is, however, found to be retarded by short-range dynamic correlations associated with the cooperative mechanism--a result of some importance from the standpoint of bimolecular reaction kinetics in solution.

A new, semi-empirical treatment for relative diffusion in liquids is developed, and is shown to reproduce the relative diffusion phenomena observed in the model fluid quite accurately. When incorporated into the standard Smoluchowski theory of diffusion-controlled reaction kinetics, the more exact treatment of relative diffusion is found to lower the predicted rate of reaction appreciably.

Finally, an entirely new approach to an understanding of the liquid state is suggested. Our experience in dealing with the simulation data--and especially, graphical displays of the simulation data--has led us to conclude that many of the more frustrating scientific problems involving the liquid state would be simplified considerably, were it possible to describe the microscopic structures characteristic of liquids in a concise and precise manner. To this end, we propose that the development of a formal language of partially-ordered structures be investigated.

I. The structure and properties of closed circular duplex DNA. II. Chemically interacting systems at equilibrium in a buoyant density gradient

I. The binding of the intercalating dye ethidium bromide to closed circular SV 40 DNA causes an unwinding of the duplex structure and a simultaneous and quantitatively equivalent unwinding of the superhelices. The buoyant densities and sedimentation velocities of both intact (I) and singly nicked (II) SV 40 DNAs were measured as a function of free dye concentration. The buoyant density data were used to determine the binding isotherms over a dye concentration range extending from 0 to 600 µg/m1 in 5.8 M CsCl. At high dye concentrations all of the binding sites in II, but not in I, are saturated. At free dye concentrations less than 5.4 µg/ml, I has a greater affinity for dye than II. At a critical amount of dye bound I and II have equal affinities, and at higher dye concentration I has a lower affinity than II. The number of superhelical turns, τ, present in I is calculated at each dye concentration using Fuller and Waring's (1964) estimate of the angle of duplex unwinding per intercalation. The results reveal that SV 40 DNA I contains about -13 superhelical turns in concentrated salt solutions.

The free energy of superhelix formation is calculated as a function of τ from a consideration of the effect of the superhelical turns upon the binding isotherm of ethidium bromide to SV 40 DNA I. The value of the free energy is about 100 kcal/mole DNA in the native molecule. The free energy estimates are used to calculate the pitch and radius of the superhelix as a function of the number of superhelical turns. The pitch and radius of the native I superhelix are 430 Å and 135 Å, respectively.

A buoyant density method for the isolation and detection of closed circular DNA is described. The method is based upon the reduced binding of the intercalating dye, ethidium bromide, by closed circular DNA. In an application of this method it is found that HeLa cells contain in addition to closed circular mitochondrial DNA of mean length 4.81 microns, a heterogeneous group of smaller DNA molecules which vary in size from 0.2 to 3.5 microns and a paucidisperse group of multiples of the mitochondrial length.

II. The general theory is presented for the sedimentation equilibrium of a macromolecule in a concentrated binary solvent in the presence of an additional reacting small molecule. Equations are derived for the calculation of the buoyant density of the complex and for the determination of the binding isotherm of the reagent to the macrospecies. The standard buoyant density, a thermodynamic function, is defined and the density gradients which characterize the four component system are derived. The theory is applied to the specific cases of the binding of ethidium bromide to SV 40 DNA and of the binding of mercury and silver to DNA.