Design of optimal RF pulses for NMR as a discrete-valued control problem
|
|
- Daniela Hicks
- 5 years ago
- Views:
Transcription
1 Design of optimal RF pulses for NMR as a discrete-valued control problem Christian Clason Faculty of Mathematics, Universität Duisburg-Essen joint work with Carla Tameling (Göttingen) and Benedikt Wirth (Münster) 9th Vienna International Conference on Mathematical Modelling Wien, February 2, 28 / 9
2 Motivation: MRT Magnetic resonance tomography (MRT): based on measurement of resonance frequency of hydrogen nuclei hydrogen nucleus (proton) acts like rotating magnet H N angular momentum : nuclear spin, axis randomly aligned S 2 / 9
3 Motivation: MRT nuclear spin in magnetic field: axis aligns with external field precesses around field direction 3 / 9
4 Motivation: MRT nuclear spin in rotating magnetic field: axis aligns with external field precesses around field direction field rotating at resonance frequency: energy is absorbed resonance frequency proportional to field strength 3 / 9
5 Motivation: MRT magnetic field off: spin relaxes......induces current in coil measurement 4 / 9
6 Mathematical model Bloch equation d dt M(t) = M(t) B(t), M() = M M(t) R 3 describes temporal evolution of spin ensemble B(t) = (u (t), u 2 (t), ω) T controlled time-dependent magnetic field ω resonance frequency control-to-state mapping S (ω) : u M(T) 5 / 9
7 Optimal control problem Goal: compute control u(t) = (u (t), u 2 (t)) such that M(T) M d M d desired magnetization state M d = M (ω) d selective to resonance frequency ( spectroscopy, slice selection) in addition: control with minimal specific absorption rate (SAR) min u L 2 2 ω S (ω) (u) M (ω) d α T u(t) 2 2 dt 2 6 / 9
8 Discrete control Technical limitation: device can only realize control from discrete set U = { u L 2 (, T; R 2 ) : u(t) {u,..., u d } a.e. } u,..., u d R 2 given (fixed amplitude, phases) non-convex discrete-valued control problem min u U 2 ω S (ω) (u) M (ω) d α T u(t) 2 2 dt 2 7 / 9
9 Multi-bang penalty convex relaxation: replace U by convex hull works only for d = 2, cf. bang-bang control (α = ) promote u(x) {u,..., u d } by convex pointwise penalty G(u) = g(u(x)) dx Ω generalize L norm: polyhedral epigraph with vertices u,..., u d not exact relaxation/penalization (in general)! 8 / 9
10 Multi-bang penalty generalize L norm: polyhedral epigraph with vertices u,..., u d 3 2 u u 2 u 3 v motivation: convex envelope of 2 u 2 + δ U multi-bang (generalized bang-bang) control non-smooth optimization in function spaces 8 / 9
11 Motivation 2 Approach Convex analysis Moreau Yosida regularization Semismooth Newton method 3 Vector-valued multi-bang 4 Numerical examples 9 / 9
12 Fenchel duality f : V R := R {+ } convex, V Banach space, V dual space subdifferential f( v) = { v V : v, v v V,V f(v) f( v) for all v V } Fenchel conjugate (always convex) f : V R, f (v ) = sup v, v V,V f(v) v V convex inverse function theorem : v f(v) v f (v ) / 9
13 Fenchel duality: application F(ū) + G(ū) = min F(u) + G(u) u Fermat principle: (F(ū) + G(ū)) 2 sum rule: F(ū) + G(ū), i.e., there is p V with { p F(ū) 3 Fenchel duality: p G(ū) { p F(ū) ū G ( p) / 9
14 Regularization G non-smooth subdifferential G set-valued regularize u, p L 2 (Ω) Hilbert space consider for γ > Proximal mapping prox γg (p) = arg min G (w) + w p 2 w 2γ single-valued, Lipschitz continuous coincides with resolvent (Id +γ G ) (p) (also required for primal-dual first-order methods) 2 / 9
15 Regularization Proximal mapping prox γg (p) = arg min G (w) + w p 2 w 2γ Complementarity formulation of u G (p) u = γ ( (p + γu) proxγg (p + γu) ) equivalent for every γ > single-valued, Lipschitz continuous, implicit 2 / 9
16 Regularization Proximal mapping prox γg (p) = arg min G (w) + w p 2 w 2γ Moreau Yosida regularization of u G (p) u = γ ( p proxγg (p) ) =: G γ (p) G γ = ( G + γ 2 2) G as γ single-valued, Lipschitz continuous, explicit nonsmooth operator equation, Newton method 2 / 9
17 Semismooth Newton method f locally Lipschitz, piecewise C : f(v) =, f : R n R Newton derivative D N f(v)δv C f(v)δv Clarke generalized gradient: convex hull of piecewise derivatives semismooth Newton method D N f(v k )δv = f(v k ), v k+ = v k + δv converges locally superlinearly 3 / 9
18 Semismooth Newton method f locally Lipschitz, piecewise C : F(u) =, F : L p (Ω) L q (Ω), [F(u)](x) = f(u(x)) Newton derivative [D N F(u)δu](x) C f(δu(x))δu(x) any measurable selection of Clarke generalized gradient semismooth Newton method D N F(u k )δu = F(u k ), u k+ = u k + δu converges locally superlinearly if p > q 3 / 9
19 Vector-valued multi-bang: penalty Here: admissible control set U of d radially distributed states, origin { ( U = ) ( ) ( )}, ω cos θ ω sin θ,..., ω cos θ d ω sin θ d fixed amplitude ω > phases θ <... < θ d < 2π multi-bang penalty g = ( δ U) convex envelope ( ( g (q) = 2 ) ) ( δ U (q) = 2 ) δ U (q) { q, u i = 2 ω2 for all i d q, u i 2 ω2 θ i +θ i 2 q θ i+θ i+ 2, q, u i 2 ω2 4 / 9
20 Vector-valued multi-bang: subdifferential Fenchel conjugate g (q) = { =: u q Q q, u i 2 ω2 q Q i Subdifferential { g {u i } q Q i i d (q) = co {u i,..., u ik } q Q i...i k i,..., i k d 5 / 9
21 Vector-valued multi-bang: subdifferential Subdifferential { g {u i } q Q i i d (q) = co {u i,..., u ik } q Q i...i k i,..., i k d Moreau Yosida regularization ( g ) γ (q) = u i ( q,ui α γω 2 2γ u i +u i+ 2 + q,u i u i+ (u i u i+ ) γ u i u i+ 2 2 q γ α γ q Q γ i ) u i q Q γ,i q Q γ i,i+ ( ω u i +u i+ 2 ) 2 (ui + u i+ ) q Q γ,i,i+ 5 / 9
22 Vector-valued multi-bang: subdifferential q 2 q 2 Q 3 Q 2 Q 2 Q 23 Q 2 Q 2 Q 23 Q Q 6 Q 3 Q Q 6 Q 4 Q 45 Q 5 Q Q 34 Q 34 Q 4 Q 56 Q 5 Q 45 Q 56 Q 6 Q 6 q Q Q γ γ 2 2 Q γ Q γ Q γ 23 2 Q γ 2 Q γ Q γ 23 Q γ Q γ 6 Q γ 3 Q γ 3 Q γ 6 Q γ q Q γ 6 Q γ Q γ Q γ 4 Q γ Q γ Q γ Q γ Q γ 56 4 Q γ Q γ 45 5 (a) subdomains for g (b) subdomains for ( g ) γ 5 / 9
23 Vector-valued multi-bang: Newton method Newton derivative q Q γ i u i u T i q Q γ D N ( g γω γ )(q) = 2,i (u i u i+ )(u i u i+ ) T q Q γ γ u i u i+ 2 i,i+ 2 γ Id q Qγ,i,i+ Superposition operator: [ DN H γ (p) ] (t) := D N ( g γ ) (p(t)) a.e. t [, T] 6 / 9
24 Vector-valued multi-bang: Newton method Newton system ( ) Id D N H γ (F (u k ))F (u k ) δu = u k + G γ (F (u k )) matrix-free Krylov method for semismooth Newton step F, F via linearized, adjoint Bloch equation discretization, adjoint from [Aigner/Clason/Rund/Stollberger 6] 6 / 9
25 Numerical examples goal: shift magnetization from M = (,, ) T at t = to M d = (,, ) T at t = T d = 3, 6 radially distributed admissible control states n =, 4 isochromats with different resonance frequencies shift all isochromats 2 shift only one isochromat α =, ω = example motivated by [Dridi/Lapert/Salomon/Glaser/Sugny 5] 7 / 9
26 Numerical examples M().5 u2 Mz.5.5 u t (a) control u(t) Mx.5 M(T) M d.5 My (b) state M(t).5 Figure: n = isochromat, d = 3 control states 8 / 9
27 Numerical examples M().5 u2 Mz.5.5 u t (a) control u(t) Mx.5 M(T) M d.5 My (b) state M(t).5 Figure: n = isochromat, d = 6 control states 8 / 9
28 Numerical examples M().5 u2 Mz.5.5 u t (a) control u(t) Mx.5 M(T) M d.5 My (b) state M(t).5 Figure: n = 4 isochromats, same target 8 / 9
29 Numerical examples M() M(T).5 u2 Mz.5.5 u t (a) control u(t) Mx.5 M(T) M d.5 My (b) state M(t).5 Figure: J = 4 isochromats, different targets 8 / 9
30 Conclusion Discrete controls in NMR: can be promoted by convex penalties linear complexity in number of parameter values efficient numerical solution (superlinear convergence) applicable to nonlinear, vector-valued problems Outlook: robust optimization, e.g., with respect to static field include signal acquisition other discrete continuous problems: switching, networks Preprint, MATLAB/Python codes: 9 / 9
ruhr.pad Vector-valued multibang control of differential equations C. Clason, C. Tameling and B. Wirth Preprint
Bild ohne Schrift: ruhr.pad UA Ruhr Zentrum für partielle Differentialgleichungen Bild mit Schrift: ruhr.pad UA Ruhr Zentrum für partielle Differentialgleichungen Vector-valued multibang control of differential
More informationAdaptive discretization and first-order methods for nonsmooth inverse problems for PDEs
Adaptive discretization and first-order methods for nonsmooth inverse problems for PDEs Christian Clason Faculty of Mathematics, Universität Duisburg-Essen joint work with Barbara Kaltenbacher, Tuomo Valkonen,
More informationThe NMR Inverse Imaging Problem
The NMR Inverse Imaging Problem Nuclear Magnetic Resonance Protons and Neutrons have intrinsic angular momentum Atoms with an odd number of proton and/or odd number of neutrons have a net magnetic moment=>
More informationFundamental MRI Principles Module 2 N. Nuclear Magnetic Resonance. X-ray. MRI Hydrogen Protons. Page 1. Electrons
Fundamental MRI Principles Module 2 N S 1 Nuclear Magnetic Resonance There are three main subatomic particles: protons positively charged neutrons no significant charge electrons negatively charged Protons
More informationMagnetic resonance imaging MRI
Magnetic resonance imaging MRI Introduction What is MRI MRI is an imaging technique used primarily in medical settings that uses a strong magnetic field and radio waves to produce very clear and detailed
More informationA Brief Introduction to Medical Imaging. Outline
A Brief Introduction to Medical Imaging Outline General Goals Linear Imaging Systems An Example, The Pin Hole Camera Radiations and Their Interactions with Matter Coherent vs. Incoherent Imaging Length
More informationA duality-based approach to elliptic control problems in non-reflexive Banach spaces
A duality-based approach to elliptic control problems in non-reflexive Banach spaces Christian Clason Karl Kunisch June 3, 29 Convex duality is a powerful framework for solving non-smooth optimal control
More informationON A CLASS OF NONSMOOTH COMPOSITE FUNCTIONS
MATHEMATICS OF OPERATIONS RESEARCH Vol. 28, No. 4, November 2003, pp. 677 692 Printed in U.S.A. ON A CLASS OF NONSMOOTH COMPOSITE FUNCTIONS ALEXANDER SHAPIRO We discuss in this paper a class of nonsmooth
More informationConvex Optimization Conjugate, Subdifferential, Proximation
1 Lecture Notes, HCI, 3.11.211 Chapter 6 Convex Optimization Conjugate, Subdifferential, Proximation Bastian Goldlücke Computer Vision Group Technical University of Munich 2 Bastian Goldlücke Overview
More informationA convex analysis approach to multi-material topology optimization
SpezialForschungsBereich F 3 Karl Franzens Universität Graz Technische Universität Graz Medizinische Universität Graz A convex analysis approach to multi-material topology optimization C. Clason K. Kunisch
More informationProximal splitting methods on convex problems with a quadratic term: Relax!
Proximal splitting methods on convex problems with a quadratic term: Relax! The slides I presented with added comments Laurent Condat GIPSA-lab, Univ. Grenoble Alpes, France Workshop BASP Frontiers, Jan.
More informationFundamental MRI Principles Module Two
Fundamental MRI Principles Module Two 1 Nuclear Magnetic Resonance There are three main subatomic particles: protons neutrons electrons positively charged no significant charge negatively charged Protons
More informationChristian Clason * Karl Kunisch
a convex analysis approach to multi-material topology optimization Christian Clason * Karl Kunisch June 5, 2015 This work is concerned with optimal control of partial differential equations where the control
More informationPrinciples of Magnetic Resonance Imaging
Principles of Magnetic Resonance Imaging Hi Klaus Scheffler, PhD Radiological Physics University of 1 Biomedical Magnetic Resonance: 1 Introduction Magnetic Resonance Imaging Contents: Hi 1 Introduction
More informationClassical behavior of magnetic dipole vector. P. J. Grandinetti
Classical behavior of magnetic dipole vector Z μ Y X Z μ Y X Quantum behavior of magnetic dipole vector Random sample of spin 1/2 nuclei measure μ z μ z = + γ h/2 group μ z = γ h/2 group Quantum behavior
More informationIntroduction to Biomedical Imaging
Alejandro Frangi, PhD Computational Imaging Lab Department of Information & Communication Technology Pompeu Fabra University www.cilab.upf.edu MRI advantages Superior soft-tissue contrast Depends on among
More informationAffine covariant Semi-smooth Newton in function space
Affine covariant Semi-smooth Newton in function space Anton Schiela March 14, 2018 These are lecture notes of my talks given for the Winter School Modern Methods in Nonsmooth Optimization that was held
More informationSpectral Broadening Mechanisms
Spectral Broadening Mechanisms Lorentzian broadening (Homogeneous) Gaussian broadening (Inhomogeneous, Inertial) Doppler broadening (special case for gas phase) The Fourier Transform NC State University
More informationPhysical fundamentals of magnetic resonance imaging
Physical fundamentals of magnetic resonance imaging Stepan Sereda University of Bonn 1 / 26 Why? Figure 1 : Full body MRI scan (Source: [4]) 2 / 26 Overview Spin angular momentum Rotating frame and interaction
More informationMRI Physics I: Spins, Excitation, Relaxation
MRI Physics I: Spins, Excitation, Relaxation Douglas C. Noll Biomedical Engineering University of Michigan Michigan Functional MRI Laboratory Outline Introduction to Nuclear Magnetic Resonance Imaging
More informationA Hands on Introduction to NMR Lecture #1 Nuclear Spin and Magnetic Resonance
A Hands on Introduction to NMR 22.920 Lecture #1 Nuclear Spin and Magnetic Resonance Introduction - The aim of this short course is to present a physical picture of the basic principles of Nuclear Magnetic
More informationON PROXIMAL POINT-TYPE ALGORITHMS FOR WEAKLY CONVEX FUNCTIONS AND THEIR CONNECTION TO THE BACKWARD EULER METHOD
ON PROXIMAL POINT-TYPE ALGORITHMS FOR WEAKLY CONVEX FUNCTIONS AND THEIR CONNECTION TO THE BACKWARD EULER METHOD TIM HOHEISEL, MAXIME LABORDE, AND ADAM OBERMAN Abstract. In this article we study the connection
More informationAlgorithms for Nonsmooth Optimization
Algorithms for Nonsmooth Optimization Frank E. Curtis, Lehigh University presented at Center for Optimization and Statistical Learning, Northwestern University 2 March 2018 Algorithms for Nonsmooth Optimization
More informationNuclear Magnetic Resonance Imaging
Nuclear Magnetic Resonance Imaging Simon Lacoste-Julien Electromagnetic Theory Project 198-562B Department of Physics McGill University April 21 2003 Abstract This paper gives an elementary introduction
More information1.The anisotropic plate model
Pré-Publicações do Departamento de Matemática Universidade de Coimbra Preprint Number 6 OPTIMAL CONTROL OF PIEZOELECTRIC ANISOTROPIC PLATES ISABEL NARRA FIGUEIREDO AND GEORG STADLER ABSTRACT: This paper
More informationBayesian inverse problems with Laplacian noise
Bayesian inverse problems with Laplacian noise Remo Kretschmann Faculty of Mathematics, University of Duisburg-Essen Applied Inverse Problems 2017, M27 Hangzhou, 1 June 2017 1 / 33 Outline 1 Inverse heat
More informationClassical Description of NMR Parameters: The Bloch Equations
Classical Description of NMR Parameters: The Bloch Equations Pascale Legault Département de Biochimie Université de Montréal 1 Outline 1) Classical Behavior of Magnetic Nuclei: The Bloch Equation 2) Precession
More informationREGULAR LAGRANGE MULTIPLIERS FOR CONTROL PROBLEMS WITH MIXED POINTWISE CONTROL-STATE CONSTRAINTS
REGULAR LAGRANGE MULTIPLIERS FOR CONTROL PROBLEMS WITH MIXED POINTWISE CONTROL-STATE CONSTRAINTS fredi tröltzsch 1 Abstract. A class of quadratic optimization problems in Hilbert spaces is considered,
More informationNuclear Magnetic Resonance Imaging
Nuclear Magnetic Resonance Imaging Jeffrey A. Fessler EECS Department The University of Michigan NSS-MIC: Fundamentals of Medical Imaging Oct. 20, 2003 NMR-0 Background Basic physics 4 magnetic fields
More informationPhysics 221A Fall 1996 Notes 13 Spins in Magnetic Fields
Physics 221A Fall 1996 Notes 13 Spins in Magnetic Fields A nice illustration of rotation operator methods which is also important physically is the problem of spins in magnetic fields. The earliest experiments
More informationc 2004 Society for Industrial and Applied Mathematics
SIAM J. OPTIM. Vol. 15, No. 1, pp. 39 62 c 2004 Society for Industrial and Applied Mathematics SEMISMOOTH NEWTON AND AUGMENTED LAGRANGIAN METHODS FOR A SIMPLIFIED FRICTION PROBLEM GEORG STADLER Abstract.
More informationBiomedical Imaging Magnetic Resonance Imaging
Biomedical Imaging Magnetic Resonance Imaging Charles A. DiMarzio & Eric Kercher EECE 4649 Northeastern University May 2018 Background and History Measurement of Nuclear Spins Widely used in physics/chemistry
More informationK. Krumbiegel I. Neitzel A. Rösch
SUFFICIENT OPTIMALITY CONDITIONS FOR THE MOREAU-YOSIDA-TYPE REGULARIZATION CONCEPT APPLIED TO SEMILINEAR ELLIPTIC OPTIMAL CONTROL PROBLEMS WITH POINTWISE STATE CONSTRAINTS K. Krumbiegel I. Neitzel A. Rösch
More informationControl of Spin Systems
Control of Spin Systems The Nuclear Spin Sensor Many Atomic Nuclei have intrinsic angular momentum called spin. The spin gives the nucleus a magnetic moment (like a small bar magnet). Magnetic moments
More informationChemistry 431. Lecture 23
Chemistry 431 Lecture 23 Introduction The Larmor Frequency The Bloch Equations Measuring T 1 : Inversion Recovery Measuring T 2 : the Spin Echo NC State University NMR spectroscopy The Nuclear Magnetic
More informationSpectral Broadening Mechanisms. Broadening mechanisms. Lineshape functions. Spectral lifetime broadening
Spectral Broadening echanisms Lorentzian broadening (Homogeneous) Gaussian broadening (Inhomogeneous, Inertial) Doppler broadening (special case for gas phase) The Fourier Transform NC State University
More informationDouble Smoothing technique for Convex Optimization Problems with Linear Constraints
1 Double Smoothing technique for Convex Optimization Problems with Linear Constraints O. Devolder (F.R.S.-FNRS Research Fellow), F. Glineur and Y. Nesterov Center for Operations Research and Econometrics
More informationSketch of the MRI Device
Outline for Today 1. 2. 3. Introduction to MRI Quantum NMR and MRI in 0D Magnetization, m(x,t), in a Voxel Proton T1 Spin Relaxation in a Voxel Proton Density MRI in 1D MRI Case Study, and Caveat Sketch
More informationShimming of a Magnet for Calibration of NMR Probes UW PHYSICS REU 2013
Shimming of a Magnet for Calibration of NMR Probes RACHEL BIELAJEW UW PHYSICS REU 2013 Outline Background The muon anomaly The g-2 Experiment NMR Design Helmholtz coils producing a gradient Results Future
More informationExact penalty decomposition method for zero-norm minimization based on MPEC formulation 1
Exact penalty decomposition method for zero-norm minimization based on MPEC formulation Shujun Bi, Xiaolan Liu and Shaohua Pan November, 2 (First revised July 5, 22) (Second revised March 2, 23) (Final
More informationPhysics of MR Image Acquisition
Physics of MR Image Acquisition HST-583, Fall 2002 Review: -MRI: Overview - MRI: Spatial Encoding MRI Contrast: Basic sequences - Gradient Echo - Spin Echo - Inversion Recovery : Functional Magnetic Resonance
More informationBasic MRI physics and Functional MRI
Basic MRI physics and Functional MRI Gregory R. Lee, Ph.D Assistant Professor, Department of Radiology June 24, 2013 Pediatric Neuroimaging Research Consortium Objectives Neuroimaging Overview MR Physics
More informationDual and primal-dual methods
ELE 538B: Large-Scale Optimization for Data Science Dual and primal-dual methods Yuxin Chen Princeton University, Spring 2018 Outline Dual proximal gradient method Primal-dual proximal gradient method
More informationChristian Clason 1 and Karl Kunisch 1
ESAIM: COCV 7 (0) 43 66 DOI: 0.05/cocv/00003 ESAIM: Control, Optimisation and Calculus of Variations www.esaim-cocv.org A DUALITY-BASED APPROACH TO ELLIPTIC CONTROL PROBLEMS IN NON-REFLEXIVE BANACH SPACES
More informationPhysical Background Of Nuclear Magnetic Resonance Spectroscopy
Physical Background Of Nuclear Magnetic Resonance Spectroscopy Michael McClellan Spring 2009 Department of Physics and Physical Oceanography University of North Carolina Wilmington What is Spectroscopy?
More informationDual methods for the minimization of the total variation
1 / 30 Dual methods for the minimization of the total variation Rémy Abergel supervisor Lionel Moisan MAP5 - CNRS UMR 8145 Different Learning Seminar, LTCI Thursday 21st April 2016 2 / 30 Plan 1 Introduction
More informationLecture #6 NMR in Hilbert Space
Lecture #6 NMR in Hilbert Space Topics Review of spin operators Single spin in a magnetic field: longitudinal and transverse magnetiation Ensemble of spins in a magnetic field RF excitation Handouts and
More informationNMR, the vector model and the relaxation
NMR, the vector model and the relaxation Reading/Books: One and two dimensional NMR spectroscopy, VCH, Friebolin Spin Dynamics, Basics of NMR, Wiley, Levitt Molecular Quantum Mechanics, Oxford Univ. Press,
More informationMagnetic Resonance Imaging. Pål Erik Goa Associate Professor in Medical Imaging Dept. of Physics
Magnetic Resonance Imaging Pål Erik Goa Associate Professor in Medical Imaging Dept. of Physics pal.e.goa@ntnu.no 1 Why MRI? X-ray/CT: Great for bone structures and high spatial resolution Not so great
More informationNMR Dynamics and Relaxation
NMR Dynamics and Relaxation Günter Hempel MLU Halle, Institut für Physik, FG Festkörper-NMR 1 Introduction: Relaxation Two basic magnetic relaxation processes: Longitudinal relaxation: T 1 Relaxation Return
More informationNUMERICAL SOLUTION OF NONSMOOTH PROBLEMS
NUMERICAL SOLUTION OF NONSMOOTH PROBLEMS Abstract. Nonsmooth minimization problems arise in problems involving constraints or models describing certain discontinuities. The existence of solutions is established
More informationLecture 12 February 11, 2016
MATH 262/CME 372: Applied Fourier Analysis and Winter 2016 Elements of Modern Signal Processing Lecture 12 February 11, 2016 Prof. Emmanuel Candes Scribe: Carlos A. Sing-Long, Edited by E. Bates 1 Outline
More informationProblem 1: Spin 1 2. particles (10 points)
Problem 1: Spin 1 particles 1 points 1 Consider a system made up of spin 1/ particles. If one measures the spin of the particles, one can only measure spin up or spin down. The general spin state of a
More informationLight and Matter. Thursday, 8/31/2006 Physics 158 Peter Beyersdorf. Document info
Light and Matter Thursday, 8/31/2006 Physics 158 Peter Beyersdorf Document info 3. 1 1 Class Outline Common materials used in optics Index of refraction absorption Classical model of light absorption Light
More informationIntroduction to MRI. Spin & Magnetic Moments. Relaxation (T1, T2) Spin Echoes. 2DFT Imaging. K-space & Spatial Resolution.
Introduction to MRI Spin & Magnetic Moments Relaxation (T1, T2) Spin Echoes 2DFT Imaging Selective excitation, phase & frequency encoding K-space & Spatial Resolution Contrast (T1, T2) Acknowledgement:
More informationClassical Description of NMR Parameters: The Bloch Equations
Classical Description of NMR Parameters: The Bloch Equations Pascale Legault Département de Biochimie Université de Montréal 1 Outline 1) Classical Behavior of Magnetic Nuclei: The Bloch Equation 2) Precession
More informationSpin Relaxation and NOEs BCMB/CHEM 8190
Spin Relaxation and NOEs BCMB/CHEM 8190 T 1, T 2 (reminder), NOE T 1 is the time constant for longitudinal relaxation - the process of re-establishing the Boltzmann distribution of the energy level populations
More informationNONSMOOTH ANALYSIS AND PARAMETRIC OPTIMIZATION. R. T. Rockafellar*
NONSMOOTH ANALYSIS AND PARAMETRIC OPTIMIZATION R. T. Rockafellar* Abstract. In an optimization problem that depends on parameters, an important issue is the effect that perturbations of the parameters
More informationJoint work with Nguyen Hoang (Univ. Concepción, Chile) Padova, Italy, May 2018
EXTENDED EULER-LAGRANGE AND HAMILTONIAN CONDITIONS IN OPTIMAL CONTROL OF SWEEPING PROCESSES WITH CONTROLLED MOVING SETS BORIS MORDUKHOVICH Wayne State University Talk given at the conference Optimization,
More informationMRI Homework. i. (0.5 pt each) Consider the following arrangements of bar magnets in a strong magnetic field.
MRI Homework 1. While x-rays are used to image bones, magnetic resonance imaging (MRI) is used to examine tissues within the body by detecting where hydrogen atoms (H atoms) are and their environment (e.g.
More informationOn the acceleration of the double smoothing technique for unconstrained convex optimization problems
On the acceleration of the double smoothing technique for unconstrained convex optimization problems Radu Ioan Boţ Christopher Hendrich October 10, 01 Abstract. In this article we investigate the possibilities
More informationIndex. p, lip, 78 8 function, 107 v, 7-8 w, 7-8 i,7-8 sine, 43 Bo,94-96
p, lip, 78 8 function, 107 v, 7-8 w, 7-8 i,7-8 sine, 43 Bo,94-96 B 1,94-96 M,94-96 B oro!' 94-96 BIro!' 94-96 I/r, 79 2D linear system, 56 2D FFT, 119 2D Fourier transform, 1, 12, 18,91 2D sinc, 107, 112
More informationOn duality theory of conic linear problems
On duality theory of conic linear problems Alexander Shapiro School of Industrial and Systems Engineering, Georgia Institute of Technology, Atlanta, Georgia 3332-25, USA e-mail: ashapiro@isye.gatech.edu
More informationA Lecture on Selective RF-pulses in MRI
A Lecture on Selective RF-pulses in MRI Charles L. Epstein February, 23 Introduction We describe the notion of selective excitation and explain how the Bloch equation is used to design selective RF-pulses.
More informationIterative regularization of nonlinear ill-posed problems in Banach space
Iterative regularization of nonlinear ill-posed problems in Banach space Barbara Kaltenbacher, University of Klagenfurt joint work with Bernd Hofmann, Technical University of Chemnitz, Frank Schöpfer and
More informationChapter 8 Magnetic Resonance
Chapter 8 Magnetic Resonance 9.1 Electron paramagnetic resonance 9.2 Ferromagnetic resonance 9.3 Nuclear magnetic resonance 9.4 Other resonance methods TCD March 2007 1 A resonance experiment involves
More informationSEVERAL PATH-FOLLOWING METHODS FOR A CLASS OF GRADIENT CONSTRAINED VARIATIONAL INEQUALITIES
SEVERAL PATH-FOLLOWING METHODS FOR A CLASS OF GRADIENT CONSTRAINED VARIATIONAL INEQUALITIES M. HINTERMÜLLER AND J. RASCH Abstract. Path-following splitting and semismooth Newton methods for solving a class
More informationFINITE-DIFFERENCE APPROXIMATIONS AND OPTIMAL CONTROL OF THE SWEEPING PROCESS. BORIS MORDUKHOVICH Wayne State University, USA
FINITE-DIFFERENCE APPROXIMATIONS AND OPTIMAL CONTROL OF THE SWEEPING PROCESS BORIS MORDUKHOVICH Wayne State University, USA International Workshop Optimization without Borders Tribute to Yurii Nesterov
More informationRelaxation times in nuclear magnetic resonance
Relaxation times in TEP Related topics Nuclear spins, atomic nuclei with a magnetic moment, precession movement of the nuclear spins, Landau-Lifshitz equation, Bloch equation, magnetisation, resonance
More informationSmoothing Proximal Gradient Method. General Structured Sparse Regression
for General Structured Sparse Regression Xi Chen, Qihang Lin, Seyoung Kim, Jaime G. Carbonell, Eric P. Xing (Annals of Applied Statistics, 2012) Gatsby Unit, Tea Talk October 25, 2013 Outline Motivation:
More informationTopics. The concept of spin Precession of magnetic spin Relaxation Bloch Equation. Bioengineering 280A Principles of Biomedical Imaging
Bioengineering 280A Principles of Biomedical Imaging Fall Quarter 2006 MRI Lecture 1 Topics The concept of spin Precession of magnetic spin Relaxation Bloch Equation 1 Spin Intrinsic angular momentum of
More informationSEMI-SMOOTH SECOND-ORDER TYPE METHODS FOR COMPOSITE CONVEX PROGRAMS
SEMI-SMOOTH SECOND-ORDER TYPE METHODS FOR COMPOSITE CONVEX PROGRAMS XIANTAO XIAO, YONGFENG LI, ZAIWEN WEN, AND LIWEI ZHANG Abstract. The goal of this paper is to study approaches to bridge the gap between
More informationSpin. Nuclear Spin Rules
Spin Bioengineering 280A Principles of Biomedical Imaging Fall Quarter 2012 MRI Lecture 1 Intrinsic angular momentum of elementary particles -- electrons, protons, neutrons. Spin is quantized. Key concept
More informationPart II: Magnetic Resonance Imaging (MRI)
Part II: Magnetic Resonance Imaging (MRI) Contents Magnetic Field Gradients Selective Excitation Spatially Resolved Reception k-space Gradient Echo Sequence Spin Echo Sequence Magnetic Resonance Imaging
More informationMeasuring Spin-Lattice Relaxation Time
WJP, PHY381 (2009) Wabash Journal of Physics v4.0, p.1 Measuring Spin-Lattice Relaxation Time L.W. Lupinski, R. Paudel, and M.J. Madsen Department of Physics, Wabash College, Crawfordsville, IN 47933 (Dated:
More informationLecture: Smoothing.
Lecture: Smoothing http://bicmr.pku.edu.cn/~wenzw/opt-2018-fall.html Acknowledgement: this slides is based on Prof. Lieven Vandenberghe s lecture notes Smoothing 2/26 introduction smoothing via conjugate
More informationNewton-type Methods for Solving the Nonsmooth Equations with Finitely Many Maximum Functions
260 Journal of Advances in Applied Mathematics, Vol. 1, No. 4, October 2016 https://dx.doi.org/10.22606/jaam.2016.14006 Newton-type Methods for Solving the Nonsmooth Equations with Finitely Many Maximum
More informationLegendre Transforms, Calculus of Varations, and Mechanics Principles
page 437 Appendix C Legendre Transforms, Calculus of Varations, and Mechanics Principles C.1 Legendre Transforms Legendre transforms map functions in a vector space to functions in the dual space. From
More informationElectron spins in nonmagnetic semiconductors
Electron spins in nonmagnetic semiconductors Yuichiro K. Kato Institute of Engineering Innovation, The University of Tokyo Physics of non-interacting spins Optical spin injection and detection Spin manipulation
More informationAUGMENTED LAGRANGIAN METHODS FOR CONVEX OPTIMIZATION
New Horizon of Mathematical Sciences RIKEN ithes - Tohoku AIMR - Tokyo IIS Joint Symposium April 28 2016 AUGMENTED LAGRANGIAN METHODS FOR CONVEX OPTIMIZATION T. Takeuchi Aihara Lab. Laboratories for Mathematics,
More informationMagnetic Resonance Imaging
http://www.qldxray.com.au/filelibrary/mri_cardiovascular_system_ca_0005.jpg Magnetic Resonance Imaging 1 Overview 1. The magnetic properties of nuclei, and how they behave in strong magnetic fields. 2.
More informationINTRODUCTION TO NMR and NMR QIP
Books (NMR): Spin dynamics: basics of nuclear magnetic resonance, M. H. Levitt, Wiley, 2001. The principles of nuclear magnetism, A. Abragam, Oxford, 1961. Principles of magnetic resonance, C. P. Slichter,
More informationNuclear Magnetic Resonance Spectroscopy
Nuclear Magnetic Resonance Spectroscopy Ecole Polytechnique Département de Chimie CHI 551 Dr. Grégory Nocton Bureau 01 30 11 A Tel: 44 02 Ecole polytechnique / CNRS Laboratoire de Chimie Moléculaire E-mail:
More informationComputation of Local ISS Lyapunov Function Via Linear Programming
Computation of Local ISS Lyapunov Function Via Linear Programming Huijuan Li Joint work with Robert Baier, Lars Grüne, Sigurdur F. Hafstein and Fabian Wirth Institut für Mathematik, Universität Würzburg
More informationSpin Dynamics Basics of Nuclear Magnetic Resonance. Malcolm H. Levitt
Spin Dynamics Basics of Nuclear Magnetic Resonance Second edition Malcolm H. Levitt The University of Southampton, UK John Wiley &. Sons, Ltd Preface xxi Preface to the First Edition xxiii Introduction
More informationSparse Regularization via Convex Analysis
Sparse Regularization via Convex Analysis Ivan Selesnick Electrical and Computer Engineering Tandon School of Engineering New York University Brooklyn, New York, USA 29 / 66 Convex or non-convex: Which
More informationMaster 2 MathBigData. 3 novembre CMAP - Ecole Polytechnique
Master 2 MathBigData S. Gaïffas 1 3 novembre 2014 1 CMAP - Ecole Polytechnique 1 Supervised learning recap Introduction Loss functions, linearity 2 Penalization Introduction Ridge Sparsity Lasso 3 Some
More informationNuclear Magnetic Resonance Spectroscopy
Nuclear Magnetic Resonance Spectroscopy Features: Used to identify products of reactions Also gives information about chemical environment, connectivity and bonding of nuclei Requirements: Pure or mostly
More informationProximal methods. S. Villa. October 7, 2014
Proximal methods S. Villa October 7, 2014 1 Review of the basics Often machine learning problems require the solution of minimization problems. For instance, the ERM algorithm requires to solve a problem
More informationUtrecht University. Radiofrequency pulse design through optimal control and model order reduction of the Bloch equation
Utrecht University Master s thesis Radiofrequency pulse design through optimal control and model order reduction of the Bloch equation Author: Willem van Valenberg Supervisors: Gerard L.G. Sleijpen Alessandro
More informationMaximal monotone operators are selfdual vector fields and vice-versa
Maximal monotone operators are selfdual vector fields and vice-versa Nassif Ghoussoub Department of Mathematics, University of British Columbia, Vancouver BC Canada V6T 1Z2 nassif@math.ubc.ca February
More informationSelf-dual Smooth Approximations of Convex Functions via the Proximal Average
Chapter Self-dual Smooth Approximations of Convex Functions via the Proximal Average Heinz H. Bauschke, Sarah M. Moffat, and Xianfu Wang Abstract The proximal average of two convex functions has proven
More informationOn the Convergence and O(1/N) Complexity of a Class of Nonlinear Proximal Point Algorithms for Monotonic Variational Inequalities
STATISTICS,OPTIMIZATION AND INFORMATION COMPUTING Stat., Optim. Inf. Comput., Vol. 2, June 204, pp 05 3. Published online in International Academic Press (www.iapress.org) On the Convergence and O(/N)
More informationOn nonexpansive and accretive operators in Banach spaces
Available online at www.isr-publications.com/jnsa J. Nonlinear Sci. Appl., 10 (2017), 3437 3446 Research Article Journal Homepage: www.tjnsa.com - www.isr-publications.com/jnsa On nonexpansive and accretive
More information1. neopentyl benzene. 4 of 6
I. 1 H NMR spectroscopy A. Theory 1. The protons and neutrons in atomic nuclei spin, as does the nucleus itself 2. The circulation of nuclear charge can generate a nuclear magnetic moment, u, along the
More informationSpectral Decompositions using One-Homogeneous Functionals
Spectral Decompositions using One-Homogeneous Functionals Martin Burger, Guy Gilboa, Michael Moeller, Lina Eckardt, and Daniel Cremers Abstract. This paper discusses the use of absolutely one-homogeneous
More informationThis can be 2 lectures! still need: Examples: non-convex problems applications for matrix factorization
This can be 2 lectures! still need: Examples: non-convex problems applications for matrix factorization x = prox_f(x)+prox_{f^*}(x) use to get prox of norms! PROXIMAL METHODS WHY PROXIMAL METHODS Smooth
More informationMagnetic Resonance Imaging (MRI)
Magnetic Resonance Imaging Introduction The Components The Technology (MRI) Physics behind MR Most slides taken from http:// www.slideworld.org/ viewslides.aspx/magnetic- Resonance-Imaging- %28MRI%29-MR-Imaging-
More informationNMR Spectroscopy: A Quantum Phenomena
NMR Spectroscopy: A Quantum Phenomena Pascale Legault Département de Biochimie Université de Montréal Outline 1) Energy Diagrams and Vector Diagrams 2) Simple 1D Spectra 3) Beyond Simple 1D Spectra 4)
More informationIntroduction to Magnetic Resonance Imaging (MRI) Pietro Gori
Introduction to Magnetic Resonance Imaging (MRI) Pietro Gori Enseignant-chercheur Equipe IMAGES - Télécom ParisTech pietro.gori@telecom-paristech.fr September 20, 2017 P. Gori BIOMED 20/09/2017 1 / 76
More information