Some monotonicity properties of symmetric pólya densities and their exponential families

Summary In this note we observe that for independent symmetric random variables X and Y, when the pdf of X is PF, the conditional distributions of ¦Y¦ given S = X + Y form a MLR family. We then show that for a function : R ⁿR that is symmetric in each coordinate and increasing on (0, )ⁿ, E((S₁,...,S_n)¦S_n = s) is even and increasing in ¦s¦. Here S₁,...,S_n are partial sums with independent symmetric PF summands. Application is made to sequential tests that minimize the maximum expected sample size when the model is a one-parameter exponential family generated by a symmetric PF density.Work supported by NSF grants MPS 72-05082 AO2 and MCS 75-23344     

Lattice-Boltzmann simulations of the dynamics of polymer solutions in periodic and confined geometries

A numerical method to simulate the dynamics of polymer solutions in confined geometries has been implemented and tested. The method combines a fluctuating lattice-Boltzmann model of the solvent [Ladd, Phys. Rev. Lett. 70, 1339 (1993)] with a point-particle model of the polymer chains. A friction term couples the monomers to the fluid [Ahlrichs and Dunweg, J. Chem. Phys. 111, 8225 (1999)], providing both the hydrodynamic interactions between the monomers and the correlated random forces. The coupled equations for particles and fluid are solved on an inertial time scale, which proves to be surprisingly simple and efficient, avoiding the costly linear algebra associated with Brownian dynamics. Complex confined geometries can be represented by a straightforward mapping of the boundary surfaces onto a regular three-dimensional grid. The hydrodynamic interactions between monomers are shown to compare well with solutions of the Stokes equations down to distances of the order of the grid spacing. Numerical results are presented for the radius of gyration, end-to-end distance, and diffusion coefficient of an isolated polymer chain, ranging from 16 to 1024 monomers in length. The simulations are in excellent agreement with renormalization group calculations for an excluded volume chain. We show that hydrodynamic interactions in large polymers can be systematically coarse-grained to substantially reduce the computational cost of the simulation. Finally, we examine the effects of confinement and flow on the polymer distribution and diffusion constant in a narrow channel. Our results support the qualitative conclusions of recent Brownian dynamics simulations of confined polymers [Jendrejack et al., J. Chem. Phys. 119, 1165 (2003) and Jendrejack et al., J. Chem. Phys. 120, 2513 (2004)].     

GROMACS: fast, flexible, and free

This article describes the software suite GROMACS (Groningen MAchine for Chemical Simulation) that was developed at the University of Groningen, The Netherlands, in the early 1990s. The software, written in ANSI C, originates from a parallel hardware project, and is well suited for parallelization on processor clusters. By careful optimization of neighbor searching and of inner loop performance, GROMACS is a very fast program for molecular dynamics simulation. It does not have a force field of its own, but is compatible with GROMOS, OPLS, AMBER, and ENCAD force fields. In addition, it can handle polarizable shell models and flexible constraints. The program is versatile, as force routines can be added by the user, tabulated functions can be specified, and analyses can be easily customized. Nonequilibrium dynamics and free energy determinations are incorporated. Interfaces with popular quantum-chemical packages (MOPAC, GAMES-UK, GAUSSIAN) are provided to perform mixed MM/QM simulations. The package includes about 100 utility and analysis programs. GROMACS is in the public domain and distributed (with source code and documentation) under the GNU General Public License. It is maintained by a group of developers from the Universities of Groningen, Uppsala, and Stockholm, and the Max Planck Institute for Polymer Research in Mainz. Its Web site is http://www.gromacs.org.     

Enamines of 1,3-dimethylbarbiturates and their symmetrical palladium(II) complexes: synthesis,characterization and biological activity

Transition Metal Chemistry - Three 1,3-dimethylbarbiturate-enamine derivatives and their symmetrical palladium(II) complexes were prepared and characterized by spectroscopic methods. In addition,...     

Model for spontaneous frequency sweeping of an Alfvén wave in a toroidal plasma

We study the frequency chirping signals arising from spontaneously excited toroidial Alfvén eigenmode (TAE) waves that are being driven by an inverted energetic particle distribution whose free energy is tapped from the generic particle/wave resonance interaction. Initially a wave is excited inside the Alfvén gap with a frequency determined from the linear tip model of Rosenbluth, Berk and Van dam (RBV) [1]. Hole/clumps structures are formed and are observed to chirp towards lower energy states. We find that the chirping signals from clump enter the Alfvén continuum which eventually produce more rapid chirping signals. The accuracy of the adiabatic approximation for the mode evolution is tested and verified by demonstrating that a WKB-like decomposition of the time response for the field phase and amplitude agree with the data. Plots of the phase space structure correlate well with the chirping dependent shape of the separatrix structure. A novel aspect of the simulation is that it performed close to the wave frame of the phase space structure, which enables the numerical time step to remain the same during the simulation, independent of the rest frame frequency.     

Supplier diversification under binomial yield

We consider supplier diversification in an EOQ type inventory setting with multiple suppliers and binomial yields. We characterize the optimal policy for the model and show that, in this case, it does not pay to diversify, in contrast to previous results in the random yield literature.     

Bayesian demand updating in the lost sales newsvendor problem: A two-moment approximation

We consider Bayesian updating of demand in a lost sales newsvendor model with censored observations. In a lost sales environment, where the arrival process is not recorded, the exact demand is not observed if it exceeds the beginning stock level, resulting in censored observations. Adopting a Bayesian approach for updating the demand distribution, we develop expressions for the exact posteriors starting with conjugate priors, for negative binomial, gamma, Poisson and normal distributions. Having shown that non-informative priors result in degenerate predictive densities except for negative binomial demand, we propose an approximation within the conjugate family by matching the first two moments of the posterior distribution. The conjugacy property of the priors also ensure analytical tractability and ease of computation in successive updates. In our numerical study, we show that the posteriors and the predictive demand distributions obtained exactly and with the approximation are very close to each other, and that the approximation works very well from both probabilistic and operational perspectives in a sequential updating setting as well.     

Biodistribution of GdCl3 and Gd-DTPA and their influence on proton magnetic relaxation in rat tissues

The biodistribution and relative molar effectiveness of the ionic (GdCl3) and chelated (Gd-DTPA) forms of gadolinium (Gd) to enhance proton relaxation rates in rat kidney, liver and spleen were evaluated. Rats were given intravenous injections of either GdCl3 (100 mumol/kg) or Gd-DTPA (178 mumol/kg). Gd-DTPA was primarily contained in the vascular compartment and was quickly accumulated in the kidney after injection with a relaxivity of 4.3 sec-1 (mumol/g kidney)-1. It was eliminated quickly from the body with only 2% of the injected dose remaining after 120 min. After GdCl3 injection, Gd was found primarily in liver and spleen. It accumulated continuously reaching 72% of the injected does in these two tissues after 120 min. Despite this continuous increase in tissue Gd concentration, the relaxation rates showed saturation in liver and spleen. The results suggest that after GdCl3 was injected it distributed either in a protein bound form that was effective at causing relaxation or in a colloid form that was not effective. The biodistribution of GdCl3 was such that it was determined by the phagocytic action of the recticuloendothelial system on a colloid. The biodistribution and tissue relaxivity of Gd-DTPA suggest it will be a useful vascular MRI contrast agent. However, the usefulness of GdCl3 as an MRI contrast agent is limited not only by its acute toxicity but also by its saturable effect on tissue relaxation rates. Consequently, GdCl3 has only a modest influence on tissue relaxivity.