Electromagnetic penguin operators and direct CP violation in K --> pi l^+ l^-

Supersymmetric extensions of the Standard Model predict a large enhancement of the Wilson coefficients of the dimension-five electromagnetic penguin operators affecting the direct CP violation in K_L --> pi^0 e^+ e^- and the charge asymmetry in K^\pm --> pi^\pm l^+ l^-. Here we compute the relevant matrix elements in the chiral quark model and compare these with the ones given by lattice calculationsComment: 12 pages, JHEP style, gluonic corrections to B_T adde

An approximate empirical Bayesian method for large-scale linear-Gaussian inverse problems

We study Bayesian inference methods for solving linear inverse problems, focusing on hierarchical formulations where the prior or the likelihood function depend on unspecified hyperparameters. In practice, these hyperparameters are often determined via an empirical Bayesian method that maximizes the marginal likelihood function, i.e., the probability density of the data conditional on the hyperparameters. Evaluating the marginal likelihood, however, is computationally challenging for large-scale problems. In this work, we present a method to approximately evaluate marginal likelihood functions, based on a low-rank approximation of the update from the prior covariance to the posterior covariance. We show that this approximation is optimal in a minimax sense. Moreover, we provide an efficient algorithm to implement the proposed method, based on a combination of the randomized SVD and a spectral approximation method to compute square roots of the prior covariance matrix. Several numerical examples demonstrate good performance of the proposed method

The anisotropy of granular materials

The effect of the anisotropy on the elastoplastic response of two dimensional packed samples of polygons is investigated here, using molecular dynamics simulation. We show a correlation between fabric coefficients, characterizing the anisotropy of the granular skeleton, and the anisotropy of the elastic response. We also study the anisotropy induced by shearing on the subnetwork of the sliding contacts. This anisotropy provides an explanation to some features of the plastic deformation of granular media.Comment: Submitted to PR

Fast Gibbs sampling for high-dimensional Bayesian inversion

Solving ill-posed inverse problems by Bayesian inference has recently attracted considerable attention. Compared to deterministic approaches, the probabilistic representation of the solution by the posterior distribution can be exploited to explore and quantify its uncertainties. In applications where the inverse solution is subject to further analysis procedures, this can be a significant advantage. Alongside theoretical progress, various new computational techniques allow to sample very high dimensional posterior distributions: In [Lucka2012], a Markov chain Monte Carlo (MCMC) posterior sampler was developed for linear inverse problems with $\ell_1$-type priors. In this article, we extend this single component Gibbs-type sampler to a wide range of priors used in Bayesian inversion, such as general $\ell_p^q$ priors with additional hard constraints. Besides a fast computation of the conditional, single component densities in an explicit, parameterized form, a fast, robust and exact sampling from these one-dimensional densities is key to obtain an efficient algorithm. We demonstrate that a generalization of slice sampling can utilize their specific structure for this task and illustrate the performance of the resulting slice-within-Gibbs samplers by different computed examples. These new samplers allow us to perform sample-based Bayesian inference in high-dimensional scenarios with certain priors for the first time, including the inversion of computed tomography (CT) data with the popular isotropic total variation (TV) prior.Comment: submitted to "Inverse Problems

Stress and Strain in Flat Piling of Disks

We have created a flat piling of disks in a numerical experiment using the Distinct Element Method (DEM) by depositing them under gravity. In the resulting pile, we then measured increments in stress and strain that were associated with a small decrease in gravity. We first describe the stress in terms of the strain using isotropic elasticity theory. Then, from a micro-mechanical view point, we calculate the relation between the stress and strain using the mean strain assumption. We compare the predicted values of Young's modulus and Poisson's ratio with those that were measured in the numerical experiment.Comment: 9 pages, 1 table, 8 figures, and 2 pages for captions of figure

Ab-initio calculation of Kerr spectra for semi-infinite systems including multiple reflections and optical interferences

Based on Luttinger's formulation the complex optical conductivity tensor is calculated within the framework of the spin-polarized relativistic screened Korringa-Kohn-Rostoker method for layered systems by means of a contour integration technique. For polar geometry and normal incidence ab-initio Kerr spectra of multilayer systems are then obtained by including via a 2x2 matrix technique all multiple reflections between layers and optical interferences in the layers. Applications to Co|Pt5 and Pt3|Co|Pt5 on the top of a semi-infinite fcc-Pt(111) bulk substrate show good qualitative agreement with the experimental spectra, but differ from those obtained by applying the commonly used two-media approach.Comment: 32 pages (LaTeX), 5 figures (Encapsulated PostScript), submitted to Phys. Rev.

Pore-scale Modeling of Viscous Flow and Induced Forces in Dense Sphere Packings

We propose a method for effectively upscaling incompressible viscous flow in large random polydispersed sphere packings: the emphasis of this method is on the determination of the forces applied on the solid particles by the fluid. Pore bodies and their connections are defined locally through a regular Delaunay triangulation of the packings. Viscous flow equations are upscaled at the pore level, and approximated with a finite volume numerical scheme. We compare numerical simulations of the proposed method to detailed finite element (FEM) simulations of the Stokes equations for assemblies of 8 to 200 spheres. A good agreement is found both in terms of forces exerted on the solid particles and effective permeability coefficients

Inverse Modeling for MEG/EEG data

We provide an overview of the state-of-the-art for mathematical methods that are used to reconstruct brain activity from neurophysiological data. After a brief introduction on the mathematics of the forward problem, we discuss standard and recently proposed regularization methods, as well as Monte Carlo techniques for Bayesian inference. We classify the inverse methods based on the underlying source model, and discuss advantages and disadvantages. Finally we describe an application to the pre-surgical evaluation of epileptic patients.Comment: 15 pages, 1 figur

Recent NA48/2 and NA62 results

The NA48/2 Collaboration at CERN has accumulated and analysed unprecedented statistics of rare kaon decays in the $K_{e4}$ modes: $K_{e4}(+-)$ ($K^\pm \to \pi^+ \pi^- e^\pm \nu$) and $K_{e4}(00)$ ($K^\pm \to \pi^0 \pi^0 e^\pm \nu$) with nearly one percent background contamination. It leads to the improved measurement of branching fractions and detailed form factor studies. New final results from the analysis of 381 $K^\pm \to \pi^\pm \gamma \gamma$ rare decay candidates collected by the NA48/2 and NA62 experiments at CERN are presented. The results include a decay rate measurement and fits to Chiral Perturbation Theory (ChPT) description.Comment: Prepared for the Proceedings of "Moriond QCD and High Energy Interactions. March 22-29 2014." conferenc