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.

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.

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.

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.

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).

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.

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.

Bioavailability of carotenoids and tocopherols from broccoli: in vivo and in vitro assessment.

Fruits and vegetables are the major sources of biologically active compounds, and carotenoids and tocopherols constitute important groups in human diets. Bioavailability is a critical feature in the assessment of the role of micronutrients in human health, and the approaches to this issue include in vitro and in vivo methods. Our aim was to evaluate the bioavailability of carotenoids and tocopherols present in broccoli and to compare in vitro and in vivo approaches. Fourteen apparently healthy volunteers consumed 200 g broccoli once a day for seven days. Blood samples were drawn at baseline and after intervention to determine changes in lutein, beta-carotene, and alpha- and gamma-tocopherol as relevant phytochemicals provided with this vegetable. Broccoli also was subjected to simulated gastrointestinal digestion to assess changes related to preabsorptive processes. Analytes in serum and at each phase of the digestion were assayed by high-performance liquid chromatography. During the intervention, the amounts supplied daily ranged from 2.4 to 3.1 mg lutein, 1.4 to 1.8 mg beta-carotene, 4.5 to 6.8 mg alpha-tocopherol, and 0.8 to 1.8 mg gamma-tocopherol. Significant changes in serum in both men and women were observed only for lutein, whereas for gamma-tocopherol a significant change was detected in women. No changes were observed for alpha-tocopherol, beta-carotene, retinol, the alpha-tocopherol-to-cholesterol ratio, or serum lipids. Using the in vitro model, more than 75% of lutein, beta-carotene, gamma-tocopherol, and alpha-tocopherol remained at the duodenal phase, whereas incorporation into the supernatants accounted for <20% of the initial content in food. Regular consumption of broccoli at dietary levels increased serum concentrations of lutein and gamma-tocopherol without affecting alpha-tocopherol or beta-carotene status in serum. The behavior of these phytochemicals under in vitro gastrointestinal conditions does not fully explain the changes observed in vivo.

[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.

Understanding the cytological diploidization mechanism of polyploid wild wheats.

The allohexaploid Aegilops species (2n = 6x = 42), Ae. neglecta 6x (UUXtXtNN), Ae. juvenalis (DcDcXcXcUU), and Ae. vavilovii (DcDcXcXcSsSs) regularly form bivalents at metaphase I. However, in Ae. crassa 6x (DcDcXcXcDD) 0.27 quadrivalents per cell were observed probably as a consequence of the partial homology displayed by the D and Dc genomes. Likewise, the synthetic amphiploid Ae. ventricosa-Secale cereale (DDNNRR) is fertile and displays a diploid-like behavior at metaphase I, despite its recent origin. The pattern of synapsis at late zygotene and pachytene in the natural and artificial allohexaploids was analyzed by whole-mount surface-spreading of synaptonemal complexes under an electron microscope. It revealed that chromosomes were mostly associated as bivalents in all cases, the mean of multivalents per nucleus ranging from 0.17 (Ae. neglecta 6x) to 1.03 (Ae. crassa 6x) in the natural species and 1.05 in the Ae. ventricosa-S. cereale amphiploid. It can be concluded that the mechanism controlling bivalent formation in these species and also in the synthetic amphiploid acts mainly at zygotene by restricting synapsis to homologous chromosomes, but also acts at pachytene by preventing chiasma formation in the homoeologous associations. These observations are discussed in relation to the origin and evolution of the mechanism of diploidization in the allopolyploid species of the Poaceae family.

Prognostic value of apoptosis in breast cancer (pT1-pT2). A TUNEL, p53, bcl-2, bag-1 and Bax immunohistochemical study.

Apoptosis or programmed cell death produces cells breaking into several fragments of nuclei, cytoplasm or both nuclei and cytoplasm, known as apoptotic bodies which can be visualized in haematoxylin-eosin staining. Some genes (promoters and suppressors) control this process and certain mutations may induce the expression of abnormal proteins, which can be detected by immunohistochemical staining. Apoptosis can be detected by the TUNEL method either identifying apoptotic bodies or cells at the initial stages of the fragmentation process. We have studied 186 cases of infiltrating ductal breast carcinoma, stages pT1-pT2, and analysed the prognostic significance of tumour recurrence and overall survival of apoptotic index (AI) through univariate and multivariate analysis. We have also studied the immunohistochemical protein expression of apoptosis promoter and suppressors gene (p53, nuclear expression; bcl-2 and Bax, cytoplasm expression; BAG-1, nuclear and cytoplasm expression). The results indicate prognostic significance of p53 and bcl-2 related to patient death and bcl-2 and tumour size to tumour recurrence, bcl-2 acting as a protector factor (apoptotic suppressor) in both situations. On the other hand, we have not found useful prognostic information of AI either to tumour recurrence or overall survival in univariate or multivariate studies. In this study, Bax expression does not provide a new prognostic role in breast carcinoma, although it contrasts to the bcl-2 action and accelerates death.

ECMO en el tratamiento de la miocarditis y en el shock séptico / No disponible

No disponible

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.