Stem-root flow effect on soil–atmosphere interactions and uncertainty assessments

Abstract. Soil water can rapidly enter deeper layers via vertical redistribution of soil water through the stem–root flow mechanism. This study develops the stem–root flow parameterization scheme and coupled this scheme with the Simplified Simple Biosphere model (SSiB) to analyze its effects on land–atmospheric interactions. The SSiB model was tested in a single column mode using the Lien Hua Chih (LHC) measurements conducted in Taiwan and HAPEX-Mobilhy (HAPEX) measurements in France. The results show that stem–root flow generally caused a decrease in the moisture content at the top soil layer and moistened the deeper soil layers. Such soil moisture redistribution results in significant changes in heat flux exchange between land and atmosphere. In the humid environment at LHC, the stem–root flow effect on transpiration was minimal, and the main influence on energy flux was through reduced soil evaporation that led to higher soil temperature and greater sensible heat flux. In the Mediterranean environment of HAPEX, the stem–root flow significantly affected plant transpiration and soil evaporation, as well as associated changes in canopy and soil temperatures. However, the effect on transpiration could either be positive or negative depending on the relative changes in the moisture content of the top soil vs. deeper soil layers due to stem–root flow and soil moisture diffusion processes

Model Calculation of Effective Three-Body Forces

We propose a scheme for extracting an effective three-body interaction originating from a two-nucleon interaction. This is based on the Q-box method of Kuo and collaborators, where folded diagrams are obtained by differentiating a sum of non-folded diagrams with respect to the starting energy. To gain insight we have studied several examples using the Lipkin model where the perturbative approach can be compared with exact results. Numerically the three-body interactions can be significant and in a matrix example good accuracy was not obtained simultaneously for both eigenvalues with two-body interactions alone.Comment: 9 pages, Revtex4, 7 figs, submitted to PR

Approximate Treatment of Hermitian Effective Interactions and a Bound on the Error

The Hermitian effective interaction can be well-approximated by (R+R^dagger)/2 if the eigenvalues of omega^dagger omega are small or state-independent(degenerate), where R is the standard non-Hermitian effective interaction and omega maps the model-space states onto the excluded space. An error bound on this approximation is given.Comment: 13 page

Three-body monopole corrections to the realistic interactions

It is shown that a very simple three-body monopole term can solve practically all the spectroscopic problems--in the $p$, $sd$ and $pf$ shells--that were hitherto assumed to need drastic revisions of the realistic potentials.Comment: 4 pages, 5figure

Rigorous treatment of electrostatics for spatially varying dielectrics based on energy minimization

A novel energy minimization formulation of electrostatics that allows computation of the electrostatic energy and forces to any desired accuracy in a system with arbitrary dielectric properties is presented. An integral equation for the scalar charge density is derived from an energy functional of the polarization vector field. This energy functional represents the true energy of the system even in non-equilibrium states. Arbitrary accuracy is achieved by solving the integral equation for the charge density via a series expansion in terms of the equation's kernel, which depends only on the geometry of the dielectrics. The streamlined formalism operates with volume charge distributions only, not resorting to introducing surface charges by hand. Therefore, it can be applied to any spatial variation of the dielectric susceptibility, which is of particular importance in applications to biomolecular systems. The simplicity of application of the formalism to real problems is shown with analytical and numerical examples.Comment: 27 pages, 5 figure

Hot new directions for quasi-Monte Carlo research in step with applications

This article provides an overview of some interfaces between the theory of quasi-Monte Carlo (QMC) methods and applications. We summarize three QMC theoretical settings: first order QMC methods in the unit cube $[0,1]^s$ and in $\mathbb{R}^s$, and higher order QMC methods in the unit cube. One important feature is that their error bounds can be independent of the dimension $s$ under appropriate conditions on the function spaces. Another important feature is that good parameters for these QMC methods can be obtained by fast efficient algorithms even when $s$ is large. We outline three different applications and explain how they can tap into the different QMC theory. We also discuss three cost saving strategies that can be combined with QMC in these applications. Many of these recent QMC theory and methods are developed not in isolation, but in close connection with applications