Three-dimensional model for multi-component reactive transport with variable density groundwater flow

PHWAT is a new model that couples a geochemical reaction model (PHREEQC-2) with a density-dependent groundwater flow and solute transport model (SEAWAT) using the split-operator approach. PHWAT was developed to simulate multi-component reactive transport in variable density groundwater flow. Fluid density in PHWAT depends not on only the concentration of a single species as in SEAWAT, but also the concentrations of other dissolved chemicals that can be subject to reactive processes. Simulation results of PHWAT and PHREEQC-2 were compared in their predictions of effluent concentration from a column experiment. Both models produced identical results, showing that PHWAT has correctly coupled the sub-packages. PHWAT was then applied to the simulation of a tank experiment in which seawater intrusion was accompanied by cation exchange. The density dependence of the intrusion and the snow-plough effect in the breakthrough curves were reflected in the model simulations, which were in good agreement with the measured breakthrough data. Comparison simulations that, in turn, excluded density effects and reactions allowed us to quantify the marked effect of ignoring these processes. Next, we explored numerical issues involved in the practical application of PHWAT using the example of a dense plume flowing into a tank containing fresh water. It was shown that PHWAT could model physically unstable flow and that numerical instabilities were suppressed. Physical instability developed in the model in accordance with the increase of the modified Rayleigh number for density-dependent flow, in agreement with previous research. (c) 2004 Elsevier Ltd. All rights reserved.

Ferromagnetism, paramagnetism, and a Curie-Weiss metal in an electron-doped Hubbard model on a triangular lattice

Motivated by the unconventional properties and rich phase diagram of NaxCoO2 we consider the electronic and magnetic properties of a two-dimensional Hubbard model on an isotropic triangular lattice doped with electrons away from half-filling. Dynamical mean-field theory (DMFT) calculations predict that for negative intersite hopping amplitudes (t < 0) and an on-site Coulomb repulsion, U, comparable to the bandwidth, the system displays properties typical of a weakly correlated metal. In contrast, for t > 0 a large enhancement of the effective mass, itinerant ferromagnetism, and a metallic phase with a Curie-Weiss magnetic susceptibility are found in a broad electron doping range. The different behavior encountered is a consequence of the larger noninteracting density of states (DOS) at the Fermi level for t > 0 than for t < 0, which effectively enhances the mass and the scattering amplitude of the quasiparticles. The shape of the DOS is crucial for the occurrence of ferromagnetism as for t > 0 the energy cost of polarizing the system is much smaller than for t < 0. Our observation of Nagaoka ferromagnetism is consistent with the A-type antiferromagnetism (i.e., ferromagnetic layers stacked antiferromagnetically) observed in neutron scattering experiments on NaxCoO2. The transport and magnetic properties measured in NaxCoO2 are consistent with DMFT predictions of a metal close to the Mott insulator and we discuss the role of Na ordering in driving the system towards the Mott transition. We propose that the Curie-Weiss metal phase observed in NaxCoO2 is a consequence of the crossover from a bad metal with incoherent quasiparticles at temperatures T > T-* and Fermi liquid behavior with enhanced parameters below T-*, where T-* is a low energy coherence scale induced by strong local Coulomb electron correlations. Our analysis also shows that the one band Hubbard model on a triangular lattice is not enough to describe the unusual properties of NaxCoO2 and is used to identify the simplest relevant model that captures the essential physics in NaxCoO2. We propose a model which allows for the Na ordering phenomena observed in the system which, we propose, drives the system close to the Mott insulating phase even at large dopings.

A new algorithm and refined bounds for extended gcd computation

Extended gcd computation is interesting itself. It also plays a fundamental role in other calculations. We present a new algorithm for solving the extended gcd problem. This algorithm has a particularly simple description and is practical. It also provides refined bounds on the size of the multipliers obtained.

Paramagnetism of caesium titanium alum

Caesium titanium alum, CsTi(SO4)(2) . 12H(2)O, is a beta alum and exhibits a large trigonal field and a dynamic Jahn-Teller effect. Exact calculations of the linear (2)T(2)xe Jahn-Teller coupling show that in the strict S-6 Site symmetry the ground multiplet consists of a Kramers doublet 2 Gamma(6) with magnetic splitting factors g(parallel to)=1.1 and g perpendicular to=0, a Gamma(4) Gamma(5) doublet at similar to 60 cm(-1) with g(parallel to)=2.51 and g(perpendicular to)=0.06 and another Gamma(4) Gamma(5) doublet at similar to 270 cm(-1) with g(parallel to)=1.67 and g(perpendicular to)=1.83. The controversial g values observed below 4.2 K, g(parallel to)=1.25 and g(perpendicular to)=1.14, are shown to arise from low symmetry distortions. These distortions couple the vibronic levels and induce into the ground state the off-diagonal axial Zeeman interaction that exists between the first excited and the ground vibronic levels. (C) 1997 American Institute of Physics.

Population versus clinical view of case fatality from acute coronary heart disease - Results from the WHO MONICA Project 1985-1990

Background The clinical view of case fatality (CF) from acute myocardial infarction (AMI) in those reaching the hospital alive is different from the population view. Registration of both hospitalized AMI cases and out-of-hospital coronary heart disease (CHD) deaths in the WHO MONICA Project allows both views to be reconciled. The WHO MONICA Project provides the largest data set worldwide to explore the relationship between CHD CF and age, sex, coronary event rate, and first versus recurrent event. Methods and Results All 79 669 events of definite AMI or possible coronary death, occurring from 1985 to 90 among 5 725 762 people, 35 to 64 years of age, in 29 MONICA populations are the basis for CF calculations. Age-adjusted CF (percentage of CHD events that were fatal) was calculated across populations, stratified for different time periods, and related to age, sex, and CHD event rate. Median 28-day population CF was 49% (range, 35% to 60%) in men and 51% (range, 34% to 70%) in women and was particularly higher in women than men in populations in which CHD event rates were low. Median 28-day CF for hospitalized events was much lower: in men 22% (range, 15% to 36%) and in women 27% (range, 19% to 46%). Among hospitalized events CF was twice as high for recurrent as for first events. Conclusions Overall 28-day CF is halved for hospitalized events compared with all events and again nearly halved for hospitalized 24-hour survivors. Because approximately two thirds of 28-day CHD deaths in men and women occurred before reaching the hospital, opportunities for reducing CF through improved care in the acute event are limited. Major emphasis should be on primary and secondary prevention.

Introduction to Discounted Cash Flow Analysis and Financial Functions in Excel

The financial and economic analysis of investment projects is typically carried out using the technique of discounted cash flow (DCF) analysis. This module introduces concepts of discounting and DCF analysis for the derivation of project performance criteria such as net present value (NPV), internal rate of return (IRR) and benefit to cost (B/C) ratios. These concepts and criteria are introduced with respect to a simple example, for which calculations using MicroSoft Excel are demonstrated.

Matrix-product codes over Fq

Codes C-1,...,C-M of length it over F-q and an M x N matrix A over F-q define a matrix-product code C = [C-1 (...) C-M] (.) A consisting of all matrix products [c(1) (...) c(M)] (.) A. This generalizes the (u/u + v)-, (u + v + w/2u + v/u)-, (a + x/b + x/a + b + x)-, (u + v/u - v)- etc. constructions. We study matrix-product codes using Linear Algebra. This provides a basis for a unified analysis of /C/, d(C), the minimum Hamming distance of C, and C-perpendicular to. It also reveals an interesting connection with MDS codes. We determine /C/ when A is non-singular. To underbound d(C), we need A to be 'non-singular by columns (NSC)'. We investigate NSC matrices. We show that Generalized Reed-Muller codes are iterative NSC matrix-product codes, generalizing the construction of Reed-Muller codes, as are the ternary 'Main Sequence codes'. We obtain a simpler proof of the minimum Hamming distance of such families of codes. If A is square and NSC, C-perpendicular to can be described using C-1(perpendicular to),...,C-M(perpendicular to) and a transformation of A. This yields d(C-perpendicular to). Finally we show that an NSC matrix-product code is a generalized concatenated code.

Colourful objects through animal eyes

To understand how bees, birds, and fish may use colour vision for food selection and mate choice, we reconstructed views of biologically important objects taking into account the receptor spectral sensitivities. Reflectance spectra a of flowers, bird plumage, and fish skin were used to calculate receptor quantum catches. The quantum catches were then coded by red, green, and blue of a computer monitor; and powers, birds, and fish were visualized in animal colours. Calculations were performed for different illumination conditions. To simulate colour constancy, we used a von Kries algorithm, i.e., the receptor quantum catches were scaled so that the colour of illumination remained invariant. We show that on land this algorithm compensates reasonably well for changes of object appearance caused by natural changes of illumination, while in water failures of von Kries colour constancy are prominent. (C) 2000 John Wiley & Sons, Inc.

Diverse solid-state and solution structures within a series of hexaamine dicopper(II) complexes

Mono- and dicopper(II) complexes of a series of potentially bridging hexaamine ligands have been prepared and characterized in the solid state by X-ray crystallography. The crystal structures of the following Cu-II complexes are reported: [Cu(HL3)](ClO4)(3), C11H31Cl3CuN6O12, monoclinic, P2(1)/n, a = 8.294(2) Angstrom, b = 18.364(3) Angstrom, c = 15.674(3) Angstrom, beta = 94.73(2)degrees, Z = 4; {[Cu-2(L-4)(CO3)](2)}(ClO4)(4). 4H(2)O, C40H100Cl4Cu4N12O26, triclinic, P (1) over bar, a = 9.4888(8) Angstrom, b=13.353(1) Angstrom,. c = 15.329(1) Angstrom, alpha = 111.250(7)degrees, beta = 90.068(8)degrees, gamma = 105.081(8)degrees, Z=1; [Cu-2(L-5)(OH2)(2)](ClO4)(4), C(13)H(36)Cl(4)Cu(2)Z(6)O(18), monoclinic, P2(1)/c, a = 7.225(2) Angstrom. b = 8.5555(5) Angstrom, c = 23.134(8) Angstrom, beta = 92.37(1)degrees, Z = 2; [Cu-2(L-6)(OH2)(2)](ClO4)(4). 3H(2)O, C14H44Cl4Cu2N6O21, monoclinic, P2(1)/a, a = 15.204(5) Angstrom, b = 7.6810(7) Angstrom, c = 29.370(1) Angstrom, beta = 100.42(2)degrees, Z = 4. Solution spectroscopic properties of the bimetallic complexes indicate that significant conformational changes occur upon dissolution, and this has been probed with EPR spectroscopy and molecular mechanics calculations.

Iminopropadienethiones, Ar-N=C=C=C=S

Aryliminopropadienethiones 9 have been generated by flash vacuum thermolysis of isoxazolones of the type 5 and characterized by mass spectrometry and matrix isolation IR spectroscopy in conjunction with DFT calculations and chemical trapping.

Analysis of acousto-ultrasonic characteristics for contact-type transducers coupled to composite laminated plates

The acousto-ultrasonic (AU) input-output characteristics for contact-type transmitting and receiving transducers coupled to composite laminated plates are considered in this paper. Combining a multiple integral transform method, an ordinary discrete layer theory for the laminates and some simplifying assumptions for the electro-mechanical transduction behaviour of the transducers, an analytical solution is developed which can deal with all the wave processes involved in the AU measurement system, i.e, wave generation, wave propagation and wave reception. The spectral response of the normal contact pressure sensed by the receiving transducer due to an arbitrary input pulse excited by the transmitting transducer is obtained. To validate the new analytical-numerical spectral technique in the low-frequency regime, the results are compared with Mindlin plate theory solutions. Based on the analytical results, numerical calculations are carried out to investigate the influence of various external parameters such as frequency content of the input pulse, transmitter/receiver spacing and transducer aperture on the output of the measurement system. The results show that the presented analytical-numerical procedure is an effective tool for understanding the input-output characteristics of the AU technique for laminated plates. (C) 2001 Elsevier Science Ltd. All rights reserved.

Interlaminar differences in the pyramidal cell phenotype in cortical areas 7m and STP (the superior temporal polysensory area) of the macaque monkey

Pyramidal neurones were injected with Lucifer Yellow in slices cut tangential to the surface of area 7m and the superior temporal polysensory area (STP) of the macaque monkey. Comparison of the basal dendritic arbors of supra- and infragranular pyramidal neurones (n=139) that were injected in the same putative modules in the different cortical areas revealed variation in their structure. Moreover, there were relative differences in dendritic morphology of supra- and infragranular pyramidal neurones in the two cortical areas. Shell analyses revealed that layer III pyramidal neurones in area STP had considerably higher peak complexity (maximum number of dendritic intersections per Shell circle) than those in layer V, whereas peak complexities were similar for supra- and infragranular pyramidal neurones in area 7m. In both cortical areas, the basal dendritic trees of layer m pyramidal neurones were characterized by a higher spine density than those in layer V. Calculations of the total number of dendritic spines in the average basal dendritic arbor revealed that layer V pyramidal neurones in area 7m had twice as many spines as cells in layer III. (4535 and 2294, respectively). A similar calculation for neurones in area STP revealed that layer III pyramidal neurones had approximately the same number of spines as cells in layer V (3585 and 3850 spines, respectively). Relative differences in the branching patterns of, and the number of spines in, the basal dendritic arbors of supra- and infragranular pyramidal neurones in the different cortical areas may allow for integration of different numbers of inputs, and different degrees of dendritic processing. These results support the thesis that intra-areal circuitry differs in different cortical areas.

Stochastic gauges in quantum dynamics for many-body simulations

Quantum dynamics simulations can be improved using novel quasiprobability distributions based on non-orthogonal Hermitian kernel operators. This introduces arbitrary functions (gauges) into the stochastic equations. which can be used to tailor them for improved calculations. A possible application to full quantum dynamic simulations of BEC's is presented. (C) 2001 Elsevier Science B.V. All rights reserved.

Quantum criticality

We investigate the theory of quantum fluctuations in non-equilibrium systems having large critical fluctuations. This allows us to treat the limits imposed by nonlinearities to quantum squeezing and noise reduction, and also to envisage future tests of quantum theory in regions of macroscopic quantum fluctuations. A long-term objective of this research is to identify suitable physical systems in which macroscopic 'Schrodinger cat'-like behaviour may be observed. We investigate two systems in particular of much current experimental interest, namely the degenerate parametric oscillator near threshold, and the evaporatively cooled (BEC). We compare the results obtained in the positive-P representation, as a fully quantum mechanical calculation, with the truncated Wigner phase space equation, also known as semi-classical theory. We show when these results agree and differ in calculations taken beyond the linearized approximation. In the region where the largest quantum fluctuations and Schrodinger cat-like behaviour might be expected, we find that the quantum predictions correspond very closely to the semi-classical theory. Nature abhors observing a Schrodinger car.

The geochemistry of late Archaean microbial carbonate: Implications for ocean chemistry and continental erosion history

Trace element concentrations and combined Sr- and Nd-isotope compositions were determined on stromatolitic carbonates (microbialites) from the 2.52 Ga Campbellrand carbonate platform (South Africa). Shale-normalised rare earth element and yttrium patterns of the ancient samples are similar to those of modern seawater in having positive La and Y anomalies and in being depleted in light rare earth elements. In contrast to modem seawater (and microbialite proxies), the 2.52 Ga samples lack a negative Ce anomaly but possess a positive Eu anomaly. These latter trace element characteristics are interpreted to reflect anoxic deep ocean waters where, unlike today, hydrothermal Fe input was not oxidised, and scavenged and rare earth elements were not coprecipitated with Fe-oxyhydroxides. The persistence of a positive Eu anomaly in relatively shallow Campbellrand platform waters indicates a dramatic reversal from hydrothermally dominated (Archaean) to continental erosion-dominated (Phanerozoic) rare earth element flux ratio. The dominant hydrothermal input is also expressed in the initial Sr- and Nd-isotope ratios. There is collinear variation in Sr-Nd systematics, which range from primitive values (Sr-87/Sr-86 of 0.702386 and epsilon (Nd) of +2.1) to more evolved crustal ratios. Mixing calculations show that the range in trace element ratios (e.g., Y/Ho) and initial isotope ratios is not a result of contamination by trapped sediment, but that the chemical band isotopic variation reflects carbonate deposition in an environment where different water masses mixed. Calculated Nd flux ratios yield a hydrothermal input into the 2.52 Ga oceans one order of magnitude larger than continental input. Such a change in flux ratio most likely required substantially reduced continental inputs, which could, in turn, reflect a plate tectonic causation (e.g., reduced topography or expansion of epicontinental seas). Copyright (C) 2001 Elsevier Science Ltd.