Statistical power calculation and sample size determination for environmental studies with data below detection limits

Power calculation and sample size determination are critical in designing environmental monitoring programs. The traditional approach based on comparing the mean values may become statistically inappropriate and even invalid when substantial proportions of the response values are below the detection limits or censored because strong distributional assumptions have to be made on the censored observations when implementing the traditional procedures. In this paper, we propose a quantile methodology that is robust to outliers and can also handle data with a substantial proportion of below-detection-limit observations without the need of imputing the censored values. As a demonstration, we applied the methods to a nutrient monitoring project, which is a part of the Perth Long-Term Ocean Outlet Monitoring Program. In this example, the sample size required by our quantile methodology is, in fact, smaller than that by the traditional t-test, illustrating the merit of our method.

Improved Algorithms for Reservoir Capacity Calculation Incorporating Storage-Dependent Losses and Reliability Norm

Two algorithms that improve upon the sequent-peak procedure for reservoir capacity calculation are presented. The first incorporates storage-dependent losses (like evaporation losses) exactly as the standard linear programming formulation does. The second extends the first so as to enable designing with less than maximum reliability even when allowable shortfall in any failure year is also specified. Together, the algorithms provide a more accurate, flexible and yet fast method of calculating the storage capacity requirement in preliminary screening and optimization models.

Error bounds for calculation of the Gittins indices

For a wide class of semi-Markov decision processes the optimal policies are expressible in terms of the Gittins indices, which have been found useful in sequential clinical trials and pharmaceutical research planning. In general, the indices can be approximated via calibration based on dynamic programming of finite horizon. This paper provides some results on the accuracy of such approximations, and, in particular, gives the error bounds for some well known processes (Bernoulli reward processes, normal reward processes and exponential target processes).

Theoretical description of time-resolved pump/probe photoemission in TaS2: a single-band DFT+DMFT(NRG) study within the quasiequilibrium approximation

In this work, we theoretically examine recent pump/probe photoemission experiments on the strongly correlated charge-density-wave insulator TaS2.We describe the general nonequilibrium many-body formulation of time-resolved photoemission in the sudden approximation, and then solve the problem using dynamical mean-field theory with the numerical renormalization group and a bare density of states calculated from density functional theory including the charge-density-wave distortion of the ion cores and spin-orbit coupling. We find a number of interesting results: (i) the bare band structure actually has more dispersion in the perpendicular direction than in the two-dimensional planes; (ii) the DMFT approach can produce upper and lower Hubbard bands that resemble those in the experiment, but the upper bands will overlap in energy with other higher energy bands; (iii) the effect of the finite width of the probe pulse is minimal on the shape of the photoemission spectra; and (iv) the quasiequilibrium approximation does not fully describe the behavior in this system.

On Calculation Of Breakdown Voltages Of Mixtures Of Electron Attaching Gases

In this paper we report the analysis of dc breakdown tests on mixtures of CC12F2, SF6, C-C4F8, 2-C4F8, N2, C02, CF4, CHF3, and 1,1,1-CH3CF3 gases on the basis of the NKH formula Vmix=k(pd)aNbUC developed by us earlier for the binary mixtures of SF6 with air, N2, N20, and CO2. It is shown that while a and c have the values 0.915 and 0.850 respectively as earlier, k and b depend on the component gases. There is a good agreement between the calculated values on the basis of the formula and measured values reported in the literature.

On Calculation Of Breakdown Voltages Of Mixtures Of Electron Attaching Gases

In this paper we report the analysis of dc breakdown tests on mixtures of CC12F2, SF6, C-C4F8, 2-C4F8, N2, C02, CF4, CHF3, and 1,1,1-CH3CF3 gases on the basis of the NKH formula Vmix=k(pd)aNbUC developed by us earlier for the binary mixtures of SF6 with air, N2, N20, and CO2. It is shown that while a and c have the values 0.915 and 0.850 respectively as earlier, k and b depend on the component gases. There is a good agreement between the calculated values on the basis of the formula and measured values reported in the literature.

DFT study of electronic and optical properties of anatase titanium dioxide tuned by nitrogen and lithium co-doping

The electronic and optical properties of anatase titanium dioxide (TiO2), co-doped by nitrogen (N) and lithium (Li), have been investigated by density functional theory plus Hubbard correction term U, namely DFT+U. It is found that Li-dopants can effectively balance the net charges brought by N-dopants and shift the local state to the top of valence band. Depending on the distribution of dopants, the adsorption edges of TiO2 may be red- or blue-shifted, being consistent with recent experimental observations.

Closo versus Hypercloso Metallaboranes: A DFT Study

The structures and electronic relationship of 9-, 10-, 11-, and 12-vertex closo and hypercloso (isocloso) etallaboranes are explored using OFT calculations. The role of the transition metal in stabilizing the hypercloso borane structures is explained using the concept of orbital compatibility. The hypercloso structures, C6H6MBn-1Hn-1 (n = 9-12; M = Fe, Ru, and Os) are taken as model complexes. Calculations on metal free polyhedral borane BnHn suggest that n vertex hypercloso structures need only n skeleton electron pairs (SEPs), but the structure will have one or more six-degree vertices, whereas the corresponding closo structures with n + 1 SEPs have only four- and five-degree vertices. This high-degree vertex of hypercloso structures can be effectively occupied by transition metal fragments with their highly diffused orbitals. Calculations further show that a heavy transition metal with more diffused orbitals prefers over a light transition metal to form hypercloso geometry, This is in accordance with the fact that there are more experimentally characterized hypercloso structures with the heavy transition metals. The size of the exohedral ligands attached to the metal atom also plays a role in deciding the stability of the hypercloso structure. The interaction between the borane and the metal fragments in the hypercloso geometry is analyzed using the fragment molecular orbital approach. The interconversion of the closo and hypercloso structures by the addition and removal of the electrons is also discussed in terms of the correlation diagrams.

Calculation of Sparking Potentials of SF6 and SF6-GAS Mixtures in Uniform and Non-Uniform Electric Fields

A simple formula is developed to predict the sparking potentials of SF6 and SF6-gas mixture in uniform and non-uniform fields. The formula has been shown to be valid over a very wide range from 1 to 1800 kPa·cm of pressure and electrode gap separation for mixtures containing 5 to 100% SF6. The calculated values are found to be in good agreement with the previously reported measurements in the literature. The formula should aid design engineers in estimating electrode-spacings and clearances in power apparatus and systems.

Tissue thickness calculation in ocular optical coherence tomography

Thickness measurements derived from optical coherence tomography (OCT) images of the eye are a fundamental clinical and research metric, since they provide valuable information regarding the eye’s anatomical and physiological characteristics, and can assist in the diagnosis and monitoring of numerous ocular conditions. Despite the importance of these measurements, limited attention has been given to the methods used to estimate thickness in OCT images of the eye. Most current studies employing OCT use an axial thickness metric, but there is evidence that axial thickness measures may be biased by tilt and curvature of the image. In this paper, standard axial thickness calculations are compared with a variety of alternative metrics for estimating tissue thickness. These methods were tested on a data set of wide-field chorio-retinal OCT scans (field of view (FOV) 60° x 25°) to examine their performance across a wide region of interest and to demonstrate the potential effect of curvature of the posterior segment of the eye on the thickness estimates. Similarly, the effect of image tilt was systematically examined with the same range of proposed metrics. The results demonstrate that image tilt and curvature of the posterior segment can affect axial tissue thickness calculations, while alternative metrics, which are not biased by these effects, should be considered. This study demonstrates the need to consider alternative methods to calculate tissue thickness in order to avoid measurement error due to image tilt and curvature.

A Methanol-Tolerant Carbon-Supported Pt-Au Alloy Cathode Catalyst for Direct Methanol Fuel Cells and Its Evaluation by DFT

A Pt-Au alloy catalyst of varying compositions is prepared by codeposition of Pt and Au nanoparticles onto a carbon support to evaluate its electrocatalytic activity toward an oxygen reduction reaction (ORR) with methanol tolerance in direct methanol fuel cells. The optimum atomic weight ratio of Pt to Au in the carbon-supported Pt-Au alloy (Pt-Au/C) as established by cell polarization, linear-sweep voltammetry (LSV), and cyclic voltammetry (CV) studies is determined to be 2:1. A direct methanol fuel cell (DMFC) comprising a carbon-supported Pt-Au (2:1) alloy as the cathode catalyst delivers a peak power density of 120 mW/cm2 at 70 °C in contrast to the peak power density value of 80 mW/cm2 delivered by the DMFC with carbon-supported Pt catalyst operating under identical conditions. Density functional theory (DFT) calculations on a small model cluster reflect electron transfer from Pt to Au within the alloy to be responsible for the synergistic promotion of the oxygen-reduction reaction on a Pt-Au electrode.

Experimental and Theoretical Charge Density Analysis of Polymorphic Structures: The Case of Coumarin 314 Dye

Experimental charge density distributions in two known conformational polymorphs (orange and yellow) of coumarin 314 dye are analyzed based on multipole modeling of X-ray diffraction data collected at 100 K. The experimental results are compared with the charge densities derived from multipole modeling of theoretical structure factors obtained from periodic quantum calculation with density functional theory (DFT) method and B3LYP/6-31G(d,p) level of theory. The presence of disorder at the carbonyl oxygen atom of ethoxycarbonyl group in the yellow form, which was not identified earlier, is addressed here. The investigationof intermolecular interactions, based on Hirshfeld surface analysis and topological properties via quantum theory of atoms in molecule and total electrostatic interaction energies, revealed significant differences between the polymorphs. The differences of electrostatic nature in these two polymorphic forms were unveiled via construction of three-dimensional deformation electrostatic potential maps plotted over the molecular surfaces. The lattice energies evaluated from ab initio calculations on the two polymorphic forms indicate that the yellow form is likely to be the most favorable thermodynamically. The dipole moments derived from experimental and theoretical charge densities and also from Lorentz tensor approach are compared with the single-molecule dipole moments. In each case, the differences of dipole moments between the polymorphs are identified.