Local and global error models to improve uncertainty quantification

In groundwater applications, Monte Carlo methods are employed to model the uncertainty on geological parameters. However, their brute-force application becomes computationally prohibitive for highly detailed geological descriptions, complex physical processes, and a large number of realizations. The Distance Kernel Method (DKM) overcomes this issue by clustering the realizations in a multidimensional space based on the flow responses obtained by means of an approximate (computationally cheaper) model; then, the uncertainty is estimated from the exact responses that are computed only for one representative realization per cluster (the medoid). Usually, DKM is employed to decrease the size of the sample of realizations that are considered to estimate the uncertainty. We propose to use the information from the approximate responses for uncertainty quantification. The subset of exact solutions provided by DKM is then employed to construct an error model and correct the potential bias of the approximate model. Two error models are devised that both employ the difference between approximate and exact medoid solutions, but differ in the way medoid errors are interpolated to correct the whole set of realizations. The Local Error Model rests upon the clustering defined by DKM and can be seen as a natural way to account for intra-cluster variability; the Global Error Model employs a linear interpolation of all medoid errors regardless of the cluster to which the single realization belongs. These error models are evaluated for an idealized pollution problem in which the uncertainty of the breakthrough curve needs to be estimated. For this numerical test case, we demonstrate that the error models improve the uncertainty quantification provided by the DKM algorithm and are effective in correcting the bias of the estimate computed solely from the MsFV results. The framework presented here is not specific to the methods considered and can be applied to other combinations of approximate models and techniques to select a subset of realizations

Impact of incorrect assumptions on the covariance structure of random effects and/or residuals in nonlinear mixed models for repeated measures data

In this paper we analyse, using Monte Carlo simulation, the possible consequences of incorrect assumptions on the true structure of the random effects covariance matrix and the true correlation pattern of residuals, over the performance of an estimation method for nonlinear mixed models. The procedure under study is the well known linearization method due to Lindstrom and Bates (1990), implemented in the nlme library of S-Plus and R. Its performance is studied in terms of bias, mean square error (MSE), and true coverage of the associated asymptotic confidence intervals. Ignoring other criteria like the convenience of avoiding over parameterised models, it seems worst to erroneously assume some structure than do not assume any structure when this would be adequate.

Vacancy-driven ordering in a two-dimensional binary alloy

Domain growth in a system with nonconserved order parameter is studied. We simulate the usual Ising model for binary alloys with concentration 0.5 on a two-dimensional square lattice by Monte Carlo techniques. Measurements of the energy, jump-acceptance ratio, and order parameters are performed. Dynamics based on the diffusion of a single vacancy in the system gives a growth law faster than the usual Allen-Cahn law. Allowing vacancy jumps to next-nearest-neighbor sites is essential to prevent vacancy trapping in the ordered regions. By measuring local order parameters we show that the vacancy prefers to be in the disordered regions (domain boundaries). This naturally concentrates the atomic jumps in the domain boundaries, accelerating the growth compared with the usual exchange mechanism that causes jumps to be homogeneously distributed on the lattice.

Recommendations to reduce hand exposure for standard nuclear medicine procedures

The optimization of the extremity dosimetry of medical staff in nuclear medicine was the aim of the Work Package 4 (WP4) of the ORAMED project, a Collaborative Project (2008-2011) supported by the European Commission within its 7th Framework Programme. Hand doses and dose distributions across the hands of medical staff working in nuclear medicine departments were evaluated through an extensive measurement program involving 32 hospitals in Europe and 139 monitored workers. The study included the most frequently used radionuclides, (99m)Tc- and (18)F-labelled radiopharmaceuticals for diagnostic and (90)Y-labelled Zevalin (R) and DOTATOC for therapy. Furthermore, Monte Carlo simulations were performed in different predefined scenarios to evaluate separately the efficacy of different radiation protection measures by comparing hand dose distributions according to various parameters. The present work gives recommendations based on results obtained with both measurements and simulations. This results in nine practical recommendations regarding the positioning of the dosemeters for an appropriate skin dose monitoring and the best protection means to reduce the personnel exposure.

Comment on "Kinetics of spinodal decomposition in the ising model with vacancy diffusion"

In a recent paper [Phys. Rev. B 50, 3477 (1994)], P. Fratzl and O. Penrose present the results of the Monte Carlo simulation of the spinodal decomposition problem (phase separation) using the vacancy dynamics mechanism. They observe that the t1/3 growth regime is reached faster than when using the standard Kawasaki dynamics. In this Comment we provide a simple explanation for the phenomenon based on the role of interface diffusion, which they claim is irrelevant for the observed behavior.

Magnetic phase separation in ordered alloys

We present a lattice model to study the equilibrium phase diagram of ordered alloys with one magnetic component that exhibits a low temperature phase separation between paramagnetic and ferromagnetic phases. The model is constructed from the experimental facts observed in Cu3-xAlMnx and it includes coupling between configurational and magnetic degrees of freedom that are appropriate for reproducing the low temperature miscibility gap. The essential ingredient for the occurrence of such a coexistence region is the development of ferromagnetic order induced by the long-range atomic order of the magnetic component. A comparative study of both mean-field and Monte Carlo solutions is presented. Moreover, the model may enable the study of the structure of ferromagnetic domains embedded in the nonmagnetic matrix. This is relevant in relation to phenomena such as magnetoresistance and paramagnetism

Two-dimensional clusters of liquid 4He

The binding energies of two-dimensional clusters (puddles) of¿4He are calculated in the framework of the diffusion Monte Carlo method. The results are well fitted by a mass formula in powers of x=N-1/2, where N is the number of particles. The analysis of the mass formula allows for the extraction of the line tension, which turns out to be 0.121 K/Å. Sizes and density profiles of the puddles are also reported.

Interfacial growth in driven diffusive systems

We have studied the growth of interfaces in driven diffusive systems well below the critical temperature by means of Monte Carlo simulations. We consider the region beyond the linear regime and of large values of the external field which has not been explored before. The simulations support the existence of interfacial traveling waves when asymmetry is introduced in the model, a result previously predicted by a linear-stability analysis. Furthermore, the generalization of the Gibbs-Thomson relation is discussed. The results provide evidence that the external field is a stabilizing effect which can be considered as effectively increasing the surface tension.

Simulations of electrophoretic collisions of DNA knots with gel obstacles

Gel electrophoresis can be used to separate nicked circular DNA molecules of equal length but forming different knot types. At low electric fields, complex knots drift faster than simpler knots. However, at high electric field the opposite is the case and simpler knots migrate faster than more complex knots. Using Monte Carlo simulations we investigate the reasons of this reversal of relative order of electrophoretic mobility of DNA molecules forming different knot types. We observe that at high electric fields the simulated knotted molecules tend to hang over the gel fibres and require passing over a substantial energy barrier to slip over the impeding gel fibre. At low electric field the interactions of drifting molecules with the gel fibres are weak and there are no significant energy barriers that oppose the detachment of knotted molecules from transverse gel fibres.

Quantum metric spaces as a model for pregeometry

A new arena for the dynamics of spacetime is proposed, in which the basic quantum variable is the two-point distance on a metric space. The scaling dimension (that is, the Kolmogorov capacity) in the neighborhood of each point then defines in a natural way a local concept of dimension. We study our model in the region of parameter space in which the resulting spacetime is not too different from a smooth manifold.

Energy and structure of dilute hard- and soft-sphere gases

The energy and structure of dilute hard- and soft-sphere Bose gases are systematically studied in the framework of several many-body approaches, such as the variational correlated theory, the Bogoliubov model, and the uniform limit approximation, valid in the weak-interaction regime. When possible, the results are compared with the exact diffusion Monte Carlo ones. Jastrow-type correlation provides a good description of the systems, both hard- and soft-spheres, if the hypernetted chain energy functional is freely minimized and the resulting Euler equation is solved. The study of the soft-sphere potentials confirms the appearance of a dependence of the energy on the shape of the potential at gas paremeter values of x~0.001. For quantities other than the energy, such as the radial distribution functions and the momentum distributions, the dependence appears at any value of x. The occurrence of a maximum in the radial distribution function, in the momentum distribution, and in the excitation spectrum is a natural effect of the correlations when x increases. The asymptotic behaviors of the functions characterizing the structure of the systems are also investigated. The uniform limit approach is very easy to implement and provides a good description of the soft-sphere gas. Its reliability improves when the interaction weakens.

High-Momentum Response of Liquid 3He

A final-state-effects formalism suitable to analyze the high-momentum response of Fermi liquids is presented and used to study the dynamic structure function of liquid 3He. The theory, developed as a natural extension of the Gersch-Rodriguez formalism, incorporates the Fermi statistics explicitly through a new additive term which depends on the semidiagonal two-body density matrix. The use of a realistic momentum distribution, calculated using the diffusion Monte Carlo method, and the inclusion of this additive correction allows for good agreement with available deep-inelastic neutron scattering data.