Physical Design Issues in Biofluidic Microchips Tamal Mukherjee MEMS Laboratory ECE Department Carnegie Mellon University Pittsburgh, PA, USA tamal@ece.cmu.edu http://www.ece.cmu.edu/ cmu edu/~mems Carnegie Mellon ISPD, April 16, 2008
Tubes to Chips: ICs Driven by Information Processing needs IBM 701 calculator (1952) Intel 4004 Calculator IC (1971) 2
Tubes to Chips: BioChips Driven by Biomolecular Analysis needs Image from Barnard College Archives Test tubes & Beakers Agilent DNA analysis (1950) Lab on a Chip (1997) 3
Portable Analysis New knowledge of molecular basis of biology e.g. Human Genome Project Massively parallel analysis infrastructure Integration and miniaturization will drive biomolecular analysis instrumentation Biomolecular Spock with Tricorder mainframes Sensor + computer Burns Science 2002 4
Typical Biological Lab Functions Synthesis Analysis A B C A+B A+B A B Mixing Reaction Separation 5
Microdevice Technology Summary BioMEMS BioChips Lab-on-a Chip Droplet Channel Pressure Electrokinetics 6
Channel-based LoC: EK drive What is Electrokinetics? Voltage driven flow Flow direction Why Electrokinetic flow? Plug velocity profile Portable kv sources EK flow can be used for electrophoresis EK flow already used in complex designs EK flow Pressure flow Serial Mixer ORNL 7
Microdevice Technologies: LoC Miniaturized Bio-chemical Lab-on-a-Chip Individual functional units demonstrated Analyzer, Reactor, Research driven by integration Design aids needed to handle complexity! Amino-Acid Analysis DNA Analysis Chemical Synthesis Immunoassay ORNL U. C. Berkeley U. Hull U. Alberta 8
Outline Introduction Motivation for Design Automation Design Hierarchy Multi-function System Simulation Multi-plex Physical Synthesis Summary 9
Multiplex Lab-on-a-Chip Same subsystem, integrated for redundancy, combinatorial i experiments integration year 10
Multifunction LoCs Example: Immunoassay 1. Load sample {Ag*, Ag} 2. Mix with reagents {Ab} 3. Rxn: Ag* + Ab Ag*-Ab Ab mixing 4. Inject sample plug 5. Separate analytes 6. Detection {[Ag*,Ag], Ab, Ag*-Ab} reaction injection loading reagent buffer separation sample waste detection buffer waste 11
Complexity Hierarchy Element Functional Subsystem System Component Increasing Integration 12
Outline Introduction Motivation for Design Automation Design Hierarchy Multi-function System Simulation Multi-plex Physical Synthesis Summary 13
Simulation Techniques Computational fluid dynamics One single turn Complimentary turns buffer Buffer Sampl e sample Flow direction Flow direction ~ 10 Hours 2~3 days ~ 10 hours Reduced order models Hierarchical decomposition and parameterization Capture geometric effects Amenable for use in design Serial Mixer (ORNL) 14
Hierarchy Example: Immunoassay Sample Buffer V+ Ag*, Ag V+ Flow Direction V- Mixing and Reaction Ag* + Ab Ag*-Ab Sample Waste Pinching Ab Ag*, Ag, Ab, V+ Ag*-Ab Ab V+ Buffer Waste 15
Synthesis Phase: Steady State Buffer V+ V+ Flow Direction V- V+ Mixing and Reaction Ag* + Ab Ag*-Ab Pinching Sample Waste Ag*, Ag, Ab, Ag*-Ab Ab V+ Buffer Waste 16
Analysis Phase: Transient V- Buffer V+ Flow Direction Sample Plug Mixture of Ag*, Ag, Ab, Ag*-Ab lv- V- Sample Waste Separatio on Channe V- Buffer Waste 17
Analysis Phase: Transient V- Buffer V+ Flow Direction V- Sample Waste Ag*-Ab Ab Ag*,Ag lv- Separatio on Channe Resolution = distance apart band broadness V- Buffer Waste 18
Component Library Library of LoC Unit Operations Compose Topology Function Type well mixer reactor injector separator splitter 19
Composition Examples Buffer Sample Sample waste A B Serpentine separation chip (ORNL) Multi-stream t mixer (M. Koch, et al.) System waste Spiral chip (ORNL) System waste Sample wa aste Buffer Sample Buffer A 1 A 2 A 3 A 4 A 5 Sample Waste-2 Was ste-1 Serial Mixing network (S.C. Jacobson, et 20al.)
Outline Introduction Motivation for Design Automation Design Hierarchy Component Models Multi-function System Simulation Multi-plex Physical Synthesis Summary 21
Simulating a Multifunction Design (Cheim, Clin. Chem., 44:3, 591-598, 1998) Real Immunoassay Chip from U. Alberta Operation Mixing/Dilution Reaction Injection Separation Detection Wang et. al. Transducers 05 22
Simulation Results 10 mm after injection Calibration curve Re elative concen ntration c Ab-Ag* Th * -Ab complex Area ratio Before turn 0.6 Th * Ag* After turn 1.0 Schematic Experimental 0.8 0.4 0.2 0.0 E i l Ag* Th * Ab-Ag* Ab-Th * 0 10 20 30 40 Antigen Theophylline (Ag) Th (mg/l) 0 5 10 15 20 25 30 35 Time (s) Electropherogram Ag unreacted Ag* Ab-Ag* Simulation matches experiment Simulation time is a few CPU seconds 23 Wang et. al. Transducers 05
Optimizing the design: NLP x 0 obj : min f s. t. g( y i ( y ) < 0 h ( y ) = 0 y i i ) x i y i y i = SIM ( x, PARAMS i ) 7.6 cm x * 10x less space Same perf injector reactor 1.22 cm mixer separation channel detector 0.75 cm 2 cm 1.73 cm 1.14 cm 7.6 cm 2.23 cm 2.33 cm Wang et. al. Transducers 05 2.5 cm 24
Outline Introduction Motivation for Design Automation Design Hierarchy Multi-function System Simulation Multi-plex Physical Synthesis Summary 25
Multiplex Physical Synthesis Input: Design Specs Overall Chip Dimensions Species/Buffer properties Operational constraints Chip fabrication Subsystem performance Family of subsystems Final Routed Layout Simultaneously determine: placement dimensions # of sections voltage Intermediate Placement Pfeiffer et. al., TCAD 06 Route subsystems to wells: single layer, planar min. length, bends 26
Placement Features Subsystem optimization: NLP System-on-Chip extensions * : * Murata, H. et al., IEEE Trans. on CAD. 1996 s e c f a d b t e c f a b d Orientation: and Well placement: e a d c Overlap constraints: Pfeiffer et. al., TCAD 06 f b E T e c a d a EL ER c f b E B Never a Penalty a d c b e f e f d b 27 b
Routing Features Routing grid graph * : e a d Expand e a d * Lengauer, T., Combinatorial Algs. for IC Layout, 1990 c f b c f b 7 8 9 Node constraints: 4 5 6 1 2 3 flows in = flows out 1 flow in/out of node (single layer, planar) +1 7 8 9 penalize bends +1 4 6 Bend reduction: 5 favor straight paths +0 +1 Pfeiffer et. al., TCAD 06 1 2 3 28
Multiplex Synthesis Example Placed and Routed Design P&R By Hand Automated Improvement Place: 20 min Time 5+ hrs. Route: 3 min > 10X faster Total: 23 min Dimensionsi 167 1.67cm x 88 8.8cm 161 1.61cm x 379 3.79cm ~ 25X 2.5X smaller Pfeiffer et. al., TCAD 06 29
Summary Lab on a Chips integrate bioanalysis functions Hierarchically decomposition used for development of fast, accurate, reusable, parameterized models Relatively few types of band profile on a chip Profile representation to simplify PDE into ODEs Models are reuseable Separation models integrated with P&R algorithms for simultaneous model based placement followed by routing Focus still at circuit -level, need to consider architecture, protocol optimization 30
Acknowledgements Collaborators Prof. James F. Hoburg (ECE) Prof. Steinar Hauan (ChE) Prof. Qiao Lin (ME, now at Columbia Univ) Students Anton Pfeiffer, Yi Wang, Ryan Magargle, Xiang He, Bikram Baidya Funding DARPA DSO SIMBIOSYS Program NSF ITR Program 31