A new ghost cell/level set method for moving boundary problems:application to tumor growth

In this paper, we present a ghost cell/level set method for the evolution of interfaces whose normal velocity depend upon the solutions of linear and nonlinear quasi-steady reaction-diffusion equations with curvature-dependent boundary conditions. Our technique includes a ghost cell method that accurately discretizes normal derivative jump boundary conditions without smearing jumps in the tangential derivative; a new iterative method for solving linear and nonlinear quasi-steady reaction-diffusion equations; an adaptive discretization to compute the curvature and normal vectors; and a new discrete approximation to the Heaviside function. We present numerical examples that demonstrate better than 1.5-order convergence for problems where traditional ghost cell methods either fail to converge or attain at best sub-linear accuracy. We apply our techniques to a model of tumor growth in complex, heterogeneous tissues that consists of a nonlinear nutrient equation and a pressure equation with geometry-dependent jump boundary conditions. We simulate the growth of glioblastoma (an aggressive brain tumor) into a large, 1 cm square of brain tissue that includes heterogeneous nutrient delivery and varied biomechanical characteristics (white matter, gray matter, cerebrospinal fluid, and bone), and we observe growth morphologies that are highly dependent upon the variations of the tissue characteristics—an effect observed in real tumor growth

Cardiolipin provides an essential activating platform for caspase-8 on mitochondria

Cardiolipin is a mitochondria-specific phospholipid known to be intimately involved with apoptosis. However, the lack of appropriate cellular models to date restricted analysis of its role in cell death. The maturation of cardiolipin requires the transacylase tafazzin, which is mutated in the human disorder Barth syndrome. Using Barth syndrome patient-derived cells and HeLa cells in which tafazzin was knocked down, we show that cardiolipin is required for apoptosis in the type II mitochondria-dependent response to Fas stimulation. Cardiolipin provides an anchor and activating platform for caspase-8 translocation to, and embedding in, the mitochondrial membrane, where it oligomerizes and is further activated, steps that are necessary for an efficient type II apoptotic response

Numerical Portrait of a Relativistic Thin Film BCS Superfluid

We present results of numerical simulations of the 2+1d Nambu - Jona-Lasinio model with a non-zero baryon chemical potential mu including the effects of a diquark source term. Diquark condensates, susceptibilities and masses are measured as functions of source strength j. The results suggest that diquark condensation does not take place in the high density phase mu>mu_c, but rather that the condensate scales non-analytically with j implying a line of critical points and long range phase coherence. Analogies are drawn with the low temperature phase of the 2d XY model. The spectrum of the spin-1/2 sector is also studied yielding the quasiparticle dispersion relation. There is no evidence for a non-zero gap; rather the results are characteristic of a normal Fermi liquid with Fermi velocity less than that of light. We conclude that the high density phase of the model describes a relativistic gapless thin film BCS superfluid.Comment: 37 pages, 16 figure

Branch-and-lift algorithm for deterministic global optimization in nonlinear optimal control

This paper presents a branch-and-lift algorithm for solving optimal control problems with smooth nonlinear dynamics and potentially nonconvex objective and constraint functionals to guaranteed global optimality. This algorithm features a direct sequential method and builds upon a generic, spatial branch-and-bound algorithm. A new operation, called lifting, is introduced, which refines the control parameterization via a Gram-Schmidt orthogonalization process, while simultaneously eliminating control subregions that are either infeasible or that provably cannot contain any global optima. Conditions are given under which the image of the control parameterization error in the state space contracts exponentially as the parameterization order is increased, thereby making the lifting operation efficient. A computational technique based on ellipsoidal calculus is also developed that satisfies these conditions. The practical applicability of branch-and-lift is illustrated in a numerical example. © 2013 Springer Science+Business Media New York

Machine Learning Framework to Identify Individuals at Risk of Rapid Progression of Coronary Atherosclerosis : From the PARADIGM Registry

Background Rapid coronary plaque progression (RPP) is associated with incident cardiovascular events. To date, no method exists for the identification of individuals at risk of RPP at a single point in time. This study integrated coronary computed tomography angiography-determined qualitative and quantitative plaque features within a machine learning (ML) framework to determine its performance for predicting RPP. Methods and Results Qualitative and quantitative coronary computed tomography angiography plaque characterization was performed in 1083 patients who underwent serial coronary computed tomography angiography from the PARADIGM (Progression of Atherosclerotic Plaque Determined by Computed Tomographic Angiography Imaging) registry. RPP was defined as an annual progression of percentage atheroma volume 651.0%. We employed the following ML models: model 1, clinical variables; model 2, model 1 plus qualitative plaque features; model 3, model 2 plus quantitative plaque features. ML models were compared with the atherosclerotic cardiovascular disease risk score, Duke coronary artery disease score, and a logistic regression statistical model. 224 patients (21%) were identified as RPP. Feature selection in ML identifies that quantitative computed tomography variables were higher-ranking features, followed by qualitative computed tomography variables and clinical/laboratory variables. ML model 3 exhibited the highest discriminatory performance to identify individuals who would experience RPP when compared with atherosclerotic cardiovascular disease risk score, the other ML models, and the statistical model (area under the receiver operating characteristic curve in ML model 3, 0.83 [95% CI 0.78-0.89], versus atherosclerotic cardiovascular disease risk score, 0.60 [0.52-0.67]; Duke coronary artery disease score, 0.74 [0.68-0.79]; ML model 1, 0.62 [0.55-0.69]; ML model 2, 0.73 [0.67-0.80]; all P<0.001; statistical model, 0.81 [0.75-0.87], P=0.128). Conclusions Based on a ML framework, quantitative atherosclerosis characterization has been shown to be the most important feature when compared with clinical, laboratory, and qualitative measures in identifying patients at risk of RPP

Quantum walks: a comprehensive review

Quantum walks, the quantum mechanical counterpart of classical random walks, is an advanced tool for building quantum algorithms that has been recently shown to constitute a universal model of quantum computation. Quantum walks is now a solid field of research of quantum computation full of exciting open problems for physicists, computer scientists, mathematicians and engineers. In this paper we review theoretical advances on the foundations of both discrete- and continuous-time quantum walks, together with the role that randomness plays in quantum walks, the connections between the mathematical models of coined discrete quantum walks and continuous quantum walks, the quantumness of quantum walks, a summary of papers published on discrete quantum walks and entanglement as well as a succinct review of experimental proposals and realizations of discrete-time quantum walks. Furthermore, we have reviewed several algorithms based on both discrete- and continuous-time quantum walks as well as a most important result: the computational universality of both continuous- and discrete- time quantum walks.Comment: Paper accepted for publication in Quantum Information Processing Journa

Edge Detection by Adaptive Splitting II. The Three-Dimensional Case

In Llanas and Lantarón, J. Sci. Comput. 46, 485–518 (2011) we proposed an algorithm (EDAS-d) to approximate the jump discontinuity set of functions defined on subsets of ℝ d . This procedure is based on adaptive splitting of the domain of the function guided by the value of an average integral. The above study was limited to the 1D and 2D versions of the algorithm. In this paper we address the three-dimensional problem. We prove an integral inequality (in the case d=3) which constitutes the basis of EDAS-3. We have performed detailed computational experiments demonstrating effective edge detection in 3D function models with different interface topologies. EDAS-1 and EDAS-2 appealing properties are extensible to the 3D cas

Chebyshev Solution of the Nearly-Singular One-Dimensional Helmholtz Equation and Related Singular Perturbation Equations: Multiple Scale Series and the Boundary Layer Rule-of-Thumb

The one-dimensional Helmholtz equation, ε 2 u xx − u = f ( x ), arises in many applications, often as a component of three-dimensional fluids codes. Unfortunately, it is difficult to solve for ε≪1 because the homogeneous solutions are exp (± x /ε), which have boundary layers of thickness O(1/ε). By analyzing the asymptotic Chebyshev coefficients of exponentials, we rederive the Orszag–Israeli rule [16] that Chebyshev polynomials are needed to obtain an accuracy of 1% or better for the homogeneous solutions. (Interestingly, this is identical with the boundary layer rule-of-thumb in [5], which was derived for singular functions like tanh([ x −1]/ε).) Two strategies for small ε are described. The first is the method of multiple scales, which is very general, and applies to variable coefficient differential equations, too. The second, when f ( x ) is a polynomial, is to compute an exact particular integral of the Helmholtz equation as a polynomial of the same degree in the form of a Chebyshev series by solving triangular pentadiagonal systems. This can be combined with the analytic homogeneous solutions to synthesize the general solution. However, the multiple scales method is more efficient than the Chebyshev algorithm when ε is very, very tiny.Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/45436/1/11075_2004_Article_2865.pd

A precise measurement of the magnetic field in the corona of the black hole binary V404 Cygni

Observations of binary stars containing an accreting black hole or neutron star often show x-ray emission extending to high energies (>10 kilo­–electron volts), which is ascribed to an accretion disk corona of energetic particles akin to those seen in the solar corona. Despite their ubiquity, the physical conditions in accretion disk coronae remain poorly constrained. Using simultaneous infrared, optical, x-ray, and radio observations of the Galactic black hole system V404 Cygni, showing a rapid synchrotron cooling event in its 2015 outburst, we present a precise 461 ± 12 gauss magnetic field measurement in the corona. This measurement is substantially lower than previous estimates for such systems, providing constraints on physical models of accretion physics in black hole and neutron star binary systems. This article has a correction. Please see: http://science.sciencemag.org/content/360/6386/eaat927

Effect of nesiritide in patients with acute decompensated heart failure

Background Nesiritide is approved in the United States for early relief of dyspnea in patients with acute heart failure. Previous meta-analyses have raised questions regarding renal toxicity and the mortality associated with this agent. Methods We randomly assigned 7141 patients who were hospitalized with acute heart failure to receive either nesiritide or placebo for 24 to 168 hours in addition to standard care. Coprimary end points were the change in dyspnea at 6 and 24 hours, as measured on a 7-point Likert scale, and the composite end point of rehospitalization for heart failure or death within 30 days. Results Patients randomly assigned to nesiritide, as compared with those assigned to placebo, more frequently reported markedly or moderately improved dyspnea at 6 hours (44.5% vs. 42.1%, P = 0.03) and 24 hours (68.2% vs. 66.1%, P = 0.007), but the prespecified level for significance (P≤0.005 for both assessments or P≤0.0025 for either) was not met. The rate of rehospitalization for heart failure or death from any cause within 30 days was 9.4% in the nesiritide group versus 10.1% in the placebo group (absolute difference, −0.7 percentage points; 95% confidence interval [CI], −2.1 to 0.7; P = 0.31). There were no significant differences in rates of death from any cause at 30 days (3.6% with nesiritide vs. 4.0% with placebo; absolute difference, −0.4 percentage points; 95% CI, −1.3 to 0.5) or rates of worsening renal function, defined by more than a 25% decrease in the estimated glomerular filtration rate (31.4% vs. 29.5%; odds ratio, 1.09; 95% CI, 0.98 to 1.21; P = 0.11). Conclusions Nesiritide was not associated with an increase or a decrease in the rate of death and rehospitalization and had a small, nonsignificant effect on dyspnea when used in combination with other therapies. It was not associated with a worsening of renal function, but it was associated with an increase in rates of hypotension. On the basis of these results, nesiritide cannot be recommended for routine use in the broad population of patients with acute heart failure. (Funded by Scios; ClinicalTrials.gov number, NCT00475852.