Coupled coarse graining and Markov Chain Monte Carlo for lattice systems

We propose an efficient Markov Chain Monte Carlo method for sampling equilibrium distributions for stochastic lattice models, capable of handling correctly long and short-range particle interactions. The proposed method is a Metropolis-type algorithm with the proposal probability transition matrix based on the coarse-grained approximating measures introduced in a series of works of M. Katsoulakis, A. Majda, D. Vlachos and P. Plechac, L. Rey-Bellet and D.Tsagkarogiannis,. We prove that the proposed algorithm reduces the computational cost due to energy differences and has comparable mixing properties with the classical microscopic Metropolis algorithm, controlled by the level of coarsening and reconstruction procedure. The properties and effectiveness of the algorithm are demonstrated with an exactly solvable example of a one dimensional Ising-type model, comparing efficiency of the single spin-flip Metropolis dynamics and the proposed coupled Metropolis algorithm.Comment: 20 pages, 4 figure

Semiparametric Multivariate Accelerated Failure Time Model with Generalized Estimating Equations

The semiparametric accelerated failure time model is not as widely used as the Cox relative risk model mainly due to computational difficulties. Recent developments in least squares estimation and induced smoothing estimating equations provide promising tools to make the accelerate failure time models more attractive in practice. For semiparametric multivariate accelerated failure time models, we propose a generalized estimating equation approach to account for the multivariate dependence through working correlation structures. The marginal error distributions can be either identical as in sequential event settings or different as in parallel event settings. Some regression coefficients can be shared across margins as needed. The initial estimator is a rank-based estimator with Gehan's weight, but obtained from an induced smoothing approach with computation ease. The resulting estimator is consistent and asymptotically normal, with a variance estimated through a multiplier resampling method. In a simulation study, our estimator was up to three times as efficient as the initial estimator, especially with stronger multivariate dependence and heavier censoring percentage. Two real examples demonstrate the utility of the proposed method

Multi-level evidence of an allelic hierarchy of USH2A variants in hearing, auditory processing and speech/language outcomes.

Language development builds upon a complex network of interacting subservient systems. It therefore follows that variations in, and subclinical disruptions of, these systems may have secondary effects on emergent language. In this paper, we consider the relationship between genetic variants, hearing, auditory processing and language development. We employ whole genome sequencing in a discovery family to target association and gene x environment interaction analyses in two large population cohorts; the Avon Longitudinal Study of Parents and Children (ALSPAC) and UK10K. These investigations indicate that USH2A variants are associated with altered low-frequency sound perception which, in turn, increases the risk of developmental language disorder. We further show that Ush2a heterozygote mice have low-level hearing impairments, persistent higher-order acoustic processing deficits and altered vocalizations. These findings provide new insights into the complexity of genetic mechanisms serving language development and disorders and the relationships between developmental auditory and neural systems

Plasmonically Enhanced Reflectance of Heat Radiation from Low-Bandgap Semiconductor Microinclusions

Increased reflectance from the inclusion of highly scattering particles at low volume fractions in an insulating dielectric offers a promising way to reduce radiative thermal losses at high temperatures. Here, we investigate plasmonic resonance driven enhanced scattering from microinclusions of low-bandgap semiconductors (InP, Si, Ge, PbS, InAs and Te) in an insulating composite to tailor its infrared reflectance for minimizing thermal losses from radiative transfer. To this end, we compute the spectral properties of the microcomposites using Monte Carlo modeling and compare them with results from Fresnel equations. The role of particle size-dependent Mie scattering and absorption efficiencies, and, scattering anisotropy are studied to identify the optimal microinclusion size and material parameters for maximizing the reflectance of the thermal radiation. For composites with Si and Ge microinclusions we obtain reflectance efficiencies of 57 - 65% for the incident blackbody radiation from sources at temperatures in the range 400 - 1600 {\deg}C. Furthermore, we observe a broadbanding of the reflectance spectra from the plasmonic resonances due to charge carriers generated from defect states within the semiconductor bandgap. Our results thus open up the possibility of developing efficient high-temperature thermal insulators through use of the low-bandgap semiconductor microinclusions in insulating dielectrics.Comment: Main article (8 Figures and 2 Tables) + Supporting Information (8 Figures

Recent advances in electronic structure theory and their influence on the accuracy of ab initio potential energy surfaces

Recent advances in electronic structure theory and the availability of high speed vector processors have substantially increased the accuracy of ab initio potential energy surfaces. The recently developed atomic natural orbital approach for basis set contraction has reduced both the basis set incompleteness and superposition errors in molecular calculations. Furthermore, full CI calculations can often be used to calibrate a CASSCF/MRCI approach that quantitatively accounts for the valence correlation energy. These computational advances also provide a vehicle for systematically improving the calculations and for estimating the residual error in the calculations. Calculations on selected diatomic and triatomic systems will be used to illustrate the accuracy that currently can be achieved for molecular systems. In particular, the F+H2 yields HF+H potential energy hypersurface is used to illustrate the impact of these computational advances on the calculation of potential energy surfaces

Beyond the required LISA free-fall performance: new LISA pathfinder results down to 20 μHz

In the months since the publication of the first results, the noise performance of LISA Pathfinder has improved because of reduced Brownian noise due to the continued decrease in pressure around the test masses, from a better correction of noninertial effects, and from a better calibration of the electrostatic force actuation. In addition, the availability of numerous long noise measurement runs, during which no perturbation is purposely applied to the test masses, has allowed the measurement of noise with good statistics down to 20  μHz. The Letter presents the measured differential acceleration noise figure, which is at (1.74±0.05)  fm s^{-2}/sqrt[Hz] above 2 mHz and (6±1)×10  fm s^{-2}/sqrt[Hz] at 20  μHz, and discusses the physical sources for the measured noise. This performance provides an experimental benchmark demonstrating the ability to realize the low-frequency science potential of the LISA mission, recently selected by the European Space Agency

Neutralino versus axion/axino cold dark matter in the 19 parameter SUGRA model

We calculate the relic abundance of thermally produced neutralino cold dark matter in the general 19 parameter supergravity (SUGRA-19) model. A scan over GUT scale parameters reveals that models with a bino-like neutralino typically give rise to a dark matter density \Omega_{\tz_1}h^2\sim 1-1000, i.e. between 1 and 4 orders of magnitude higher than the measured value. Models with higgsino or wino cold dark matter can yield the correct relic density, but mainly for neutralino masses around 700-1300 GeV. Models with mixed bino-wino or bino-higgsino CDM, or models with dominant co-annihilation or A-resonance annihilation can yield the correct abundance, but such cases are extremely hard to generate using a general scan over GUT scale parameters; this is indicative of high fine-tuning of the relic abundance in these cases. Requiring that m_{\tz_1}\alt 500 GeV (as a rough naturalness requirement) gives rise to a minimal probably dip in parameter space at the measured CDM abundance. For comparison, we also scan over mSUGRA space with four free parameters. Finally, we investigate the Peccei-Quinn augmented MSSM with mixed axion/axino cold dark matter. In this case, the relic abundance agrees more naturally with the measured value. In light of our cumulative results, we conclude that future axion searches should probe much more broadly in axion mass, and deeper into the axion coupling.Comment: 23 pages including 17 .eps figure

The role of input noise in transcriptional regulation

Even under constant external conditions, the expression levels of genes fluctuate. Much emphasis has been placed on the components of this noise that are due to randomness in transcription and translation; here we analyze the role of noise associated with the inputs to transcriptional regulation, the random arrival and binding of transcription factors to their target sites along the genome. This noise sets a fundamental physical limit to the reliability of genetic control, and has clear signatures, but we show that these are easily obscured by experimental limitations and even by conventional methods for plotting the variance vs. mean expression level. We argue that simple, global models of noise dominated by transcription and translation are inconsistent with the embedding of gene expression in a network of regulatory interactions. Analysis of recent experiments on transcriptional control in the early Drosophila embryo shows that these results are quantitatively consistent with the predicted signatures of input noise, and we discuss the experiments needed to test the importance of input noise more generally.Comment: 11 pages, 5 figures minor correction

The nuclear immune receptor RPS4 is required for RRS1SLH1-dependent constitutive defense activation in Arabidopsis thaliana

Plant nucleotide-binding leucine-rich repeat (NB-LRR) disease resistance (R) proteins recognize specific ‘‘avirulent’’ pathogen effectors and activate immune responses. NB-LRR proteins structurally and functionally resemble mammalian Nod-like receptors (NLRs). How NB-LRR and NLR proteins activate defense is poorly understood. The divergently transcribed Arabidopsis R genes, RPS4 (resistance to Pseudomonas syringae 4) and RRS1 (resistance to Ralstonia solanacearum 1), function together to confer recognition of Pseudomonas AvrRps4 and Ralstonia PopP2. RRS1 is the only known recessive NBLRR R gene and encodes a WRKY DNA binding domain, prompting suggestions that it acts downstream of RPS4 for transcriptional activation of defense genes. We define here the early RRS1-dependent transcriptional changes upon delivery of PopP2 via Pseudomonas type III secretion. The Arabidopsis slh1 (sensitive to low humidity 1) mutant encodes an RRS1 allele (RRS1SLH1) with a single amino acid (leucine) insertion in the WRKY DNA-binding domain. Its poor growth due to constitutive defense activation is rescued at higher temperature. Transcription profiling data indicate that RRS1SLH1-mediated defense activation overlaps substantially with AvrRps4- and PopP2-regulated responses. To better understand the genetic basis of RPS4/RRS1-dependent immunity, we performed a genetic screen to identify suppressor of slh1 immunity (sushi) mutants. We show that many sushi mutants carry mutations in RPS4, suggesting that RPS4 acts downstream or in a complex with RRS1. Interestingly, several mutations were identified in a domain C-terminal to the RPS4 LRR domain. Using an Agrobacterium-mediated transient assay system, we demonstrate that the P-loop motif of RPS4 but not of RRS1SLH1 is required for RRS1SLH1 function. We also recapitulate the dominant suppression of RRS1SLH1 defense activation by wild type RRS1 and show this suppression requires an intact RRS1 P-loop. These analyses of RRS1SLH1 shed new light on mechanisms by which NB-LRR protein pairs activate defense signaling, or are held inactive in the absence of a pathogen effector