Microscopic study of He2-SF6 trimers

The He2-SF6 trimers, in their different He isotopic combinations, are studied in the framework of both the correlated Jastrow approach and the correlated hyperspherical harmonics (CHH) expansion method. The energetics and structure of the He-SF6 dimers are analyzed, and the existence of a characteristic rotational band in the excitation spectrum is discussed, as well as the isotopic differences. The binding energies and the spatial properties of the trimers, in their ground and lowest lying excited states, obtained by the Jastrow ansatz are in excellent agreement with the results of the converged CHH expansion. The introduction of the He-He correlation makes all trimers bound by largely suppressing the short range He-He repulsion. The structural properties of the trimers are qualitatively explained in terms of the shape of the interactions, Pauli principle, and masses of the constituents.

Microscopic description of the twist mode in normal and superfluid trapped Fermi gases

We investigate the "twist" mode (rotation of the upper against the lower hemisphere) of a dilute atomic Fermi gas in a spherical trap. The normal and superfluid phases are considered. The linear response to this external perturbation is calculated within the microscopic Hartree-Fock-Bogoliubov approach. In the normal phase the excitation spectrum is concentrated in a rather narrow peak very close to the trapping frequency. In the superfluid phase the strength starts to be damped and fragmented and the collectivity of the mode is progressively lost when the temperature decreases. In the weak-pairing regime some reminiscence of the collective motion still exists, whereas in the strong-pairing regime the twist mode is completely washed out. The disappearance of the twist mode in the strong-pairing regime with decreasing temperature is interpreted in the framework of the two-fluid model.

Sex differences in urinary levels of several biological indicators of exposure: a simulation study using a compartmental based toxicokinetic model

Toxicokinetic modeling is a useful tool to describe or predict the behavior of a chemical agent in the human or animal organism. A general model based on four compartments was developed in a previous study in order to quantify the effect of human variability on a wide range of biological exposure indicators. The aim of this study was to adapt this existing general toxicokinetic model to three organic solvents, which were methyl ethyl ketone, 1-methoxy-2-propanol and 1,1,1,-trichloroethane, and to take into account sex differences. We assessed in a previous human volunteer study the impact of sex on different biomarkers of exposure corresponding to the three organic solvents mentioned above. Results from that study suggested that not only physiological differences between men and women but also differences due to sex hormones levels could influence the toxicokinetics of the solvents. In fact the use of hormonal contraceptive had an effect on the urinary levels of several biomarkers, suggesting that exogenous sex hormones could influence CYP2E1 enzyme activity. These experimental data were used to calibrate the toxicokinetic models developed in this study. Our results showed that it was possible to use an existing general toxicokinetic model for other compounds. In fact, most of the simulation results showed good agreement with the experimental data obtained for the studied solvents, with a percentage of model predictions that lies within the 95% confidence interval varying from 44.4 to 90%. Results pointed out that for same exposure conditions, men and women can show important differences in urinary levels of biological indicators of exposure. Moreover, when running the models by simulating industrial working conditions, these differences could even be more pronounced. In conclusion, a general and simple toxicokinetic model, adapted for three well known organic solvents, allowed us to show that metabolic parameters can have an important impact on the urinary levels of the corresponding biomarkers. These observations give evidence of an interindividual variablity, an aspect that should have its place in the approaches for setting limits of occupational exposure.

Numerical simulation of gel electrophoresis of DNA knots in weak and strong electric fields.

Gel electrophoresis allows one to separate knotted DNA (nicked circular) of equal length according to the knot type. At low electric fields, complex knots, being more compact, drift faster than simpler knots. Recent experiments have shown that the drift velocity dependence on the knot type is inverted when changing from low to high electric fields. We present a computer simulation on a lattice of a closed, knotted, charged DNA chain drifting in an external electric field in a topologically restricted medium. Using a Monte Carlo algorithm, the dependence of the electrophoretic migration of the DNA molecules on the knot type and on the electric field intensity is investigated. The results are in qualitative and quantitative agreement with electrophoretic experiments done under conditions of low and high electric fields.

A simple method for the simulation of multiple elastic scattering of electrons

A screened Rutherford cross section is modified by means of a correction factor to obtain the proper transport cross section computed by partial¿wave analysis. The correction factor is tabulated for electron energies in the range 0¿100 keV and for elements in the range from Z=4 to 82. The modified screened Rutherford cross section is shown to be useful as an approximation for the simulation of plural and multiple scattering. Its performance and limitations are exemplified for electrons scattered in Al and Au.

Monte Carlo simulation of kilovolt electron transport in solids

A Monte Carlo procedure to simulate the penetration and energy loss of low¿energy electron beams through solids is presented. Elastic collisions are described by using the method of partial waves for the screened Coulomb field of the nucleus. The atomic charge density is approximated by an analytical expression with parameters determined from the Dirac¿Hartree¿Fock¿Slater self¿consistent density obtained under Wigner¿Seitz boundary conditions in order to account for solid¿state effects; exchange effects are also accounted for by an energy¿dependent local correction. Elastic differential cross sections are then easily computed by combining the WKB and Born approximations to evaluate the phase shifts. Inelastic collisions are treated on the basis of a generalized oscillator strength model which gives inelastic mean free paths and stopping powers in good agreement with experimental data. This scattering model is accurate in the energy range from a few hundred eV up to about 50 keV. The reliability of the simulation method is analyzed by comparing simulation results and experimental data from backscattering and transmission measurements.

Simple method for the simulation of multiple elastic scattering of electrons

A screened Rutherford cross section is modified by means of a correction factor to obtain the proper transport cross section computed by partial¿wave analysis. The correction factor is tabulated for electron energies in the range 0¿100 keV and for elements in the range from Z=4 to 82. The modified screened Rutherford cross section is shown to be useful as an approximation for the simulation of plural and multiple scattering. Its performance and limitations are exemplified for electrons scattered in Al and Au.

Prediction of Infant Drug Exposure Through Breastfeeding: Population PK Modeling and Simulation of Fluoxetine Exposure.

The likelihood of significant exposure to drugs in infants through breast milk is poorly defined, given the difficulties of conducting pharmacokinetics (PK) studies. Using fluoxetine (FX) as an example, we conducted a proof-of-principle study applying population PK (popPK) modeling and simulation to estimate drug exposure in infants through breast milk. We simulated data for 1,000 mother-infant pairs, assuming conservatively that the FX clearance in an infant is 20% of the allometrically adjusted value in adults. The model-generated estimate of the milk-to-plasma ratio for FX (mean: 0.59) was consistent with those reported in other studies. The median infant-to-mother ratio of FX steady-state plasma concentrations predicted by the simulation was 8.5%. Although the disposition of the active metabolite, norfluoxetine, could not be modeled, popPK-informed simulation may be valid for other drugs, particularly those without active metabolites, thereby providing a practical alternative to conventional PK studies for exposure risk assessment in this population.

Microscopic calculations of the excitation spectrum of one 3He impurity in liquid 4He

We calculate the chemical potential ¿0 and the effective mass m*/m3 of one 3He impurity in liquid 4He. First a variational wave function including two- and three-particle dynamical correlations is adopted. Triplet correlations bring the computed values of ¿0 very close to the experimental results. The variational estimate of m*/m3 includes also backflow correlations between the 3He atom and the particles in the medium. Different approximations for the three-particle distribution function give almost the same values for m*/m3. The variational approach underestimates m*/m3 by ~10% at all of the considered densities. Correlated-basis perturbation theory is then used to improve the wave function to include backflow around the particles of the medium. The perturbative series built up with one-phonon states only is summed up to infinite order and gives results very close to the variational ones. All the perturbative diagrams with two independent phonons have then been summed to compute m*/m3. Their contribution depends to some extent on the form used for the three-particle distribution function. When the scaling approximation is adopted, a reasonable agreement with the experimental results is achieved.

Seismic attenuation due to wave-induced fluid flow at microscopic and macroscopic scales

Wave-induced fluid flow at microscopic and mesoscopic scales arguably constitutes the major cause of intrinsic seismic attenuation throughout the exploration seismic and sonic frequency ranges. The quantitative analysis of these phenomena is, however, complicated by the fact that the governing physical processes may be dependent. The reason for this is that the presence of microscopic heterogeneities, such as micro-cracks or broken grain contacts, causes the stiffness of the so-called modified dry frame to be complex-valued and frequency-dependent, which in turn may affect the viscoelastic behaviour in response to fluid flow at mesoscopic scales. In this work, we propose a simple but effective procedure to estimate the seismic attenuation and velocity dispersion behaviour associated with wave-induced fluid flow due to both microscopic and mesoscopic heterogeneities and discuss the results obtained for a range of pertinent scenarios.

Endocytic and secretory traffic in Arabidopsis merge in the trans-Golgi network/early endosome, an independent and highly dynamic organelle.

Plants constantly adjust their repertoire of plasma membrane proteins that mediates transduction of environmental and developmental signals as well as transport of ions, nutrients, and hormones. The importance of regulated secretory and endocytic trafficking is becoming increasingly clear; however, our knowledge of the compartments and molecular machinery involved is still fragmentary. We used immunogold electron microscopy and confocal laser scanning microscopy to trace the route of cargo molecules, including the BRASSINOSTEROID INSENSITIVE1 receptor and the REQUIRES HIGH BORON1 boron exporter, throughout the plant endomembrane system. Our results provide evidence that both endocytic and secretory cargo pass through the trans-Golgi network/early endosome (TGN/EE) and demonstrate that cargo in late endosomes/multivesicular bodies is destined for vacuolar degradation. Moreover, using spinning disc microscopy, we show that TGN/EEs move independently and are only transiently associated with an individual Golgi stack.

Effects of traffic control on the soil physical quality and the cultivation of sugarcane

The cultivation of sugarcane with intensive use of machinery, especially for harvest, induces soil compaction, affecting the crop development. The control of agricultural traffic is an alternative of management in the sector, with a view to preserve the soil physical quality, resulting in increased sugarcane root growth, productivity and technological quality. The objective of this study was to evaluate the physical quality of an Oxisol with and without control traffic and the resulting effects on sugarcane root development, productivity and technological quality. The following managements were tested: no traffic control (NTC), traffic control consisting of an adjustment of the track width of the tractor and sugarcane trailer (TC1) and traffic control consisting of an adjustment of the track width of the tractor and trailer and use of an autopilot (TC2). Soil samples were collected (layers 0.00-0.10; 0.10-0.20 and 0.20-0.30 m) in the plant rows, inter-row center and seedbed region, 0.30 m away from the plant row. The productivity was measured with a specific weighing scale. The technological variables of sugarcane were measured in each plot. Soil cores were collected to analyze the root system. In TC2, the soil bulk density and compaction degree were lowest and total porosity and macroporosity highest in the plant row. Soil penetration resistance in the plant row, was less than 2 MPa in TC1 and TC2. Soil aggregation and total organic carbon did not differ between the management systems. The root surface and volume were increased in TC1 and TC2, with higher productivity and sugar yield than under NTC. The sugarcane variables did not differ between the managements. The soil physical quality in the plant row was preserved under management TC1 and TC2, with an improved root development and increases of 18.72 and 20.29 % in productivity and sugar yield, respectively.

Accurate and efficient simulation of multiphase flow in a heterogeneous reservoir by using error estimate and control in the multiscale finite-volume framework

The multiscale finite-volume (MSFV) method is designed to reduce the computational cost of elliptic and parabolic problems with highly heterogeneous anisotropic coefficients. The reduction is achieved by splitting the original global problem into a set of local problems (with approximate local boundary conditions) coupled by a coarse global problem. It has been shown recently that the numerical errors in MSFV results can be reduced systematically with an iterative procedure that provides a conservative velocity field after any iteration step. The iterative MSFV (i-MSFV) method can be obtained with an improved (smoothed) multiscale solution to enhance the localization conditions, with a Krylov subspace method [e.g., the generalized-minimal-residual (GMRES) algorithm] preconditioned by the MSFV system, or with a combination of both. In a multiphase-flow system, a balance between accuracy and computational efficiency should be achieved by finding a minimum number of i-MSFV iterations (on pressure), which is necessary to achieve the desired accuracy in the saturation solution. In this work, we extend the i-MSFV method to sequential implicit simulation of time-dependent problems. To control the error of the coupled saturation/pressure system, we analyze the transport error caused by an approximate velocity field. We then propose an error-control strategy on the basis of the residual of the pressure equation. At the beginning of simulation, the pressure solution is iterated until a specified accuracy is achieved. To minimize the number of iterations in a multiphase-flow problem, the solution at the previous timestep is used to improve the localization assumption at the current timestep. Additional iterations are used only when the residual becomes larger than a specified threshold value. Numerical results show that only a few iterations on average are necessary to improve the MSFV results significantly, even for very challenging problems. Therefore, the proposed adaptive strategy yields efficient and accurate simulation of multiphase flow in heterogeneous porous media.