Methodology for the calculation of the potential of mean force for a cation-pi complex in water.

We report potential of mean force (PMF) calculations on the interaction between the p-sulfonatocalix[4]arene and a monovalent cation (Cs(+)). It has been recently shown from microcalorimetry and (133)Cs NMR experiments that the association with Cs(+) is governed by favourable cation-pi interactions and is characterized by the insertion of the cation into the cavity of the macrocycle. We show that the PMF calculation based upon a classical model is not able to reproduce both the thermodynamic properties of association and the insertion of the cation. In order to take into account the different contributions of the cation-pi interactions, we develop a new methodology consisting of changing the standard PMF by an additional contribution resulting from quantum calculations. The calculated thermodynamic properties of association are thus in line with the microcalorimetry and (133)Cs NMR experiments and the structure of the complex at the Gibbs free-energy minimum shows the insertion of the cation into the cavity of the calixarene.     

Structure and dynamics of hydrated NH(4) (+): an ab initio QM/MM molecular dynamics simulation

A combined ab initio quantum mechanical/molecular mechanical (QM/MM) molecular dynamics simulation has been performed to investigate solvation structure and dynamics of NH(4) (+) in water. The most interesting region, the sphere includes an ammonium ion and its first hydration shell, was treated at the Hartree-Fock level using DZV basis set, while the rest of the system was described by classical pair potentials. On the basis of detailed QM/MM simulation results, the solvation structure of NH(4) (+) is rather flexible, in which many water molecules are cooperatively involved in the solvation shell of the ion. Of particular interest, the QM/MM results show fast translation and rotation of NH(4) (+) in water. This phenomenon has resulted from multiple coordination, which drives the NH(4) (+) to translate and rotate quite freely within its surrounding water molecules. In addition, a "structure-breaking" behavior of the NH(4) (+) is well reflected by the detailed analysis on the water exchange process and the mean residence times of water molecules surrounding the ion.     

The self-consistent mean spherical approximation for the one-component plasma

The consequences of choosing the adjustable hard-core diameter in the mean spherical approximation for the one-component plasma so as to achieve thermodynamic consistency between the energy and compressibility equations are investigated. Such a choice is found to be possible only for >8.5 and, although the resulting correlation functions are discontinuous, the height of the main peak in the static structure factor is remarkably accurate. Two especially noteworthy aspects of the thermodynamic results are that the compressibility equation is much more accurate than in any previous approximation free of input from computer simulations and that the nonstatic part of the internal energy has a ^1/4 dependence in the strong coupling limit in agreement with Monte Carlo data.     

Ion Association in Lanthanide Chloride Solutions

A better understanding of the solution chemistry of the lanthanide (Ln) salts in water would have wide ranging implications in materials processing, waste management, element tracing, medicine and many more fields. This is particularly true for minerals processing, given governmental concerns about lanthanide security of supply and the drive to identify environmentally sustainable processing routes. Despite much effort, even in simple systems, the mechanisms and thermodynamics of Ln^III association with small anions remain unclear. In the present study, molecular dynamics (MD), using a newly developed force field, provide new insights into LnCl₃(aq) solutions. The force field accurately reproduces the structure and dynamics of Nd³⁺, Gd³⁺ and Er³⁺ in water when compared to calculations using density functional theory (DFT). Adaptive-bias MD simulations show that the mechanisms for ion pairing change from dissociative to associative exchange depending upon cation size. Thermodynamics of association reveal that whereas ion pairing is favourable, the equilibrium distribution of species at low concentration is dominated by weakly bound solvent-shared and solvent-separated ion pairs, rather than contact ion pairs, reconciling a number of contrasting observations of Ln^III–Cl association in the literature. In addition, we show that the thermodynamic stabilities of a range of inner sphere and outer sphere coordination complexes are comparable and that the kinetics of anion binding to cations may control solution speciation distributions beyond ion pairs. The techniques adopted in this work provide a framework with which to investigate more complex solution chemistries of cations in water.     

通过脱碳控制全球平均温度

Establishing a reliable method to predict the global mean temperature (T_e) is of great importance because CO₂ reduction activities require political and global cooperation and significant financial resources. The current climate models all seem to predict that the earth's temperature will continue to increase, mainly based on the assumption that CO₂ emissions cannot be lowered significantly in the foreseeable future. Given the earth's multifactor climate system, attributing atmospheric CO₂ as the only cause for the observed temperature anomaly is most likely an oversimplification; the presence of water (H₂O) in the atmosphere should at least be considered. As such, T_e is determined by atmospheric water content controlled by solar activity, along with anthropogenic CO₂ activities. It is possible that the anthropogenic CO₂ activities can be reduced in the future. Based on temperature measurements and thermodynamic data, a new model for predicting T_e has been developed. Using this model, past, current, and future CO₂ and H₂O data can be analyzed and the associated T_e calculated. This new, esoteric approach is more accurate than various other models, but has not been reported in the open literature. According to this model, by 2050, T_e may increase to 15.5 ℃ under "business-as-usual" emissions. By applying a reasonable green technology activity scenario, T_e may be reduced to approximately 14.2 ℃. To achieve CO₂ reductions, the scenario described herein predicts a CO₂ reduction potential of 513 gigatons in 30 years. This proposed scenario includes various CO₂ reduction activities, carbon capturing technology, mineralization, and bio-char production; the most important CO₂ reductions by 2050 are expected to be achieved mainly in the electricity, agriculture, and transportation sectors. Other more aggressive and plausible drawdown scenarios have been analyzed as well, yielding CO₂ reduction potentials of 1051 and 1747 gigatons, respectively, in 30 years, but they may reduce global food production. It is emphasized that the causes and predictions of the global warming trend should be regarded as open scientific questions because several details concerning the physical processes associated with global warming remain uncertain. For example, the role of solar activities coupled with Milankovitch cycles are not yet fully understood. In addition, other factors, such as ocean CO₂ uptake and volcanic activity, may not be negligible.     

Hypersurfaces of prescribed mean curvature enclosing a given body

Given a smooth domain Ω in ℝ ^m+1 with compact closure and a smooth integrable functionh: ℝ ^m+1→ℝ satisfyingh(x)≥H _∂Ω (x) on ∂Ω whereH _∂ω denotes the mean curvature of ∂Ω calculated w.r.t. the interior unit normal we show that there is a setA⊂ℝ ^m+1 with the properties andH _∂A=h on ∂A.     

Optimal convergence properties of kernel density estimators without differentiability conditions

Let X ₁, X ₂, ..., X _n be independent observations from an (unknown) absolutely continuous univariate distribution with density f and let % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqefm0B1jxALjhiov2D% aebbfv3ySLgzGueE0jxyaibaiiYdd9qrFfea0dXdf9vqai-hEir8Ve% ea0de9qq-hbrpepeea0db9q8as0-LqLs-Jirpepeea0-as0Fb9pgea% 0lrP0xe9Fve9Fve9qapdbaqaaeGacaGaaiaabeqaamaabaabcaGcba% GabmOzayaajaGaaiikaiaadIhacaGGPaGaeyypa0Jaaiikaiaad6ga% caWGObGaaiykamaaCaaaleqabaGaeyOeI0IaaGymaaaakmaaqadaba% Gaam4saiaacUfadaWcgaqaaiaacIcacaWG4bGaeyOeI0Iaamiwamaa% BaaaleaacaWGPbaabeaakiaacMcaaeaacaWGObGaaiyxaaaaaSqaai% aadMgacqGH9aqpcaaIXaaabaGaamOBaaqdcqGHris5aaaa!5356!\[\hat f(x) = (nh)^{ - 1} \sum\nolimits_{i = 1}^n {K[{{(x - X_i )} \mathord{\left/ {\vphantom {{(x - X_i )} {h]}}} \right. \kern-\nulldelimiterspace} {h]}}} \] be a kernel estimator of f(x) at the point x, \s-<x<, with h=h _n (h _n O and nh _n , as n) the bandwidth and K a kernel function of order r. Optimal rates of convergence to zero for the bias and mean square error of such estimators have been studied and established by several authors under varying conditions on K and f. These conditions, however, have invariably included the assumption of existence of the r-th order derivative for f at the point x. It is shown in this paper that these rates of convergence remain valid without any differentiability assumptions on f at x. Instead some simple regularity conditions are imposed on the density f at the point of interest. Our methods are based on certain results in the theory of semi-groups of linear operators and the notions and relations of calculus of finite differences.This research was supported in part by grants from the Natural Sciences and Engineering Research Council of Canada and the University of Alberta Central Research Fund.     

Proof of some conjectures on the mean-value of Titchmarsh series — III

With some applications in view, the following problem is solved in some special case which is not too special. LetF(s) =Σ _n ^∞ =1a_n λ _n ^−s be a generalized Dirichlet series with 1 =λ ₁ <λ ₂ < …,λ _n≤Dn, andλ _n+1 -λ _n ≥D ^{− 1} λ _n+1 ^{− a} where α>0 andD(≥ 1) are constants. Then subject to analytic continuation and some growth conditions, a lower bound is obtained for . These results will be applied in other papers to appear later.     

Partition function zeros for the one-dimensional ordered plasma in Dirichlet boundary conditions

We consider the grand canonical partition function for the ordered one-dimensional, two-component plasma at fugacity in an applied electric fieldE with Dirichlet boundary conditions. The system has a phase transition from a low-coupling phase with equally spaced particles to a high-coupling phase with particles clustered into dipolar pairs. An exact expression for the partition function is developed. In zero applied field the zeros in the plane occupy the imaginary axis from –i to –i_c and i_c to i for some _c. They also occupy the diamond shape of four straight lines from ±i_c to _c and from ±i_c to –_c. The fugacity acts like a temperature or coupling variable. The symmetry-breaking field is the applied electric fieldE. A finite-size scaling representation for the partition in scaled coupling and scaled electric field is developed. It has standard mean field form. When the scaled coupling is real, the zeros in the scaled field lie on the imaginary axis and pinch the real scaled field axis as the scaled coupling increases. The scaled partition function considered as a function of two complex variables, scaled coupling and scaled field, has zeros on a two-dimensional surface in a domain of four real variables. A numerical discussion of some of the properties of this surface is presented.     

Negative ion formation in compounds relevant to SF6 decomposition in electrical discharges

Dissociative and nondissociative electron attachment in the electron impact energy range 0–14 eV are reported for SOF₂ SOF₄, SO₂F₂, SF₄, SO₂, and SiF₄ compounds which can be formed by electrical discharges in SF₆. The electron energy dependences of the mass-identified negative ions were determined in a time-of-flight mass spectrometer. The ions studied include F^– and SOF ₂ ^–* from SOF₂; SOF ₃ ^– and F^– from SOF₄; SO₂F ₂ ^–* , SO₂F^–, F ₂ ^– , and F^– from SO₂F₂; SF ₄ ^–* and F^– from SF₄; O^–, SO^–, and S^– from SO₂; and SiF ₃ ^– and F^– from SiF₄. Thermochemical data have been determined from the threshold energies of some of the fragment negative ions. Lifetimes of the anions SOF ₂ ^–* , SO₂F ₂ ^–* , and SF ₄ ^–* are also reported.