Synaptic transmission with the glia.

For decades, scientists thought that all of the missing secrets of brain function resided in neurons. However, a wave of new findings indicates that glial cells, formerly considered mere supporters and subordinate to neurons, participate actively in synaptic integration and processing of information in the brain.

Citrate- vs. acetate-based dialysate in bicarbonate haemodialysis: consequences on haemodynamics, coagulation, acid-base status, and electrolytes.

BACKGROUND: A concentrate for bicarbonate haemodialysis acidified with citrate instead of acetate has been marketed in recent years. The small amount of citrate used (one-fifth of the concentration adopted in regional anticoagulation) protects against intradialyser clotting while minimally affecting the calcium concentration. The aim of this study was to compare the impact of citrate- and acetate-based dialysates on systemic haemodynamics, coagulation, acid-base status, calcium balance and dialysis efficiency. METHODS: In 25 patients who underwent a total of 375 dialysis sessions, an acetate dialysate (A) was compared with a citrate dialysate with (C+) or without (C) calcium supplementation (0.25 mmol/L) in a randomised single-blind cross-over study. Systemic haemodynamics were evaluated using pulse-wave analysis. Coagulation, acid-base status, calcium balance and dialysis efficiency were assessed using standard biochemical markers. RESULTS: Patients receiving the citrate dialysate had significantly lower systolic blood pressure (BP) (-4.3 mmHg, p < 0.01) and peripheral resistances (PR) (-51 dyne.sec.cm-5, p < 0.001) while stroke volume was not increased. In hypertensive patients there was a substantial reduction in BP (-7.8 mmHg, p < 0.01). With the C+ dialysate the BP gap was less pronounced but the reduction in PR was even greater (-226 dyne.sec.cm-5, p < 0.001). Analyses of the fluctuations in PR and of subjective tolerance suggested improved haemodynamic stability with the citrate dialysate. Furthermore, an increase in pre-dialysis bicarbonate and a decrease in pre-dialysis BUN, post-dialysis phosphate and ionised calcium were noted. Systemic coagulation activation was not influenced by citrate. CONCLUSION: The positive impact on dialysis efficiency, acid-base status and haemodynamics, as well as the subjective tolerance, together indicate that citrate dialysate can significantly contribute to improving haemodialysis in selected patients.

Serum calcium levels are associated with novel cardiometabolic risk factors in the population-based CoLaus study.

BACKGROUND: Associations of serum calcium levels with the metabolic syndrome and other novel cardio-metabolic risk factors not classically included in the metabolic syndrome, such as those involved in oxidative stress, are largely unexplored. We analyzed the association of albumin-corrected serum calcium levels with conventional and non-conventional cardio-metabolic risk factors in a general adult population. METHODOLOGY/PRINCIPAL FINDINGS: The CoLaus study is a population-based study including Caucasians from Lausanne, Switzerland. The metabolic syndrome was defined using the Adult Treatment Panel III criteria. Non-conventional cardio-metabolic risk factors considered included: fat mass, leptin, LDL particle size, apolipoprotein B, fasting insulin, adiponectin, ultrasensitive CRP, serum uric acid, homocysteine, and gamma-glutamyltransferase. We used adjusted standardized multivariable regression to compare the association of each cardio-metabolic risk factor with albumin-corrected serum calcium. We assessed associations of albumin-corrected serum calcium with the cumulative number of non-conventional cardio-metabolic risk factors. We analyzed 4,231 subjects aged 35 to 75 years. Corrected serum calcium increased with both the number of the metabolic syndrome components and the number of non-conventional cardio-metabolic risk factors, independently of the metabolic syndrome and BMI. Among conventional and non-conventional cardio-metabolic risk factors, the strongest positive associations were found for factors related to oxidative stress (uric acid, homocysteine and gamma-glutamyltransferase). Adiponectin had the strongest negative association with corrected serum calcium. CONCLUSIONS/SIGNIFICANCE: Serum calcium was associated with the metabolic syndrome and with non-conventional cardio-metabolic risk factors independently of the metabolic syndrome. Associations with uric acid, homocysteine and gamma-glutamyltransferase were the strongest. These novel findings suggest that serum calcium levels may be associated with cardiovascular risk via oxidative stress.

PTH and 1.25 vitamin D response to a low-calcium diet is associated with bone mineral density in renal stone formers.

BACKGROUND: Renal calcium stones and hypercalciuria are associated with a reduced bone mineral density (BMD). Therefore, the effect of changes in calcium homeostasis is of interest for both stones and bones. We hypothesized that the response of calciuria, parathyroid hormone (PTH) and 1.25 vitamin D to changes in dietary calcium might be related to BMD. METHODS: A single-centre prospective interventional study of 94 hyper- and non-hypercalciuric calcium stone formers consecutively retrieved from our stone clinic. The patients were investigated on a free-choice diet, a low-calcium diet, while fasting and after an oral calcium load. Patient groups were defined according to lumbar BMD (z-score) obtained by dual X-ray absorptiometry (group 1: z-score <-0.5, n = 30; group 2: z-score -0.5-0.5, n = 36; group 3: z-score >0.5, n = 28). The effect of the dietary interventions on calciuria, 1.25 vitamin D and PTH in relation to BMD was measured. RESULTS: An inverse relationship between BMD and calciuria was observed on all four calcium intakes (P = 0.009). On a free-choice diet, 1.25 vitamin D and PTH levels were identical in the three patient groups. However, the relative responses of 1.25 vitamin D and PTH to the low-calcium diet were opposite in the three groups with the highest increase of 1.25 vitamin D in group 1 and the lowest in group 3, whereas PTH increase was most pronounced in group 3 and least in group 1. CONCLUSION: Calcium stone formers with a low lumbar BMD exhibit a blunted response of PTH release and an apparently overshooting production of 1.25 vitamin D following a low-calcium diet.

Calcium-sensing receptors modulate renin release in vivo and in vitro in the rat.

OBJECTIVES: Calcium-sensing receptors (CaSRs) have been localized in the juxtaglomerular apparatus where they may contribute to the regulation of renin release. In the present study, we investigated the in-vitro and in-vivo effects of the calcimimetic R-568 on renin release. METHODS: In vitro, the effect of calcimimetics on renin release was assessed by incubating freshly isolated rat juxtaglomerular cells with or without R-568 (1 and 10 mumol/l) in serum-free medium in the presence or absence of forskolin or CaCl2. In vivo, we measured the impact of R-568 (20 ng/min intravenously) on the acute changes in plasma renin activity (PRA) induced by either a 90 min infusion of the angiotensin-converting enzyme inhibitor captopril, or the beta-receptor agonist isoproterenol, or of a vehicle in or after a furosemide challenge in conscious Wistar rats. RESULTS: In vitro, R-568 dose-dependently blunted renin release, but also reduced the increase in renin due to forskolin (P < 0.01). Both isoproterenol and enalapril increased in vivo PRA to 3.1 +/- 0.3 and 3.7 +/- 0.5 ng Ang I/ml per h, respectively (P < 0.01), compared with vehicle (1.5 +/- 0.2 ng Ang I/ml per h). R-568 significantly reduced PRA to 2.1 +/- 0.1 ng/ml per h in isoproterenol-treated rats and to 1.6 +/- 0.2 ng/ml per h in enalapril-treated rats (P < 0.05). In low-salt treated animals, acute infusion of furosemide increased PRA from 8.7 +/- 3.2 to 18.6 +/- 2.3, whereas R-568 partially blunted this rise to 11.2 +/- 1.5 (P = 0.02). In vivo, R-568 significantly lowered serum calcium and PTH1-84, but the drug-induced changes in PRA were independent of the changes in calcium and parathyroid hormone. CONCLUSION: After the recent discovery of CaSRs in juxtaglomerular cells of mice, our results confirm the presence of such receptors in rats and demonstrate that these receptors modulate renin release both in vitro and in vivo. This suggests that CaSRs play a role as a regulatory pathway of renin release.

Influence of nutrition on urinary oxalate and calcium in preterm and term infants.

Few data for normal urinary oxalate (Ox) and calcium (Ca) excretion related both to gestational age and nutritional factors have been reported in preterm or term infants. We therefore determined the molar Ox and Ca to creatinine (Cr) ratios in spot urines from 64 preterm and 37 term infants aged 1-60 days, either fed formula or human milk (HM). Only vitamin D was supplemented; renal or metabolic diseases were excluded. Urinary Ox/Cr ratio was higher in preterm than in term infants, both when formula fed (1st month 253 vs. 180 mmol/mol and 2nd month 306 vs. 212 mmol/mol; P&lt;0.05) or HM fed (206 vs. 169 mmol/ mol and 283* vs. 232 mmol/mol; *P&lt;0.05). Ox/Cr was also higher in formula- than HM-fed preterm infants. The ratio increased during the first 2 months of life irrespective of nutrition. Urinary Ca/Cr ratio was comparable in all groups during the 1st month of life, except for a lower (P &lt; 0.05) value in term infants fed HM (0.10 mol/mol). It increased in all groups during the 2nd month of life, being highest in HM-fed preterm infants (1.86 mol/mol). In conclusion, urinary Ox and Ca excretion is influenced by both gestational age and nutrient intake in preterm and term infants.

High sensitivity of immature GABAergic neurons to blockers of voltage-gated calcium channels.

The involvement of voltage-gated calcium channels in the survival of immature CNS neurons was studied in aggregating brain cell cultures by examining cell type-specific effects of various channel blockers. Nifedipine (10 microM), a specific blocker of L-type calcium channels, caused a pronounced and irreversible decrease of glutamic acid decarboxylase activity, whereas the activity of choline acetyltransferase was significantly less affected. Flunarizine (1-10 microM, a relatively unspecific ion channel blocker) elicited similar effects, that were attenuated by NMDA. The glia-specific marker enzymes, glutamine synthetase and 2',3'-cyclic nucleotide 3'-phosphohydrolase, were affected only after treatment with high concentrations of nifedipine (50 microM) or NiCl2 (100 microM, shown to block T-type calcium channels). Nifedipine (50 microM), NiCl2 (100 microM), and flunarizine (5 microM) also caused a significant increase in the soluble nucleosome concentration, indicating increased apoptotic cell death. This effect was prevented by cycloheximide (1 microM). Furthermore, the combined treatment with calcicludine (10 nM, blocking L-type calcium channels) and funnel-web spider toxin-3.3 (100 nM, blocking T-type channels) also caused a significant increase in free nucleosomes as well as a decrease in glutamic acid decarboxylase activity. In contrast, cell viability was not affected by peptide blockers specific for N-, P-, and/or Q-type calcium channels. Highly differentiated cultures showed diminished susceptibility to nifedipine and flunarizine. The present data suggest that the survival of immature neurons, and particularly that of immature GABAergic neurons, requires the sustained entry of Ca2+ through voltage-gated calcium channels.

Chest pulse-wave velocity: a novel approach to assess arterial stiffness.

Pulse-wave velocity (PWV) is considered as the gold-standard method to assess arterial stiffness, an independent predictor of cardiovascular morbidity and mortality. Current available devices that measure PWV need to be operated by skilled medical staff, thus, reducing the potential use of PWV in the ambulatory setting. In this paper, we present a new technique allowing continuous, unsupervised measurements of pulse transit times (PTT) in central arteries by means of a chest sensor. This technique relies on measuring the propagation time of pressure pulses from their genesis in the left ventricle to their later arrival at the cutaneous vasculature on the sternum. Combined thoracic impedance cardiography and phonocardiography are used to detect the opening of the aortic valve, from which a pre-ejection period (PEP) value is estimated. Multichannel reflective photoplethysmography at the sternum is used to detect the distal pulse-arrival time (PAT). A PTT value is then calculated as PTT = PAT - PEP. After optimizing the parameters of the chest PTT calculation algorithm on a nine-subject cohort, a prospective validation study involving 31 normo- and hypertensive subjects was performed. 1/chest PTT correlated very well with the COMPLIOR carotid to femoral PWV (r = 0.88, p < 10 (-9)). Finally, an empirical method to map chest PTT values onto chest PWV values is explored.

A pseudo-spectral method for simulating of poro-elastic seismic wave propagation in complex borehole environments

We present a novel numerical approach for the comprehensive, flexible, and accurate simulation of poro-elastic wave propagation in cylindrical coordinates. An important application of this method is the modeling of complex seismic wave phenomena in fluid-filled boreholes, which represents a major, and as of yet largely unresolved, computational problem in exploration geophysics. In view of this, we consider a numerical mesh consisting of three concentric domains representing the borehole fluid in the center, the borehole casing and the surrounding porous formation. The spatial discretization is based on a Chebyshev expansion in the radial direction, Fourier expansions in the other directions, and a Runge-Kutta integration scheme for the time evolution. A domain decomposition method based on the method of characteristics is used to match the boundary conditions at the fluid/porous-solid and porous-solid/porous-solid interfaces. The viability and accuracy of the proposed method has been tested and verified in 2D polar coordinates through comparisons with analytical solutions as well as with the results obtained with a corresponding, previously published, and independently benchmarked solution for 2D Cartesian coordinates. The proposed numerical solution also satisfies the reciprocity theorem, which indicates that the inherent singularity associated with the origin of the polar coordinate system is handled adequately.

Time-lapse and probabilistic inversion strategies for plane-wave electromagnetic data.

L'utilisation efficace des systèmes géothermaux, la séquestration du CO2 pour limiter le changement climatique et la prévention de l'intrusion d'eau salée dans les aquifères costaux ne sont que quelques exemples qui démontrent notre besoin en technologies nouvelles pour suivre l'évolution des processus souterrains à partir de la surface. Un défi majeur est d'assurer la caractérisation et l'optimisation des performances de ces technologies à différentes échelles spatiales et temporelles. Les méthodes électromagnétiques (EM) d'ondes planes sont sensibles à la conductivité électrique du sous-sol et, par conséquent, à la conductivité électrique des fluides saturant la roche, à la présence de fractures connectées, à la température et aux matériaux géologiques. Ces méthodes sont régies par des équations valides sur de larges gammes de fréquences, permettant détudier de manières analogues des processus allant de quelques mètres sous la surface jusqu'à plusieurs kilomètres de profondeur. Néanmoins, ces méthodes sont soumises à une perte de résolution avec la profondeur à cause des propriétés diffusives du champ électromagnétique. Pour cette raison, l'estimation des modèles du sous-sol par ces méthodes doit prendre en compte des informations a priori afin de contraindre les modèles autant que possible et de permettre la quantification des incertitudes de ces modèles de façon appropriée. Dans la présente thèse, je développe des approches permettant la caractérisation statique et dynamique du sous-sol à l'aide d'ondes EM planes. Dans une première partie, je présente une approche déterministe permettant de réaliser des inversions répétées dans le temps (time-lapse) de données d'ondes EM planes en deux dimensions. Cette stratégie est basée sur l'incorporation dans l'algorithme d'informations a priori en fonction des changements du modèle de conductivité électrique attendus. Ceci est réalisé en intégrant une régularisation stochastique et des contraintes flexibles par rapport à la gamme des changements attendus en utilisant les multiplicateurs de Lagrange. J'utilise des normes différentes de la norme l2 pour contraindre la structure du modèle et obtenir des transitions abruptes entre les régions du model qui subissent des changements dans le temps et celles qui n'en subissent pas. Aussi, j'incorpore une stratégie afin d'éliminer les erreurs systématiques de données time-lapse. Ce travail a mis en évidence l'amélioration de la caractérisation des changements temporels par rapport aux approches classiques qui réalisent des inversions indépendantes à chaque pas de temps et comparent les modèles. Dans la seconde partie de cette thèse, j'adopte un formalisme bayésien et je teste la possibilité de quantifier les incertitudes sur les paramètres du modèle dans l'inversion d'ondes EM planes. Pour ce faire, je présente une stratégie d'inversion probabiliste basée sur des pixels à deux dimensions pour des inversions de données d'ondes EM planes et de tomographies de résistivité électrique (ERT) séparées et jointes. Je compare les incertitudes des paramètres du modèle en considérant différents types d'information a priori sur la structure du modèle et différentes fonctions de vraisemblance pour décrire les erreurs sur les données. Les résultats indiquent que la régularisation du modèle est nécessaire lorsqu'on a à faire à un large nombre de paramètres car cela permet d'accélérer la convergence des chaînes et d'obtenir des modèles plus réalistes. Cependent, ces contraintes mènent à des incertitudes d'estimations plus faibles, ce qui implique des distributions a posteriori qui ne contiennent pas le vrai modèledans les régions ou` la méthode présente une sensibilité limitée. Cette situation peut être améliorée en combinant des méthodes d'ondes EM planes avec d'autres méthodes complémentaires telles que l'ERT. De plus, je montre que le poids de régularisation des paramètres et l'écart-type des erreurs sur les données peuvent être retrouvés par une inversion probabiliste. Finalement, j'évalue la possibilité de caractériser une distribution tridimensionnelle d'un panache de traceur salin injecté dans le sous-sol en réalisant une inversion probabiliste time-lapse tridimensionnelle d'ondes EM planes. Etant donné que les inversions probabilistes sont très coûteuses en temps de calcul lorsque l'espace des paramètres présente une grande dimension, je propose une stratégie de réduction du modèle ou` les coefficients de décomposition des moments de Legendre du panache de traceur injecté ainsi que sa position sont estimés. Pour ce faire, un modèle de résistivité de base est nécessaire. Il peut être obtenu avant l'expérience time-lapse. Un test synthétique montre que la méthodologie marche bien quand le modèle de résistivité de base est caractérisé correctement. Cette méthodologie est aussi appliquée à un test de trac¸age par injection d'une solution saline et d'acides réalisé dans un système géothermal en Australie, puis comparée à une inversion time-lapse tridimensionnelle réalisée selon une approche déterministe. L'inversion probabiliste permet de mieux contraindre le panache du traceur salin gr^ace à la grande quantité d'informations a priori incluse dans l'algorithme. Néanmoins, les changements de conductivités nécessaires pour expliquer les changements observés dans les données sont plus grands que ce qu'expliquent notre connaissance actuelle des phénomenès physiques. Ce problème peut être lié à la qualité limitée du modèle de résistivité de base utilisé, indiquant ainsi que des efforts plus grands devront être fournis dans le futur pour obtenir des modèles de base de bonne qualité avant de réaliser des expériences dynamiques. Les études décrites dans cette thèse montrent que les méthodes d'ondes EM planes sont très utiles pour caractériser et suivre les variations temporelles du sous-sol sur de larges échelles. Les présentes approches améliorent l'évaluation des modèles obtenus, autant en termes d'incorporation d'informations a priori, qu'en termes de quantification d'incertitudes a posteriori. De plus, les stratégies développées peuvent être appliquées à d'autres méthodes géophysiques, et offrent une grande flexibilité pour l'incorporation d'informations additionnelles lorsqu'elles sont disponibles. -- The efficient use of geothermal systems, the sequestration of CO2 to mitigate climate change, and the prevention of seawater intrusion in coastal aquifers are only some examples that demonstrate the need for novel technologies to monitor subsurface processes from the surface. A main challenge is to assure optimal performance of such technologies at different temporal and spatial scales. Plane-wave electromagnetic (EM) methods are sensitive to subsurface electrical conductivity and consequently to fluid conductivity, fracture connectivity, temperature, and rock mineralogy. These methods have governing equations that are the same over a large range of frequencies, thus allowing to study in an analogous manner processes on scales ranging from few meters close to the surface down to several hundreds of kilometers depth. Unfortunately, they suffer from a significant resolution loss with depth due to the diffusive nature of the electromagnetic fields. Therefore, estimations of subsurface models that use these methods should incorporate a priori information to better constrain the models, and provide appropriate measures of model uncertainty. During my thesis, I have developed approaches to improve the static and dynamic characterization of the subsurface with plane-wave EM methods. In the first part of this thesis, I present a two-dimensional deterministic approach to perform time-lapse inversion of plane-wave EM data. The strategy is based on the incorporation of prior information into the inversion algorithm regarding the expected temporal changes in electrical conductivity. This is done by incorporating a flexible stochastic regularization and constraints regarding the expected ranges of the changes by using Lagrange multipliers. I use non-l2 norms to penalize the model update in order to obtain sharp transitions between regions that experience temporal changes and regions that do not. I also incorporate a time-lapse differencing strategy to remove systematic errors in the time-lapse inversion. This work presents improvements in the characterization of temporal changes with respect to the classical approach of performing separate inversions and computing differences between the models. In the second part of this thesis, I adopt a Bayesian framework and use Markov chain Monte Carlo (MCMC) simulations to quantify model parameter uncertainty in plane-wave EM inversion. For this purpose, I present a two-dimensional pixel-based probabilistic inversion strategy for separate and joint inversions of plane-wave EM and electrical resistivity tomography (ERT) data. I compare the uncertainties of the model parameters when considering different types of prior information on the model structure and different likelihood functions to describe the data errors. The results indicate that model regularization is necessary when dealing with a large number of model parameters because it helps to accelerate the convergence of the chains and leads to more realistic models. These constraints also lead to smaller uncertainty estimates, which imply posterior distributions that do not include the true underlying model in regions where the method has limited sensitivity. This situation can be improved by combining planewave EM methods with complimentary geophysical methods such as ERT. In addition, I show that an appropriate regularization weight and the standard deviation of the data errors can be retrieved by the MCMC inversion. Finally, I evaluate the possibility of characterizing the three-dimensional distribution of an injected water plume by performing three-dimensional time-lapse MCMC inversion of planewave EM data. Since MCMC inversion involves a significant computational burden in high parameter dimensions, I propose a model reduction strategy where the coefficients of a Legendre moment decomposition of the injected water plume and its location are estimated. For this purpose, a base resistivity model is needed which is obtained prior to the time-lapse experiment. A synthetic test shows that the methodology works well when the base resistivity model is correctly characterized. The methodology is also applied to an injection experiment performed in a geothermal system in Australia, and compared to a three-dimensional time-lapse inversion performed within a deterministic framework. The MCMC inversion better constrains the water plumes due to the larger amount of prior information that is included in the algorithm. The conductivity changes needed to explain the time-lapse data are much larger than what is physically possible based on present day understandings. This issue may be related to the base resistivity model used, therefore indicating that more efforts should be given to obtain high-quality base models prior to dynamic experiments. The studies described herein give clear evidence that plane-wave EM methods are useful to characterize and monitor the subsurface at a wide range of scales. The presented approaches contribute to an improved appraisal of the obtained models, both in terms of the incorporation of prior information in the algorithms and the posterior uncertainty quantification. In addition, the developed strategies can be applied to other geophysical methods, and offer great flexibility to incorporate additional information when available.

A link between eumelanism and calcium physiology in the barn owl.

In many animals, melanin-based coloration is strongly heritable and is largely insensitive to the environment and body condition. According to the handicap principle, such a trait may not reveal individual quality because the production of different melanin-based colorations often entails similar costs. However, a recent study showed that the production of eumelanin pigments requires relatively large amounts of calcium, potentially implying that melanin-based coloration is associated with physiological processes requiring calcium. If this is the case, eumelanism may be traded-off against other metabolic processes that require the same elements. We used a correlative approach to examine, for the first time, this proposition in the barn owl, a species in which individuals vary in the amount, size, and blackness of eumelanic spots. For this purpose, we measured calcium concentration in the left humerus of 85 dead owls. Results showed that the humeri of heavily spotted individuals had a higher concentration of calcium. This suggests either that plumage spottiness signals the ability to absorb calcium from the diet for both eumelanin production and storage in bones, or that lightly spotted individuals use more calcium for metabolic processes at the expense of calcium storage in bones. Our study supports the idea that eumelanin-based coloration is associated with a number of physiological processes requiring calcium.

Numerical analysis of wave-induced fluid flow effects on seismic data: Application to monitoring of CO2 storage at the Sleipner Field

In this work we analyze how patchy distributions of CO2 and brine within sand reservoirs may lead to significant attenuation and velocity dispersion effects, which in turn may have a profound impact on surface seismic data. The ultimate goal of this paper is to contribute to the understanding of these processes within the framework of the seismic monitoring of CO2 sequestration, a key strategy to mitigate global warming. We first carry out a Monte Carlo analysis to study the statistical behavior of attenuation and velocity dispersion of compressional waves traveling through rocks with properties similar to those at the Utsira Sand, Sleipner field, containing quasi-fractal patchy distributions of CO2 and brine. These results show that the mean patch size and CO2 saturation play key roles in the observed wave-induced fluid flow effects. The latter can be remarkably important when CO2 concentrations are low and mean patch sizes are relatively large. To analyze these effects on the corresponding surface seismic data, we perform numerical simulations of wave propagation considering reservoir models and CO2 accumulation patterns similar to the CO2 injection site in the Sleipner field. These numerical experiments suggest that wave-induced fluid flow effects may produce changes in the reservoir's seismic response, modifying significantly the main seismic attributes usually employed in the characterization of these environments. Consequently, the determination of the nature of the fluid distributions as well as the proper modeling of the seismic data constitute important aspects that should not be ignored in the seismic monitoring of CO2 sequestration problems.

Accurate calculations of intermolecular interaction energies using explicitly correlated coupled cluster wave functions and a dispersion-weighted MP2 method.

Explicitly correlated coupled-cluster calculations of intermolecular interaction energies for the S22 benchmark set of Jurecka, Sponer, Cerny, and Hobza (Chem. Phys. Phys. Chem. 2006, 8, 1985) are presented. Results obtained with the recently proposed CCSD(T)-F12a method and augmented double-zeta basis sets are found to be in very close agreement with basis set extrapolated conventional CCSD(T) results. Furthermore, we propose a dispersion-weighted MP2 (DW-MP2) approximation that combines the good accuracy of MP2 for complexes with predominately electrostatic bonding and SCS-MP2 for dispersion-dominated ones. The MP2-F12 and SCS-MP2-F12 correlation energies are weighted by a switching function that depends on the relative HF and correlation contributions to the interaction energy. For the S22 set, this yields a mean absolute deviation of 0.2 kcal/mol from the CCSD(T)-F12a results. The method, which allows obtaining accurate results at low cost, is also tested for a number of dimers that are not in the training set.