Molecular dynamics simulation of methane in sodium montmorillonite clay hydrates at elevated pressures and temperatures

Computer simulation has been used to study the structure and dynamics of methane in hydrated sodium montmorillonite clays under conditions encountered in sedimentary basins. Systems containing approximately one, two, three and four molecular layers of water have followed gradients of 150 bar km-1 and 30Kkm-1, to a maximum burial depth of 6 km (900 bar and 460 K). Methane is coordinated to approximately 19 oxygen atoms, of which typically 6 are provided by the clay surface. Only in the three- and four-layer hydrates is methane able to leave the clay surface. Diffusion depends strongly on the porosity (water content) and burial depth: self-diffusion coefficients are in the range 0.12 × 10-9m2s-1 for water and 0.04 × 10−9m2s−1 < D < 8.64 × 10−9m2s−1 for methane. Bearing in mind that porosity decreases with burial depth, it is estimated that maximum diffusion occurs at around 3 km. This is in good agreement with the known location of methane reservoirs in sedimentary basins.

MRI pressure and stress measurement in novel homogeneous soft solids

This paper presents MRI measurements of a novel semi solid MR contrast agent to pressure. The agent is comprised of potassium chloride cross linked carageenan gum at a concentration of 2% w/v, with micron size lipid coated bubbles of air at a concentration of 3% v/v. The choice for an optimum suspending medium, the methods of production and the preliminary MRI results are presented herein. The carageenan gum is shown to be ideally elastic for compressions relating to volume changes less than 15%, in contrast to the inelastic gellan gum also tested. Although slightly lower than that of gellan gum, carageenan has a water diffusion coefficient of 1.72×10-9 m2.s-1 indicating its suitability to this purpose. RARE imaging is performed whilst simultaneously compressing test and control samples and a maximum sensitivity of 1.6% MR signal change per % volume change is found which is shown to be independent of proton density variations due to the presence of microbubbles and compression. This contrast agent could prove useful for numerous applications, and particularly in chemical engineering. More generally the method allows the user to non-invasively image with MRI any process that causes, within the solid, local changes either in bubble size or bubble shape. © 2008 American Institute of Physics.

Fish population dynamics in a seasonally varying wetland

Small fishes in seasonally flooded environments such as the Everglades are capable of spreading into newly flooded areas and building up substantial biomass. Passive drift cannot account for the rapidity of observed population expansions. To test the ‘reaction–diffusion’ mechanism for spread of the fish, we estimated their diffusion coefficient and applied a reaction–diffusion model. This mechanism was also too weak to account for the spatial dynamics. Two other hypotheses were tested through modeling. The first—the ‘refuge mechanism’—hypothesizes that small remnant populations of small fishes survive the dry season in small permanent bodies of water (refugia), sites where the water level is otherwise below the surface. The second mechanism, which we call the ‘dynamic ideal free distribution mechanism’ is that consumption by the fish creates a prey density gradient and that fish taxis along this gradient can lead to rapid population expansion in space. We examined the two alternatives and concluded that although refugia may play an important role in recolonization by the fish population during reflooding, only the second, taxis in the direction of the flooding front, seems capable of matching empirical observations. This study has important implications for management of wetlands, as fish biomass is an essential support of higher trophic levels.

Mechanistic Models of Anti-HIV Microbicide Drug Delivery

A new modality for preventing HIV transmission is emerging in the form of topical microbicides. Some clinical trials have shown some promising results of these methods of protection while other trials have failed to show efficacy. Due to the relatively novel nature of microbicide drug transport, a rigorous, deterministic analysis of that transport can help improve the design of microbicide vehicles and understand results from clinical trials. This type of analysis can aid microbicide product design by helping understand and organize the determinants of drug transport and the potential efficacies of candidate microbicide products.

Microbicide drug transport is modeled as a diffusion process with convection and reaction effects in appropriate compartments. This is applied here to vaginal gels and rings and a rectal enema, all delivering the microbicide drug Tenofovir. Although the focus here is on Tenofovir, the methods established in this dissertation can readily be adapted to other drugs, given knowledge of their physical and chemical properties, such as the diffusion coefficient, partition coefficient, and reaction kinetics. Other dosage forms such as tablets and fiber meshes can also be modeled using the perspective and methods developed here.

The analyses here include convective details of intravaginal flows by both ambient fluid and spreading gels with different rheological properties and applied volumes. These are input to the overall conservation equations for drug mass transport in different compartments. The results are Tenofovir concentration distributions in time and space for a variety of microbicide products and conditions. The Tenofovir concentrations in the vaginal and rectal mucosal stroma are converted, via a coupled reaction equation, to concentrations of Tenofovir diphosphate, which is the active form of the drug that functions as a reverse transcriptase inhibitor against HIV. Key model outputs are related to concentrations measured in experimental pharmacokinetic (PK) studies, e.g. concentrations in biopsies and blood. A new measure of microbicide prophylactic functionality, the Percent Protected, is calculated. This is the time dependent volume of the entire stroma (and thus fraction of host cells therein) in which Tenofovir diphosphate concentrations equal or exceed a target prophylactic value, e.g. an EC50.

Results show the prophylactic potentials of the studied microbicide vehicles against HIV infections. Key design parameters for each are addressed in application of the models. For a vaginal gel, fast spreading at small volume is more effective than slower spreading at high volume. Vaginal rings are shown to be most effective if inserted and retained as close to the fornix as possible. Because of the long half-life of Tenofovir diphosphate, temporary removal of the vaginal ring (after achieving steady state) for up to 24h does not appreciably diminish Percent Protected. However, full steady state (for the entire stromal volume) is not achieved until several days after ring insertion. Delivery of Tenofovir to the rectal mucosa by an enema is dominated by surface area of coated mucosa and whether the interiors of rectal crypts are filled with the enema fluid. For the enema 100% Percent Protected is achieved much more rapidly than for vaginal products, primarily because of the much thinner epithelial layer of the mucosa. For example, 100% Percent Protected can be achieved with a one minute enema application, and 15 minute wait time.

Results of these models have good agreement with experimental pharmacokinetic data, in animals and clinical trials. They also improve upon traditional, empirical PK modeling, and this is illustrated here. Our deterministic approach can inform design of sampling in clinical trials by indicating time periods during which significant changes in drug concentrations occur in different compartments. More fundamentally, the work here helps delineate the determinants of microbicide drug delivery. This information can be the key to improved, rational design of microbicide products and their dosage regimens.

Low valence cation doping of bulk Cr2O3: Charge compensation and oxygen vacancy formation

The different oxidation states of chromium allow its bulk oxide form to be reducible, facilitating the oxygen vacancy formation process, which is a key property in applications such as catalysis. Similar to other useful oxides such as TiO2, and CeO2, the effect of substitutional metal dopants in bulk Cr2O3 and its effect on the electronic structure and oxygen vacancy formation are of interest, particularly in enhancing the latter. In this paper, density functional theory (DFT) calculations with a Hubbard + U correction (DFT+U) applied to the Cr 3d and O 2p states, are carried out on pure and metal-doped bulk Cr2O3 to examine the effect of doping on the electronic and geometric structure. The role of dopants in enhancing the reducibility of Cr2O3 is examined to promote oxygen vacancy formation. The dopants are Mg, Cu, Ni, and Zn, which have a formal +2 oxidation state in their bulk oxides. Given this difference in host and, dopant oxidation states, we show that to predict the correct ground state two metal dopants charge compensated with an oxygen vacancy are required. The second oxygen atom removed is termed "the active" oxygen vacancy and it is the energy required to remove this atom that is related to the reduction process. In all cases, we find that substitutional doping improves the oxygen vacancy formation of bulk Cr2O3 by lowering the energy cost.

(Table T1) Density, porosity, void ratio, diffusivity and permeability of ODP Site 185-1149 sediments

The subduction of oceanic plates regulates crustal growth, influences arc volcanism, and refertilizes the mantle. Continental growth occurs by subduction of crustal material (seawater components, marine sediments, and basaltic crust). The geochemical and physical evolution of the Earth's crust depends, in large part, on the fate of subducted material at convergent margins (Armstrong, 1968, doi:10.1029/RG006i002p00175; Karig and Kay, 1981, 10.1098/rsta.1981.0108). The crustal material on the downgoing plate is recycled to various levels in the subduction zone. The recycling process that takes place in the "Subduction Factory" is difficult to observe directly but is clearly illuminated using chemical tracers. Von Huene and Scholl (1991, doi:10.1029/91RG00969) and Plank and Langmuir (1993, doi:10.1038/362739a0) preliminarily calculated a large flux of subducted materials. By mass balancing the chemical tracers and measuring the fractionations that occur between them, the Subduction Factory work and the effect on the Earth's evolution can be estimated. In order to elucidate this mass balance, Ocean Drilling Program Leg 185 drilled two deepwater shales into the oceanic crust situated in the Mariana-Izu Trenches and recovered core samples of incoming oceanic crust. The calculations of mass circulation in the subduction zone, however, did not take into account the mass transfer properties within subducted oceanic crust, although the dewatering fluid and diffused ions may play an important role in various activities such as seismogeneity, serpentine diapiring, and arc volcanism. Thus, this paper focuses on the quantitative measurements of the physical and mass transfer properties of subducted oceanic crust.

C1-C8 hydrocarbons in marine sediments

Light hydrocarbon (C1-C8) profiles are compared for three wells of varying maturities: two immature DSDP wells (Site 397 near the Canary Islands and Site 530A near the Walvis Ridge in the south-east Atlantic) and a mature well, the East Cameron well in the Texas Gulf Coast. Primary migration of C1 and C2 appears to be occurring in all of the sedimentary rocks examined. Primary migration of C3+ components becomes important only as fine-grained sedimentary rocks enter the catagenetic hydrocarbon generation zone or over short distances in more permeable sections. Lateral migration along bedding planes was more important than vertical migration in sedimentary rocks of all maturities. The lightest (methane, ethane and propane gases) hydrocarbon show greater fractionation than do the C4-C8 alkanes which generally show minimal fractionation during the migrational process. Subsurface diffusion coefficients for these p.p.b. quantities of C2-C5 alkanes from immature sediments from DSDP Site 530 are estimated to be several orders of magnitude less than values reported in the literature for diffusion of much larger amounts of these compounds from mature water wet sediments into air or sandstones. Since our calculations suggest light hydrocarbons are present in amounts less than their reported solubilities in pure water at 25°C, we postulate that the sediment organic matter has a substantial effect on retarding the movement of these light hydrocarbons.

Refractometry Studies of the Optical Properties of Polymer Films and the Development of Polymer Coated Refractive Index Sensors

Sensors for real-time monitoring of environmental contaminants are essential for protecting ecosystems and human health. Refractive index sensing is a non-selective technique that can be used to measure almost any analyte. Miniaturized refractive index sensors, such as silicon-on-insulator (SOI) microring resonators are one possible platform, but require coatings selective to the analytes of interest. A homemade prism refractometer is reported and used to characterize the interactions between polymer films and liquid or vapour-phase analytes. A camera was used to capture both Fresnel reflection and total internal reflection within the prism. For thin-films (d = 10 μm - 100 μm), interference fringes were also observed. Fourier analysis of the interferogram allowed for simultaneous extraction of the average refractive index and film thickness with accuracies of ∆n = 1-7 ×10-4 and ∆d < 3-5%. The refractive indices of 29 common organic solvents as well as aqueous solutions of sodium chloride, sucrose, ethylene glycol, glycerol, and dimethylsulfoxide were measured at λ = 1550 nm. These measurements will be useful for future calibrations of near-infrared refractive index sensors. A mathematical model is presented, where the concentration of analyte adsorbed in a film can be calculated from the refractive index and thickness changes during uptake. This model can be used with Fickian diffusion models to measure the diffusion coefficients through the bulk film and at the film-substrate interface. The diffusion of water and other organic solvents into SU-8 epoxy was explored using refractometry and the diffusion coefficient of water into SU-8 is presented. Exposure of soft baked SU-8 films to acetone, acetonitrile and methanol resulted in rapid delamination. The diffusion of volatile organic compound (VOC) vapours into polydimethylsiloxane and polydimethyl-co-polydiphenylsiloxane polymers was also studied using refractometry. Diffusion and partition coefficients are reported for several analytes. As a model system, polydimethyl-co-diphenylsiloxane films were coated onto SOI microring resonators. After the development of data acquisition software, coated devices were exposed to VOCs and the refractive index response was assessed. More studies with other polymers are required to test the viability of this platform for environmental sensing applications.

Moisture Desorption Studies on Polymer Hydrated and Vacuum Extruded Bentonite Clay Mat

Moisture desorption observations from two bentonite clay mats subjected to ten environmental zones with individually different combinations of laboratory-controlled constant temperatures (between 20 °C and 40 °C) and relative humidity (between 15% and 70%) are presented. These laboratory observations are compared with predictions from mathematical models, such as thin-layer drying equations and kinetic drying models proposed by Page, Wang and Singh, and Henderson and Pabis. The quality of fit of these models is assessed using standard error (SE) of estimate, relative percent of error, and coefficient of correlation. The Page model was found to better predict the drying kinetics of the bentonite clay mats for the simulated tropical climates. Critical study on the drying constant and moisture diffusion coefficient helps to assess the efficacy of a polymer to retain moisture and control desorption through water molecule bonding. This is further substantiated with the Guggenheim–Anderson–De Boer (GAB) desorption isotherm model which is presented.

Shapes and Dynamics of Blood Cells in Poiseuille and Shear Flows

The dynamics, shape, deformation, and orientation of red blood cells in microcirculation affect the rheology, flow resistance and transport properties of whole blood. This leads to important correlations of cellular and continuum scales. Furthermore, the dynamics of RBCs subject to different flow conditions and vessel geometries is relevant for both fundamental research and biomedical applications (e.g drug delivery). In this thesis, the behaviour of RBCs is investigated for different flow conditions via computer simulations. We use a combination of two mesoscopic particle-based simulation techniques, dissipative particle dynamics and smoothed dissipative particle dynamics. We focus on the microcapillary scale of several μm. At this scale, blood cannot be considered at the continuum but has to be studied at the cellular level. The connection between cellular motion and overall blood rheology will be investigated. Red blood cells are modelled as viscoelastic objects interacting hydrodynamically with a viscous fluid environment. The properties of the membrane, such as resistance against bending or shearing, are set to correspond to experimental values. Furthermore, thermal fluctuations are considered via random forces. Analyses corresponding to light scattering measurements are performed in order to compare to experiments and suggest for which situations this method is suitable. Static light scattering by red blood cells characterises their shape and allows comparison to objects such as spheres or cylinders, whose scattering signals have analytical solutions, in contrast to those of red blood cells. Dynamic light scattering by red blood cells is studied concerning its suitability to detect and analyse motion, deformation and membrane fluctuations. Dynamic light scattering analysis is performed for both diffusing and flowing cells. We find that scattering signals depend on various cell properties, thus allowing to distinguish different cells. The scattering of diffusing cells allows to draw conclusions on their bending rigidity via the effective diffusion coefficient. The scattering of flowing cells allows to draw conclusions on the shear rate via the scattering amplitude correlation. In flow, a RBC shows different shapes and dynamic states, depending on conditions such as confinement, physiological/pathological state and cell age. Here, two essential flow conditions are studied: simple shear flow and tube flow. Simple shear flow as a basic flow condition is part of any more complex flow. The velocity profile is linear and shear stress is homogeneous. In simple shear flow, we find a sequence of different cell shapes by increasing the shear rate. With increasing shear rate, we find rolling cells with cup shapes, trilobe shapes and quadrulobe shapes. This agrees with recent experiments. Furthermore, the impact of the initial orientation on the dynamics is studied. To study crowding and collective effects, systems with higher haematocrit are set up. Tube flow is an idealised model for the flow through cylindric microvessels. Without cell, a parabolic flow profile prevails. A single red blood cell is placed into the tube and subject to a Poiseuille profile. In tube flow, we find different cell shapes and dynamics depending on confinement, shear rate and cell properties. For strong confinements and high shear rates, we find parachute-like shapes. Although not perfectly symmetric, they are adjusted to the flow profile and maintain a stationary shape and orientation. For weak confinements and low shear rates, we find tumbling slippers that rotate and moderately change their shape. For weak confinements and high shear rates, we find tank-treading slippers that oscillate in a limited range of inclination angles and strongly change their shape. For the lowest shear rates, we find cells performing a snaking motion. Due to cell properties and resultant deformations, all shapes differ from hitherto descriptions, such as steady tank-treading or symmetric parachutes. We introduce phase diagrams to identify flow regimes for the different shapes and dynamics. Changing cell properties, the regime borders in the phase diagrams change. In both flow types, both the viscosity contrast and the choice of stress-free shape are important. For in vitro experiments, the solvent viscosity has often been higher than the cytosol viscosity, leading to a different pattern of dynamics, such as steady tank-treading. The stress-free state of a RBC, which is the state at zero shear stress, is still controversial, and computer simulations enable direct comparisons of possible candidates in equivalent flow conditions.

Thermodynamics and Heat Transport of doped Spin-Ice Materials

In this thesis the low-temperature magnetism of the spin-ice systems Dy2Ti2O7 and Ho2Ti2O7 is investigated. In general, a clear experimental evidence for a sizable magnetic contribution kappa_{mag} to the low-temperature, zero-field heat transport of both spin-ice materials is observed. This kappa_{mag} can be attributed to the magnetic monopole excitations, which are highly mobile in zero field and are suppressed by a rather small external field resulting in a drop of kappa(H). Towards higher magnetic fields, significant field dependencies of the phononic heat conductivities kappa_{ph}(H) of Ho2Ti2O7 and Dy2Ti2O7 are found, which are, however, of opposite signs, as it is also found for the highly dilute reference materials (Ho0.5Y0.5)2Ti2O7 and (Dy0.5Y0.5)2Ti2O7. The dominant effect in the Ho-based materials is the scattering of phonons by spin flips which appears to be significantly stronger than in the Dy-based materials. Here, the thermal conductivity is suppressed due to enhanced lattice distortions observed in the magnetostriction. Furthermore, the thermal conductivity of Dy2Ti2O7 has been investigated concerning strong hysteresis effects and slow-relaxation processes towards equilibrium states in the low-temperature and low-field regime. The thermal conductivity in the hysteretic regions slowly relaxes towards larger values suggesting that there is an additional suppression of the heat transport by disorder in the non-equilibrium states. The equilibration can even be governed by the heat current for particular configurations. A special focus was put on the dilution series Dy2Ti2O7x. From specific heat measurements, it was found that the ultra-slow thermal equilibration in pure spin ice Dy2Ti2O7 is rapidly suppressed upon dilution with non-magnetic yttrium and vanishes completely for x>=0.2 down to the lowest accessible temperatures. In general, the low-temperature entropy of (Dy1-xYx)2Ti2O7, considerably decreases with increasing x, whereas its temperature-dependence drastically increases. Thus, it could be clarified that there is no experimental evidence for a finite zero-temperature entropy in (Dy1-xYx)2Ti2O7 above x>=0.2, in clear contrast to the finite residual entropy S_{P}(x) expected from a generalized Pauling approximation. A similar discrepancy is also present between S_{P}(x) and the low-temperature entropy obtained by Monte Carlo simulations, which reproduce the experimental data from 25 K down to 0.7 K, whereas the data at 0.4 K are overestimated. A straightforward description of the field-dependence kappa(H) of the dilution series with qualitative models justifies the extraction of kappa_{mag}. It was observed that kappa_{mag} systematically scales with the degree of dilution and its low-field decrease is related to the monopole excitation energy. The diffusion coefficient D_{mag} for the monopole excitations was calculated by means of c_{mag} and kappa_{mag}. It exhibits a broad maximum around 1.6 K and is suppressed for T<=0.5 K, indicating a non-degenerate ground state in the long-time limit, and in the high-temperature range for T>=4 K where spin-ice physics is eliminated. A mean-free path of 0.3 mum is obtained for Dy2Ti2O7 at about 1 K within the kinetic gas theory.

Intravoxel Incoherent Motion Metrics as Potential Biomarkers for Survival in Glioblastoma.

OBJECTIVE: Intravoxel incoherent motion (IVIM) is an MRI technique with potential applications in measuring brain tumor perfusion, but its clinical impact remains to be determined. We assessed the usefulness of IVIM-metrics in predicting survival in newly diagnosed glioblastoma. METHODS: Fifteen patients with glioblastoma underwent MRI including spin-echo echo-planar DWI using 13 b-values ranging from 0 to 1000 s/mm2. Parametric maps for diffusion coefficient (D), pseudodiffusion coefficient (D*), and perfusion fraction (f) were generated for contrast-enhancing regions (CER) and non-enhancing regions (NCER). Regions of interest were manually drawn in regions of maximum f and on the corresponding dynamic susceptibility contrast images. Prognostic factors were evaluated by Kaplan-Meier survival and Cox proportional hazards analyses. RESULTS: We found that fCER and D*CER correlated with rCBFCER. The best cutoffs for 6-month survival were fCER>9.86% and D*CER>21.712 x10-3mm2/s (100% sensitivity, 71.4% specificity, 100% and 80% positive predictive values, and 80% and 100% negative predictive values; AUC:0.893 and 0.857, respectively). Treatment yielded the highest hazard ratio (5.484; 95% CI: 1.162-25.88; AUC: 0.723; P = 0.031); fCER combined with treatment predicted survival with 100% accuracy. CONCLUSIONS: The IVIM-metrics fCER and D*CER are promising biomarkers of 6-month survival in newly diagnosed glioblastoma.

Existence and wandering of bumps in a spiking neural network model

We study spatially localized states of a spiking neuronal network populated by a pulse coupled phase oscillator known as the lighthouse model. We show that in the limit of slow synaptic interactions in the continuum limit the dynamics reduce to those of the standard Amari model. For non-slow synaptic connections we are able to go beyond the standard firing rate analysis of localized solutions allowing us to explicitly construct a family of co-existing one-bump solutions, and then track bump width and firing pattern as a function of system parameters. We also present an analysis of the model on a discrete lattice. We show that multiple width bump states can co-exist and uncover a mechanism for bump wandering linked to the speed of synaptic processing. Moreover, beyond a wandering transition point we show that the bump undergoes an effective random walk with a diffusion coefficient that scales exponentially with the rate of synaptic processing and linearly with the lattice spacing.

Solução semi-analítica para modelagem das condições inicial e de contorno aplicáveis ao ensaio de difusão pura

Dissertação (mestrado)—Universidade de Brasília, Faculdade de Tecnologia, Departamento de Engenharia Civil e Ambiental, 2015.

High protein flexibility and reduced hydration water dynamics are key pressure adaptive strategies in prokaryotes

Water and protein dynamics on a nanometer scale were measured by quasi-elastic neutron scattering in the piezophile archaeon Thermococcus barophilus and the closely related pressure-sensitive Thermococcus kodakarensis, at 0.1 and 40 MPa. We show that cells of the pressure sensitive organism exhibit higher intrinsic stability. Both the hydration water dynamics and the fast protein and lipid dynamics are reduced under pressure. In contrast, the proteome of T. barophilus is more pressure sensitive than that of T. kodakarensis. The diffusion coefficient of hydration water is reduced, while the fast protein and lipid dynamics are slightly enhanced with increasing pressure. These findings show that the coupling between hydration water and cellular constituents might not be simply a master-slave relationship. We propose that the high flexibility of the T. barophilus proteome associated with reduced hydration water may be the keys to the molecular adaptation of the cells to high hydrostatic pressure.