Accuracy limit of non-polarizable four-point water models: TIP4P/2005 vs OPC. Should water models reproduce the experimental dielectric constant?

The last generation of four center non-polarizable models of water can be divided into two groups: those reproducing the dielectric constant of water, as OPC, and those significantly underestimating its value, as TIP4P/2005. To evaluate the global performance of OPC and TIP4P/2005, we shall follow the test proposed by Vega and Abascal in 2011 evaluating about 40 properties to fairly address this comparison. The liquid-vapor and liquid-solid equilibria are computed, as well as the heat capacities, isothermal compressibilities, surface tensions, densities of different ice polymorphs, the density maximum, equations of state at high pressures, and transport properties. General aspects of the phase diagram are considered by comparing the ratios of different temperatures (namely, the temperature of maximum density, the melting temperature of hexagonal ice, and the critical temperature). The final scores are 7.2 for TIP4P/2005 and 6.3 for OPC. The results of this work strongly suggest that we have reached the limit of what can be achieved with non-polarizable models of water and that the attempt to reproduce the experimental dielectric constant deteriorates the global performance of the water force field. The reason is that the dielectric constant depends on two surfaces (potential energy and dipole moment surfaces), whereas in the absence of an electric field, all properties can be determined simply from just one surface (the potential energy surface). The consequences of the choice of the water model in the modeling of electrolytes in water are also discussed.

The temperature of maximum in density of aqueous solutions of nitrate and ammonium salts: Testing the Madrid-2019 force field.

The shift in the temperature of maximum in density (TMD) at room pressure of aqueous solutions of a set of five salts containing NO3- and/or NH4+ groups is studied both through experiments and through molecular dynamics simulations using the Madrid-2019 force field for ions and the TIP4P/2005 model for water. The experiments demonstrate the potential transferability and limitations of the Madrid-2019 force field for nitrate and ammonium ions recently developed by our group at different temperatures and add updated information to the reported datasets of TMDs for strong electrolytes. By using the Despretz law, individual ion contributions are extracted for predictive purposes from the experimental values of the shift in the TMD. Interesting findings for the behavior of the shift in the TMD in nitrate salts expose that this property might be particularly challenging for modelization approaches when dealing with polyatomic species.

Three-phase equilibria of hydrates from computer simulation. I. Finite-size effects in the methane hydrate.

Clathrate hydrates are vital in energy research and environmental applications. Understanding their stability is crucial for harnessing their potential. In this work, we employ direct coexistence simulations to study finite-size effects in the determination of the three-phase equilibrium temperature (T3) for methane hydrates. Two popular water models, TIP4P/Ice and TIP4P/2005, are employed, exploring various system sizes by varying the number of molecules in the hydrate, liquid, and gas phases. The results reveal that finite-size effects play a crucial role in determining T3. The study includes nine configurations with varying system sizes, demonstrating that smaller systems, particularly those leading to stoichiometric conditions and bubble formation, may yield inaccurate T3 values. The emergence of methane bubbles within the liquid phase, observed in smaller configurations, significantly influences the behavior of the system and can lead to erroneous temperature estimations. Our findings reveal finite-size effects on the calculation of T3 by direct coexistence simulations and clarify the system size convergence for both models, shedding light on discrepancies found in the literature. The results contribute to a deeper understanding of the phase equilibrium of gas hydrates and offer valuable information for future research in this field.

Three-phase equilibria of hydrates from computer simulation. II. Finite-size effects in the carbon dioxide hydrate.

In this work, the effects of finite size on the determination of the three-phase coexistence temperature (T3) of the carbon dioxide (CO2) hydrate have been studied by molecular dynamic simulations and using the direct coexistence technique. According to this technique, the three phases involved (hydrate-aqueous solution-liquid CO2) are placed together in the same simulation box. By varying the number of molecules of each phase, it is possible to analyze the effect of simulation size and stoichiometry on the T3 determination. In this work, we have determined the T3 value at 8 different pressures (from 100 to 6000 bar) and using 6 different simulation boxes with different numbers of molecules and sizes. In two of these configurations, the ratio of the number of water and CO2 molecules in the aqueous solution and the liquid CO2 phase is the same as in the hydrate (stoichiometric configuration). In both stoichiometric configurations, the formation of a liquid drop of CO2 in the aqueous phase is observed. This drop, which has a cylindrical geometry, increases the amount of CO2 available in the aqueous solution and can in some cases lead to the crystallization of the hydrate at temperatures above T3, overestimating the T3 value obtained from direct coexistence simulations. The simulation results obtained for the CO2 hydrate confirm the sensitivity of T3 depending on the size and composition of the system, explaining the discrepancies observed in the original work by Míguez et al. [J. Chem Phys. 142, 124505 (2015)]. Non-stoichiometric configurations with larger unit cells show a convergence of T3 values, suggesting that finite-size effects for these system sizes, regardless of drop formation, can be safely neglected. The results obtained in this work highlight that the choice of a correct initial configuration is essential to accurately estimate the three-phase coexistence temperature of hydrates by direct coexistence simulations.

Three-phase equilibria of hydrates from computer simulation. III. Effect of dispersive interactions in the methane and carbon dioxide hydrates.

In this work, the effect of the range of dispersive interactions in determining the three-phase coexistence line of the CO2 and CH4 hydrates has been studied. In particular, the temperature (T3) at which solid hydrate, water, and liquid CO2/gas CH4 coexist has been determined through molecular dynamics simulations using different cutoff values (from 0.9 to 1.6 nm) for dispersive interactions. The T3 of both hydrates has been determined using the direct coexistence simulation technique. Following this method, the three phases in equilibrium are put together in the same simulation box, the pressure is fixed, and simulations are performed at different temperatures T. If the hydrate melts, then T > T3. Conversely, if the hydrate grows, then T < T3. The effect of the cutoff distance on the dissociation temperature has been analyzed at three different pressures for CO2 hydrate: 100, 400, and 1000 bar. Then, we have changed the guest and studied the effect of the cutoff distance on the dissociation temperature of the CH4 hydrate at 400 bar. Moreover, the effect of long-range corrections for dispersive interactions has been analyzed by running simulations with homo- and inhomogeneous corrections and a cutoff value of 0.9 nm. The results obtained in this work highlight that the cutoff distance for the dispersive interactions affects the stability conditions of these hydrates. This effect is enhanced when the pressure is decreased, displacing the T3 about 2-4 K depending on the system and the pressure.

Madrid-2019 force field: An extension to divalent cations Sr2+ and Ba2.

In this work, we present a parameterization of Sr2+ and Ba2+ cations, which expands the alkali earth set of cations of the Madrid-2019 force field. We have tested the model against the experimental densities of eight different salts, namely, SrCl2, SrBr2, SrI2, Sr(NO3)2, BaCl2, BaBr2, BaI2, and Ba(NO3)2. The force field is able to reproduce the experimental densities of all these salts up to their solubility limit. Furthermore, we have computed the viscosities for two selected salts, finding that the experimental values are overestimated, but the predictions are still reasonable. Finally, the structural properties for all the salts have been calculated with this model and align remarkably well with experimental observations.

Growth rate of CO2 and CH4 hydrates by means of molecular dynamics simulations.

CO2 and CH4 hydrates are of great importance both from an energetic and from an environmental point of view. It is therefore highly relevant to quantify and understand the rate with which they grow. We use molecular dynamics simulations to shed light on the growth rate of these hydrates. We put the solid hydrate phase in contact with a guest aqueous solution in equilibrium with the pure guest phase and study the growth of both hydrates at 400 bars with temperature. We compare our results with previous calculations of the ice growth rate. We find a growth rate maximum as a function of the supercooling in all cases. The incorporation of guest molecules into the solid structure strongly decelerates hydrate growth. Consistently, ice grows faster than either hydrate and the CO2 hydrate grows faster than the CH4 one because of the higher solubility of CO2. We also quantify the molecular motion required to build the solids under study and find that the distance traveled by liquid molecules exceeds by orders of magnitude that advanced by any solid. Less molecular motion is needed in order for ice to grow as compared to the hydrates. Moreover, when temperature increases, more motion is needed for solid growth. Finally, we find a good agreement between our growth rate calculations and experiments of hydrate growth along the guest-solution interface. However, more work is needed to reconcile experiments of hydrate growth toward the solution among each other and with simulations.

Homogeneous nucleation rate of methane hydrate formation under experimental conditions from seeding simulations.

In this work, we shall estimate via computer simulations the homogeneous nucleation rate for the methane hydrate at 400 bars for a supercooling of about 35 K. The TIP4P/ICE model and a Lennard-Jones center were used for water and methane, respectively. To estimate the nucleation rate, the seeding technique was employed. Clusters of the methane hydrate of different sizes were inserted into the aqueous phase of a two-phase gas-liquid equilibrium system at 260 K and 400 bars. Using these systems, we determined the size at which the cluster of the hydrate is critical (i.e., it has 50% probability of either growing or melting). Since nucleation rates estimated from the seeding technique are sensitive to the choice of the order parameter used to determine the size of the cluster of the solid, we considered several possibilities. We performed brute force simulations of an aqueous solution of methane in water in which the concentration of methane was several times higher than the equilibrium concentration (i.e., the solution was supersaturated). From brute force runs, we infer the value of the nucleation rate for this system rigorously. Subsequently, seeding runs were carried out for this system, and it was found that only two of the considered order parameters were able to reproduce the value of the nucleation rate obtained from brute force simulations. By using these two order parameters, we estimated the nucleation rate under experimental conditions (400 bars and 260 K) to be of the order of log10 (J/(m3 s)) = -7(5).

Scaled charges for ions: An improvement but not the final word for modeling electrolytes in water.

In this work, we discuss the use of scaled charges when developing force fields for NaCl in water. We shall develop force fields for Na+ and Cl- using the following values for the scaled charge (in electron units): ±0.75, ±0.80, ±0.85, and ±0.92 along with the TIP4P/2005 model of water (for which previous force fields were proposed for q = ±0.85 and q = ±1). The properties considered in this work are densities, structural properties, transport properties, surface tension, freezing point depression, and maximum in density. All the developed models were able to describe quite well the experimental values of the densities. Structural properties were well described by models with charges equal to or larger than ±0.85, surface tension by the charge ±0.92, maximum in density by the charge ±0.85, and transport properties by the charge ±0.75. The use of a scaled charge of ±0.75 is able to reproduce with high accuracy the viscosities and diffusion coefficients of NaCl solutions for the first time. We have also considered the case of KCl in water, and the results obtained were fully consistent with those of NaCl. There is no value of the scaled charge able to reproduce all the properties considered in this work. Although certainly scaled charges are not the final word in the development of force fields for electrolytes in water, its use may have some practical advantages. Certain values of the scaled charge could be the best option when the interest is to describe certain experimental properties.

Melting points of water models: Current situation.

By using the direct coexistence method, we have calculated the melting points of ice Ih at normal pressure for three recently proposed water models, namely, TIP3P-FB, TIP4P-FB, and TIP4P-D. We obtained Tm = 216 K for TIP3P-FB, Tm = 242 K for TIP4P-FB, and Tm = 247 K for TIP4P-D. We revisited the melting point of TIP4P/2005 and TIP5P obtaining Tm = 250 and 274 K, respectively. We summarize the current situation of the melting point of ice Ih for a number of water models and conclude that no model is yet able to simultaneously reproduce the melting temperature of ice Ih and the temperature of the maximum in density at room pressure. This probably points toward our both still incomplete knowledge of the potential energy surface of water and the necessity of incorporating nuclear quantum effects to describe both properties simultaneously.

Maximum in density of electrolyte solutions: Learning about ion-water interactions and testing the Madrid-2019 force field.

In this work, we studied the effect of Li+, Na+, K+, Mg2+, and Ca2+ chlorides and sulfates on the temperature of maximum density (TMD) of aqueous solutions at room pressure. Experiments at 1 molal salt concentration were carried out to determine the TMD of these solutions. We also performed molecular dynamics simulations to estimate the TMD at 1 and 2 m with the Madrid-2019 force field, which uses the TIP4P/2005 water model and scaled charges for the ions, finding an excellent agreement between experiment and simulation. All the salts studied in this work shift the TMD of the solution to lower temperatures and flatten the density vs temperature curves (when compared to pure water) with increasing salt concentration. The shift in the TMD depends strongly on the nature of the electrolyte. In order to explore this dependence, we have evaluated the contribution of each ion to the shift in the TMD concluding that Na+, Ca2+, and SO4 2- seem to induce the largest changes among the studied ions. The volume of the system has been analyzed for salts with the same anion and different cations. These curves provide insight into the effect of different ions upon the structure of water. We claim that the TMD of electrolyte solutions entails interesting physics regarding ion-water and water-water interactions and should, therefore, be considered as a test property when developing force fields for electrolytes. This matter has been rather unnoticed for almost a century now and we believe it is time to revisit it.

The Madrid-2019 force field for electrolytes in water using TIP4P/2005 and scaled charges: Extension to the ions F-, Br-, I-, Rb+, and Cs.

In this work, an extension of the Madrid-2019 force field is presented. We have added the cations Rb+ and Cs+ and the anions F-, Br-, and I-. These ions were the remaining alkaline and halogen ions, not previously considered in the Madrid-2019 force field. The force field, denoted as Madrid-2019-Extended, does not include polarizability and uses the TIP4P/2005 model of water and scaled charges for the ions. A charge of ±0.85e is assigned to monovalent ions. The force field developed provides an accurate description of aqueous solution densities over a wide range of concentrations up to the solubility limit of each salt studied. Good predictions of viscosity and diffusion coefficients are obtained for concentrations below 2 m. Structural properties obtained with this force field are also in reasonable agreement with the experiment. The number of contact ion pairs has been controlled to be low so as to avoid precipitation of the system at concentrations close to the experimental solubility limit. A comprehensive comparison of the performance for aqueous solutions of alkaline halides of force fields of electrolytes using scaled and integer charges is now possible. This comparison will help in the future to learn about the benefits and limitations of the use of scaled charges to describe electrolyte solutions.

The impact of water vapor on the OH reactivity toward CH3CHO at ultra-low temperatures (21.7-135.0 K): Experiments and theory.

The role of water vapor (H2O) and its hydrogen-bonded complexes in the gas-phase reactivity of organic compounds with hydroxyl (OH) radicals has been the subject of many recent studies. Contradictory effects have been reported at temperatures between 200 and 400 K. For the OH + acetaldehyde reaction, a slight catalytic effect of H2O was previously reported at temperatures between 60 and 118 K. In this work, we used Laval nozzle expansions to reinvestigate the impact of H2O on the OH-reactivity with acetaldehyde between 21.7 and 135.0 K. The results of this comprehensive study demonstrate that water, instead, slows down the reaction by factors of ∼3 (21.7 K) and ∼2 (36.2-89.5 K), and almost no effect of added H2O was observed at 135.0 K.

Gas phase kinetics of the OH + CH3CH2OH reaction at temperatures of the interstellar medium (T = 21-107 K).

Ethanol, CH3CH2OH, has been unveiled in the interstellar medium (ISM) by radioastronomy and it is thought to be released into the gas phase after the warm-up phase of the grain surface, where it is formed. Once in the gas phase, it can be destroyed by different reactions with atomic and radical species, such as hydroxyl (OH) radicals. The knowledge of the rate coefficients of all these processes at temperatures of the ISM is essential in the accurate interpretation of the observed abundances. In this work, we have determined the rate coefficient for the reaction of OH with CH3CH2OH (k(T)) between 21 and 107 K by employing the pulsed and continuous CRESU (Cinétique de Réaction en Ecoulement Supersonique Uniforme, which means Reaction Kinetics in a Uniform Supersonic Flow) technique. The pulsed laser photolysis technique was used for generating OH radicals, whose time evolution was monitored by laser induced fluorescence. An increase of approximately 4 times was observed for k(21 K) with respect to k(107 K). With respect to k(300 K), the OH-reactivity at 21 K is enhanced by two orders of magnitude. The obtained T-expression in the investigated temperature range is k(T) = (2.1 ± 0.5) × 10-11 (T/300 K)-(0.71±0.10) cm3 molecule-1 s-1. In addition, the pressure dependence of k(T) has been investigated at several temperatures between 21 K and 90 K. No pressure dependence of k(T) was observed in the investigated ranges. This may imply that this reaction is purely bimolecular or that the high-pressure limit is reached at the lowest total pressure experimentally accessible in our system. From our results, k(T) at usual IS temperatures (∼10-100 K) is confirmed to be very fast. Typical rate coefficients can be considered to range within about 4 × 10-11 cm3 molecule-1 s-1 at 100 K and around 1 × 10-10 cm3 molecule-1 s-1 at 20 K. The extrapolation of k at the lowest temperatures of the dense molecular clouds of ISM is also discussed in this paper.

[Analysis of tracheostomies in a pediatric intensive care unit during the period 2003-2013]. / Análisis de las traqueotomías en cuidados intensivos pediátricos durante el periodo 2003-2013.

INTRODUCTION: Tracheotomy in pediatric patients is a rare procedure. In this pediatric series, perioperative complications, mortality related to surgical procedure and overall mortality are analyzed. PATIENTS AND METHODS: This is a retrospective study conducted from January 2003 to December 2013. Data were retrieved from patients who were tracheotomized and admitted to our PICU in the postoperative period. RESULTS: Data were collected from 25 tracheotomized patients admitted during the study period. The mean age was 3.3 months (median 14 months, range 1-144 months), and PICU length of stay was 53 days (median 37 days, range 1-338 days). Most patients (68%) had comorbidities before their admission, with a higher prevalene of craniofacial anomalies/polymalformative syndromes (32%) and prematurity related disorders (12%) being obserevd. The most common etiologies related to the procedure were congenital airway obstruction (16%) and several types of spinal cord injury (16%), followed by tracheobronchomalacia (12%) and subglottic stenosis (12%). Some kind of complication was detected in 40% of patients, with accidental decannulation being the most frequent. Accidental or unexpected decannulation was present in a percentage as high as 20% of our patients, mainly in the first 24 hours after surgery. One of the patients died as a result of this. CONCLUSIONS: The postoperative course of a tracheotomy is associated with a high rate of complications, some of them related to life-threatening events.

Behavior of natural radionuclides in wastewater treatment plants.

56 samples, including influent, primary effluent, secondary effluent and final effluent wastewater from two Spanish municipal wastewater treatment plants (WWTPs), were analyzed to assess both the occurrence and behavior of natural radioactivity during 12 sampling campaigns carried out over the period 2007-2010. Influent and final effluent wastewaters were sampled by taking into account the hydraulic residence time within the WWTP. A wide range of gross alpha activities (15-129 mBq/L) and gross beta activities (477-983 mBq/L) in liquid samples were obtained. A correlation analysis between radioactivity in liquid samples and the performance characteristics of the WWTPs was performed. The results in liquid samples showed that gross beta activities were not influenced by treatment in the studied WWTPs. However, gross alpha activities behave differently and an increase was detected in the effluent values compared with influent wastewater. This behavior was due to the increase in the total dissolved uranium produced during secondary treatment. The results indicate that the radiological characteristics of the effluents do not present a significant radiological risk and make them suitable for future applications.

Distribution of uranium isotopes in surface water of the Llobregat river basin (Northeast Spain).

A study is presented on the distribution of (234)U, (238)U, (235)U isotopes in surface water of the Llobregat river basin (Northeast Spain), from 2001 to 2006. Sixty-six superficial water samples were collected at 16 points distributed throughout the Llobregat river basin. Uranium isotopes were measured by alpha spectrometry (PIPS detectors). The test procedure was validated according to the quality requirements of the ISO17025 standard. The activity concentration for the total dissolved uranium ranges from 20 to 261mBqL(-1). The highest concentrations of uranium were detected in an area with formations of sedimentary rock, limestone and lignite. A high degree of radioactive disequilibrium was noted among the uranium isotopes. The (234)U/(238)U activity ratio varied between 1.1 and 1.9 and the waters with the lowest uranium activity registered the highest level of (234)U/(238)U activity ratio. Correlations between uranium activity in the tested water and chemical and physical characteristics of the aquifer were found.

Natural and anthropogenic radionuclides in airborne particulate samples collected in Barcelona (Spain).

Results for naturally occurring (7)Be, (210)Pb, (40)K, (214)Bi, (214)Pb, (212)Pb, (228)Ac and (208)Tl and anthropogenic (137)Cs in airborne particulate matter in the Barcelona area during the period from January 2001 to December 2005 are presented and discussed. The (212)Pb and (208)Tl, (214)Bi and (214)Pb, (7)Be and (210)Pb radionuclide levels showed a significant correlation with each other, with correlation coefficients of 0.99, 0.78 and 0.69, respectively, suggesting similar origin/behaviour of these radionuclides in the air. Caessium-137 and Potassium-40 were transported to the air as resuspended particle from the soil. The (7)Be and (210)Pb concentrations showed similar seasonal variations, with a tendency for maximum concentrations during the summer months. An inverse relationship was observed between the (7)Be, (210)Pb, (40)K and (137)Cs concentrations and weekly rainfall, indicating washout of atmospheric aerosols carrying these radionuclides.

Dermoscopic pattern of intermediate stage in seborrhoeic keratosis regressing to lichenoid keratosis: report of 24 cases.

BACKGROUND: Lichenoid keratosis (LK) is a well-described entity which has been proposed to represent an immunological or regressive response to pre-existing epidermal lesions such as solar lentigines or seborrhoeic keratoses. OBJECTIVES: To evaluate the dermoscopic criteria of a series of cases of LK with remaining areas of seborrhoeic keratosis which were both dermoscopically and histologically diagnosed. METHODS: Pigmented lesions with dermoscopic areas of seborrhoeic keratosis and LK in the same tumour were consecutively diagnosed and prospectively included in the study. All pigmented lesions were examined and registered using DermLite Foto equipment (3Gen, LLC, Dana Point, CA, U.S.A.), at 10-fold magnification, at the Dermatology Department of Hospital de Sant Pau i Santa Tecla (Tarragona, Spain), between 1 January 2003 and 31 December 2005. RESULTS: In total, 24 cases of lesions with dermoscopic areas of seborrhoeic keratosis and LK were collected. In four lesions (17%), the clinical differential diagnosis without dermoscopy included malignant melanoma and in seven lesions (29%), basal cell carcinoma. The diagnosis of LK was clinically considered without dermoscopy in only six cases (25%). A granular pattern was observed to be distributed throughout the LK areas of the lesions. This pattern consisted of the presence of brownish-grey, bluish-grey or whitish-grey coarse granules that formed, in 11 cases (46%), globules and/or short lines. In one lesion, located on the face, these short lines produced annular or rhomboid structures as seen in lentigo maligna melanoma. CONCLUSIONS: Dermoscopy is a useful tool which assists in the correct clinical recognition of LK, which may also potentially illuminate the pathogenesis of these tumours, showing the intermediate stage of regressing epidermal lesions in an LK.

Changes in carotenoid intake from fruit and vegetables in the Spanish population over the period 1964-2004.

OBJECTIVE: To assess changes in carotenoid intake based on the variations in the consumption of fresh fruit and vegetables in the Spanish population over the period 1964-2004. DESIGN: Consumption data of fresh fruit and vegetables from Family Budget Surveys carried out in 1964, 1980, 1990 and 2004. Consumption data (g per person per day) accounted for >90% of fruit and vegetable consumption at each time point. Quality controlled high-performance liquid chromatography analysis of the carotenoid composition of Spanish fruit and vegetable was used. SUBJECTS: Randomly selected, private households throughout Spain (20,800 households in 1964, 30,311 households in 1980, 21,155 households in 1990 and 6000 households in 2004). Twelve vegetables and 16 fruits representing 89-96% of the total consumption of fresh fruit and vegetables were used. RESULTS: Individual consumption of fruit and vegetables has changed over this period, altering the total and individual intake of carotenoids. Total carotenoid intake increased from 2.5 mg per person per day in 1964 to 4.1 mg per person per day in 1990, with a decrease to 3.3 mg per person per day in 2004. These increments are due to an increase in lycopene, alpha- and beta-carotene, while a decrease in lutein and zeaxanthin is observed during the last decade. A continuous and consistent decrease in the relative contribution of lutein in the diet is observed over the period studied. CONCLUSION: Although the consumption of fruit and vegetables is still consistent with a Mediterranean-type pattern, modifications in the consumption of individual fruits and vegetables have provoked changes in total and specific carotenoid intake with potential relevance in human health.