BEYOND THE “LEAST LIMITING WATER RANGE”: RETHINKING SOIL PHYSICS RESEARCH IN BRAZIL

As opposed to objective definitions in soil physics, the subjective term “soil physical quality” is increasingly found in publications in the soil physics area. A supposed indicator of soil physical quality that has been the focus of attention, especially in the Brazilian literature, is the Least Limiting Water Range (RLL), translated in Portuguese as "Intervalo Hídrico Ótimo" or IHO. In this paper the four limiting water contents that define RLLare discussed in the light of objectively determinable soil physical properties, pointing to inconsistencies in the RLLdefinition and calculation. It also discusses the interpretation of RLL as an indicator of crop productivity or soil physical quality, showing its inability to consider common phenological and pedological boundary conditions. It is shown that so-called “critical densities” found by the RLL through a commonly applied calculation method are questionable. Considering the availability of robust models for agronomy, ecology, hydrology, meteorology and other related areas, the attractiveness of RLL as an indicator to Brazilian soil physicists is not related to its (never proven) effectiveness, but rather to the simplicity with which it is dealt. Determining the respective limiting contents in a simplified manner, relegating the study or concern on the actual functioning of the system to a lower priority, goes against scientific construction and systemic understanding. This study suggests a realignment of the research in soil physics in Brazil with scientific precepts, towards mechanistic soil physics, to replace the currently predominant search for empirical correlations below the state of the art of soil physics.

Random mixtures with orientational order, and the anisotropic resistivity tensor of high Tc superconductors

By generalizing effective-medium theory to the case of orientationally ordered but positionally disordered two component mixtures, it is shown that the anisotropic dielectric tensor of oxide superconductors can be extracted from microwave measurements on oriented crystallites of YBa2Cu3O7¿x embedded in epoxy. Surprisingly, this technique appears to be the only one which can access the resistivity perpendicular to the copper¿oxide planes in crystallites that are too small for depositing electrodes. This possibility arises in part because the real part of the dielectric constant of oxide superconductors has a large magnitude. The validity of the effective-medium approach for orientationally ordered mixtures is corroborated by simulations on two¿dimensional anisotropic random resistor networks. Analysis of the experimental data suggests that the zero-temperature limit of the finite frequency resistivity does not vanish along the c axis, a result which would simply the existence of states at the Fermi surface, even in the superconducting state

Finite-size scaling investigation of the liquid-liquid critical point in ST2 water and its stability with respect to crystallization.

The liquid-liquid critical point scenario of water hypothesizes the existence of two metastable liq- uid phases low-density liquid (LDL) and high-density liquid (HDL) deep within the supercooled region. The hypothesis originates from computer simulations of the ST2 water model, but the stabil- ity of the LDL phase with respect to the crystal is still being debated. We simulate supercooled ST2 water at constant pressure, constant temperature, and constant number of molecules N for N ≤ 729 and times up to 1 μs. We observe clear differences between the two liquids, both structural and dynamical. Using several methods, including finite-size scaling, we confirm the presence of a liquid-liquid phase transition ending in a critical point. We find that the LDL is stable with respect to the crystal in 98% of our runs (we perform 372 runs for LDL or LDL-like states), and in 100% of our runs for the two largest system sizes (N = 512 and 729, for which we perform 136 runs for LDL or LDL-like states). In all these runs, tiny crystallites grow and then melt within 1 μs. Only for N ≤ 343 we observe six events (over 236 runs for LDL or LDL-like states) of spontaneous crystal- lization after crystallites reach an estimated critical size of about 70 ± 10 molecules.

Experimental Method to Determine Some Physical Properties in Physics Classes

ABSTRACT Particle density, gravimetric and volumetric water contents and porosity are important basic concepts to characterize porous systems such as soils. This paper presents a proposal of an experimental method to measure these physical properties, applicable in experimental physics classes, in porous media samples consisting of spheres with the same diameter (monodisperse medium) and with different diameters (polydisperse medium). Soil samples are not used given the difficulty of working with this porous medium in laboratories dedicated to teaching basic experimental physics. The paper describes the method to be followed and results of two case studies, one in monodisperse medium and the other in polydisperse medium. The particle density results were very close to theoretical values for lead spheres, whose relative deviation (RD) was -2.9 % and +0.1 % RD for the iron spheres. The RD of porosity was also low: -3.6 % for lead spheres and -1.2 % for iron spheres, in the comparison of procedures – using particle and porous medium densities and saturated volumetric water content – and monodisperse and polydisperse media.

Diffusive transport properties in monovalent and divalent metal-ion halide melts: A computer simulation study

Self- and cross-velocity correlation functions and related transport coefficients of molten salts are studied by molecular-dynamics simulation. Six representative systems are considered, i.e., NaCl and KCl alkali halides, CuCl and CuBr noble-metal halides, and SrCl2 and ZnCl2 divalent metal-ion halides. Computer simulation results are compared with experimental self-diffusion coefficients and electrical conductivities. Special attention is paid to dynamic cross correlations and their dependence on the Coulomb interactions as well as on the size and mass differences between anions and cations.

Monte Carlo simulation of a whole-body counter using IGOR phantoms.

Whole-body counting is a technique of choice for assessing the intake of gamma-emitting radionuclides. An appropriate calibration is necessary, which is done either by experimental measurement or by Monte Carlo (MC) calculation. The aim of this work was to validate a MC model for calibrating whole-body counters (WBCs) by comparing the results of computations with measurements performed on an anthropomorphic phantom and to investigate the effect of a change in phantom's position on the WBC counting sensitivity. GEANT MC code was used for the calculations, and an IGOR phantom loaded with several types of radionuclides was used for the experimental measurements. The results show a reasonable agreement between measurements and MC computation. A 1-cm error in phantom positioning changes the activity estimation by >2%. Considering that a 5-cm deviation of the positioning of the phantom may occur in a realistic counting scenario, this implies that the uncertainty of the activity measured by a WBC is ∼10-20%.

Evaporation and coarsening dynamics with open boundaries

We present a study of the evaporation dynamics of a substance undergoing a coarsening process. The system is modeled by the Cahn-Hilliard equation with absorbing boundaries. We have found that the dynamics, although of a diffusive nature, is much slower than the usual one without coarsening. Analytical and simulation results are in reasonable agreement.

Stochastic generation of homogeneous isotropic turbulence with well-defined-spectra

A precise and simple computational model to generate well-behaved two-dimensional turbulent flows is presented. The whole approach rests on the use of stochastic differential equations and is general enough to reproduce a variety of energy spectra and spatiotemporal correlation functions. Analytical expressions for both the continuous and the discrete versions, together with simulation algorithms, are derived. Results for two relevant spectra, covering distinct ranges of wave numbers, are given.

Morphology and segregation in two-component diffusion-limited aggregation

A diffusion-limited-aggregation (DLA) model with two components (A and B species) is presented to investigate the structure of the composite deposits. The sticking probability PAB (=PBA) between the different species is introduced into the original DLA model. By using computer simulation it is shown that various patterns are produced with varying the sticking probabilities PAB (=PBA) and PAA (= PBB), where PAA (=PBB) is the sticking probability between the same species. Segregated patterns can be analyzed under the condition PAB < PAA, assumed throughout the paper. With decreasing sticking probability PAB, a clustering of the same species occurs. With sufficiently small values of both sticking probabilities PAB and PAA, the deposit becomes dense and the segregated patterns of the composite deposit show a striped structure. The effect of the concentration on the pattern morphology is also shown.

Spatiotemporal order out of noise

Experimental observations of self-organized behavior arising out of noise are also described, and details on the numerical algorithms needed in the computer simulation of these problems are given.

A model of cellular dosimetry for macroscopic tumors in radiopharmaceutical therapy.

PURPOSE: In the radiopharmaceutical therapy approach to the fight against cancer, in particular when it comes to translating laboratory results to the clinical setting, modeling has served as an invaluable tool for guidance and for understanding the processes operating at the cellular level and how these relate to macroscopic observables. Tumor control probability (TCP) is the dosimetric end point quantity of choice which relates to experimental and clinical data: it requires knowledge of individual cellular absorbed doses since it depends on the assessment of the treatment's ability to kill each and every cell. Macroscopic tumors, seen in both clinical and experimental studies, contain too many cells to be modeled individually in Monte Carlo simulation; yet, in particular for low ratios of decays to cells, a cell-based model that does not smooth away statistical considerations associated with low activity is a necessity. The authors present here an adaptation of the simple sphere-based model from which cellular level dosimetry for macroscopic tumors and their end point quantities, such as TCP, may be extrapolated more reliably. METHODS: Ten homogenous spheres representing tumors of different sizes were constructed in GEANT4. The radionuclide 131I was randomly allowed to decay for each model size and for seven different ratios of number of decays to number of cells, N(r): 1000, 500, 200, 100, 50, 20, and 10 decays per cell. The deposited energy was collected in radial bins and divided by the bin mass to obtain the average bin absorbed dose. To simulate a cellular model, the number of cells present in each bin was calculated and an absorbed dose attributed to each cell equal to the bin average absorbed dose with a randomly determined adjustment based on a Gaussian probability distribution with a width equal to the statistical uncertainty consistent with the ratio of decays to cells, i.e., equal to Nr-1/2. From dose volume histograms the surviving fraction of cells, equivalent uniform dose (EUD), and TCP for the different scenarios were calculated. Comparably sized spherical models containing individual spherical cells (15 microm diameter) in hexagonal lattices were constructed, and Monte Carlo simulations were executed for all the same previous scenarios. The dosimetric quantities were calculated and compared to the adjusted simple sphere model results. The model was then applied to the Bortezomib-induced enzyme-targeted radiotherapy (BETR) strategy of targeting Epstein-Barr virus (EBV)-expressing cancers. RESULTS: The TCP values were comparable to within 2% between the adjusted simple sphere and full cellular models. Additionally, models were generated for a nonuniform distribution of activity, and results were compared between the adjusted spherical and cellular models with similar comparability. The TCP values from the experimental macroscopic tumor results were consistent with the experimental observations for BETR-treated 1 g EBV-expressing lymphoma tumors in mice. CONCLUSIONS: The adjusted spherical model presented here provides more accurate TCP values than simple spheres, on par with full cellular Monte Carlo simulations while maintaining the simplicity of the simple sphere model. This model provides a basis for complementing and understanding laboratory and clinical results pertaining to radiopharmaceutical therapy.

Simulation of interface dynamics across pore discontinuities

We study the dynamics of a water-oil meniscus moving from a smaller to a larger pore. The process is characterised by an abrupt change in the configuration, yielding a sudden energy release. A theoretic study for static conditions provides analytical solutions of the surface energy content of the system. Although the configuration after the sudden energy release is energetically more convenient, an energy barrier must be overcome before the process can happen spontaneously. The energy barrier depends on the system geometry and on the flow parameters. The analytical results are compared to numerical simulations that solve the full Navier-Stokes equation in the pore space and employ the Volume Of Fluid (VOF) method to track the evolution of the interface. First, the numerical simulations of a quasi-static process are validated by comparison with the analytical solutions for a static meniscus, then numerical simulations with varying injection velocity are used to investigate dynamic effects on the configuration change. During the sudden energy jump the system exhibits an oscillatory behaviour. Extension to more complex geometries might elucidate the mechanisms leading to a dynamic capillary pressure and to bifurcations in final distributions of fluid phases in porous

Ion-beam induced oxidation of GaAs and AlGaAs

The oxidation of GaAs and AlxGa1−xAs targets by oxygen irradiation has been studied in detail. It was found that the oxidation process is characterized by the strong preferential oxidation of Al as compared to Ga, and of Ga as compared to As. This experimental observation, which has been accurately quantified by using x‐ray photoelectron spectroscopy, is connected to the different heats of formation of the corresponding oxides. The oxide grown by ion beam oxidation shows a strong depletion in As and relatively low oxidation of As as well. The depletion can be associated with the preferential sputtering of the As oxide in respect to other compounds whereas the low oxidation is due to the low heat of formation. In contrast Al is rapidly and fully oxidized, turning the outermost layer of the altered layer to a single Al2O3 overlayer, as observed by transmission electron microscopy. The radiation enhanced diffusion of oxygen and aluminum in the altered layer explains the large thickness of these altered layers and the formation of Al oxides on top of the layers. For the case of ion‐beam oxidation of GaAs a simulation program has been developed which describes adequately the various growth mechanisms experimentally observed