The electrical double layer for a fully asymmetric electrolyte around a spherical colloid: an integral equation study

The hypernetted chain/mean spherical approximation (HNC/MSA) integral equation is obtained and solved numerically for a totally asymmetric primitive model electrolyte around a spherical macroparticle. The ensuing radial distribution functions show a very good agreement when compared to our Monte Carlo and molecular dynamics simulations for spherical geometry and with respect to previous anisotropic reference HNC calculations in the planar limit. We report an analysis of the potential vs charge relationship, radial distribution functions, mean electrostatic potential and cumulative reduced charge for representative cases of 1:1 and 2:2 salts with a size asymmetry ratio of 2. Our results are collated with those of the Modified Gouy-Chapman (MGC) and unequal radius Modified Gouy-Chapman (URMGC) theories and with those of HNC/MSA in the restricted primitive model (RPM) to assess the importance of size asymmetry effects. One of the most striking characteristics found is that,\textit{contrary to the general belief}, away from the point of zero charge the properties of an asymmetric electrical double layer (EDL) are not those corresponding to a symmetric electrolyte with the size and charge of the counterion, i.e. \textit{counterions do not always dominate}. This behavior suggests the existence of a new phenomenology in the EDL that genuinely belongs to a more realistic size-asymmetric model where steric correlations are taken into account consistently. Such novel features can not be described by traditional mean field theories like MGC, URMGC or even by enhanced formalisms, like HNC/MSA, if they are based on the RPM.Comment: 29 pages, 13 figure

Ion pairing in model electrolytes: A study via three particle correlation functions

A novel integral equations approach is applied for studying ion pairing in the restricted primitive model (RPM) electrolyte, i. e., the three point extension (TPE) to the Ornstein-Zernike integral equations. In the TPE approach, the three-particle correlation functions $g^{[3]}({\bf r}_{1},{\bf r}_{2},{\bf r}_{3})$ are obtained. The TPE results are compared to molecular dynamics (MD) simulations and other theories. Good agreement between TPE and MD is observed for a wide range of parameters, particularly where standard integral equations theories fail, i. e., low salt concentration and high ionic valence. Our results support the formation of ion pairs and aligned ion complexes.Comment: 43 pages (including 18 EPS figs) - RevTeX 4 - J. Chem. Phys. (in press

Electrophoresis of a rod macroion under polyelectrolyte salt: Is mobility reversed for DNA?

By molecular dynamics simulation, we study the charge inversion phenomenon of a rod macroion in the presence of polyelectrolyte counterions. We simulate electrophoresis of the macroion under an applied electric field. When both counterions and coions are polyelectrolytes, charge inversion occurs if the line charge density of the counterions is larger than that of the coions. For the macroion of surface charge density equal to that of the DNA, the reversed mobility is realized either with adsorption of the multivalent counterion polyelectrolyte or the combination of electrostatics and other mechanisms including the short-range attraction potential or the mechanical twining of polyelectrolyte around the rod axis.Comment: 8 pages, 5 figures, Applied Statistical Physics of Molecular Engineering (Mexico, 2003). Journal of Physics: Condensed Matters, in press (2004). Journal of Physics: Condensed Matters, in press (2004

Charge Fluctuation Forces Between Stiff Polyelectrolytes in Salt Solution: Pairwise Summability Re-examined

We formulate low-frequency charge-fluctuation forces between charged cylinders - parallel or skewed - in salt solution: forces from dipolar van der Waals fluctuations and those from the correlated monopolar fluctuations of mobile ions. At high salt concentrations forces are exponentially screened. In low-salt solutions dipolar energies go as $R^{-5}$ or $R^{-4}$; monopolar energies vary as $R^{-1}$ or $\ln{R}$, where $R$ is the minimal separation between cylinders. However, pairwise summability of rod-rod forces is easily violated in low-salt conditions. Perhaps the most important result is not the derivation of pair potentials but rather the demonstration that some of these expressions may not be used for the very problems that originally motivated their derivation.Comment: 8 pages and 1 fig in ps forma

Characterizing, modelling and understanding the climate variability of the deep water formation in the North-Western Mediterranean Sea

Observing, modelling and understanding the climate-scale variability of the deep water formation (DWF) in the North-Western Mediterranean Sea remains today very challenging. In this study, we first characterize the interannual variability of this phenomenon by a thorough reanalysis of observations in order to establish reference time series. These quantitative indicators include 31 observed years for the yearly maximum mixed layer depth over the period 1980–2013 and a detailed multi-indicator description of the period 2007–2013. Then a 1980–2013 hindcast simulation is performed with a fully-coupled regional climate system model including the high-resolution representation of the regional atmosphere, ocean, land-surface and rivers. The simulation reproduces quantitatively well the mean behaviour and the large interannual variability of the DWF phenomenon. The model shows convection deeper than 1000 m in 2/3 of the modelled winters, a mean DWF rate equal to 0.35 Sv with maximum values of 1.7 (resp. 1.6) Sv in 2013 (resp. 2005). Using the model results, the winter-integrated buoyancy loss over the Gulf of Lions is identified as the primary driving factor of the DWF interannual variability and explains, alone, around 50 % of its variance. It is itself explained by the occurrence of few stormy days during winter. At daily scale, the Atlantic ridge weather regime is identified as favourable to strong buoyancy losses and therefore DWF, whereas the positive phase of the North Atlantic oscillation is unfavourable. The driving role of the vertical stratification in autumn, a measure of the water column inhibition to mixing, has also been analyzed. Combining both driving factors allows to explain more than 70 % of the interannual variance of the phenomenon and in particular the occurrence of the five strongest convective years of the model (1981, 1999, 2005, 2009, 2013). The model simulates qualitatively well the trends in the deep waters (warming, saltening, increase in the dense water volume, increase in the bottom water density) despite an underestimation of the salinity and density trends. These deep trends come from a heat and salt accumulation during the 1980s and the 1990s in the surface and intermediate layers of the Gulf of Lions before being transferred stepwise towards the deep layers when very convective years occur in 1999 and later. The salinity increase in the near Atlantic Ocean surface layers seems to be the external forcing that finally leads to these deep trends. In the future, our results may allow to better understand the behaviour of the DWF phenomenon in Mediterranean Sea simulations in hindcast, forecast, reanalysis or future climate change scenario modes. The robustness of the obtained results must be however confirmed in multi-model studies

Spectral analysis and zeta determinant on the deformed spheres

We consider a class of singular Riemannian manifolds, the deformed spheres $S^N_k$, defined as the classical spheres with a one parameter family $g[k]$ of singular Riemannian structures, that reduces for $k=1$ to the classical metric. After giving explicit formulas for the eigenvalues and eigenfunctions of the metric Laplacian $\Delta_{S^N_k}$, we study the associated zeta functions $\zeta(s,\Delta_{S^N_k})$. We introduce a general method to deal with some classes of simple and double abstract zeta functions, generalizing the ones appearing in $\zeta(s,\Delta_{S^N_k})$. An application of this method allows to obtain the main zeta invariants for these zeta functions in all dimensions, and in particular $\zeta(0,\Delta_{S^N_k})$ and $\zeta'(0,\Delta_{S^N_k})$. We give explicit formulas for the zeta regularized determinant in the low dimensional cases, $N=2,3$, thus generalizing a result of Dowker \cite{Dow1}, and we compute the first coefficients in the expansion of these determinants in powers of the deformation parameter $k$.Comment: 1 figur

Attraction between DNA molecules mediated by multivalent ions

The effective force between two parallel DNA molecules is calculated as a function of their mutual separation for different valencies of counter- and salt ions and different salt concentrations. Computer simulations of the primitive model are used and the shape of the DNA molecules is accurately modelled using different geometrical shapes. We find that multivalent ions induce a significant attraction between the DNA molecules whose strength can be tuned by the averaged valency of the ions. The physical origin of the attraction is traced back either to electrostatics or to entropic contributions. For multivalent counter- and monovalent salt ions, we find a salt-induced stabilization effect: the force is first attractive but gets repulsive for increasing salt concentration. Furthermore, we show that the multivalent-ion-induced attraction does not necessarily correlate with DNA overcharging.Comment: 51 pages and 13 figure

Population inversion of a NAHS mixture adsorbed into a cylindrical pore

A cylindrical nanopore immersed in a non-additive hard sphere binary fluid is studied by means of integral equation theories and Monte Carlo simulations. It is found that at low and intermediate values of the bulk total number density the more concentrated bulk species is preferentially absorbed by the pore, as expected. However, further increments of the bulk number density lead to an abrupt population inversion in the confined fluid and an entropy driven prewetting transition at the outside wall of the pore. These phenomena are a function of the pore size, the non-additivity parameter, the bulk number density, and particles relative number fraction. We discuss our results in relation to the phase separation in the bulk.Comment: 7 pages, 8 Figure

Origin of Variability in Northern Hemisphere Winter Blocking on Interannual to Decadal Time Scales

Variability of mid-latitude blocking in the boreal winter northern hemisphere is investigated for the period 1960/61 to 2001/02 by means of relaxation experiments with the model of the European Center for Medium-Range Weather Forecasts. It is shown that there is pronounced interannual and decadal variability in blocking, especially over the Eurasian continent, consistent with previous studies. The relaxation experiments show that realistic variability in the tropics can account for a significant part of observed interannual blocking variability, but also that about half of the observed variability can only be explained by extratropical tropospheric variability. On the quasi-decadal time scale, extratropical sea surface temperature and sea-ice, in addition to tropical variability, play a more important role. The stratosphere, which has been shown to influence interannual variability of the North Atlantic Oscillation in previous studies, has no significant influence on blocking according to our analysis

Mechanisms of decadal variability in the Labrador Sea and the wider North Atlantic in a high-resolution climate model

A necessary step before assessing the performance of decadal predictions is the evaluation of the processes that bring memory to the climate system, both in climate models and observations. These mechanisms are particularly relevant in the North Atlantic, where the ocean circulation, related to both the Subpolar Gyre and the Meridional Overturning Circulation (AMOC), is thought to be important for driving significant heat content anomalies. Recently, a rapid decline in observed densities in the deep Labrador Sea has pointed to an ongoing slowdown of the AMOC strength taking place since the mid 90s, a decline also hinted by in-situ observations from the RAPID array. This study explores the use of Labrador Sea densities as a precursor of the ocean circulation changes, by analysing a 300-year long simulation with the state-of-the-art coupled model HadGEM3-GC2. The major drivers of Labrador Sea density variability are investigated, and are characterised by three major contributions. First, the integrated effect of local surface heat fluxes, mainly driven by year-to-year changes in the North Atlantic Oscillation, which accounts for 62% of the total variance. Additionally, two multidecadal-to-centennial contributions from the Greenland-Scotland Ridge outflows are quantified; the first associated with freshwater exports via the East Greenland Current, and the second with density changes in the Denmark Strait Overflow. Finally, evidence is shown that decadal trends in Labrador Sea densities are followed by important atmospheric impacts. In particular, a negative winter NAO response appears to follow the positive Labrador Sea density trends, and provides a phase reversal mechanism