Continuum, Discrete, and Explicit Solvation Models for Describing the Low-Lying Absorption Spectrum of the Pterin Acid in Aqueous Environment

The absorption spectrum of the acid form of pterin in water was investigated theoretically. Different procedures using continuum, discrete, and explicit models were used to include the solvation effect on the absorption spectrum, characterized by two bands. The discrete and explicit models used Monte Carlo simulation to generate the liquid structure and time-dependent density functional theory (B3LYP/6-31G+(d)) to obtain the excitation energies. The discrete model failed to give the correct qualitative effect on the second absorption band. The continuum model, in turn, has given a correct qualitative picture and a semiquantitative description. The explicit use of 29 solvent molecules, forming a hydration shell of 6 angstrom, embedded in the electrostatic field of the remaining solvent molecules, gives absorption transitions at 3.67 and 4.59 eV in excellent agreement with the S(0)-S(1) and S(0)-S(2) absorption bands at of 3.66 and 4.59 eV, respectively, that characterize the experimental spectrum of pterin in water environment. (C) 2010 Wiley Periodicals, Inc. Int J Quantum Chem 110: 2371-2377, 2010

Solvent effects on the electronic absorption spectrum of camphor using continuum, discrete or explicit approaches

We address the effect of solvation on the lowest electronic excitation energy of camphor. The solvents considered represent a large variation in-solvent polarity. We consider three conceptually different ways of accounting for the solvent using either an implicit, a discrete or an explicit solvation model. The solvatochromic shifts in polar solvents are found to be in good agreement with the experimental data for all three solvent models. However, both the implicit and discrete solvation models are less successful in predicting solvatochromic shifts for solvents of low polarity. The results presented suggest the importance of using explicit solvent molecules in the case of nonpolar solvents. (C) 2009 Elsevier B.V. All rights reserved.

Extending the use of canonical and microcanonical Monte Carlo algorithms to spin models

We describe the canonical and microcanonical Monte Carlo algorithms for different systems that can be described by spin models. Sites of the lattice, chosen at random, interchange their spin values, provided they are different. The canonical ensemble is generated by performing exchanges according to the Metropolis prescription whereas in the microcanonical ensemble, exchanges are performed as long as the total energy remains constant. A systematic finite size analysis of intensive quantities and a comparison with results obtained from distinct ensembles are performed and the quality of results reveal that the present approach may be an useful tool for the study of phase transitions, specially first-order transitions. (C) 2009 Elsevier B.V. All rights reserved.

Dynamic effects of breakup on fusion reactions of weakly bound nuclei

The traditional reduction methods to represent the fusion cross sections of different systems are flawed when attempting to completely eliminate the geometrical aspects, such as the heights and radii of the barriers, and the static effects associated with the excess neutrons or protons in weakly bound nuclei. We remedy this by introducing a new dimensionless universal function, which allows the separation and disentanglement of the static and dynamic aspects of the breakup coupling effects connected with the excess nucleons. Applying this new reduction procedure to fusion data of several weakly bound systems, we find a systematic suppression of complete fusion above the Coulomb barrier and enhancement below it. Different behaviors are found for the total fusion cross sections. They are appreciably suppressed in collisions of neutron-halo nuclei, while they are practically not affected by the breakup coupling in cases of stable weakly bound nuclei. (C) 2009 Elsevier B.V. All rights reserved.

Could saturation effects be visible in a future electron-ion collider?

We expect to observe parton saturation in a future electron-ion collider. In this Letter we discuss this expectation in more detail considering two different models which are in good agreement with the existing experimental data on nuclear structure functions. In particular, we study the predictions of saturation effects in electron-ion collisions at high energies, using a generalization for nuclear targets of the b-CGC model, which describes the ep HERA quite well. We estimate the total. longitudinal and charm structure functions in the dipole picture and compare them with the predictions obtained using collinear factorization and modern sets of nuclear parton distributions. Our results show that inclusive observables are not very useful in the search for saturation effects. In the small x region they are very difficult to disentangle from the predictions of the collinear approaches. This happens mainly because of the large uncertainties in the determination of the nuclear parton distribution functions. On the other hand, our results indicate that the contribution of diffractive processes to the total cross section is about 20% at large A and small Q(2), allowing for a detailed study of diffractive observables. The study of diffractive processes becomes essential to observe parton Saturation. (C) 2008 Elsevier B.V. All rights reserved.

Sensitivity of complex networks measurements

Complex networks obtained from real-world networks are often characterized by incompleteness and noise, consequences of imperfect sampling as well as artifacts in the acquisition process. Because the characterization, analysis and modeling of complex systems underlain by complex networks are critically affected by the quality and completeness of the respective initial structures, it becomes imperative to devise methodologies for identifying and quantifying the effects of the sampling on the network structure. One way to evaluate these effects is through an analysis of the sensitivity of complex network measurements to perturbations in the topology of the network. In this paper, measurement sensibility is quantified in terms of the relative entropy of the respective distributions. Three particularly important kinds of progressive perturbations to the network are considered, namely, edge suppression, addition and rewiring. The measurements allowing the best balance of stability (smaller sensitivity to perturbations) and discriminability (separation between different network topologies) are identified with respect to each type of perturbation. Such an analysis includes eight different measurements applied on six different complex networks models and three real-world networks. This approach allows one to choose the appropriate measurements in order to obtain accurate results for networks where sampling bias cannot be avoided-a very frequent situation in research on complex networks.

Interaction of oligonucleotide-based amphiphilic block copolymers with cell membrane models

Oligonucleotides have unique molecular recognition properties, being involved in biological mechanisms such as cell-surface receptor recognition or gene silencing. For their use in human therapy for drug or gene delivery, the cell membrane remains a barrier, but this can be obviated by grafting a hydrophobic tail to the oligonucleotide. Here we demonstrate that two oligonucleotides, one consisting of 12 guanosine units (G(12)), and the other one consisting of five adenosine and seven guanosine (A(5)G(7)) units, when functionalized with poly(butadiene), namely PB-G(12) and PB-A(5)G(7), can be inserted into Langmuir monolayers of dipalmitoyl phosphatidyl choline (DPPC), which served as a cell membrane model. PB-G(12) and PB-A(5)G(7) were found to affect the DPPC monolayer even at high surface pressures. The effects from PB-G(12) were consistently stronger, particularly in reducing the elasticity of the DPPC monolayers, which may have important biological implications. Multilayers of DPPC and nucleotide-based copolymers could be adsorbed onto solid supports, in the form of Y-type LB films, in which the molecular-level interaction led to lower energies in the vibrational spectra of the nucleotide-based copolymers. This successful deposition of solid films opens the way for devices to be produced which exploit the molecular recognition properties of the nucleotides. (C) 2010 Elsevier Inc. All rights reserved.

Effects of electric fields on ultracold Rydberg atom interactions

The behaviour of interacting ultracold Rydberg atoms in both constant electric fields and laser fields is important for designing experiments and constructing realistic models of them. In this paper, we briefly review our prior work and present new results on how electric fields affect interacting ultracold Rydberg atoms. Specifically, we address the topics of constant background electric fields on Rydberg atom pair excitation and laser-induced Stark shifts on pair excitation.

Sampling, WLS, and Mixed Models

Mixed models may be defined with or without reference to sampling, and can be used to predict realized random effects, as when estimating the latent values of study subjects measured with response error. When the model is specified without reference to sampling, a simple mixed model includes two random variables, one stemming from an exchangeable distribution of latent values of study subjects and the other, from the study subjects` response error distributions. Positive probabilities are assigned to both potentially realizable responses and artificial responses that are not potentially realizable, resulting in artificial latent values. In contrast, finite population mixed models represent the two-stage process of sampling subjects and measuring their responses, where positive probabilities are only assigned to potentially realizable responses. A comparison of the estimators over the same potentially realizable responses indicates that the optimal linear mixed model estimator (the usual best linear unbiased predictor, BLUP) is often (but not always) more accurate than the comparable finite population mixed model estimator (the FPMM BLUP). We examine a simple example and provide the basis for a broader discussion of the role of conditioning, sampling, and model assumptions in developing inference.

Improved testing inference in mixed linear models

Mixed linear models are commonly used in repeated measures studies. They account for the dependence amongst observations obtained from the same experimental unit. Often, the number of observations is small, and it is thus important to use inference strategies that incorporate small sample corrections. In this paper, we develop modified versions of the likelihood ratio test for fixed effects inference in mixed linear models. In particular, we derive a Bartlett correction to such a test, and also to a test obtained from a modified profile likelihood function. Our results generalize those in [Zucker, D.M., Lieberman, O., Manor, O., 2000. Improved small sample inference in the mixed linear model: Bartlett correction and adjusted likelihood. Journal of the Royal Statistical Society B, 62,827-838] by allowing the parameter of interest to be vector-valued. Additionally, our Bartlett corrections allow for random effects nonlinear covariance matrix structure. We report simulation results which show that the proposed tests display superior finite sample behavior relative to the standard likelihood ratio test. An application is also presented and discussed. (C) 2008 Elsevier B.V. All rights reserved.

Likelihood ratio tests for variance components in linear mixed models

Although the asymptotic distributions of the likelihood ratio for testing hypotheses of null variance components in linear mixed models derived by Stram and Lee [1994. Variance components testing in longitudinal mixed effects model. Biometrics 50, 1171-1177] are valid, their proof is based on the work of Self and Liang [1987. Asymptotic properties of maximum likelihood estimators and likelihood tests under nonstandard conditions. J. Amer. Statist. Assoc. 82, 605-610] which requires identically distributed random variables, an assumption not always valid in longitudinal data problems. We use the less restrictive results of Vu and Zhou [1997. Generalization of likelihood ratio tests under nonstandard conditions. Ann. Statist. 25, 897-916] to prove that the proposed mixture of chi-squared distributions is the actual asymptotic distribution of such likelihood ratios used as test statistics for null variance components in models with one or two random effects. We also consider a limited simulation study to evaluate the appropriateness of the asymptotic distribution of such likelihood ratios in moderately sized samples. (C) 2008 Elsevier B.V. All rights reserved.

Predicting random effects with an expanded finite population mixed model

Prediction of random effects is an important problem with expanding applications. In the simplest context, the problem corresponds to prediction of the latent value (the mean) of a realized cluster selected via two-stage sampling. Recently, Stanek and Singer [Predicting random effects from finite population clustered samples with response error. J. Amer. Statist. Assoc. 99, 119-130] developed best linear unbiased predictors (BLUP) under a finite population mixed model that outperform BLUPs from mixed models and superpopulation models. Their setup, however, does not allow for unequally sized clusters. To overcome this drawback, we consider an expanded finite population mixed model based on a larger set of random variables that span a higher dimensional space than those typically applied to such problems. We show that BLUPs for linear combinations of the realized cluster means derived under such a model have considerably smaller mean squared error (MSE) than those obtained from mixed models, superpopulation models, and finite population mixed models. We motivate our general approach by an example developed for two-stage cluster sampling and show that it faithfully captures the stochastic aspects of sampling in the problem. We also consider simulation studies to illustrate the increased accuracy of the BLUP obtained under the expanded finite population mixed model. (C) 2007 Elsevier B.V. All rights reserved.

Beta-binomial/Poisson regression models for repeated bivariate counts

We analyze data obtained from a study designed to evaluate training effects on the performance of certain motor activities of Parkinson`s disease patients. Maximum likelihood methods were used to fit beta-binomial/Poisson regression models tailored to evaluate the effects of training on the numbers of attempted and successful specified manual movements in 1 min periods, controlling for disease stage and use of the preferred hand. We extend models previously considered by other authors in univariate settings to account for the repeated measures nature of the data. The results suggest that the expected number of attempts and successes increase with training, except for patients with advanced stages of the disease using the non-preferred hand. Copyright (c) 2008 John Wiley & Sons, Ltd.

Effects of Micelles and Vesicles on the Oximolysis of p-Nitrophenyl Diphenyl Phosphate: A Model System for Surfactant-Based Skin-defensive Formulations against Organophosphates

The rates of oximolysis of p-nitrophenyl diphenyl phosphate (PNPDPP) by Acetophenoxime; 10-phenyl-10-hydi-oxyiminodecanoic acid; 4-(9-carboxynonanyl)-1-(9-carboxy-1-hydroyiminononanyl) benzene; 1-dodecyl-2-[(hydroxyimino)methyl]-pyridinium chloride (IV) and N-methylpyridinium-2-aldoxime chloride were determined in micelles of N-hexadecyl-N,N,N-trimethylammonium chloride (CTAC), N-hexadecyl-N,N-dimethylammonium propanesulfonate and dioctadecyldimethylammonium chloride (DODAC) vesicles. The effects of CTAC micelles and DODAC vesicles on the rates of oxymolysis of O,O-Diethyl O-(4-nitrophenyl) phosphate (paraoxon) by oxime IV were also determined. Analysis of micellar and vesicular effects on oximolysis of PNPDPP, using pseudophase or pseudophase with explicit consideration of ion exchange models, required the determination of the aggregate`s effects on the pK(a), of oximes and on the rates of PNPDPP hydrolysis. All aggregates increased the rate of oximolysis of PNPDPP and the results were analyzed quantitatively. In particular, DODAC vesicles catalyzed the reaction and increased the rate of oximolysis of PNPDPP by IV several million fold at pH`s compatible with pharmaceutical formulations. The rate increase produced by DODAC vesicles on the rate of oximolysis paraoxon by IV demonstrates the pharmaceutical potential of this system, since the substrate is used as an agricultural defensive agent and the surfactant is extensively employed in cosmetic formulations. (C) 2008 Wiley-Liss, Inc. and the American Pharmacists Association J Pharm Sci 98:1040-1052, 2009

Solvatochromic and lonochromic Effects of Iron(II)bis-(1,10-phenanthroline)dicyano: a Theoretical Study

Solvatochromic and ionochromic effects of the iron(II)bis(1,10-phenanthroline)dicyano (Fe(phen)(2)(CN)(2)) complex were investigated by means of combined DFT/TDDFT calculations using the PBE and B3LYP functionals. Extended solvation models of Fe(phen)(2)(CN)(2) in acetonitrile and aqueous solution, as well as including interaction with Mg(2+), were constructed. The calculated vertical excitation energies reproduce well the observed solvatochromism in acetonitrile and aqueous solutions, the ionochromism in acetonitrile in the presence of Mg(2+), and the absence of ionochromic effect in aqueous solution. The vertical excitation energies and the nature of the transitions were reliably predicted after inclusion of geometry relaxation upon aqueous micro- and global solvation and solvent polarization effect in the TDDFT calculations. The two intense UV-vis absorption bands occurring for all systems studied are interpreted as transitions from a hybrid Fe(II)(d)/cyano N(p) orbital to a phenanthroline pi* orbital rather than a pure metal-to-ligand-charge transfer (MLCT). The solvatochromic and ionochromic blue band shifts of Fe(phen)(2)(CN)(2) were explained with preferential stabilization of the highest occupied Fe(II)(d)/cyano N(p) orbitals as a result of specific interactions with water solvent molecules or Mg(2+) ions in solution. Such interactions occur through the CN(-) groups in the complex, and they have a decisive role for the observed blue shifts of UV-vis absorption bands.