Enhanced Ion Anisotropy by Nonconventional Coordination Geometry: Single-Chain Magnet Behavior for a [{(FeL)-L-II}(2){Nb-IV(CN)(8)}] Helical Chain Compound Designed with Heptacoordinate Fe-II

Nonconventional heptacoordination in combination with efficient magnetic exchange coupling is shown to yield a 1-D heteronuclear {(FeNbIV)-Nb-II} compound with remarkable magnetic features when compared to other Fe(II)-based single chain magnets (SCM). Cyano-bridged heterometallic {3d-4d} and {3d-5d} chains are formed upon assembling Fe(II) bearing a pentadentate macrocycle as the blocking ligand with octacyano metallates, [M(CN)(8)](4-) (M = Nb-IV, Mo-IV, W-IV.) X-ray diffraction (single-crystal and powder) measurements reveal that the [{(H2O)Fe(L-1)}{M(CN)(8)}{Fe(L-1)}](infinity) architectures consist of isomorphous 1-D polymeric structures based on the alternation of {Fe(L-1)}(2+) and {M(CN)(8)}(4-) units (L-1 stands for the pentadentate macrocycle). Analysis of the magnetic susceptibility behavior revealed cyano-bridged {Fe-Nb} exchange interaction to be antiferromagnetic with J = -20 cm(-1) deduced from fitting an Ising model taking into account the noncollinear spin arrangement. For this ferrimagnetic chain a slow relaxation of its magnetization is observed at low temperature revealing a SCM behavior with Delta/k(B) = 74 K and tau(0) = 4.6 x 10(-11) s. The M versus H behavior exhibits a hysteresis loop with a coercive field of 4 kOe at 1 K and reveals at 380 mK magnetic avalanche processes, i.e., abrupt reversals in magnetization as H is varied. The origin of these characteristics is attributed to the combination of efficient {Fe-Nb} exchange interaction and significant anisotropy of the {Fe(L-1)) unit. High field EPR and magnetization experiments have revealed for the parent compound [Fe(L-1)(H2O)(2)]Cl-2 a negative zero field splitting parameter of D approximate to -17 cm(-1). The crystal structure, magnetic behavior, and Mossbauer data for [Fe(L-1)(H2O)(2)]Cl-2 are also reported.

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.

Effect of coordinated ions on structure and flexibility of parallel G-quadruplexes: A molecular dynamics study

ingle tract guanine residues can associate to form stable parallel quadruplex structures in the presence of certain cations. Nanosecond scale molecular dynamics simulations have been performed on fully solvated fibre model of parallel d(G(7)) quadruplex structures with Na+ or K+ ions coordinated in the cavity formed by the O6 atoms of the guanine bases. The AMBER 4.1 force field and Particle Mesh Ewald technique for electrostatic interactions have been used in all simulations. There quadruplex structures are stable during the simulation, with the middle four base tetrads showing root mean square deviation values between 0.5 to 0.8 Angstrom from the initial structure as well the high resolution crystal structure. Even in the absence of any coordinated ion in the initial structure, the G-quadruplex structure remains intact throughout the simulation. During the 1.1 ns MD simulation, one Nai counter ion from the solvent as well as several water molecules enter the central cavity to occupy the empty coordination sites within the parallel quadruplex and help stabilize the structure. Hydrogen bonding pattern depends on the nature of the coordinated ion, with the G-tetrad undergoing local structural variation to accommodate cations of different sizes. in the absence of any coordinated ion. due to strong mutual repulsion, O6 atoms within G-tetrad are forced farther apart from each other, which leads to a considerably different hydrogen bonding scheme within the G-tetrads and very favourable interaction energy between the guanine bases constituting a G-tetrad. However, a coordinated ion between G-tetrads provides extra stacking energy for the G-tetrads and makes the quadruplex structure more rigid. Na+ ions, within the quadruplex cavity, are more mobile than coordinated K+ ions. A number of hydrogen bonded water molecules are observed within the grooves of all quadruplex structures.

Molecular theory of solvation and solvation dynamics of a classical ion in a dipolar liquid

A microscopic theory of equilibrium solvation and solvation dynamics of a classical, polar, solute molecule in dipolar solvent is presented. Density functional theory is used to explicitly calculate the polarization structure around a solvated ion. The calculated solvent polarization structure is different from the continuum model prediction in several respects. The value of the polarization at the surface of the ion is less than the continuum value. The solvent polarization also exhibits small oscillations in space near the ion. We show that, under certain approximations, our linear equilibrium theory reduces to the nonlocal electrostatic theory, with the dielectric function (c(k)) of the liquid now wave vector (k) dependent. It is further shown that the nonlocal electrostatic estimate of solvation energy, with a microscopic c(k), is close to the estimate of linearized equilibrium theories of polar liquids. The study of solvation dynamics is based on a generalized Smoluchowski equation with a mean-field force term to take into account the effects of intermolecular interactions. This study incorporates the local distortion of the solvent structure near the ion and also the effects of the translational modes of the solvent molecules.The latter contribution, if significant, can considerably accelerate the relaxation of solvent polarization and can even give rise to a long time decay that agrees with the continuum model prediction. The significance of these results is discussed.

Correlation of Oxygen Storage Capacity and Structural Distortion in Transition-Metal-, Noble-Metal-, and Rare-Earth-Ion-Substituted CeO2 from First Principles Calculation

Oxygen storage/release (OSC) capacity is an important feature common to all three-way catalysts to combat harmful exhaust emissions. To understand the mechanism of improved OSC for doped CeO2, we undertook the structural investigation by X-ray diffraction (XRD), X-ray photoelectron spectroscopy (XPS), H-2-TPR (temperature-programmed hydrogen reduction) and density functional theoretical (DFT) calculations of transition-metal-, noble-metal-, and rare-earth (RE)-ion-substituted ceria. In this report, we present the relationship between the OSC and structural changes induced by the dopant ion in CeO2. Transition metal and noble metal ion substitution in ceria greatly enhances the reducibility of Ce1-xMxO2-delta (M = Mn, Fe, Co, Ni, Cu, Pd, Pt, Ru), whereas rare-earth-ion-substituted Ce(1-x)A(x)O(2-delta) (A = La, Y) have very little effect in improving the OSC. Our simulated optimized structure shows deviation in cation oxygen bond length from ideal bond length of 2.34 angstrom (for CeO2). For example, our theoretical calculation for Ce28Mn4O62 structure shows that Mn-O bonds are in 4 + 2 coordination with average bond lengths of 2.0 and 3.06 angstrom respectively. Although the four short Mn-O bond lengths spans the bond distance region of Mn2O3, the other two Mn-O bonds are moved to longer distances. The dopant transition and noble metal ions also affects Ce coordination shell and results in the formation of longer Ce-O bonds as well. Thus longer cation oxygen bonds for both dopant and host ions results in enhanced synergistic reduction of the solid solution. With Pd ion substitution in Ce1-xMxO2-delta (M = Mn, Fe, Co, Ni, Cu) further enhancement in OSC is observed in H-2-TPR. This effect is reflected in our model calculations by the presence of still longer bonds compared to the model without Pd ion doping. The synergistic effect is therefore due to enhanced reducibility of both dopant and host ion induced due to structural distortion of fluorite lattice in presence of dopant ion. For RE ions (RE = Y, La), our calculations show very little deviation of bonds lengths from ideal fluorite structure. The absence of longer Y-O/La-O and Ce-O bonds make the structure much less susceptible to reduction.

Shell-model calculation of defect energies in alkali halides employing crystal-independent interionic potential parameters

Shell model calculation of defect energies in alkali halides have been carried out using the ion-dependent, crystal-independent potential parameters of Sangster and Atwood (1978). Results indicate that appreciable differences exist between barrier heights for migration of cations and anions. While barrier heights for cations are generally lower than for anions in alkali halides of NaCl structure, the opposite is true in alkali halides of CsCl structure.

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.

Ion-acoustic solutions and shock waves in multicomponent plasmas

The present investigation of ion-acoustic waves is based on the study of the nonlinearity of plasma waves in a dispersive medium. Here the authors study ion-acoustic solitary waves in a warm ion plasma with non-isothermal electrons and then the results for solitary waves in a plasma with isothermal electrons are obtained. Incorporating the previous results obtained from the solitary wave solutions, the authors generalize the effect of negative ions on ion-acoustic waves in plasmas consisting of either a warm or cold ion species. A reflection phenomenon of ions in these waves is also studied. These results can be generalized, but the discussion is limited to a particular model of the plasma.

Effect of coordinated ions on structure and flexiblity of parallel G-quandruplexes: a molecular dynamics study

Single tract guanine residues can associate to form stable parallel quadruplex structures in the presence of certain cations. Nanosecond scale molecular dynamics simulations have been performed on fully solvated fibre model of parallel d(G7) quadruplex structures with Na+ or K+ ions coordinated in the cavity formed by the 06 atoms of the guanine bases. The AMBER 4.1 force field and Particle Mesh Ewald technique for electrostatic interactions have been used in all simulations. These quadruplex structures are stable during the simulation, with the middle four base tetrads showing root mean square deviation values between 0.5 to 0.8 A from the initial structure as well the high resolution crystal structure. Even in the absence of any coordinated ion in the initial structure, the G-quadruplex structure remains intact throughout the simulation. During the 1.1 ns MD simulation, one Na+ counter ion from the solvent as well as several water molecules enter the central cavity to occupy the empty coordination sites within the parallel quadruplex and help stabilize the structure. Hydrogen bonding pattern depends on the nature of the coordinated ion, with the G-tetrad undergoing local structural variation to accommodate cations of different sizes. In the absence of any coordinated ion, due to strong mutual repulsion, 06 atoms within G-tetrad are forced farther apart from each other, which leads to a considerably different hydrogen bonding scheme within the G-tetrads and very favourable interaction energy between the guanine bases constituting a G-tetrad. However, a coordinated ion between G-tetrads provides extra stacking energy for the G-tetrads and makes the quadruplex structure more rigid. Na+ ions, within the quadruplex cavity, are more mobile than coordinated K+ ions. A number of hydrogen bonded water molecules are observed within the grooves of all quadruplex structures

A simple computational model for predicting .pi.-facial selectivity in reductions of sterically unbiased ketones. Relative importance of electrostatic and orbital interactions

Various factore controlling the preferred facial selectivity in the reductions of a number of sterically unbiased ketones have been evaluated using a semiempirical MO procedure. MNDO optimized geometries do not reveal any significant ground-state distortions which can be correlated with the observed face selectivities. Electrostatic effecta due to an approaching reagent were modeled by placing a test negative charge at a fixed distance from the carbonyl carbon on each of the two faces. A second series of calculations was carried out using the hydride ion as a test nucleophile. The latter calculations effectively include orbital interactions involving the u and u* orbitals of the newly formed bond in the reaction. The computed energy differences with the charge model are generally much larger compared to those with the hydride ion. However, both models lead to predictions which are qualitatively consistent with the experimentally determined facial preferences for most of the systems. Thus, electrostatic interactions between the nucleophile and the substrate seem to effectively determine the face selectivities in these molecules. However, there are a few exceptions in which orbital interactions are found to contribute significantly and occasionally reverse the preference dictated by electrostatic effecta. The remarkable succew of the hydride model calculations, in spite of retaining the unperturbed geometries of the substrates, points to the unimportance of torsional effeds and orbital distortions associated with the pyramidalized carbonyl unit in the transition state in most of the substrates considered. Additional experimental results are reported which provide useful calibration for the present computational approach.

The scratch hardness and friction of a soft rigid-plastic solid

The paper describes an experimental and analytical study of the normal and scratch hardnesses of a model soft rigid-plastic solid. The material known as ‘Plasticine’, a mixture of dry particles and a mineral oil, has been deformed with a range of rigid conical indentors with included angles of between 30° and 170°. The sliding velocity dependence of the computed scratch hardness and friction has been examined in the velocity range 0.19 mm/s to 7.3 m/s. Data are also described for the time dependence of the normal hardness and also the estimated rate dependence of the intrinsic flow stress. The latter values were estimated from data obtained during the upsetting of right cylinders. Three major conclusions are drawn from these data and the associated analysis. (1) A first-order account of the scratching force may be provided by adopting a model which sums the computed plastic deformation and interfacial sliding contributions to the total sliding work. This is tantamount to the adoption of the two-term non-interacting model of friction. (2) For this system during sliding, at high sliding velocities at least, the interface shear stress which defines the boundary condition is not directly related to the bulk shear stress. The interface rheological characteristics indicate an appreciable dependence on the imposed strain or strain rate. In particular, the relative contributions of the slip and stick boundary conditions appear to be a function of the imposed sliding velocity. (3) The computed normal and scratch hardness values are not simply interrelated primarily because of the evolving boundary conditions which appear to exist in the scratching experiments.

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.

Chain folding and A:T pairing in human telomeric DNA: a model-building and molecular dynamics study

The various types of chain folding and possible intraloop as well as interloop base pairing in human telomeric DNA containing d(TTAG(3)) repeats have been investigated by model-building, molecular mechanics, and molecular dynamics techniques. Model-building and molecular mechanics studies indicate that it is possible to build a variety of energetically favorable folded-back structures with the two TTA loops on same side and the 5' end thymines in the two loops forming TATA tetrads involving a number of different intraloop as well as interloop A:T pairing schemes. In these folded-back structures, although both intraloop and interloop Watson-Crick pairing is feasible, no structure is possible with interloop Hoogsteen pairing. MD studies of representative structures indicate that the guanine-tetraplex stem is very rigid and, while the loop regions are relatively much more flexible, most of the hydrogen bonds remain intact throughout the 350-ps in vacuo simulation. The various possible TTA loop structures, although they are energetically similar, have characteristic inter proton distances, which could give rise to unique cross-peaks in two-dimensional nuclear Overhauser effect spectroscopy (NOESY) experiments. These folded-back structures with A:T pairings in the loop region help in rationalizing the data from chemical probing and other biochemical studies on human telomeric DNA.

Non-Marcus energy gap dependence of electron transfer rate in contact ion pairs. Novel interplay between relaxation and reaction in solution

Electron transfer reactions between donor-acceptor pairs in solution and in organized media exhibit diverse behaviour. Recent experiments have indicated an interesting breakdown of the Marcus parabolic energy gap dependence in the normal regime for back electron transfer from contact ion pairs. A novel explanation of this breakdown has recently been proposed (M. Tachiya and S. Murata, J. Am. Chem. Sec., 116(1994) 2434) which attributes the breakdown to the interplay between the relaxation in the reactant well and the reaction. A particularly interesting aspect of the model is that it envisages the electron transfer in the normal regime to take place from a completely non-equilibrium condition. In this article a time dependent solution of the model is presented for the first time, after generalizing it to include a realistic initial population distribution. The decay of the contact ion pair population is completely non-exponential. This can be used to check the validity of the Tachiya-Murata model. The dynamics of electron transfer from the solvent separated ion pair, which seem to obey the Marcus relation, is exponential.

Influence of carbide particles on the dynamic recrystallization of as-cast 304 stainless steel: a study using secondary ion mass spectrometric (SIMS) analysis

The domain of dynamic recrystallization (DRX) in as-cast 304 stainless steel material occurs at higher temperatures (1250 degrees C) and lower strain rates (0.001 s(-1)) than in wrought 304 stainless steel (1100 degrees C and 0.01 s(-1)). The above result has been explained earlier on the basis of a simple theoretical DRX model involving the rate of nucleation versus rate of grain boundary migration. The present investigation is aimed at examining experimentally the influence of carbide particles on the DRX of ascast 304 using secondary ion mass spectrometric (SIMS) analysis. Isothermal compression tests at a constant true strain rate have been performed on wrought 304 and as-cast 304 materials in the temperature and strain rate ranges of 1000 to 1250 degrees C and 0.001 to 1 s(-1) respectively. The SIMS analysis carried out on the deformed samples revealed that the large carbides present in the as-cast 304 material strongly influence the DRX process. In as-cast 304 material, the presence of large carbide particles in the microstructure shifts the DRX domain to higher temperature and lower strain rate in comparison with wrought 304 material.