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

A computational model integrating brain electrophysiology and metabolism highlights the key role of extracellular potassium and oxygen

The human brain is a small organ which uses a disproportional amount of the total metabolic energy pro- duction in the body. While it is well understood that the most significant energy sink is the maintenance of the neuronal membrane potential during the brain signaling activity, the role of astrocytes in the energy balance continues to be the topic of a lot of research. A key function of astrocytes, besides clearing glutamate from the synaptic clefts, is the potassium clearing after neuronal activation. Extracellular potassium plays a significant role in triggering neuronal firing, and elevated concentration of potassium may lead to abnormal firing pattern, e.g., seizures, thus emphasizing the importance of the glial K+ buffering role. The predictive mathematical model proposed in this paper elucidates the role of glial potassium clearing in brain energy metabolism, integrating a detailed model of the ion dynamics which regulates neuronal firing with a three compartment metabolic model. Because of the very different characteristic time scales of electrophysiology and metabolism, care must be taken when coupling the two models to ensure that the predictions, e.g., neuronal firing frequencies and the oxygen- glucose index (OGI) of the brain during activation and rest, are in agreement with empirical observations. The temporal multi-scale nature of the problem requires the design of new computational tools to ensure a stable and accurate numerical treatment of the problem. The model predictions for different protocols, including combinations of elevated activation and ischemic episodes, are in good agreement with experimental observations reported in the literature.This work was supported by the Bizkaia Talent and European Commission through CO- FUND under the grant CIPAS: Computational Inverse Problems Across Scales (AYD-000-278, 2015), by the Basque Government through the BERC 2014-2017 program, and by the Spanish Ministry of Economics and Competitive- ness MINECO through the BCAM Severo Ochoa excellence accreditation SEV-2013-0323 and the Spanish ”Plan Estatal de Investigacio ́n, Desarrollo e Innovacio ́n Orientada a los Retos de la Sociedad” under Grant BELEMET - Brain ELEctro-METabolic modeling and numerical approximation (MTM2015-69992-R). The work of Daniela Cal- vetti was partly supported by Grant Number 246665 from the Simons Foundation, and the work of Erkki Somersalo was partly supported by NSF Grant DMS 1016183. Daniela Calvetti and Erkki Somersalo were partly supported by NIH, grant 1U01GM111251-01

Brain energetics plays a key role in the coordination of electrophysiology, metabolism and hemodynamics: evidence from an integrated computational model

The energetic needs of brain cells at rest and during elevated neuronal activation has been the topic of many investigations where mathematical models have played a significant role providing a context for the interpretation of experimental findings. A recently proposed mathematical model, comprising a double feedback between cellular metabolism and electrophysiology, sheds light on the interconnections between the electrophysiological details associated with changes in the frequency of neuronal firing and the corresponding metabolic activity. We propose a new extended mathematical model comprising a three-way feedback connecting metabolism, electrophysiology and hemodynamics. Upon specifying the time intervals of higher neuronal activation, the model generates a potassium based signal leading to the concomitant increase in cerebral blood flow with associated vasodilation and metabolic changes needed to sustain the increased energy demand. The predictions of the model are in good qualitative and quantitative agreement with experimental findings reported in the literature, even predicting a slow after-hyperpolarization of a duration of approximately 16 s matching experimental observations.The work of Daniela Calvetti was partly support by NSF grants DMS-1522334 and NIH grant 1U01 GM111251-01. The work of Erkki Somersalo was partly support by NSF grants DMS 1714617 and NIH grant 1U01GM111251-01

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

VARIAN CLINAC 6 MeV Photon Spectra Unfolding using a Monte Carlo Meshed Model

[EN] Energy spectrum is the best descriptive function to determine photon beam quality of a Medical Linear Accelerator (LinAc). The use of realistic photon spectra in Monte Carlo simulations has a great importance to obtain precise dose calculations in Radiotherapy Treatment Planning (RTP). Reconstruction of photon spectra emitted by medical accelerators from measured depth dose distributions in a water cube is an important tool for commissioning a Monte Carlo treatment planning system. Regarding this, the reconstruction problem is an inverse radiation transport function which is ill conditioned and its solution may become unstable due to small perturbations in the input data. This paper presents a more stable spectral reconstruction method which can be used to provide an independent confirmation of source models for a given machine without any prior knowledge of the spectral distribution. Monte Carlo models used in this work are built with unstructured meshes to simulate with realism the linear accelerator head geometry.Ha recibido financiación de DESARROLLO E IMPLEMENTACION DE UN SISTEMA DE PLANIFICACION RADIOTERAPÉUTICA (TELETERAPIA Y BRAQUITERAPIA) BASADO EN EL MÉTODO MONTECARLO PARA EL HUPLAFE (UNIVERSIDAD POLITECNICA DE VALENCIA), DESARROLLO Y VALIDACION EXPERIMENTAL DE UN SISTEMA DE PLANIFICACION DE TRATAMIENTOS RADIOTERAPEUTICOS UTILIZANDO TECNICAS DE MONTE CARLO (UNIVERSIDAD POLITECNICA DE VALENCIA), Dosimetría in vivo en radioterapia sin filtro aplanador (FFF) (UNIVERSIDAD POLITECNICA DE VALENCIA), MONTE CARLO TREATMENT PLANNING SYSTEM: SOFTWARE PARA EL CALCULO DOSIMETRICO DE ALTA PRECISION EN RADIOTERAPIA. (UNIVERSIDAD POLITECNICA DE VALENCIA) y SISTEMA DE PLANIFICACION DE TRATAMIENTOS EN RADIOTERAPIA BASADO EN EL METODO DE MOTE CARLO Y EL PROCESADO 3D DE IMAGENES MEDICAS (GENERALITAT VALENCIANA)Morató-Rafet, S.; Juste Vidal, BJ.; Miró Herrero, R.; Verdú Martín, GJ. (2017). VARIAN CLINAC 6 MeV Photon Spectra Unfolding using a Monte Carlo Meshed Model. EPJ Web of Conferences. 153(04012):1-7. https://doi.org/10.1051/epjconf/201715304012S171530401

Nab-paclitaxel in substitution of oxaliplatin and irinotecan in folfirinox schedule as first-line therapy in patients with metastatic pancreatic cancer: results of phase I dose finding of NabucCO study by GOIRC

n/

Block matrix formulations for evolving networks

Many types of pairwise interaction take the form of a fixed set of nodes with edges that appear and disappear over time. In the case of discrete-time evolution, the resulting evolving network may be represented by a time-ordered sequence of adjacency matrices. We consider here the issue of representing the system as a single, higher dimensional block matrix, built from the individual time-slices. We focus on the task of computing network centrality measures. From a modeling perspective, we show that there is a suitable block formulation that allows us to recover dynamic centrality measures respecting time's arrow. From a computational perspective, we show that the new block formulation leads to the design of more effective numerical algorithms.Comment: 18 pages, 2 figure

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