Murray Z. Frank, Vidhan K. Goyal

合集下载

Michael Gratzel简介

Michael Gratzel简介

Michael GrätzelFrom Wikipedia, the free encyclopediaMichael Grätzel (born 11 May 1944, in Dorfchemnitz, Saxony, Germany)[1] is a professor at the École Polytechnique Fédérale de Lausanne where he directs the Laboratory of Photonics and Interfaces. He pioneered research on energy and electron transfer reactions in mesoscopic-materials and their optoelectronic applications. He discovered a new type of solar cell based on dye sensitized mesoscopic oxide particles and pioneered the use of nanomaterials in lithium ion batteries.[2][3]Author of over 900 publications, two books and inventor or co-inventor of over 50 patents, he has been the Mary Upton Visiting Professor at Cornell University and a Distinguished VisitingProfessor at the National University of Singapore. He was an Invited Professor at the University of California, Berkeley, the École normale supérieure de Cachan(Paris) and Delft University of Technology. His work has been cited over 107’000 times (h-index 151) making him one of the 10 most highly cited chemists in the world.[4]He was a frequent guest scientist at the National Renewable Energy Laboratory (NREL) in Golden, Colorado, was a fellow of the Japanese Society for the Promotion of Science. In 2009 he was named Distinguished Honorary Professor by the Chinese Academy of Science (Changchun) and the Huazhong University of Science and Technology.He has received numerous awards including the Millennium 2000 European innovation prize, the 2001 Faraday Medal of the British Royal Society, the 2001 Dutch Havinga Award, the 2004 Italgas Prize, two McKinsey Venture awards in 1998 and 2002 and the 2005 Gerischer Prize. In 2007 he was awarded the Harvey Prize of Technion for pioneered research on energy and electron transfer reactions in mesoscopic-materials and their optoelectronic applications. In 2009 he was awarded the Balzan Prize for the Science of New Materials. His most recent awards include the 2012 Albert Einstein World Award of Science,[5]2011 Gutenberg Research Award, 2011 Paul Karrer Gold Medal and the 2010 Millennium Technology Grand Prize. He holds a doctorate from the Technical University of Berlin and honorary doctorates from the Universities of Uppsala, Turin and Nova Gorica. He was elected honorary member of the Société Vaudoise des Sciences Naturelles. Dr. Grätzel is a member of the Scientific Advisory Committee at the IMDEA Nanoscience Institute.。

BOR2011.12

BOR2011.12

WORLD OF REPRODUCTIVE BIOLOGY ..................................................1089Charlotte SchubertCommentary.Dissecting the Phthalate-Induced Sertoli Cell Injury:The Fragile Balance of Proteases and Their Inhibitors .........................................................................1091Se ´verine Mazaud-GuittotIn the current issue of Biology of Reproduction ,Yao et al.explore the transcriptional mechanism of Timp2down-regulation that follows phthalate exposure.Minireview.Uterine Histotroph and Conceptus Development:Select Nutrients and Secreted Phosphoprotein 1Affect Mechanistic Target of Rapamycin Cell Signaling in Ewes ..................................1094Fuller W.Bazer,Guoyao Wu,Greg A.Johnson,Jinyoung Kim,and Gwonhwa SongSelect nutrients and secreted phosphoprotein 1in uterine histotroph activate the mechanistic target of rapamycin cell signaling pathway required for conceptus development.Seasonal Androgen Cycles in Adult Male American Alligators (Alligator mississippiensis )from a Barrier Island Population ........................................................................1108Heather J.Hamlin,Russell H.Lowers,and Louis J.Guillette,Jr.Plasma testosterone and dehydroepiandrosterone levels in adult male alligators reflect reproductive activity,but patterns differ between large and small adults.Spermatogonial Stem Cell Self-Renewal Requires ETV5-Mediated Downstream Activation of Brachyury inMice ............................................................................1114Xin Wu,Shaun M.Goodyear,John W.Tobias,Mary R.Avarbock,and Ralph L.BrinsterMicroarray profiling of ETV5-associated gene expression identified Brachyury as a novel target gene expressed in spermatogonial stem cells that is critical for self-renewal.CDK7and CCNH Are Components of CDK-Activating Kinase and Are Required for Meiotic Progression of Pig Oocytes ..........................................................................1124Wataru Fujii,Takanori Nishimura,Kiyoshi Kano,Koji Sugiura,and Kunihiko NaitoCDK7and CCNH are components of CDK-activating kinase (CAK),which activates CDC2by T161phosphorylation and are required for meiotic progression of porcine oocytes.Oxytocin Increases Invasive Properties of Endometrial Cancer Cells Through Phosphatidylinositol 3-Kinase/AKT-Dependent Up-Regulation of Cyclooxygenase-1,-2,and X-Linked Inhibitor of Apoptosis Protein ...........1133Marie-Claude De ´ry,Parvesh Chaudhry,Vale ´rie Leblanc,Sophie Parent,Anne-Marie Fortier,and Eric Asselin Oxytocin regulates endometrial cancer cell invasion through phosphatidylinositol 3-kinase/AKT-dependent up-regulation of cyclooxygenase-1,-2,and X-linked inhibitor of apoptosis protein.RNA Sequencing Reveals Novel Gene Clusters in Bovine Conceptuses Associated with Maternal Recognition of Pregnancy and Implantation ...........................................................1143Solomon Mamo,Jai P.Mehta,Paul McGettigan,Trudee Fair,Thomas E.Spencer,Fuller W.Bazer,and Patrick LonerganBovine conceptus transcript discovery reveals dynamics consisting of known and novel genes unique to key stages of pre-and peri-implantation embryo development.Effects of Gonadotropin-Releasing Hormone Immunization on Reproductive Function and Behavior in Captive Female Rocky Mountain Elk (Cervus elaphus nelsoni )..........................................1152Jenny G.Powers,Dan L.Baker,Tracy L.Davis,Mary M.Conner,Anneke H.Lothridge,and Terry t A single GnRH immunization with GonaCon-B decreased fertility for three subsequent reproductive seasons in female elk.Stem Leydig Cell Differentiation:Gene Expression During Development of the Adult Rat Population of Leydig Cells ............................................................................1161Erin L.Stanley,Daniel S.Johnston,Jinjiang Fan,Vassilios Papadopoulos,Haolin Chen,Ren-Shan Ge,Barry R.Zirkin,and Scott A.JelinskyUsing microarray gene expression analysis,changes during the transition of stem to adult Leydig cells in the rat can be identified,and differences between stem Leydig cells and bone marrow stem cells can be demonstrated.BIOLOGY of REPRODUCTIONOfficial Journal of the Society forthe Study of ReproductionDecember 2011VOLUME 85NUMBER 6iiiRole of Angiotensin II in the Periovulatory Epidermal Growth Factor-Like Cascade in Bovine Granulosa Cells In Vitro (1167)Vale´rio M.Portela,Gustavo Zamberlam,Paulo B.D.Gonc¸alves,Joa˜o F.C.de Oliveira,and Christopher A.PriceAngiotensin II is a cofactor with LH in the early induction of ADAM17in periovulatory granulosa cells.Loss of Gremlin Delays Primordial Follicle Assembly but Does Not Affect Female Fertility in Mice (1175)Michelle Myers,Swamy K.Tripurani,Brooke Middlebrook,Aris N.Economides,Ernesto Canalis,and Stephanie A.PangasThe Grem1knockout has novel defects in early postnatal follicle development,but the conditional knockout in theovary has normal fertility,due to compensation by other BMP antagonists expressed in the ovary.Intracytoplasmic Sperm Injection with Mouse Spermatozoa Preserved Without Freezing for Six Months CanLead to Full-Term Development (1183)Chong Li,Eiji Mizutani,Tetsuo Ono,Yukari Terashita,Xiao-feng Jia,Hui-juan Shi,and Teruhiko WakayamaMouse spermatozoa preserved up to six months in a BSA-or Ficoll-added media without freezing can support full-term development via ICSI.Importin Alpha2-Interacting Proteins with Nuclear Roles During Mammalian Spermatogenesis (1191)Jennifer D.Ly-Huynh,Kim G.Lieu,Andrew T.Major,Penelope A.F.Whiley,Janet E.Holt,Kate L.Loveland,andDavid A.JansThe role of importin alpha2in mediating nuclear import of Hop2,Arip3,and a novel Chrp variant,is critical forhomologous chromosome pairing and the progression of spermatogenesis.Transcriptional Suppression of Sertoli Cell Timp2in Rodents Following Mono-(2-ethylhexyl)Phthalate ExposureIs Regulated by CEBPA and MYC (1203)Pei-Li Yao,Yi-Chen Lin,and John H.RichburgMono-(2-ethylhexyl)phthalate decreases Timp2expression in Sertoli cells via deactivation of CEBPA andactivation of MYC.New Insights on the Morphology of Adult Mouse Penis (1216)Esequiel Rodriguez,Jr.,Dana A.Weiss,Jennifer H.Yang,Julia Menshenina,Max Ferretti,Tristan J.Cunha,DaleBarcellos,Lok Yun Chan,Gail Risbridger,Gerald R.Cunha,and Laurence S.BaskinComplete formation of the adult mouse penis is the end point of masculine sex differentiation of the externalgenitalia.Mediators of the JAK/STAT Signaling Pathway in Human Spermatozoa (1222)Catherine Lachance and Pierre LeclercMany mediators of the JAK/STAT signaling pathway are expressed in cytoskeletal structures from the head andflagellum of human spermatozoa,suggesting an involvement in functions other than their well-known transcriptionfactor activity.Adrenomedullin2/Intermedin Regulates HLA-G in Human Trophoblasts (1232)Madhu Chauhan,Meena Balakrishnan,Uma Yallampalli,Janice Endsley,Gary D.V.Hankins,Regan Theiler,andChandra YallampalliAdrenomedullin2enhances trophoblast cell invasion via stimulation of MAPK3/1signaling and regulation of theHLA-G molecule.The Proximal Promoter Region of the Zebrafish gsdf Gene Is Sufficient to Mimic the Spatio-Temporal Expression Pattern of the Endogenous Gene in Sertoli and Granulosa Cells (1240)Aude Gautier,Fre´de´ric Sohm,Jean-Ste´phane Joly,Florence Le Gac,and Jean-Jacques LareyreThe proximal promoter region of the zebrafish gsdf gene is sufficient to drive the spatio-temporal expression patternof the endogenous genes in Sertoli and granulosa cells in vivo in transgenic zebrafish lines.Gap Junction-Mediated Communications Regulate Chromatin Remodeling During Bovine Oocyte Growth and Differentiation Through cAMP-Dependent Mechanism(s) (1252)Alberto M.Luciano,Federica Franciosi,Silvia C.Modina,and Valentina LoddeGap junction-mediated communications regulate chromatin remodeling,transcriptional activity,and developmentalcompetence acquisition during growth and differentiation of bovine oocytes from small antral follicles via a cAMP-mediated mechanism(s).Avian SERPINB11Gene:Characteristics,Tissue-Specific Expression,and Regulation of Expression by Estrogen (1260)Whasun Lim,Ji-Hye Kim,Suzie E.Ahn,Wooyoung Jeong,Jinyoung Kim,Fuller W.Bazer,Jae Yong Han,andGwonhwa SongSERPINB11is a novel estrogen-induced gene expressed only in the epithelial cells of the chicken oviduct thatregulates oviductal development and differentiated functions.iii Aging Results in Differential Regulation of DNA Repair Pathways in Pachytene Spermatocytes in the BrownNorway Rat (1269)Catriona Paul,Makoto Nagano,and Bernard RobaireTesticular aging coincides with a decrease in the efficiency of the base excision repair pathway and increased DNAdamage in germ cells.Age-Dependent Susceptibility of Chromosome Cohesion to Premature Separase Activation in Mouse Oocytes (1279)Teresa Chiang,Richard M.Schultz,and Michael mpsonAge increases susceptibility of chromosome cohesion to premature separase activation,and cohesion is protectedindependently by securin and an inhibitory phosphorylation of separase.Additions and Corrections (1284)Author Index (1288)Volume Contents (1291)Contents by CategoryCommentary1091Dissecting the Phthalate-Induced Sertoli Cell Injury:The Fragile Balance of Proteases and Their Inhibitors.Se´verine Mazaud-Guittot Embryo1143RNA Sequencing Reveals Novel Gene Clusters in Bovine Conceptuses Associated with Maternal Recognition of Pregnancyand Implantation.Solomon Mamo,Jai P.Mehta,Paul McGettigan,Trudee Fair,Thomas E.Spencer,Fuller W.Bazer,and PatrickLonerganFemale Reproductive Tract1133Oxytocin Increases Invasive Properties of Endometrial Cancer Cells Through Phosphatidylinositol3-Kinase/AKT-Dependent Up-Regulation of Cyclooxygenase-1,-2,and X-Linked Inhibitor ofApoptosis Protein.Marie-Claude De´ry,Parvesh Chaudhry,Vale´rieLeblanc,Sophie Parent,Anne-Marie Fortier,and Eric Asselin1232Adrenomedullin2/Intermedin Regulates HLA-G in Human Trophoblasts.Madhu Chauhan,Meena Balakrishnan,UmaYallampalli,Janice Endsley,Gary D.V.Hankins,Regan Theiler,and Chandra Yallampalli1260Avian SERPINB11Gene:Characteristics,Tissue-Specific Expression,and Regulation of Expression by Estrogen.Whasun Lim, Ji-Hye Kim,Suzie E.Ahn,Wooyoung Jeong,Jinyoung Kim,Fuller W.Bazer,Jae Yong Han,and Gwonhwa SongGamete Biology1124CDK7and CCNH Are Components of CDK-Activating Kinase and Are Required for Meiotic Progression of Pig Oocytes.Wataru Fujii,Takanori Nishimura,Kiyoshi Kano,Koji Sugiura,and Kunihiko Naito 1222Mediators of the JAK/STAT Signaling Pathway in Human Spermatozoa.Catherine Lachance and Pierre Leclerc1252Gap Junction-Mediated Communications Regulate Chromatin Remodeling During Bovine Oocyte Growth and DifferentiationThrough cAMP-Dependent Mechanism(s).Alberto M.Luciano,Federica Franciosi,Silvia C.Modina,and Valentina Lodde1279Age-Dependent Susceptibility of Chromosome Cohesion to Premature Separase Activation in Mouse Oocytes.Teresa Chiang,Richard M.Schultz,and Michael mpsonMale Reproductive Tract1108Seasonal Androgen Cycles in Adult Male American Alligators (Alligator mississippiensis)from a Barrier Island Population.Heather J.Hamlin,Russell H.Lowers,and Louis J.Guillette,Jr.1216New Insights on the Morphology of Adult Mouse Penis.Esequiel Rodriguez,Jr.,Dana A.Weiss,Jennifer H.Yang,Julia Menshenina, Max Ferretti,Tristan J.Cunha,Dale Barcellos,Lok Yun Chan,GailRisbridger,Gerald R.Cunha,and Laurence S.BaskinMinireview1094Uterine Histotroph and Conceptus Development:Select Nutrients and Secreted Phosphoprotein1Affect Mechanistic Target of Rapamycin Cell Signaling in Ewes.Fuller W.Bazer,Guoyao Wu,Greg A.Johnson,Jinyoung Kim,and Gwonhwa Song Ovary1167Role of Angiotensin II in the Periovulatory Epidermal Growth Factor-Like Cascade in Bovine Granulosa Cells In Vitro.Vale´rio M.Portela, Gustavo Zamberlam,Paulo B.D.Gonc¸alves,Joa˜o F.C.de Oliveira,and Christopher A.Price1175Loss of Gremlin Delays Primordial Follicle Assembly but Does Not Affect Female Fertility in Mice.Michelle Myers,Swamy K.Tripurani, Brooke Middlebrook,Aris N.Economides,Ernesto Canalis,andStephanie A.PangasReproductive Technology1152Effects of Gonadotropin-Releasing Hormone Immunization on Reproductive Function and Behavior in Captive Female RockyMountain Elk(Cervus elaphus nelsoni).Jenny G.Powers,Dan L.Baker,Tracy L.Davis,Mary M.Conner,Anneke H.Lothridge,andTerry t1183Intracytoplasmic Sperm Injection with Mouse Spermatozoa Preserved Without Freezing for Six Months Can Lead to Full-TermDevelopment.Chong Li,Eiji Mizutani,Tetsuo Ono,Yukari Terashita, Xiao-feng Jia,Hui-juan Shi,and Teruhiko WakayamaTestis1114Spermatogonial Stem Cell Self-Renewal Requires ETV5-Mediated Downstream Activation of Brachyury in Mice.Xin Wu,Shaun M.Goodyear,John W.Tobias,Mary R.Avarbock,and Ralph L.Brinster 1161Stem Leydig Cell Differentiation:Gene Expression During Development of the Adult Rat Population of Leydig Cells.Erin L.Stanley,Daniel S.Johnston,Jinjiang Fan,Vassilios Papadopoulos,Haolin Chen,Ren-Shan Ge,Barry R.Zirkin,and Scott A.Jelinsky 1191Importin Alpha2-Interacting Proteins with Nuclear Roles During Mammalian Spermatogenesis.Jennifer D.Ly-Huynh,Kim G.Lieu,Andrew T.Major,Penelope A.F.Whiley,Janet E.Holt,Kate L.Loveland,and David A.Jans1240The Proximal Promoter Region of the Zebrafish gsdf Gene Is Sufficient to Mimic the Spatio-Temporal Expression Pattern of theEndogenous Gene in Sertoli and Granulosa Cells.Aude Gautier,Fre´de´ric Sohm,Jean-Ste´phane Joly,Florence Le Gac,and Jean-Jacques Lareyre1269Aging Results in Differential Regulation of DNA Repair Pathways in Pachytene Spermatocytes in the Brown Norway Rat.Catriona Paul, Makoto Nagano,and Bernard RobaireToxicology1203Transcriptional Suppression of Sertoli Cell Timp2in Rodents Following Mono-(2-ethylhexyl)Phthalate Exposure Is Regulated byCEBPA and MYC.Pei-Li Yao,Yi-Chen Lin,and John H.Richburgiv。

翻译原文

翻译原文

The primary goal of paper is to evaluate bridge management system in reducing maintenance and repair costs and studying impact of these systems on the structure’s life cycle. The paper a combination of management, engineering and economy has been proposed that: 1) increase safe traffics in roads network. 2) the extension of infrastructure life cycle using high benefit cost ratio maintenance programs 3)improvement of bridge condition by supplying decision makers with a specific tool for analysis of condition status and prediction of remedy needed. The paper addresses issues regarding with maintenance and repair policies and decisions for replaced bridge elements that from studies of previous experiences and interviewing with structural engineers in different companies. Bridge companies throughout the country will benefit from results of this research by comparing different maintenance methods. Findings show that using the proposed bridge management system has high influence on reducing the maintenance costs and increasing life cycle. However, it has also been revealed that, in special conditions, using engineers’ experience and decision-making capabilities has high impact on solving problems of bridg bridge management system for deciding of bridge maintenance and repair is very important, few researches have been performed in this topic. Keywords: Bridge Management; Bridge maintenance; Repair; Rehabilitation; Costs 1. INTRODUCTION Bridges cannot last forever. Whatever form of construction is used and whatever materials are adopted, sooner or later the effects of degradation begin to appear. A large percentage of existing infrastructure assets are deteriorating due to age, harsh environmental conditions, and insufficient capacity (Bordogna, 1995). Bridges are affected by their environment no less than people are, and the wear and tear of traffic, pollution, abuse, neglect, and just plain old age take their toll. It is implicit, and often made quite explicit, in the design of every product of engineering that there are limits to its health and strength, and therefore limits to Bridge Management what it can be subjected to. Recognition of those limits and regular check-ups and inspections of the artifact are required, as is a certain amount of preventive maintenance and repair. Bridges are considered to be vital links in any roadway network. Complete or partial failure to maintain these links paralyses the overall performance of the roadway network and causes excessive public and private losses. Therefore, bridge networks need to be managed in a way that ensures their uninterrupted performance throughout their design life. 2. LITERATURE VIEW Bridge management is the means by which a bridge stock is cared for from conception to the end of its useful life. Unfortunately many politicians and bridge authorities throughout the world whilst acknowledging the need for regular inspection and maintenance during the service life of their bridges, failed to appreciate the need for forward planning at the conception and design stages to ensure that sound principles were applied which would maximize their long term durability. The mechanism by which the coordination and implementation is achieved is the Bridge Management System (BMS) with the specific aims of assisting bridge managers and managing agencies: 1- To have a clear picture of all the bridges being managed and to prioritize them in terms of importance relative to the overall road and rail traffic infrastructure. 2To understand the maintenance needs of a particular bridge and by considering a number of intervention strategies to optimize the cost-benefit ratio. 3-To initiate and control the chosen maintenance action. 4To assess the value of the bridges on a periodic basis by the inclusion of performance indicators. At the heart of a BMS is the database built up using information obtained from the regular inspection and maintenance activities and containing a register of the

小分子芯片-SPR

小分子芯片-SPR

ORIGINAL ARTICLESmall molecule microarray screening methodology based on surface plasmon resonance imagingVikramjeet Singh a ,b ,*,Kuldeep Singh c ,Amita Nand a ,b ,f ,Huanqin Dai d ,Jianguo Wang e ,Lixin Zhang d ,Alejandro Merino a ,b ,Jingsong Zhu a ,baNational Center for Nanoscience and Technology,Beijing 100190,People’s Republic of China b University of Chinese Academy of Sciences,100049Beijing,People’s Republic of China cDepartment of chemistry,Maharishi Markandeshwar University,133207Ambala,India dChinese Academy of Sciences Key Laboratory of Pathogenic Microbiology &Immunology,Institute of Microbiology,CAS,Beijing 100190,People’s Republic of China eState-Key Laboratory and Institute of Elemento-Organic Chemistry,Nankai University,Tianjin 300073,People’s Republic of China fGuangzhou Xinren Biotechnology Co.,Ltd.,Guangzhou 510663,People’s Republic of ChinaReceived 16September 2014;accepted 13December 2014KEYWORDSSmall molecule microarray;Surface plasmon resonance;14-3-3f protein;Isatin and ligand–protein interactionAbstract In order to increase the scope and utility of small molecule microarrays (SMMs)we have combined SMMs and SPRi to screen small molecule antagonists against protein targets.Several small molecules,including immunosuppressive drugs (rapamycin and FK506)and reported inhib-itors (FOBISIN and Blapsin)of 14-3-3f proteins have been used to validate this technology.Fur-thermore,a small library of isatin derivatives have been synthesized and screened on developed platform against 14-3-3f protein.Three molecules,derived from the endogenous intermediate isatin termed,FZIB-35,FZIB-36and FZIB-38were identified as novel inhibitors which shows significant interaction with 14-3-3f .A mutation in the binding groove of 14-3-3f ,(K49E),almost abolishes the binding of these compounds to 14-3-3f protein.To exclude the probability of false positives,two more purified proteins (PtpA and BirA)were also tested.Furthermore,in order to confirm the bind-ing pocket specificity,competition assay against R18peptide was also carried out on presented plat-form.We show that SMMs in combination with SPRi are a powerful method to identify lead compounds in high throughput manner without the need to develop an activity based assay.ª2015The Authors.Production and hosting by Elsevier B.V.on behalf of King Saud University.This is an open access article under the CC BY-NC-ND license (/licenses/by-nc-nd/4.0/).1.IntroductionSmall molecule microarrays represent valuable tools for high throughput screening (HTS)in drug discovery (Kuruvilla et al.,2002)and enable the discovery of important and unknown protein–ligand interaction resulting in modulation of protein function (Koehler et al.,2003).SMMs in integration with cell based assay and confocal laser scanning microscopy*Corresponding author at:National Center for Nanoscience and Technology,Beijing 100190,People’s Republic of China.Tel.:+918901474914.E-mail address:kasana.chem@ (V.Singh).Peer review under responsibility of King SaudUniversity./10.1016/j.arabjc.2014.12.0201878-5352ª2015The Authors.Production and hosting by Elsevier B.V.on behalf of King Saud University.This is an open access article under the CC BY-NC-ND license (/licenses/by-nc-nd/4.0/).(CLSM)have been also described(Darvas et al.,2004;Molna r et al.,2013).To date,a number of elegant methods have been described for screening of small molecule inhibitors against protein targets.Conventional HTS methods such as TR-FRET,fluorescence polarization and ALPHAscreen face daunting challenges due to a number of limitations such as fluorescence interference,protein labeling,small molecule sol-ubility,and lengthy analysis times.Therefore,an alternative label free detection technology can be significantly advanta-geous.A great advantage of SPRi over classical SPR technique (Redman,2007)is throughput,allowing the parallel evaluation of hundreds or thousands of compounds simultaneously(Pillet et al.,2010).Moreover it provides a rapid identification of bio-molecular interaction along with their kinetic parameters in real time(Mcdonnell,2001).A variety of small molecules have been reported on SPRi for measuring protein–ligand interac-tion and protein–protein inhibition(Jung et al.,2005;Pillet et al.,2011).In this article,a combination of SMMs and SPRi has been used to detect ligand–protein interaction.Different strategies have been described for developing diverse linker systems on solid supports capable of anchoring small mole-cules(Hackler et al.,2003).A selective immobilization strategy was used for the fabrication of the SMMs through,either amino or hydroxy functional group of small -pounds were covalently captured on gold chip through simple EDC/NHS chemistry linked via PEG chains.Three different types of experiments were carried out to check the specificity of the ligands to the related target and to exclude false posi-tives.We validate this technology by using the interaction between FKBP12-Rapa-FK506and some known inhibitors of14-3-3f including the compounds FOBISIN(Zhao et al., 2011)and Blapsin(Yan et al.,2012).14-3-3proteins are a fam-ily of eukaryotic proteins that can bind to many phosphoser-ine/phospho-threonine containing signaling proteins such as kinases,phosphatases,and trans-membrane receptors (Aitken,2006).Hundreds of signaling and disease associated proteins including p53(Rajagopalan et al.,2010),C-Raf-1 (Molzan et al.,2010),BAD(Jiping et al.,1996),and histone deacetylases(Wang et al.,2000)have been documented to bind to14-3-3proteins.The dimeric14-3-3f isoform(Liu et al., 1996),in particular,is one of the most widely expressed and plays a major role in apoptosis.Additionally,a recent investi-gation identified the f isoform as a biomarker with high spec-ificity and sensitivity for the diagnosis and prognosis of head and neck cancer(Macha et al.,2010).Due to the involvement of14-3-3proteins in major cellular processes and diseases,cur-rent research has shifted toward the discovery of small mole-cule inhibitors which can provide good therapeutic opportunities.Over the last decade,a number of small mole-cule antagonists for14-3-3proteins have been studied(Yan et al.,2012)including some non-peptidic antagonists which act as inhibitors as well as stabilizers(Milroy et al.,2012).Cur-rently,there is no reported use of SMMs and SPRi in the dis-covery of new14-3-3proteins inhibitors.The main purpose of this research is to evaluate the SPRi technology for the screen-ing of small molecule inhibitors against14-3-3f.Further,a small library of compounds derived from isatin,which con-tained at least one NH2or OH group were immobilized and generate small molecule microarrays.Isatin is an endogenous Indole widely distributed in mammalian brain,peripheral tis-sue,and bodyfluids(Medvedev et al.,1996).14-3-3f represents one of these targets having specific and comparatively high interaction with isatin(Buneeva et al.,2010).Recently,an isat-in derivative has been reported(ID45)against coxsackievirus B3(CVB3)replication(Zhang et al.,2014).The primary screening of all isatin derivatives results3potential hits against 14-3-3f.Four different purified proteins,FKBP12,PtpA and BirA including K49E mutant of14-3-3f were tested against screened hits followed by competition approach against R18 peptide(Wang et al.,1999)shows promising inhibitory activity on SPR assay of identified compounds.In order to validate, these compounds further tested in ELISA and able to disrupt 14-3-3f interaction with its binding partner PRAS40protein. Combination of these two advanced technologies,SMMs and SPRi provides rapid screening and kinetics parameters of the tested inhibitors.We believe that this method can be applied for large scale primary screenings at low cost and with-out the need to develop an activity assay.2.Material and methods2.1.ReagentsUnless otherwise noted,material and solvents were obtained from commercial suppliers and used without further purifica-tion.Gold coated slides(Plexera),SH-(PEG)n-COOH(M.W. 1000)and SH-(PEG)n-OH(M.W.346)(Shanghai Yan Yi bio-tech.).EDC-HCl(1-(3-Dimethylaminopropyl)-3-ethylcarbodi-imide hydrochloride)and NHS(N-hydroxy succinimide), DMAP(N,N-dimethyl amino pyridine)(Aladdin Chemistry). DMSO,ethanol and ACN(Aldrich).Superblock solution was ordered from Thermo Scientific.FOBISIN101and FOBI-SIN106were purchased from Sigma.FKBP12protein was purchased from Sinobiological Inc.R18peptide,Blapsin inhibitors,isatin library(34compounds)and proteins such as14-3-3f,14-3-3f K49E mutant,PtpA BirA,were obtained from Prof.Lixin Zhang’s laboratory(Institute of Microbiol-ogy,Chinese Academy of Sciences).Synthesis procedure and NMR of identified inhibitors are presented in supplementary information.2.2.SMMs protocolA schematic representation for the screening process of SMMs is provided in Fig.1.Freshly deposited gold coated standard SPRi chips were cleaned with piranha solution(70%H2SO4/ 30%H2O2)for10min.The chips were extensively rinsed with Millipore water for30min.The chips were then immersed in ethanol containing1mM solution of SH-(PEG)n-COOH and SH-(PEG)n-OH(1:10)at4°C overnight and washed(shaker) in pure ethanol for30min before drying with nitrogen.Here we used the standard EDC/NHS chemistry for covalent immo-bilization of the small molecules on the surface of the chips.The carboxylic group(–COOH)from the SH-PEG-COOH was modified with a1:1mixture of EDC(0.39M)/NHS(0.1M). N-hydroxy succinimide ester is a robust chemistry widely uti-lized and able to attack amine and hydroxyl nucleophile groups (Ma dler et al.,2009)of small molecules and form stable amide and ester bonds pounds at10mM concentra-tion in100%DMSO were spotted into duplicate using a Genet-ix Qarray mini printer(contact mode printing)produces 250l M features,covalently immobilized on the sensor chip and blocked by superblock solution to minimize non-specificadsorption of proteins on the surface.A typical array image on PlexArrayÒHT system(Plexera)is shown in supplementary Fig.1.N,N-dimethyl amino pyridine(1uM)aq.solution wasadded to the printing solutions to facilitate nucleophile attack to form the desired ester bond.The slides were subsequently washed with DMSO,CAN,DMF,ethanol,PBS andfinally with distilled water for30min respectively to remove non-spe-cifically adsorbed compounds.2.3.SPRi methodAll the experiments were carried out using the PlexArrayÒHT system which is based on surface plasmon resonance imaging (Guan and Cong,2007).Small molecules containing at least one amino or hydroxy functional group are suitable to be immobilized using this strategy.Purified recombinant proteins, FKBP12,14-3-3f,14-3-3f(K49E),PtpA and BirA were in PBS pH7.4containing tween20(0.05%)and10%glycerol.Differ-ent concentrations of proteins were used as analyte.A solution of NaOH(10mM)was used to regenerate the surface and remove bound proteins from the SMMs enabling the sensor chip to be reused for additional analyte injections.All pre-sented data were repeated three times to derive the standard deviations.2.4.Binding experiments and data analysisAll the stock solutions of small molecules were stored in100% dimethyl sulphoxide(DMSO)atÀ20°C.Protein samples were stored in PBS with10%glycerol atÀ20°C.PBS was used as both assay and running buffer.A typical sample injection cycle consists of200s association phase with analyte solution and 300s dissociation phase with running buffer at3ul/sflow rate. Multiple concentrations of each protein14-3-3f(200,400and 600nM)and FKBP12(25,50and100nM)wereflowed on the SPRi instrument as analyte to get accurate kinetic parameters. Other purified proteins such as14-3-3f(K49E),PtpA and BirA were tested to confirm binding pocket specificity.The highest concentration tested for each protein was600nM.For data analysis,we used two software packages:data were analyzed according to our previous work(Singh et al.,2014).The spe-cific binding of protein to the immobilized small molecules was determined by subtracting the nonspecific physical adsorp-tion on reference spots using the Plexera SPR Data Analysis Module.3.Results3.1.High throughput screening of inhibitors by SPR imaging assayThe microarrays were then blocked and washed before expos-ing them to the purified recombinant proteins FKBP12and14-3-3f.As shown in Fig.2A,the Rapamycin and FK506spots bound the FKBP12protein specifically.Conversely,FOBISIN and Blapsin showed specific binding to14-3-3f,(Fig.2B).The resultant arrays can be regenerated with10mM aqueous NaOH solution and reused several times showing a great reproducibility.Unrelated compounds and surface back-Figure1Schematic representation of small molecule microarray.Figure2Identification of inhibitors by SPRi(A)SPRi graph showing interaction of Rapamycin and FK506with FKBP12 protein with FOBISIN as a negative control and(B)identification of FOBISIN and Blapsin inhibitors against14-3-3f protein(Rapa was taken as negative control)on SMMs platform.Figure3Identification and structure of inhibitors(A)SPRigraph showing interaction of three identified inhibitors,FZIB-38,FZIB-35and FZIB-36including R18as a positive control andrapamycin as a negative control and(B)chemical structure ofidentified inhibitors.Figure4Screening results against mutant and other unrelated proteins.(A)SPR response of all protein targets inhibitors and(B)response of new identified inhibitor toward all target proteins.(C)injection of14-3-3f protein followed injection shows complete abolishment of binding with known inhibitors and(D)new isatin inhibitors which the specific pocket of14-3-3f.Structural analysis of14-3-3f has determined that the amphipathic groove is the primary ligand binding site.The amphipathic groove lines up with the surface residue which is conserved between all isoforms of14-3-3proteins.Lys-49 is located in the conserved ligand binding site and plays a crit-ical role in ligand interaction((Zhang et al.,1997).Charge reversal mutation K49E in14-3-3f has shown to decrease its interaction with Raf-1kinase and thus with R18peptide (Wang et al.,1998).In order to demonstrate that the interac-tion of14-3-3f with the aforementioned compounds was via the specific binding pocket,we tested the14-3-3f(K49E) mutant.Two subsequent injections of14-3-3f and14-3-3f (K49E)separated by single regeneration wereflowed on a sin-gle chip.As shown in Fig.4C and D the binding of the14-3-3f (K49E)mutant to each inhibitor was dramatically reduced to negligible.This again strongly suggests that,known inhibitors including novel hits represent bonafide inhibitors that bind to the primary ligand binding site.petition assay on SPR imagingTo further confirm that the SMMs combined with SPRi can detect specific binding events of14-3-3f toward their inhibitors, a competition assay based on SPR imaging was developed.R18 is a high affinity peptide antagonist of14-3-3f protein which has strong interaction in the range of70nM.We used the R18peptide as a competitive inhibitor for the immobilized FOBISIN101,FOBISIN106,BLAP1,BLAP2,and BLAP3. 14-3-3f was injected either alone,or in a mixture with two con-centrations of the R18peptide(Zeta+R18_300nM and Zeta+R18_600nM).In all of three injections(Fig.5A and B),the concentrations of14-3-3f were constant(600nM). The mixture containing300nM R18peptide shows a dramatic reduction in the signal.The binding signal was almost negligi-ble when the concentration of R18was increased to600nM (Fig.5C and D).These data together with the lack of binding of the14-3-3f(K49E)mutant to the each inhibitor spots strongly support the ability of these compounds to disrupt functional interactions with relevant physiological partners. 3.4.Verification by ELISATo validate and see whether new hits screened from SPRi assay has some inhibition activity in solution,compounds were tested in ELISA.ELISA was performed in the same conditions used in the identification of the FOBISIN inhibitor of14-3-3 protein(see supplementary info.)by Dr.Haian Fu(Zhao et al.,2011)As a whole,ELISA analysis provides further evi-dence that these inhibitors can interrupt the interaction of14-3-3with PRAS40protein(Fig.6).However their IC50values in the low micromolar range,are3.92,5.44and5.47for FZIB-38,FZIB35and FZIB36respectively.It is important to note that the KD values determined by SMM-SPR method are in general lower than the corresponding IC50values reported in the literature for known inhibitors also.This could be due to either the enhanced affinity of the immobilized inhibitors on sensor surface or the relatively high concentrations required for protein–protein in vitro inhibition.Competition assay of all14-3-3f inhibitors(A)sensorgram showing competition assay against R18peptide.by two injections of same concentrations in addition to300nM and600nM of R18peptideknown inhibitors and(B)new identified isatin inhibitors to further confirm specific pocket phenomenon.known inhibitors and(D)identified inhibitors in completion assay.3.5.Kinetics analysis from SPR imagingDespite the fact that the kinetic parameters can change signifi-cantly upon the immobilization of the compounds,we mea-sured the kinetic parameters for all known compounds that bind FKBP12and 14-3-3f (Table 1).Here we used global fitting of a kinetic model in which a 1:1complex forms between inhib-itors and target proteins in data analysis module software.The data fit very well to this model;however,our values for kinetic rate constants determined from our SPRi experiments for Rap-amycin and FK506molecules are significantly different from the ones reported in the literature.This could be due to steric hindrance caused by the immobilization strategy Kinetics for known 14-3-3f inhibitors were not available in the literature.For all 14-3-3f inhibitors,only IC50values for protein–protein inhibition have been reported which is based on in vitro (FRET between KD and IC50is may be due to that IC50was obtained for protein–protein inhibition instead of direct measurement ligands affinity toward the target proteins.4.DiscussionWe have demonstrated here that small molecule microarray technology is quite useful in combination with surface plasmon resonance imaging for screening of small molecules modula-tors against targets of interest.Identification of three novel specific isatin derived compounds that showed potential utility as 14-3-3f inhibitors support this methodology.Furthermore,when these compounds were used in ELISA based 14-3-3f -PRAS40binding assay,all three compounds show promising activity suggested that presented methodology has the poten-tial to be used in high throughput manner without the need of development of an activity based assay that in some cases could be difficult to implement.However,during the course of this work,we realize that there is still a lot of room for improvement.Uniformity of spots and signal strength can be increased by trying different length of PEG linker.Photo-cross-linkers that bind randomly to any chemical group in a compound have proved to work well (unpublished data).This will allow the functional immobilization of larger sets of com-pounds that lack OH or NH 2groups or for which these groups arenecessary for binding to their targets.Another approach that facilitates the creation and functionality of SMMs is the use of 3dimensional surface chemistries instead of the 2dimensional surfaces utilized in this work.Although,this plat-form has some drawbacks at present,it has proved to be suit-able for screening of FKBP12and 14-3-3f ligands.Although this approach can also be used in conjunction with other exist-ing detection platforms including the use of fluorescence and microscopic readouts,we believe that the real time kinetics information gives this methodology a significant advantage.Low reagent requirements and rapid screening time make SMM technology particularly useful to academic and indus-trial discovery programs.The specificity and affinity obtained on this SMM platform can avoid long,laborious and costly efforts of primary screening in this field.Further developments on this technology are in progress in our laboratory.6Inhibition of 14-3-3f -PRAS40(PPIs)interaction identified inhibitors in ELISA.Table 2Kinetic parameters and IC50values of new identified inhibitors from SPRi and ELISA respectively.CompoundsProtein Ka (1/Ms)Kd (1/s)KA (1/M)KD (nM)IC50(l M)FZIB-3814-3-3f 5.07·103 2.8·10À4 1.81·10755.3±2.2 3.92FZIB-3514-3-3f 2.94·103 2.24·10À4 1.31·10776.6±3.8 5.44FZIB-3614-3-3f1.64·1032.93·10À41.24·10779.6±4.15.47Table 1Kinetic parameters of known inhibitors from SPRi.Compounds Protein Ka (1/Ms)Kd (1/s)KA (1/M)KD (nM)Rapamycin FKBP12 6.6·1041.87e À3 3.53·10728.2±2.3FK506FKBP12 4.35·104 2.35e À3 5.73·10754.1±2.44FOBISIN 10114-3-3f 1.18·104 5.64·10À4 2.09·10747.8±2.81FOBISIN 10614-3-3f 1.38·104 4.94·10À4 2.8·10735.8±2.1BLAP114-3-3f 1.39·104 5.54·10À4 2.5·10740±3.92BLAP214-3-3f 1.2·104 6.08·10À4 1.97·10750.8±3.76BLAP314-3-3f7.8·1038.54·10À49.14·106109±3.74AcknowledgmentsThis work wasfinancially supported by the following Grants: National Natural Science Foundation of China(Nos. 61077064/60921001)and National Major Scientific Instru-ments and Equipments Development Project(No. 2011YQ03012405).Appendix A.Supplementary materialSupplementary data associated with this article can be found, in the online version,at /10.1016/j.arabjc. 2014.12.020.ReferencesAitken,A.,2006.14-3-3proteins:a historic overview.Semin.Cancer Biol.16,162–172.Banaszynski,A.L.,Liu,W.C.,Wandless,J.T.,2005.Characterization of the FKBP–Rapamycin–FRB ternary complex.J.Am.Chem.Soc.127,4715–4721.Buneeva,O.,Gnedenko,O.,Zgoda,V.,Kopylov,A.,Glover,V., Lvanov,A.,Medvedev,A.,Archakov,A.,2010.Isatin-binding proteins of rat and mouse brain:proteomic identification and optical biosensor validation.Proteomics10,23–37.Darvas, F.,Dorma n,G.,Krajcsi,P.,Puska s,L.G.,Kova ri,Z., Lo rincz,Z.,Urge,L.,2004.Recent advances in chemical genomics.Curr.Med.Chem.11,3119–3145.Guan,H.W.,Cong,S.,2007.Surface Plasmon Resonance Spectrom-eter with an Actuator Driven Angle Scanning Patent,US20070222996A1.Hackler Jr.,L.,Dorma n,G.,Kele,Z.,Urge,L.,Darvas,F.,Puska s, L.G.,2003.Development of chemically modified glass surfaces for nucleic acid,protein and small molecule microarrays.Mol.Divers.7,25–36.Jiping,Z.,Hisashi,H.,Elizabeth,Y.,Jennifer,J.,Stanley,K.,1996.Serine phosphorylation of death agonist BAD in response to survival factor results in binding to14-3-3not BCL-XL.Cell87, 619–628.Jung,O.S.,Ro,H.S.,Kho,H.S.,Shin,Y.B.,Kim,M.G.,Chung,B., 2005.Surface plasmon resonance imaging-based protein arrays for high-throughput screening of protein–protein interaction inhibi-tors.Proteomics5,4427–4431.Koehler,A.N.,Shamji,A.F.,Schreiber,S.L.,2003.Discovery of an inhibitor of a transcription factor using small molecule microarrays and diversity-oriented synthesis.J.Am.Chem.Soc.125,8420–8842. Kuruvilla,F.G.,Shamji,A.F.,Sternson,S.M.,Hergenrother,P.J., Schreiber,S.L.,2002.Dissecting glucose signalling with diversity-oriented synthesis and small-molecule microarrays.Nature416, 653–657.Liu,D.,Bienkowska,J.,Petosa,C.,Collier,J.R.,Fu,H.,Liddington, R.,1996.Crystal structure of the zeta isoform of the14-3-3 proteins.Nature376,191–194.Macha,M.A.,Matta,A.,Chauhan,S.,Siu,K.M.,Ralhan,R.,2010.14-3-3zeta is a molecular target in guggulsterone induced apoptosis in Head and Neck cancer cells.BMC Cancer10,655.Ma dler,S.,Bich,C.,Touboul,D.,Zenobi,R.,2009.Chemical cross-linking with NHS esters:a systematic study on amino acid reactivities.J.Mass Spectrom.44,694–706.Mcdonnell,M.J.,2001.Surface plasmon resonance:toward an understanding of the mechanisms of biological molecular recogni-tion.Curr.Opin.Chem.Biol.5,572–577.Medvedev,A.E.,Clow,A.,Sandler,M.,Glover,V.,1996.Isatin:a link between natriuretic peptides and monoamines.Biochem.Pharma-col.52,385–391.Milroy,L.G.,Brunsveld,L.,Ottmann,C.,2012.Stabilization and inhibition of protein–protein interactions:the14-3-3case study.ACS Chem.Biol.8,27–35.Molna r,E.,Kuntam,S.,Cingaram,P.K.R.,Peksel,B.,Suresh,B., Fa bia n,G.,Fehe r,L.Z.,Bokros,A.,Medgyesi,A.,Ayaydin,F., Puska s,L.G.,bination of small molecule microarray and confocal microscopy techniques for live cell stainingfluorescent dye discovery.Molecules18,9999–10013.Molzan,M.,Schumacher,B.,Baljuls,A.,Polzien,L.,Weyand,M., Thiel,P.,Rose,M.,Kuhenne,P.,Kaiser,M.,Rapp,U.R., Kuhlmann,J.,Ottmann,C.,2010.Impaired binding of14-3-3to C-RAF in Noonan syndrome suggests new approaches in diseases with increased Ras signaling.Mol.Cell Biol.30,4698–4711. Pillet,F.,Thibault,C.,Bellon,S.,Maillart,E.,Tre visiol,E.,Vieu,C., Francois,J.M.,Anton Leberre,V.,2010.Simple surface chemistry to immobilize DNA probes that significantly increases sensitivity and spots density of surface plasmon resonance imaging based microarray systems.Sens.Actuators,B147,87–92.Pillet,F.,Romera,C.,Tre visiol,E.,Bellond,S.,Teulade-Fichoug, M.P.,Franc oisa,J.M.,Pratviele,G.,Leberre,V.A.,2011.Surface plasmon resonance imaging(SPRi)as an alternative technique for rapid and quantitative screening of small molecules,useful in drug discovery.Sens.Actuators,B157,304–309.Rajagopalan,S.,Sade,R.S.,Townsley,F.M.,Fersht,A.R.,2010.Mechanistic differences in the transcriptional activation of p53by 14-3-3isoforms.Nucleic Acids Res.238,893–906.Redman,J.E.,2007.Surface plasmon resonance for probing quadru-plex folding and interactions with proteins and small molecules.Methods43,302–312.Singh,V.,Nand,A.,Cheng,Z.,Yang,M.,Zhu,J.,2014.3D small molecule microarray with enhanced sensitivity and immobilization capacity monitored by surface plasmon resonance imaging.RSC Adv./10.1039/C4RA07306A.Wang,H.,Zhang,L.,Liddington,R.,Fu,H.,1998.Mutations in the hydrophobic surface of an amphi-pathic groove of14-3-3disrupt its interaction with Raf-1kinase.J.Biol.Chem.273,16297–16304. Wang,B.,Yang,H.,Liu,C.Y.,Jelinek,T.,Zhang,L.,Ruoslahti,E., Fu,H.,1999.Isolation of high-affinity peptide antagonists of14-3-3proteins by phage display.Biochemistry38,12499–1254. Wang, A.H.,Kruhlak,M.J.,Wu,J.,Bertos,N.R.,Vezmar,M., Posner,B.I.,Bazett-Jones,D.P.,Yang,X.J.,2000.Regulation of histone deacetylase4by binding of14-3-3proteins.Mol.Cell.Biol.20,6904–6912.Yan,Y.,Dai,H.,Du,Y.B.,Guo,H.,Li,D.,Zhang,L.,Fu,H.,Dong, X.,Cheng,Y.,2012.Identification of Blapsins A and B as potent small-molecule14-3-3inhibitors from the insect Blaps Japanensis.Bioorg.Med.Chem.Lett.22,4179–4181.Zhang,L.,Wang,H.,Liu,D.,Liddington,R.,Fu,H.,1997.Raf-1 kinase and exoenzyme S interact with14-3-3f through a common site involving lysine49.J.Biol.Chem.272,13717–13724. Zhang,H.M.,Dai,H.,Hanson,P.J.,Li,H.,Guo,H.,Ye,X.,Hemida, M.G.,Wang,L.,Tong,Y.,Qiu,Y.,Liu,S.,Wang,F.,Song,F., Zhang,B.,Wang,J.G.,Zhang,L.X.,Yang,D.,2014.Antiviral activity of an isatin derivative via induction of PERK-Nrf2-mediated suppression of cap-independent translation.ACS Chem.Biol.9,1015–1024.Zhao,J.,Du,Y.,Horton,J.,Upadhyay,A.,Lou,B.,Bai,Y.,Zhang, X.,Du,L.,Li,M.,Wang,B.,Zhang,L.,Barbieri,J.,Khuri,F., Cheng,X.,Fu,H.,2011.Discovery and structural characterization of a small molecule14-3-3protein–protein interaction inhibitor.Proc.Natl.Acad.Sci.108,16212–16216.。

管理科学奖历届名单

管理科学奖历届名单

管理科学奖历届名单管理科学奖(The Management Science Award)是美国科学界最高荣誉之一,颁发给在管理科学领域做出卓越贡献的人士。

自1957年设立以来,已经颁发了37次,以下是其历届得主名单:1957年:H. Irving Grousbeck1958年:E. P. Hollander1959年:Paul E. Green、Franklin S. Forsythe1960年:George Dantzig、William W. Cooper、Tjalling Koopmans、Richard Bellman1961年:Harold Hotelling1962年:No Winner1963年:Leonid Hurwicz、Kenneth Arrow、Paul Samuelson1964年:J. Richard Nisbett1965年:Donald F. Morrison1966年:Edward J. Green1967年:Melvin J. Hinich、Kenneth J. Arrow1968年:C. R. Rao1969年:A. E. Lawson、David McNeill1970年:William S. Cooper、Abraham Charnes1971年:Harold W. Kuhn、Martin Shubik1972年:Carl R. Kordesch、Arnold O. Allen1973年:Franco Modigliani1974年:No Winner1975年:Robert M. Thrall1976年:Fred W. Glover1977年:Oliver E. Williamson1978年:Leonard D. White、Paul J. Zarembka1979年:Lester B. Lave、Jay W. Forrester1980年:George P. Richardson、Edwin R. Crowell、Jay W. Forrester1981年:No Winner1982年:Richard M. Cyert、James G. March1983年:Harry M. Markowitz、Merton H. Miller、William F. Sharpe1984年:Ralph L. Keeney1985年:Howard Raiffa1986年:Arnold Ross、David S. Sibley、Ralph E. Gomory1987年:Kenneth J. Arrow、John C. Harsanyi、John F. Nash Jr. 1988年:Richard C. Larson、John A. Howard1989年:Andrew Whinston、Benjamin M. Blau1990年:Robert A. Muhlnickel、Bernard W. Taylor III、Ralph P. Sawyer1991年:Richard M. Cyert、Robert J. Weber1992年:Marshall W. Meyer、John H. McArthur1993年:J. Edward Russo、Paul J. H. Schoemaker1994年:John R. Hauser1995年:Danny Samson、Ravi Aron1996年:David L. Olson、Jie Wu1997年:George S. Day、David J. Reibstein1998年:Henry Mintzberg、John C. Wood1999年:W. Chan Kim、Renée Mauborgne2000年:No Winner2001年:Michael I. Jordan、David E. Rumelhart、Yoshua Bengio 2002年:No Winner2003年:Richard H. Thaler2004年:No Winner2005年:D. Clay Whybark、Jayashankar M. Swaminathan2006年:Rajeev K. Tyagi、Surendra M. Gupta、Rajesh K. Tyagi 2007年:Robert S. Kaplan2008年:Edward H. Bowman、Tom J. Brown Jr. 2009年:Robert J. Gordon2010年:No Winner2011年:Robert B. Wilson2012年:Allen S. Lee、Michael J. Shaw 2013年:F. Thomas Juster、Richard Suzman 2014年:No Winner2015年:No Winner2016年:No Winner以上是管理科学奖历届得主名单。

The Decision Reliability of MAP, Log-MAP,

The Decision Reliability of MAP, Log-MAP,

The Decision Reliability of MAP, Log-MAP, Max-Log-MAP and SOV A Algorithmsfor Turbo CodesAbstract —In this paper, we study the reliability of decisions ofe Codes, Channel Reliability,e N comm llular, satellite and we also consider two improved versions, named Log-MAP two different or identicalRecursi s, connectedin pFig. 1. The turbo encoder with rate 1/3.The first encoder operat ed b e u , i ond encoderp Lucian Andrei Peri şoar ă, and Rodica Stoianth MAP, Log-MAP, Max-Log-MAP and SOVA decoding algorithms for turbo codes, in terms of the a priori information, a posteriori information, extrinsic information and channel reliability. We also analyze how important an accurate estimate of channel reliability factor is to the good performances of the iterative turbo decoder. The simulations are made for parallel concatenation of two recursive systematic convolutional codes with a block interleaver at the transmitter, AWGN channel and iterative decoding with mentioned algorithms at the receiver.Keywords —Convolutional Turbo D cision Reliability, Extrinsic Information, Iterative Decoding.I. I NTRODUCTIONunication systems, like ce computer fields, the information is represented as a sequence of binary digits. The binary message is modulated to an analog signal and transmitted over a communication channel affected by noise that corrupt the transmitted signal.The channel coding is used to protect the information fromnoise and to reduce the number of error bits.One of the most used channel codes are convolutional codes, with the decoding strategy based on the Viterbialgorithm. The advantages of convolutional codes are used inTurbo Codes (TC), which can achieve performances within a2 dB of channel capacity [1]. These codes are parallelconcatenation of two Recursive Systematic Convolutional (RSC) codes separated by an interleaver. The performances of the turbo codes are due to parallel concatenation ofcomponent codes, the interleaver schemes and the iterative decoding using the Soft Input Soft Output (SISO) algorithms [2], [3].In this paper we study the decision reliability problem for turbo coding schemes in the case of two different decodingstrategies: Maximum A Posteriori (MAP) algorithm and Soft Output Viterbi Algorithm (SOVA). For the MAP algorithmand Max-Log-MAP algorithms. The first one is a simplified algorithm which offers the same optimal performance with a reasonable complexity. The second one and the SOVA are less complex again, but give a slightly degraded performance. The paper is organized as follows. In Section II, the turbo encoder is presented. In Section III, the turbo decoder is ex Manuscript received December 10, 2008. This work was supported in part by the Romanian National University Research Council (CNCSIS) under theGrant type TD (young doctoral students), no. 24.L. A. Peri şoar ă is with the Applied Electronics and InformationEngineering Department, Politehnica University of Bucharest, Romania (e-mail: lucian@orfeu.pub.ro, lperisoara@, www.orfeu.pub.ro).R. Stoian is with the Applied Electronics and Information Engineering Department, Politehnica University of Bucharest, Romania (e-mail: rodica@orfeu.pub.ro, rodicastoian2004@, www.orfeu.pub.ro).plained in detail, presenting firstly the iterative decoding principle (turbo principle), specifying the concepts of a priori information, a posteriori information, extrinsic information, channel reliability and source reliability. Then, we review the MAP, Log-MAP, Max-Log-MAP and SOVA decoding algorithms for which we discuss the decision reliability. In Section IV is analyzed the influence of channel reliability factor on decoding performances for the mentioned decoding algorithms. Section V presents some simulation results, which we obtained.II. T HE T URBO C ODING S CHEME The turbo encoder can use ve Systematic Convolutional (RSC) code arallel, see Fig. 1.es on the input bits represent n their original order, while the sec y the fram o erates on the input bits which are permuted by the interleaver, frame u ’, [4]. The output of the turbo encoder is represented by the frame: I2)()()1211,12,121,22,21,2,,,,,,,,,...,,,k k k u c c u c c u c c ==v u c c /R k n = to b , (1)is less likely where frame c1 is the output of the first RSC and frame c2 is the output of the second RSC. If the input frame u is of length k and the output frame x is of length n , then the encoder rate is .For block encoding data is segmented into non-overlapping blocks of length k with each block encoded (and decoded)independently. This scheme imposes the use of a blockinterleaver with the constraint that the RSC’s must begin in the same state for each new block. This requires either trellis termination or trellis truncation. Trellis termination need appending extra symbols (usually named tail bits) to the inputframe to ensure that the shift registers of the constituent RSC encoders starts and ends at the same zero state. If the encoder has code rate 1/3, then it maps k data bits into 3k coded bits plus 3m tail bits. Trellis truncation simply involves resettingthe state of the RSC’s for each new block.The interleaver used for parallel concatenation is a device that permutes coordinates either on a block basis (a generalized “block” interleaver) or on a sliding window basis(a generalized “convolutional” interleaver). The interleaver ensures that the set of code sequences generated by the turbo code has nice weight properties, which reduces the probabilitythat the decoder will mistake one codeword for another.The output codeword is then modulated, for example with Binary Phase Shift Keying (BPSK), resulting the sequence , which is transmitted over an Additive White Gaussian Noise (AWGN) channel.(12,,=v u c c 12,)p x x )(,s p =x x e a low weight codeword due to the interleaver in front of it. The interleaver shuffles the inputsequence It is known that turbo codes are the best practical codes due to their performance at low SNR. One reason for their better performance is that turbo codes produce high weight code words [4]. For example, if the input sequence u is originally low weight, the systematic u and parity c 1 outputs mayproduce a low weight codeword. However, the parity output c 2 is less likely to be a low weight codeword due to the u , in such a way that when introduced to the second encoder, it is more likely probable to produce a high weight codeword. This is ideal for the code because high weight code words result in better decoder performance. III. T HE T URBO D ECODING S CHEME Let be the received sequence of length n , 12(,,)s p p =y y y y where the vector y s is formed only by the received informationsymbols s y 222222(,,...,)p p p p n y y y =y p 1 and y p 2and . These three streams are applied to the input of the turbo decoder presented in Fig. 2. 11112(,,...,)p p p p n y y y 1=y y At time j , decoder 1 using partial received information 1,s p j j y y (), makes its decision and outputs the a posterioriinformation s j L x +()()()e s s s s . Then, the extrinsic information is computed j j j c jL x L x L x L y +−=−−. Decoder 2 makes itsdecision based on the extrinsic information ()e sj L x 2 from decoder 1 and the received information ',s p j jy y . The term(')s j L x + is the a posteriori information derived from decoder 2 and used by decoder 1 as a priori information about thereceived sequence, noted with (')sj L x −(). Now, the second iteration can begin, and the first decoder decodes the same channel symbols, but now with additional information about the value of the input symbols provided by the second decoder in the first iteration. After some iterations, the algorithm converges and the extrinsic information values remains the same. Now the decision about the message bits u j is made based on the a posteriori values s j L x +.e s y p 2y p 1y sFig. 2. The turbo decoder.Each constituent decoder operates based on the Logarithm Likelihood Ratio (LLR).A. The Decision Reliability of MAP DecoderBahl, Cocke, Jelinek and Raviv proposed the Maximum APosteriori (MAP) decoding algorithm for convolutional codesin 1974 [1]. The iterative decoder developed by Berrou et al.[5] in 1993 has a greatly increased attention. In their paper,they considered the iterative decoding of two RSC codesconcatenated in parallel through a non-uniform interleaver and the MAP algorithm was modified to minimize the sequence error probability instead of bit error probability.Because of its increased complexity, the MAP algorithm was simplified in [6] and the optimal MAP algorithm calledthe Log-MAP algorithm was developed. The LLR of a transmitted bit is defined as [2]:(1)()log ()(1)s Wenoted def j s sj j s j P x L x L x P x −⎛⎞=+==⎜⎟⎜⎟=−⎝⎠where the sign of the LLR ()s j L x indicate whether the bit s j xis more likely to be +1 or -1 and the magnitude of the LLRgives an indication of the correct value of s j x . The term()sj L x − is defined as the a priori information about s j x .In channel coding theory we are interested in theprobability that , based or conditioned on some received sequence 1s j x =±s j y . Hence, we use the conditional LLR: ()()()1||log (1|s s We noted def j j s s s j j j s s j j P x y L x y L x P x y +⎛⎞=+⎜⎟=⎜⎟=−⎝⎠=) The conditional probabilities (1|s sj j P x y =± are the a posteriori probabilities of the decoded bit s j x and ()s j L x + is thea posteriori information about sj x , which is the information that the decoder gives us, including the received frame, the a priori information for the systematic symbols y s j and the apriori information for symbol x s j . It is the output of the MAPalgorithm. In addition, we will use the conditional LLR ()|s s j j L y x based on the probability that the receiver’s output would be s j y when the transmitted bit s j x was either +1 or -1:()()()|1|log |1s s defj j s s jjs s j j P y x L y x P y x ⎛⎞=+⎜=⎜=−⎝⎠⎟⎟. (3)For AWGN channel using BPSK modulation, we can write the conditional probability density functions, [7]:()()20|12s s b j j j EP y x y a N ⎡⎤=±=−⎢⎣⎦m ⎥, (4)where is the transmitted energy per bit, a is the fadingamplitude and is the noise variance.b E 0/2N We can rewrite the (3) as follows: ()()()2200|4,s s s s b j j j j Noteds s b j c j E L y x y a y a N E a y L y N ⎡⎤=−−−+⎢⎥⎣⎦== (5) the fading amplitude and is the noise power. For nonfading AWGN channels a = 1 and 0N /204c b L E N =. Theratio is defined as the Signal to Noise Ration (SNR) of thechannel.0/b E N The extrinsic information can be computed as [1], [2], [9]: ()()()()()()1|()log 1|1|log log 1|()().s s j j e sj s s jj s s j j s sj j s s s j j c j P x y L x P x y P x P y x P x P y x L x L x L y +−⎛⎞=+⎜⎟=⎜⎟=−⎝⎠⎛⎞⎛=+=+⎜⎟⎜−−⎜⎟⎜=−=−⎝⎠⎝=−−11s j s j ⎞⎟⎟⎠ (6)The a posteriori information defined in (2), can be written asthe following [1], [10]:11(')()(',)()log (')()(',)e j j j s j e j j j s s s s L x s s s s −++−−α⋅β⋅γ=α⋅β⋅γ∑∑, (7)where +∑is the summation over all possible transition branch pairs (s ’,s ) in the trellis, at time j , given the transmittedsymbol x s j = +1. Analog, −∑is for transmitted symbol x s j =-1.The forward and backward terms, represented in Fig. 3 as transitions between two consecutive states from the trellis, can be computed recursively as following [7], [10], [11]:1'()(')(',)j j j s s s s s −α=αγ∑, (8)1(')()(',)j j j ss s s s −β=βγ∑. (9)For systematic codes, which is our case, the branch transition probabilities (',)js s γ are given by the relation:11(',)exp ()(',)22s s s s e j j j c jj j s s L x x L x y s −⎡γ=+⋅γ⎢⎣⎦s ⎤⎥, (10) where:112211(',)exp 22e p p j c j j c p p j j s s L x y L x ⎡⎤γ=+⎢⎥⎣⎦y .(11)At each iteration and for each frame y, ()s j L x + is computedat the output of the second decoder and the decision is done,symbol by symbol j = 1…k , based on the sign of ()sj L x +, original information bit u j being estimated as [2], [3]: {ˆ()sj usign L x +=}j . (12) In the iterative decoding procedure, the extrinsicinformation ()e s j L x is permuted by the interleaver andbecomes the a priori information ()sj L x − for the next decoder. influence on ()s j L x + is insignificant.B. The Decision Reliability of Max-Log-MAP DecoderThe MAP algorithm as described in previous section is much more complex than the Viterbi algorithm and with hard decision outputs performs almost identically to it. Therefore for almost 20 years it was largely ignored. However, its application in turbo codes renewed interest in this algorithm. Its complexity can be dramatically reduced without affecting its performance by using the sub-optimal Max-Log-MAP algorithm, proposed in [12]. This technique simplifies the MAP algorithm by transferring the recursions into the log domain and invoking the approximation: ln max()i x i ii e x ⎛⎞≈⎜⎟⎝⎠∑. (13)where max()i i x means the maximum value of x i . If we note:()()ln ()j j A s =αs , (14)()()ln ()j j B s s =β, (15)and:()(',)ln (',)j j G s s s s =γ, (16)then the equations (8), (9) and (10) can be written as: ()(()1'1'1'()ln ()ln (')(',)ln exp (')(',)max (')(',),j j j j s j j s j j s )A s s s s A s G s s A s G s s −−−⎛⎞=α=αγ⎜⎟⎝⎠⎛=+⎜⎝⎠≈+∑∑s ⎞⎟⎞⎟(17) ()()()11(')ln (')ln ()(',)ln exp ()(',)max ()(',),j j j j s j j s j j s B s s s s s B s G s s B s G s s −−⎛⎞=β=βγ⎜⎟⎝⎠⎛=+⎜⎝⎠≈+∑∑ (18) 11(',)()22s s s s jj j c j G s s C x L x L x y −=++j , (19) term ()s s j j x L x −.Finally, the a posteriori LLR ()s j L x + which the Max-Log-MAP algorithm calculates is:Fig. 3. Trellis states transitions.for ()j s αfor 1(')j s −β((1(',)11(',)1()max(')(',)()max (')(',)().j j s j j j j s s for u j j j s s for u L x As G s s B s ))A s G s s B s +−=+−=−≈++−++ (20)In [12] and [13] the authors shows that the complexity of Max-Log-MAP algorithm is bigger than two times that of a classical Viterbi algorithm Unfortunately, the storage requirements are much greater for Max-Log-MAP algorithm, due to the need to store both the forward and backward recursively calculated metrics and before the ()j A s ()j B s ()s j L x + values can be calculated.C. The Decision Reliability of Log-MAP DecoderThe Max-Log-MAP algorithm gives a slight degradation in performance compared to the MAP algorithm due to the approximation of (13). When used for the iterative decodingof turbo codes, Robertson found this degradation to result in a drop in performance of about 0.35 dB, [12]. However, the approximation of (13) can be made exact by using the Jacobian logarithm:()(()121212121212ln()max(,)ln 1exp ||max(,)||(,),x x e e x x x x )x x f x x g x x +=++−−=+−= (21)where ()f δ can be thought of as a correction term. However,the maximization in (17) and (18) is completed by the correction term ()f δ in (21). This means that the exact ratherthan approximate values of and are calculated. For binary trellises, the maximization will be done only for two terms. Therefore we can correct the approximations in (17) and (18) by adding the term ()j A s ()j B s ()f δ, where δ is the magnitude of the difference between the metrics of the twomerging paths. This is the basis of the Log-MAP algorithmproposed by Robertson, Villebrun and Hoeher in [12]. Thus we must generalize the previous equation for more than two 1x terms, by nesting the 12(,)g x x operations as follows: (((13211ln ,,,(,)i n x n n i e g x g x g x g x x −=⎛⎞=⎜⎟⎝⎠∑K ))), (22)The correction term ()f δδ need not to be computed for every value of , but instead can be stored in a look-up table. In [12], Robertson found that such a look-up table need containonly eight values for , ranging between 0 and 5. This meansthat the Log-MAP algorithm is only slightly more complexthan the Max-Log-MAP algorithm, but it gives exactly the same performance as the MAP algorithm. Therefore, it is a very attractive algorithm to use in the component decoders of an iterative turbo decoder. δD. The Decision Reliability of SOVA DecoderThe MAP algorithm has a high computational complexityfor providing the Soft Input Soft Output (SISO) decoding. However, we obtain easily the optimal a posteriori probabilities for each decoded symbol. The Viterbi algorithm provides the Maximum Likelihood (ML) decoding for convolutional codes, with optimalsequence estimation. The conventional Viterbi decoder has two main drawbacks for a serial decoding scheme: the inner Viterbi decoder produces bursts of error bits and hard decision output, which degrades the performance of the outer Viterbi decoder [3]. Hagenauer and Hoeher modified the classical Viterbi algorithms and they provided a substantially less complex and suboptimal alternative in their Soft OutputViterbi Algorithm (SOVA). The performance improvement is obtained if the Viterbi decoders are able to produce reliability values or soft outputs by using a modified metric [14]. These reliability values are passed on to the subsequent Viterbi decoders as a priori information .In soft decision decoding, the receiver doesn’t assign a zero or a one to each received symbol from the AWGN channel, but uses multi-bit quantized values for the received sequence y , because the channel alphabet is greater than the sourcealphabet [3]. In this case, the metric derived from Maximum Likelihood principle, is used instead of Hamming distance. For an AWGN channel, the soft decision decoding produces again of 2÷3 dB over hard decision decoding, and an eight-level quantization offers enough performance in comparison with an infinite bit quantization [7].The original Viterbi algorithm searches for an informationsequence u that maximizes the a posteriori probability, s being the states sequence generated by the message u . Using the Bayes theorem and taking into account that thereceived sequence y is fixed for the metric computation and it can be discarded, the maximization of is: (|)P s y (|)P s y {}{max (|)max (|)()P P =u us }P y y s s . (23)For a systematic code, this relation can be expanded to:(1211max (,,)|,()k s p p j j j j j j j P y y y s s P s −=)⎧⎫⎪⎪⎨⎬⎪⎪⎩⎭∏u. (24) Taking into account that:()()()(1211122(,,)|,|||s p p j j j j j s s p p p p j j j j j j P y y y s s P y x P y x P y x −==⋅⋅), (25)where 1(,)j j s s − denotes the transitions between the states attime j -1 and the states at time j , the SOVA metric is obtained from (24) as [15]:()()***1***|1(1)log log ,(0)|1j j j j j j j j j jP y x P u M M x u P u P y x −⎛⎞=+⎛⎞=⎜⎟=++⎜⎟⎜⎟⎜⎟==−⎝⎠⎝⎠∑ (26)where *1,2,(,,)j j j j x u c c = is the RSC output code word at timej , at channel input and *1(,,)s p p j j j j 2y y y y = is the channeloutput. The summation is made for each pair of information symbols (,s j j u y ) and for each pair of parity symbols (11,,p j j c y )and (2,2,p j j y 1*c ).According [14] and [7], the relation (26) can be reduced as: **c j j ()j j j j M M L −=+∑x y u L u +(), (27)where the source reliability j L u , defined in (26), is the log-likelihood ratio of the binary symbol u j . The sign of ()j L u ) is the hard decision of u j and the magnitude of (j L u is the decision reliability .According [10], the SOVA metric includes values from the past metric M j -1, the channel reliability L c and the source reliability ()j L u (, as an a priori value. If the channel is very good, the second term in (27) is greater than the third term andthe decoding relies on the received channel values. If thechannel is very bad, the decoding relies on the a priori information )j L u . If M 1j , M 2j are two metrics of the survivor path and concurrent path in the trellis, at time j , then the metric difference is defined as [7]:01212j j j M M −)(s m Δ=. (28)The probability of path m , at time j , is related as:()/2mjM (path )exp m j P m P ==. (29) where j s is a states vector and mj M is the metric. The probability of choosing the survivor path is: 001)(path (correc ath 1)(path 2)1jjP e P P P e ΔΔ==++t)(p . (30)The reliability of this path decision is calculated as:(correct)orrect)log 1-(c j P P =Δ. (31) The reliability values along the survivor paths, for aparticular node and time j , are denoted as d j Δ, where d is the distance from the current node at time j . If the survivor path bit for is the same with the associated bit on the competing path, then there would be no error if the competing path is chosen. The reliability value remains unchanged.d j =To improve the reliability values an updating process must be used, so the “soft” values of a decision symbol are:(')'di j d j di L u u −−=j=Δ∑, (32)which can be approximated as:{}0...(')'min i j d j d i d L u u −−=j =⋅Δ. (33)The SOVA algorithm described in this section is the least complex of all the SISO decoders discussed in this section. In [12], Robertson shows that the SOVA algorithm is about halfas complex as the Max-Log-MAP algorithm. However, theSOVA algorithm is also the least accurate of the algorithmsdescribed in this section and, when used in an iterative turbo decoder, performs about 0.6 dB worse than a decoder using the MAP algorithm. If we represent the outputs of the SOVA algorithm they will be significantly more noisy than thosefrom the MAP algorithm, so an increased number of decodingiterations must be used for SOVA to obtain the sameperformances as for MAP algorithm.The same results are reported also for the iterative decoding (turbo decoding) of the turbo product codes, which are basedon two concatenated Hamming block codes not on convolutional codes [19]. IV. T HE INFLUENCE OF L C ON DECODING PERFORMANCE In this section we analyze the importance of an accurate estimate of the channel reliability factor L c is to the good performance of an iterative turbo decoder which uses the MAP, SOVA, Max-Log-MAP and Log-MAP algorithms. In the MAP algorithm the channel inputs and the a priori information are used to calculate the transition probabilities from one state to another, that are then used to calculate theforward and backward recursion terms [2], [8]. Finally, the aposteriori information ()s j L x + is computed and the decision about the original message is made based on it. In the iterative decoding with MAP algorithm, the channelreliability is calculated from the received channel values. At first iteration, the decoder 1 has no a priori information available (the ()s j L x − is zero) and the output from thealgorithm is calculated based on channel values. If an incorrect value of L c is used the decoder will make more decision errors and the extrinsic information from the output of the first decoder will have incorrect values, for the softchannel inputs [16].In the SOVA algorithm the channel values are used torecursively calculate the metric *c j L y j M for the current state s along a path from the metric 1j M − for the previous state along that path added to an a priori information term and to a cross-correlation term between the transmitted and the receivedchannel values, *j x and *j y , using (27). The channel reliability factor is used to scale this cross-correlation. When we usec Lan incorrect value of , e.g. , we are scaling the channel values applied to the inputs of component decoders by a factor of one instead of the correct value of . This has the effect of scaling all the metrics by the same factor, see (8), and the metric differences are also scaled by the same factor, see (9). This scaling of the metrics do not affect the path chosen by the algorithm as a survivor path or as a Maximum Likelihood (ML) path, so the hard decisions given by the algorithm are not affected by using an incorrect value of L c [16]-[18].c L ()j B s 1c L =c L c In the iterative decoding with SOVA algorithm, in the first iteration we assume that no a-priori information about the transmitted bits is available to the decoder (the a-priori information is zero), the first component decoder takes only the channel values. If channel reliability factor is incorrect, the channel values are scaled, the extrinsic information will be also scaled by the same factor and the a-priori information for the second decoder will also be scaled. Because of the linearity of the SOVA, the effect of using an incorrect value of the channel reliability factor is that the output LLR from the decoder is scaled by a constant factor. The relative importance of the two inputs to the decoder, the a priori information and the channel information, will not change, since the LLRs for both these sources of information will be scaled by the same factor. In the final iteration, the soft outputs from the final component decoder will have the same sign as those that would have been calculated using the correct value of . So, the hard outputs from the turbo decoder using the SOVA algorithm are not affected by the channel reliability factor [16].L c L The Max-Log-MAP algorithm has the same linearity that is found in the SOVA algorithm. Instead of one metric, now two metrics and are calculated, for forward andbackward recursions, see (17), (18) and (19), were are used only simple additions of the cross-correlation of the transmitted and received symbols. But, if an incorrect value of the channel reliability value is used, all the metrics are simply scaled by a factor as in the SOVA algorithm. The soft outputs given by the differences in metrics between different paths will also be scaled by the same factor, with the sign unchanged and the final hard decisions given by the turbo decoder will not be affected.()j A s The Log-MAP algorithm is identical to the Max-Log-MAP algorithm, except for a correction term ()()ln exp()f δ=−δ1+, used in the calculation of the forward and backward metrics and ()j A s ()j B s , and the soft output LLRs. The function()f δ is not a linear function, it decreases asymptoticallytowards zero as δ increases. Hence the linearity that is present in the Max-Log-MAP and SOVA algorithms is not present in the Log-MAP algorithm. This effect of non-linearity determines more hard decision errors of thecomponent decoders if the channel reliability factor is incorrect, and the extrinsic information derived from the first component decoder have incorrect amplitudes, which become the a-priori information for the second decoder in the first iteration. Both decoders in subsequent iterations will have incorrect amplitudes relative to the soft channel inputs.c L In the iterative decoding with Log-MAP algorithm, the extrinsic information exchange from one component decoder to another, determines a rapid decrease in the BER as the number of iterations increases. When the incorrect value of is used, no such rapid fall in the BER occurs due to the incorrect scaling of the a priori information relative to the channel inputs. In fact, the performance of the decoder is largely unaffected by the number of iterations used.c L For wireless communications, some of them modeled as Multiple Input Multiple Output (MIMO) systems [23], the channel is considered to be Rayleigh or Rician fading channel. If the Channel State Information (CSI) is not known at the receiver, a natural approach is to estimate the channel impulse response and to use the estimated values to compute the channel reliability factor required by the MAP algorithm to calculate the correct decoding metric.c L In [20], the degradation in the performance of a turbo decoder using the MAP algorithm is studied when the channel SNR is not correctly estimated. The authors propose a method for blind estimation of the channel SNR, using the ratio of the average squared received channel value to the square of the average of the magnitudes of the received channel values. In addition, they showed that using these estimates for SNR gives almost identical performances for the turbo decoder to that given using the true SNR.In [8], the authors proposes a simple estimation scheme for from the statistical computation on the block observation of matched filter outputs. The channel estimator includes the error variance of the channel estimates. In [24], is used the Minimum Mean Squared Error (MMSE) estimation criterion and is studied an iterative joint channel MMSE estimation and MAP decoding.c L None of above works requires a training sequence with pilot symbols to estimate the channel reliability factor. Other studies used pilot symbols to estimate the channel parameters, like [22] and [25].In [22] it is shown that it is not necessary to estimate the channel SNR for a turbo decoder with Max-Log-MAP or SOVA algorithms. If the MAP or the Log-MAP algorithm is used then the value of does not have to be very close to the true value for a good BER performance to be obtained. c LV. S IMULATION RESULTSThis section presents some simulation results for the turbo codes ensembles, with MAP, Max-Log-MAP, Log-MAP and SOVA decoding algorithms. The turbo encoder is the same for the four decoding algorithms and is described by two identical RSC codes with constraint length 3 and the generator polynomials and . No tail bitsand no puncturing are performed. The two constituent encoders are parallel concatenated by a classical block interleaver, with dimensions variable according to the frame21f G =+D D 21b G D =++。

Bibliography

Bibliography[1]Alamir,M.&Bornard,G.(1994)On the stability of receding horizon control of nonlinear discrete-time systems.Systems&Control Letters23,291-296[2]Etay,J.,M.Alamir&Y.Fautrelle(1994)Dynamics of electromagnetically driven spheres inparison between numerical simulations and experiments.Magnetohydrodynamics, 30,No3,389-400[3]Alamir,M.&Bornard,G.(1995)Stability of a truncated infinite constrained receding horizonscheme:The general nonlinear case.Automatica31,No9,1353-1356[4]Alamir,M.&Bornard,G.(1996)Globally stabilizing discontinuous feedback law for a class ofnonlinear systems.European Journal of Control2,44-50[5]Alamir,M.(1999)Optimization Based Nonlinear Observers Revisited.International Journalof Control72,No13,1204-1217[6]Alamir,M.&Balloul,I.(1999)Robust Constrained Control Algorithm for General Batch Pro-cesses.International Journal of Control72,No14,1271-1287[7]Alamir,M.&Marchand,N.(1999)Numerical Stabilization of Nonlinear Systems-Exact Theoryand Approximate Numerical Implementation.European Journal of Control5,87-97[8]Alamir,M.(2001)Nonlinear Receding Horizon sub-optimal Guidance Law for minimum intercep-tion time problem.Control Engineering Practice.Volume9,Issue1,107-116.[9]Alamir,M.(2001)Solutions of Nonlinear Optimal and Robust Control Problems via a MixedCollocation/DAE’s Based Algorithm Automatica.Volume37,Issue7,pp.1109-1115[10]Alamir,M.,L.Calvillo-Corona(2002)Further results on nonlinear receding-horizon observersIEEE Trans.on Automatic Control,Volume47,No7,pp.1188-1193[11]Alamir,M.(2002)Sensitivity analysis in simltaneous state/parameter estimation for inductionmotors International Journal of Control.Vol75,No10,pp.753-758[12]Alamir,M.&Balloul,I.&Marchand,N.(2002)Robust stabilization of nonlinear systems bydiscontinuous dynamic state feedback International Journal Control.Vol75,No6,pp 421-433[13]Alamir,M.and N.Marchand(2003)Constrained minimum-time-oriented feedback control for thestabilization of nonholonomic systems in chained form Journal of Optimization theory and applications.Vol118,No2,pp229-244[14]Alamir,M.(2004)New path generation based receding horizon formulation for constrained stabil-isation of nonlinear systems Automatica.Volume40,Issue4,pages647-652(2004)1[15]Marchand,N.&Alamir,M.(2003)Discontinuous Exponential Stabilisation of Chained FormSystems.Automatica.Vol39,pp343-348[16]Marchand,N.&Alamir,M.(2003)Numerical stabilization of a rigid spacecraft with two actuators.Journal of Dynamic Systems,Measurement,and Control.vol.125,no.3,pp.489-491.[17]Alamir,M.&Corriou,J.P.(2003)Nonlinear Receding-horizon state estimation for dispersiveadsorption columns with nonlinear isotherm.Journal of Process Control.Vol13,Number6, pp517-523[18]Alamir,M.&F.Boyer(2003)Fast generation of attractive trajectories for an under-actuatedsatellite:Application to feedback control design.Journal of Optimization&Engineering Vol4.pp215-230[19]Attia,S.A.,Alamir,M.&C.Canudas de Wit(2006)Voltage Collapse Avoidance in PowerSystems:A Receding Horizon Approach International Journal on intelligent automation and soft computing.Vol12.Special issue on”intelligent automation in power systems. [20]Chemori,A.&Alamir,M.(2005)Limit cycle generation for a class of nonlinear systems withjumps using a low dimensional predictive control International Journal of Control,Vol78, Issue15,pp1206-1217[21]Chemori,A.&Alamir,M.(2006)Multi-step limit cycle generation for Rabbit’s walking based on anonlinear low dimensional predictive control scheme International Journal of Mechatronics, Vol16/5pp259-277(2006)[22]Alamir,M.and Ibrahim,F.and Corriou,J.P.(2005)Aflexible nonlinear model predictive con-trol scheme for quality/performance handling in nonlinear SMB chromatography Journal of process control Volume16,Number4,pp.333-344[23]Corriou,J.P.and Alamir,M.(2005)A hybrid nonlinear state observer for concentration pro-files reconstruction in nonlinear simulated moving bed Journal of process control Volume16, Number4,pp.345-353[24]Boyer,F.and Alamir,M.(2006)Further results on the accessibility of a satellite with two reactionwheels Journal of Guidance,Control,and Dynamics.Vol30(2)2007[25]Alamir,M.(2006)A benchmark for optimal control problem’s solvers for hybrid nonlinear systemsAutomatica,Volume42,Issue9,pp1593-1598[26]Alamir,M.;Sheibat-Othman,N.and Othman,S.(2007)Constrained Nonlinear Predictive Con-trol for maximizing Production in Polymerization Processes.IEEE Transactions on Control Systems Technology,Vol15,No2,March,2007.[27]Commault,C.and Alamir,M.(2007)On the reachability in arbitrary time for positive continuouslinear systems Systems&Control Letters.Vol56,Issue4,pp272-276(2007).[28]Wittrant,E.and Canudas-de-wit,C.and Georges,D.and Alamir,M.(2007)Remote stabilizationvia communication network with a distributed control law IEEE Transactions on Automatic Control(To appear)[29]Boyer,F.Alamir,M.Chablat,D.Khalil,W.Leroyer,A.and Lemoine,Ph.(2006)Robot Anguillesous marin en3D Revue”techniques de l’ing’nieur”,Vol.S7856,pp1-16(2006). [30]Alamir,M.and Chareyron,S.(2007)State Constrained Optimal Control Applied to Cell-Cycle-Specific cancer Chemotherapy”Optimal Control Appliactions and Methods”,Vol28,Is-sue3,pp175-190(2007).[31]Alamir,M.and Murilo,A.(2007)Swing-up and stabilization of a Twin-Pendulum under stateand control constraints by fast NMPC scheme”To appear in Automatica”(2007).[32]Alamir,M.and Sheibat-Othman,N.and Othman,S.(2007)On the use of Nonlinear RecedingHorizon Observers in Batch Terpolymerization Reactors.”To appear in the International Journal of Modeling,Identification and Control”(2007).[33]Alamir,M.&Bornard,G.(1994)Sur la stabilit´e de la commande optimale`a horizon fuyant dessyst`e mes discrets non lin´e aires.C.R.Acad.Sci.,s’rie I,318,pp769-773Papiers invit’s et papiers de discussion(revues internationales)[34]Alamir,M.&S.A.Attia(2003)Discussion on the paper”An Optimal Control Approach forHybrid Systems.”European Journal of Control.Vol.9,pp459-460Brevets[35]Thomas,J.L,tteux,J.C.Alacoque,G.Bornard&M.Alamir(1998)Proc´e d´e d’estimation dela vitesse d’un v´e hicule ou d’une association de v´e hicules.Brevet´e tendu sur l’Europe et l’Extrˆe me Orient98402960.3-2213.Alsthom Transport.[36]Alamir,M&Alacoque,J.C(2000)Proc´e d´e et dispositif de localisation d’un v´e hicule sur une voie.Brevet Europe,USA,Canada,Chine.ALSTOM.AT275,A30253Livres&chapitres dans des livres[37]Alamir,M.(2000)Estimateurs d’´e tat.horizon glissant.Dans:Automatique et proc´e d´e s chim-iques.Ed Herm`e s.ISBN2-7462-0207-7.Coordinateur:J.P.Corriou[38]Marchand,N.,M.Alamir&I.Balloul(2000)Stabilization of nonlinear systems by dynamic statefeedback.Dans:NCN Springer-Verlag-Series.Ed Springer-Veralg[39]Marchand,N.&M.Alamir(2000)Asymptotic controllability implies stabilizability:a more generalresult.Dans:NCN Springer-Verlag-Series.[40]Balloul,I.&M.Alamir(2000)Optimal Control with harmonic reduction and induction machines.Dans:NCN Springer-Verlag-Series.[41]M.Alamir&S.Chareyron(2006)State Constrained Optimal Control Applied to Cell-Cycle-Specific Cancer Chemotherapy.Dans:Positive Systems Theory and Applications.Lecture Notes in Control and Information Science.Springer-Verlag2006.[42]M.Alamir(2006)A low dimensional contractive NMPC scheme for nonlinear systems stabi-lization:Theoretical framework and numerical investigation on relatively fast systems.In:R.Findeisen,R.Allgower,L.Biegler(Eds.),Assessment and Future Directions of Nonlinear Model Predictive Control,Lecture Notes in Control and Information Sciences,Springer.[43]M.Alamir&Boyer,F.(2006)Re-injecting the structure in NMPC schemes-Application tothe constrained stabilization of a snakeboard.In:Fast Motions in BioMechanics and Robotics, Optimization and Feedback Control,Lecture Notes in Control and Information Sciences,Springer.[44]Alamir,M.(2006)La commande Pr’dictive Non lin’aire.Dans:La Commande Pr’dictive.EdHermSs.(Ouvrage collectif).Coordinateur:D.Dumur[45]Witrant,E.and C.Canudas de Wit and F.Georges and and M.Alamir(2006)On the Use ofState Predictors in networked Control Systems.In:LNCIS Springer-Verlag-Series.[46]Alamir,M.(2006)Stabilization of Nonlinear Systems Using Receding-Horizon Control Schemes:A Parametrized Approach for Fast Systems.Lecture Notes in Control and Information Sciences,Springer,London,ISBN1-84628-470-8(2006)[47]Alamir,M.(2007)Nonlinear Moving Horizon Observers:Theory&Real-Time Implementation.InNonlinear Observers and Applications,Gildas Besan¸c on(Ed).Lecture Notes on Communication and Information Science.Springer-Verlag-Series.Rapports de contrats industriels ou Europ´e ens(LAG)[48]Alamir,M.(1995)A Friction Based Transient Model for Slug Flow in Two-Phase Pipelines.Rapport de Recherche.RSI,Montbonnot,France[49]Alamir,M.(1995)Simulation of Two-Phase Flow in Pipelines.Rapport de Recherche.RSI,Montbonnot,France[50]Alamir,M.(1997)Sur l’estimation de la vitesse d’un train en pr´e sence de glissement au niveaudu contact roue/rail.Rapport de G/ALCATEL-ALSTHOM-RECHERCHE [51]Guey,V.&M.Alamir(1997)Mod`e le de simulation en vue de la commande des trains pendulaires.Rapport de Contrat-Phase G/GEC-ALSTHOM9LQ4[52]Alamir,M.(1997)Sur la pr´e diction des courbes pour trains pendulaires.Rapport de Contrat-Phase G/GEC-ALSTHOM9LQ4[53]Alamir,M.,V.Guey&G.Bornard(1998)Sur la commande des trains pendulaires.Rapport deContrat-Phase G/GEC-ALSTHOM9LQ4[54]Alamir,M.(1998)Un observateur non lin´e aire pour la localisation d’un train sur une carte demesures.Rapport de Contrat-Phase G/GEC-ALSTHOM9LQ4[55]M.Alamir&S.A.Attia&C.Canudas de Wit(2003)Nonlinear Predictive Controller for thesimplified ABB Test Power System Stabilization Problem.Rapport de Contrat Europ´e en CC-Project(LAG-VERIMAG)[56]M.Alamir&S.A.Attia(2003)On solving optimal control problems for switched hybrid nonlinearsystems by strong variations algorithms.Rapport de Contrat Europ´e en CC-Project(LAG-VERIMAG)[57]Besan¸c on,G.,M.Alamir&G.Bornard(1998)Commande des machines asynhcronesavec stabilisation dufiltre de tˆe te.Rapport de G/ALCATEL-ALSTHOM-RECHERCHE[58]Bornard,G.and Alamir M.(2007)Am´e lioration et validation d’un simulateur d’un g´e n´e rateurde vapeur:Simulation de la loi de commande Calor.Rapport de G/CALOR-Tˆa che1[59]Alamir M.&Guy Bornard(2007)Estimation du niveau d’eau`a partir des mesures de temp´e ratureRapport de G/CALOR-Tˆa che2.a[60]Guy Bornard&Alamir M.(2007)Am´e lioration des performances d’un g´e n´e rateur de vapeur de fer`a repasser.Conception de lois de commande.Elaboration d’une loi et´e tude comparative.Rapport de G/CALOR-Tˆa che2.b[61]Alamir M.&G.Bornard(2007)Analyse en aveugle d’un g´e n´e rateur de vapeur Rapport deContrat.Gipsa-lab/CALOR-Rapport Phase1:Test de faisabilit´e sur donn´e es disponibles.[62]Alamir M.&G.Bornard(2007)Analyse en aveugle d’un g´e n´e rateur de vapeur Rapport deContrat.Gipsa-lab/CALOR-Rapport Phase2:Test sur les donn´e es exp´e rimentales.Rapports dans le cadre du conseil scientifique aux industriels[63]Alamir,M.(2001)On simultaneous estimation of vehicle’s load and road inclination Rapportde consulting.Sherpa Engineering,Nanterre,FranceArticles dans des conf´e rences internationales avec actes[64]Alamir,M.&G.Bornard(1994)New sufficient conditions for global stability of receding horizoncontrol for discrete-time nonlinear systems.in Advances in model predictive control,Oxford Science Publications,173-181[65]Alamir,M.&H.Khennouf(1995)Discontinuous receding horizon control-based stabilizing feedbackfor nonholonomic systems in power form.Proceedings of the IEEE Conference on Decision and Control,New Orleans[66]Alamir,M.&G.Bornard(1995)Minimum energy-based discontinuous globally stabilizing strategyfor a class of nonlinear systems.Proceedings of the European Control Conference,Rome, Italy[67]Etay,J.M.Alamir&Y.Fauterelle(1993)Dynamics of electromagnetically driven spheres inliquid,comparaison between numerical simulations and experiments.Proceedings of the7th Beer-Sheva seminar on MHDflows and turbulence,Israel.[68]Marchand,N.&M.Alamir(1998)From Open Loop Trajectories to Stabilizing State Feedback-Application to a CSTR-.Proceedings of the IF AC Conference on Systems Structures and Control,Nantes,France.[69]Alamir,M.&G.Bornard(1995)Optimization based stabilizing strategy for nonlinear discretetime systems with unmatched uncertainties.Proceedings of the Second International Sym-posium on Methods and Models in Automation and Robotics,Miedzyzdroje,Poland.193-198[70]Marchand,N.&M.Alamir(1998)Numerical Stabilization of a Rigid Spacecraft with Two Ac-tuators.Proceedings of the IF AC Conference on Systems Structures and Control, Nantes,France.[71]Benthabet,M.M.Alamir,M.Bailly&J.P.Corriou(1997)Nonlinear control of a simulatedmoving bed.Proceedings of the AIChE’s1997Annual Meeting,Los Angeles,USA.16-21 [72]Alamir,M.(1999)Optimization Based Nonlinear Observers Revisited.Proceedings of theIF AC World Congress,Beijin,China.[73]Alamir,M.(1999)On the Feasibility of Finite Difference Homotopy Based Algorithm for NonlinearOptimal Control Problems.Proceedings of the IF AC World Congress,Beijin,China. [74]Alamir,M.&I.Balloul(1999)Robust Constrained Control Algorithm for General Batch Processes.Proceedings of the European Control Conference,Karlsruhe,Germany[75]Balloul,I.&M.Alamir(1999)Analysis Of H Infinity Feedback Systems Using The Gap Metric.Proceedings of the12th International Conference on Control Systems and Computer Science Bucharest,Romania[76]Alamir,M.(2000)Nonlinear Robust Controller for Rotating Stall and Surge in Axial Flow Com-pressors Proceedings of IEEE39th Conference on decision and Control Sydney,Australia [77]Alamir,M.(2000)Solution of nonlinear optimal and robust control problems via a mixed col-location/DAE’s based algorithm Proceedings of IEEE39th Conference on decision and Control Sydney,Australia[78]Alamir,M.(2001)New path generation based receding horizon formulation for constrained stabili-sation of nonlinear systems Proceedings of IEEE40th Conference on decision and Control Orlando,Florida[79]Alamir,N;L.Calvillo(2001)Further results on nonlinear receding horizon observers Proceedingsof IEEE40th Conference on decision and Control Orlando,Florida[80]Alamir,M.;H.Khennouf(2001)A Matlab/Gui c Case-Study Environment for Nonlinear Con-trol Learning Proceedings of IEEE40th Conference on decision and Control Orlando, Florida[81]L.Calvillo-Coronna,M.Alamir(2002)A Towards the construction of a robust optimal observerIF AC World Congress,Barcelona,Spain.[82]Alamir,M.;F.Boyer(2001)Fast generation of attractive trajectories for a deficient satellite.Application to feedback control design Proceedings of the IF AC Symposium on s Systems Structures and Control.Prague.[83]Alamir,M.(2002)On Friction Compensation without Friction Model.Proceedings of the IF ACWorld Congress.Bareclona,Spain.[84]Alamir,M.(2002)Sensitivity Analysis in Simultaneous State/Parameter Estimation for InductionMotors.Proceedings of the IF AC World Congress.Bareclona,Spain.[85]Alamir,M.;J.P.Corriou(2002)Nonlinear Receding-Horizon State Estimation for Dispersive Ad-sorption Columns with Nonlinear Isotherm Proceedings of the41st Conference on Decision and Control,Las Vegas,USA.[86]Alamir,M.;N.Marchand(2002)Constrained sub-optimal feedback control for minimum-timestabilization of nonholonomic systems in chained form Proceedings of the41st Conference on Decision and Control,Las Vegas,USA.[87]N.Marchand&M.Alamir(2002)Discontinuous Exponential Stabilization of Chained Form Sys-tems Proceedings of the41st Conference on Decision and Control,Las Vegas,USA.[88]B.Youssef&M.Alamir(2003)Generic signature based tool for diagnosis and parametric esti-mation for multi-variable dynamical nonlinear systems.Proceedings of the42th Conference on Decision and Control,Hawa¡,USA.[89]Alamir,M.;S.A.Attia and Carlos Canudas-de-Wit(2004)Nonlinear Predictive Controller for aTest Power Systems Stabilization Problem.(Invited session)Proceedings of World Automa-tion Congress,June2004,Spain.[90]Chemori,A;M.Alamir(2004)Nonlinear Predictive Control of Under-actuated Mechanical Sys-tems Application:the ECP505inverted pendulum Proceedings of the16th International Symposium on Mathematica Theory of Networks and Systems,Leuven,Belgium. [91]Chemori,A;M.Alamir(2004)Low dimensional predictive control scheme for limit cycle gen-eration in nonlinear hybrid control system CCCT’04–2nd International Conference on Computing,Communication and Control Technologies–Austin,Texas,USA[92]Chemori,A;M.Alamir(2004)Generation of multi-steps limit cycles for Rabbit using a lowdimensional(scalar)nonlinear predictive control scheme IROS2004IEEE/RSJ international Conference on Intelligent Robots and Systems,Sendai,Japan[93]Witrant,E.;C.Canudas de Wit;D.Georges and M.Alamir(2004)Remote Stabilization viaTime-Varying Communication Network Delays:Application to TCP networks Proceedings of the2004IEEE Conference on Control Applications,Taipei,Taiwan.(Invited paper)[94]Alamir,M.;S.A.Attia(2004)On solving optimal control problems for switched nonlinear sys-tems by strong variations algorithms NOLCOS2004Symposium on Nonlinear Control Systems.Stuttgart,Germany.[95]Alamir,M.; mault(2004)A New Equal-Partition Measurement Encoding Scheme forNetworked Control Systems NOLCOS2004Symposium on Nonlinear Control Systems.Stuttgart,Germany.[96]Alamir,M.;F.Ibrahim and J.P.Corriou(2005)Aflexible nonlinear MPC Scheme for Qual-ity/Performance Handling in Nonlinear SMB Chromatography Proceedings of the IF AC World Congress.Praha,Czech Republic.[97]Dejardin,N.Alamir,M.and J.P.Corriou(2005)Model Predictive Control of a Simulated MovingBed Proceedings of the IF AC World Congress.Praha,Czech Republic.[98]Corriou,J.P.;M.Alamir(2005)A nonlinear observer for concentration profiles in simulatedmoving bed Proceedings of the IF AC World Congress.Praha,Czech Republic.[99]Youssef,B.;M.Alamir(2005)Diagnosis and on-line parametric estimation of automotive elec-tronic throttle control system Proceedings of the IF AC World Congress.Praha,Czech Republic.[100]Youssef,B.;M.Alamir and F.Ibrahim(2005)Diagnosis and on-line parametric estimation of Simulated Moving Bed Proceedings of the44th IEEE Conference on Decision and Con-trol.[101]Attia,S.A.;M.Alamir and C.Canudas de Wit(2005)Suboptimal Control of switched nonlinear systems under location and switching constraints Proceedings of the IF AC World Congress.Praha,Czech Republic.[102]Houdebine,M.;S.Dedieu;M.Alamir and O.Sename(2005)A new fractional frequency syn-thesizer architecture with stability and robustness analysis Proceedings of the IF AC World Congress.Praha,Czech Republic.[103]Houdebine,M.;Dedieu,S.;Sename,O.and Alamir,M.(2006)A sampled spur free fractional frequency synthesizer and its noise analysis Proceedings of ESSCIRC2006,32nd European Solid-State Circuits Conference,18-22September2006,Montreux,Switzerland[104]Alamir,M.(2005)A low dimensional contractive NMPC scheme for nonlinear systems stabi-lization:Theoretical framework and numerical investigation on relatively fast systems Proceed-ings of the international workshop on assessment and future directions of NMPC, Freudenstadt,Germany,August2005(Invited paper)[105]Alamir,M.and Boyer,F.(2005)Re-injecting the strucrture in NMPC schemes:Application to the constrained stabilization of a snake board Proceedings of the international Symposium on Fast Motion in Biomechanics and Robotics,September2005(Invited paper) [106]Hably,A.;Marchand,N.and Alamir,M.(2005)Constrained Minimum-time-oriented Stabilization of Extended Chained Form Systems Proceedings of44th IEEE Conference on Decision and Control and European Control Conference(CDC-ECC’05),Seville,Spain,2005.[107]Attia,S.A.and Alamir,M.(2005)On Complexity Reduction of Voltage Stabilization MPC Schemes by Partial Explicit Feedbacks Proceedings of44th IEEE Conference on Deci-sion and Control and European Control Conference(CDC-ECC’05),Seville,Spain, 2005.(Invited Session)[108]Alamir,M.;Sheibat-Othman,N.and Othman,S.(2006)Constrained nonlinear receding-horizon control for maximizing production in polymerization processes Proceedings of thefirst IF AC workshop on nonlinear model predictive control for fast systems,NMPC-FS’06, Grenoble,France.[109]Attia,S.A.and Alamir,M.(2006)A predictive switching strategy for a class of hybrid systems: wave suppression influid dynamic systems Proceedings of thefirst IF AC workshop on nonlinear model predictive control for fast systems,NMPC-FS’06,Grenoble,France.[110]Alamir,M.;Hafidi,G.;Marchand,N.;Elrafei,M.;Porez,M.and Boyer,F.(2007)Feedback design for3D movement of an Eel-like robot Proceedings of the2007IEEE International Conference on Robotics and Automation,Roma,Italy.[111]Alamir,M.(2007)On state-dependant sampling for nonlinear controlled systems sharing limited computational resources Proceedings of the2007Europen Control Conference,Greece, 2007.[112]Houdebine,M.;Alamir,M.;Sename,O.and Dedieu,S.(2007)Nonlinear Observer for Frequency Meter in Digital Phase Locked Loops Proceedings of the2007Europen Control Confer-ence,Greece,2007.[113]Alamir,M.;Sheinat-Othman N.and Othman S.(2007)On The Use of Nonlinear Receding-Horizon Observers in Batch Terpolymerization Processes with Partially Unmodelled Dynamics: Formulation and Experimental Validation.Proceedings of the7th IF AC Symposium on Nonlinear Control Systems(NOLCOS2007),Pretoria,South Africa,August2007 (Invited session on nonlinear observers in bioreactors).[114]Murilo,A.P.and Alamir,M.(2007)Output Feedback Design of a Twin Pendulum System in Presence of Sensor Bias Proceedings of the7th IF AC Symposium on Nonlinear Control Systems(NOLCOS2007),Pretoria,South Africa,August(2007)[115]Amari R.,Alamir,M.and Tona P.(2008)Unified MPC strategy for idle-speed control,vehicle start-up and gearing applied to an Automated Manual Transmission Proceedings of the IF AC World Congress,South Korea,(2008)[116]El-Rafei M.,Alamir,M.,Marchand N.,Porez M.and Boyer F.(2008)Motion Control of a Three-Dimensional Eel-like Robot Without Pectoral Fins Proceedings of the IF AC World Congress,South Korea,(2008)[117]Alamir,M.(2008)On Useful Redundancy in Dynamic Inverse Problems Related Optimization Proceedings of the IF AC World Congress,South Korea,(2008)Th`e ses[118]Alamir,M.(1995)Contributions`a l’´e tude de commande optimale.horizon fuyant des systSmes discrets non lin´e aires ThSse de doctorat,Laboratoire d’Automatique de Grenoble,Institut National Polytechnique de Grenoble.[119]Alamir,M.(2000)Promenades non lin´e aires Th`e se d’habilitation`a diriger les recherches (HDR),Laboratoire d’Automatique de Grenoble,Institut National Polytechnique de Grenoble.。

CURRENTGRADUATESTUDENTSFORMERGRADUATE…

George NagyCURRENT GRADUATE STUDENTSFORMER GRADUATE STUDENTSABDALI, Kamal NSF, Stanford U MS 1968 U. Montreal ABU TARIF, Asad GE Medical Image Syracuse MS May 1998ABU TARIF, Asad GE Medical Image Syracuse PhD Dec 2002 advisor ADAMS, Marseta (now M. Dill) FAA ME Dec 2003AHMED, Zubair X6032 274 8236 PhD May 1991 committee ALI, Michael (Prof. Stephanou) PhD July 1999 committee AL KHOFAHI, Khalid (Badri Roysam) Thompson R&D, MN PhD 2001 comittee ALLEN, David Lincoln Labs MS May 1986 ANAGNOSTOPOULOS, Tasso MS May 1996ANDRA, Srinivas Soros NYC PhD 2006 advisor ANSON, Ed Tulip MS 1975 UNL BATTACHARYYA, Anoop Epson Research PhD Dec 1994 committee Belik, David (Prof. Nelson) PhD May 1983 committee BARNEY SMITH, Elisa Boise State U.PhD Dec 1998 advisor BHASKAR, Kasi MS 1977 UNLBOHN, Jan (Prof. Wozny) VPI PhD Aug 1993 committee CARPINELLI, John Stevens Institute? PhD Aug 1987 committee CHANDRASEKHAR, Narayanaswami (Pf. Franklin)IBM Watson PhD Dec 1990 committee CHEN, Ying (B. Roysam) PhD Dec 2009 committee CHRISTENSEN, Neil PhD May 1986 committee CHUGANI, Mahesh (Prof. Savic) National Instruments PhD May 1996 committee COGDELL, David U. Miss. PhD Dec 1987 committee DEFFENBAUGH, Grant CMU PhD Dec 2003 advisorEL-NASAN, Adnan U. Yarmouk PhD Aug 2003 advisor DOUGLASS, Barry UT Dallas PhD Dec 1989 committee FALSAFI, Aram Digital MS May 1989GATTANI, Abhishek Stryker Endoscopy MS August 2005GODA, Brian (Jack McDonald) US Army PhD May 2001 committee GUO, Hwei Chi PhD May 1997 committee GREEN, Ed (Prof. Moorthy) Union College PhD May 1996 committee GUERRIERI, Ernesto (Prof. Freeman?) PhD May 1989 committee GUR-ALI, Ozden (Dec’n Science) GE CRD PhD Aug 1994 committee Hallquist, Roy (Prof. Nelson, UNL) PhD May 1973 committee HARMSEN Jeremiah (Prof. Pearlman) PhD Aug 2005 committee Jeong Yongwon (Prof. Radke) PhD Dec 2006 committee HILAIRE, Thierry (ME) Renault R&D PhD Dec 1993 committee HIRAOGLU, Muzaffer Analogic Corp Peabody ? PhD May 1992 advisor HIJAZI, Nabil (Prof. Savic) PhD May 2000 committee HOMBAL, vadi (Prof. Sanderosn) PhD Dec 2009 committee HONG, Shen (Badri Roysam) Siemens NJ PhD 2000 committee INANC, Metin (WRF) Sync NJ PhD 2008 Committee Jandhyala, Chakradhar, Ramana Blomberg NYC MS mAY 2010 advisor JEANNE, Philippe France MS Aug 1990Jha, Piyushee BEA Endicott MS May 2008 advisor JIANG, Ming (Prof. Ji) PhD Dec 2005 committee, JOGLEKAR, Indrajeet MS 2008 advisor JORGE, Joachim (Prof. Glinert) U.Lisbon PhD May 1995 committee JOSHI, Ashutosh Fair Isaac San Francisco PhD May 2005 advisor JOSHI, Raviv (Prof. Sanderson) PhD Dec 1996 committee JUNG, Dz-Mou Yahoo San Diego MS Dec 1989" PhD May 1995 advisor KANAI, Junichi RPI PhD Aug 1990 advisor KALAFALA, Kerim IBM Fishkill MS Dec 1997KANG, Hang Bong Samsung PhD Dec 1993 committee KANJILAL, Shuvayu Oracle MS May 1997 KANKAHALLI, Mohan (Prof. Franklin) Nat'l U Singapore PhD Dec 1990 committee KIM, Jaesok (Prof. E. Walker) Bell Labs PhD 1988 committee KIM, Beong-Jo (Prof. Pearlman) Korean Army PhD Dec 1997 committee KLINE, Jaclyn Moorny NYS Labor Dept MS 2009KODEIH, Mohamad SUNY, New Paltz PhD 1988 committee Kucharen, Promote (Prof. Modestino) ME May 1996KURADA, Lakshmi Vijaya MS Aug 1996KWAK, William Digital PhD Dec 1988 committee Leibovich, Ilya ME May 1996LIAO, Wenhui Thomson Reuters PhD 2006?? committee??LIMNER, Joel (ME) GMR PhD Aug 1997 committeeLU, Zhitao (Bill Pearlman) PhD Dec 2000 committeeLu, Renzhi (Radke) PhD` 2007 committeeLui, Roy MS 2010 advisorLYON, Doug Fairfield U. PhD Dec 1991 advisor MACULOTTI, Marina U. Genova Laurea July 1988MALLOCH, Chuck (Prof. Gerhardt) Pratt&whitney – Atrey PhD Dec 1989 committee MANJUNATH, D (Prof. Pearlman) PhD July 1993 committee MARTIN, Kenneth Morris (Prof. Wozny) GE PhD May 1998 committee MARTINO, Peter (C. Stewart?) Digital PhD May 1988 committee MAULIK, Amitava Connectiva, Kolkota PhD May 1992 advisorMEHTA, Shashank IIT Kanpur PhD Aug 1985 advisorMILLER, Anne AOL NYC (Roger Grice) MS Dec 2008 co-advisor MILLER, Jimmy (Prof. Stewart) GE CRD PhD Aug 1997 committee MITHAL, Sam Digital MS May 1989MUKHERJEE, Maharaj IBM Fishkill PhD 1992 advisorNAIR, Hari X2896 I MS May 1988NARENDRA, N.C. PhD May 1991 advisor NICEWARNER, Keith PhD Feb 1995 committee PADFIELD, Dirk GE PhD May 2009 committee PADMANABHAN, Raghav RPI PhD Program MS 2009PERARA, Amitha (Prof. Stewart) GE Research PhD May 2003,committee PEDRINO, Helio (Prof. Franklin) PhD May 2000 committee PHILHOWER, Robert (Prof. McDonald) IBM Watson PhD 1993 committee PRUETT, David off:914-385-6190 IBM MS Dec 1988RAJAIDEHKORDY, Barry MS 1983 UNLRAY, Clark (Prof. Franklin) West Point MC PhD May 1994 committee ROUSSEL, Nicolas PhD 2007 committeeSALLA, Trevor, Media Tech, Philly MS May 1995SATHYAMURTHY, Radhika Dupont MS 1992SARKAR, Prateek PARC, Palo Alto MS Dec 1994SARKAR, Prateek Palo Alto Research Ctr (PARC)PhD May 2000 advisor SHAPIRA, Andrew (Prof. Moorthy) OneXero Seattle MS 1998?SHAPIRA, Andrew (Prof. Moorthy) OneXero Seattle PhD Dec 2000 committeeSHEN, Weicheng U. New Hamp. PhD Dec 1987 committee SORENSEN, Jeff (Prof. Savic) IBM Watson PhD 1993 committee? SRIDHARAN, K. IIT Madras) PhD May 1996 committee SIVASUBRAMANIAM, Suthaharan Oracle ME Dec 1998SUMMERS, Jason PostDoc, NEC Japan PhD Aug 2003SWANN, Jonquil METHOMAS, Mathews Digital MS May 1988TONG, Yan (Prof. Q. Ji) GE Research PhD Dec 2007 Committee TSERKEZOU, Polly Zurich MS May 1988VISWANATHAN, Mahesh IBM Watson Yorktown NY PhD Dec 1991 advisor VEERAMACHANENI, Harsha Thompson R&D Eagan MN PhD Dec 2002 advisorVIZCAYA, Jose(Prof. Gerhardt)U. Javeriana Colombia PhD May 1998 committee VOUGIKIAS, Stavros PhD Dec 1995 committee WACLAWIK, Jim Boston MS Dec 1991WAGLE, Sharad retired PhD May 1978 advisorWANG, Xiaoyin Qualcom MS Dec 1995WANG, Peng (Prof. Ji) U. Penn PhD Dec 2005 committeeXU, Yihong EMC PhD Aug 1998 advisor Yanamadala, Bhavani Shankar Atlanta MS 2007YU, Jim Bell Labs MS Dec 1986ZHANG, Tong Brontes Tech, Woburn MA PhD May 2004 advisorZ NLM, NIH Bathesda MD ZOU, Jie NLM, NIH Bathesda MD PhD May 2004 advisorEXTERNAL READER OR EXAMNER FOR:SKUCE, D.R. McGill University MS Jun 1971DYDYK, R.B. U. Waterloo PhD Mar 1972 HUSSAIN, A.B.S. UBC PhD Jun 1972 POULSEN, R.S. McGill University PhD Apr 1973 NAGARAJA, G. IISc Bangalore PhD Apr 1975BANSAL, Veena IIT Kanpur PhD Dec 1997PAL, Umada ISI Calcutta PhD Dec 1997RICE, Stephen UNLV PhD 1996DE JESUS, Edison, U. Campinas PhD 1997PABST, Frederic U. Fribourg PhD Dec 1998WALTER, Fredrik SLU, Uppsala PhD Oct 1999WIMMER, Zsolt ENST, Paris PhD Dec 1998BOLDO, Didier Sorbonne, Paris PhD July 2002 BAGDANOV, Andrew U. Amsterdam PhD June 2002MURALI, S. U. Mysore PhD Nov 2002 OLIVETTI, Emanuele U. Trento PhD May 2008GARDES, Joel U. Lyon PhD 2009LONG, Vanessa Mcquerrie, Australia PhD 2010Zhu, Yaoyao Prof. Huang Lehigh PhD 2012Chen, Jin Prof. Lopresti Lehign PhD 2013 ?RENSSELAER UNDERGRADUATESBARGHAVA, Anuba Cervitor SURP 2010 BAUSEWEIN, Jason X 7928 (518) 863-4811 URP 1990 BERG, Andrew Ballow Camera NYS CATS 2009 CAMPOFOIERE, Kyle RFID, URP for credit 2009 CELENTANO, Kathryn Cervitor RCOS 2009 CHAN, Hing Lun 1992-93? CHENG, Greenie CAVIAR Summer project 2003 CHOW, Man Chit Francis 1992? CLIFFORD, Bryan Cybertrust REU, RCOS 2008-2009 DERBY, Laura CAVIAR URP 2004 DING, Mike Cervitor RCOS URP 2009 FLIZARDO, Daniel Ballots UTP 2011 FELKAMP, Amanda Cervitor URP 2011-12 GAGOSKI, Borjan CAVIAR SUMMER URP 2003 GREEN, Matthew SUMMER 2004 HILDEBRAND, Dan CAVIAR RCOS 2009 HUBER,John independent research 1996 HUNTER Travis TANGO URP 2009 ISLAM, Ashfaqul Calligraphy URP 2011 JIANG, Haimei CAVIAR independent research URP 2001 KELLEY, Sean TANGO RCOS 2009 KIM, Sung Hun URP 1994 KONG, Jialiang, Jason ASR MS CMU Qualcom 09 URP 2005 KORDON, Mark X 8897 304 William URP 1990 KYRIAZIS, George Athens 30 1 2026780 272-6197 Tables Senior project 1990 MCAVOY, Dave Ballot camera RCOS 2009 MCCAUGHRIN,Eric,******************X4599URP 1990 MOSHER, Doug SURP 2010 MURPHY, Luke TANGO URA 2007 MUTALATHU, Max TANGO, Cybertrust REU 2009 NGUYEN, Tram Model car? URP 1997 POLYAK, Ilya URP 1991 POPLAWSKI Seven TANGO SURP 2010 ROBERTS, Sam Ballot Images RCOS 2009 ROTHCHILD, Russ CAI Student monitor Senior project 1989TRANLONG, Luke Model car Senior Project 1991 SAJJAD, Syed Senior Project 1993 SHRIVASTAVA, Vivek X 7800 Senior Project 1991 SILVA, Gregory R-dropping URP 2011 SILVERSMITH, William TANGO REU, NSF 2009 STEVENS, Robert RFID, URP for credit Lutron 2009 SUH, Ria URP 1994-95 SWEIS, Slameh Line wrap Design Project 1994 TAMHANKAR, Mangesh TANGO URP 2011 VERNIKOWSKY, Makim TANGO URP 2011 VULIN, Lillian URP 1996 WARREN, Jeff Cervix URP 2011 WONG, Tyler Chen IBM URP 1997 WONG, Lance 273-8281 URP 1992 YU, Chang, OCR Summer Project 2005 YU, Desong GeoWeb Independent Research URP 2001 ZHANG, Qianyi “Landy” Cervitor SURP 2010 WU, Ziyan Ballots, Grad Student in PicProc 2011?INCOMPLETESLANGER, Jefferey URP 1999 JOHNSON, Kurt MS 1990? ??? LEU, She-Wan 1994SHIRALI, Nagesh 273-9249 PhD May 1990? supervisor Cadence Design 408 987-5221ZHONG, D. 1995STUEBEN Gregg (Moorthy) 2001MANTRI, Anup PhD 2010 advisor。

GaN- Processing, defects, and devices

APPLIED PHYSICS REVIEWSGaN:Processing,defects,and devicesS.J.Pearton a)Department of Materials Science and Engineering,University of Florida,Gainesville,Florida32611J.C.ZolperOffice of Naval Research,Arlington,Virginia22217R.J.ShulSandia National Laboratories,Albuquerque,New Mexico87185F.RenDepartment of Chemical Engineering,University of Florida,Gainesville,Florida32611͑Received16October1998;accepted for publication3March1999͒The role of extended and point defects,and key impurities such as C,O,and H,on the electrical and optical properties of GaN is reviewed.Recent progress in the development of high reliability contacts,thermal processing,dry and wet etching techniques,implantation doping and isolation,and gate insulator technology is detailed.Finally,the performance of GaN-based electronic and photonic devices such asfield effect transistors,UV detectors,laser diodes,and light-emitting diodes is covered,along with the influence of process-induced or grown-in defects and impurities on the device physics.͓S0021-8979͑99͒00613-1͔TABLE OF CONTENTSI.Introduction (1)II.Processing (2)A.Ohmic and Schottky contacts (2)B.Rapid thermal processing (6)C.Gate dielectrics (22)D.Wet etching (24)E.Dry etching (28)F.Implant isolation (39)III.Role of Impurities (40)A.Hydrogen (40)B.Oxygen (54)C.Carbon (56)IV.Devices (57)A.AlGaN/GaN Electronics (57)B.Ultrahigh power switches (64)ser diodes (66)D.Light-emitting diodes (67)E.UV Photodetectors (68)V.Summary (69)References (69)I.INTRODUCTIONCurrent GaN-based device technologies include light-emitting diodes͑LEDs͒,laser diodes,and UV detectors on the photonic side and microwave power and ultrahigh power switches on the electronics side.1The LED technology is by now relatively mature,with lifetimes of blue and green emit-ters apparently determined mostly by light-induced degrada-tion of the polymer package that encapsulates the devices.2 The main trends in this technology appear to be optimization of optical output efficiency and solving the polymer package degradation issue.For the laser diodes one of the main life-time limiters is p-ohmic contact metal migration along dis-locations which short out the GaN p-contact layer by spiking all the way to the n side of the junction.3,4This is exacer-bated by the generally high specific contact resistance(R C) of the p-ohmic contact and the associated heating of this area during device operation.The advent of lower threshold de-vices and dislocation-free GaN overgrowth of SiO2-masked regions has allowed achievement of laser lifetimes over 10000h.5Facet formation on the laser has been achieved by dry etching,cleaving,polishing,and selective/crystallo-graphic growth.In structures grown on Al2O3both contacts must be made on the top of the device and hence dry etching is necessary to expose the n side of the junction.Fabrication of UV detectors is relatively straightforward and the main issue seems to be one of improving material purity and quality.With respect to electronic devices for microwave power applications,the main process improvements needed are in the areas of low R C n-ohmic contacts͑the requirements are more stringent than for photonic devices,with R C р10Ϫ7⍀cmϪ2being desirable͒,stable and reproducible Schottky contacts,and low damage dry etching that main-tains surface stoichiometry.For the proposed high power switches͑capable of25kA with3kV-voltage standoff͒there are a number of possible device structures,including thyris-a͒Electronic mail:spear@mse.ufl.eduJOURNAL OF APPLIED PHYSICS VOLUME86,NUMBER11JULY19991tors and several types of power metal-oxide-semiconductor field-effect transistor ͑MOSFET ͒.A schematic of a lateral GaN MOSFET is shown in Fig.1.In this case,critical tech-nologies include high implant activation efficiency,gate in-sulator,trench etching for capacitor formation,and stable high temperature/high current stable ohmic contacts.Recent progress in the development of dry and wet etch-ing techniques,implant doping and isolation,thermal pro-cessing,gate insulator technology,and high reliability con-tacts is first reviewed.Etch selectivities up to 10for InN over AlN are possible in inductively coupled plasmas using a Cl 2/Ar chemistry,but in general selectivities for each binary nitride relative to each other are low ͑р2͒because of the high ion energies required to initiate etching.Improved n -type ohmic contact resistances are obtained by selective area Si ϩimplantation followed by very high temperature ͑Ͼ1300°C ͒anneals to minimize the thermal budget and AlN encapsulation which prevents GaN surface decomposition.Implant isolation is effective in GaN,AlGaN,and AlInN,but marginal in InGaN.Candidate gate insulators for GaN in-clude AlN,AlON,and Ga ͑Gd ͒O x ,but interface state densi-ties must still be decreased to realize state-of-the-art metal-insulator-semiconductor ͑MIS ͒devices.Many outstanding reviews on GaN materials and devices have appeared previously,6–14so we will focus on processing and the influence of defects and impurities on devices.II.PROCESSINGA.Ohmic and Schottky contacts1.Schottky contactsThere are still large variations in barrier heights reported by different workers for standard metals on GaN.Pt appears to produce the highest consistent values ͑ϳ1.0–1.1eV ͒with Ti producing the lowest ͑0.1–0.6eV ͒.The variability ap-pears to result from the presence of several transport mecha-nisms,and to materials and process factors such as defects present in these films,the effectiveness of surface cleans prior to metal deposition,local stoichiometry variations,and variations in surface roughness which could affect unifor-mity of the results.New work on silicides showspromise.15,16For Schottky contacts Pt appeared to be stable to approximately 400°C for 1h,while PtSi is somewhat more stable ͑500°C,1h ͒,and also has barrier heights of ϳ0.8eV.Recent reviews and studies of Schottky contact proper-ties on GaN have appeared.15–28The measured barrier heights in most cases are a function of the difference be-tween the metal work function and the electron affinity of GaN.Some typical values for barrier height for different metals are 1.1eV for Pt,190.91–1.15eV for Au,20,210.6eV for Ti,22and 0.94eV for Pd.23For Ni there is a fairly large discrepancy in reported values,ranging from 0.6624to 0.99eV.25For deposition onto n -GaN (ϳ1017cm Ϫ3),rectifying behavior was observed for Pt,Ni,Pd,Au,Co,Cu,Ag,ohmic behavior for So,Hf,Zn,Al,and V,while intermediate be-havior ͑slightly rectifying ͒was obtained for Nb,Ti,Cr,W,and Mo.16Schmitz et al.16calculated from their data that the density of surface states on GaN was ϳ1.8ϫ1013cm Ϫ2eV Ϫ1,suggesting the degree of pinning of the barrier height is less than GaAs where the surface state den-sity is roughly an order of magnitude higher.A comparison of barrier height data from various sources is shown in Fig.2.16,29–33In the early days of form-ing rectifying contacts on GaN it was often believed the Fermi level at the surface and at the metal–nitrides interface was unpinned.The data of Fig.2shows that indeed the bar-rier height does vary with metal work function.The strategy is then to use a metal with a large work function on GaN ͑such as Pt ͒to form a Schottky barrier,while a metal with a low work function ͑such as Ti ͒should be selected for ohmic contacts.The influence of the surface cleanliness is obviously most important in determining the quality of the Schottky contact.In situ deposition of Ga,followed by thermal desorption under ultrahigh vacuum conditions is found to produce clean GaN surfaces,34,35while in situ N 2ϩionsputtering can also remove native oxides.36,37Liu and Lau 15,FIG. 1.Schematic of an ultra high breakdown voltage GaN powerMOSFET.pilation of published results for Schottky barrier heights on GaN ͑after Refs.18,28,and 68͒.Liu et al.,38and Mohney and Lau18have reviewed surfacecleaning processes for GaN.A number of different acid so-lutions,including HNO3/HCl,HCl/H2O,and HF/H2O,havebeen examined for removing the native oxide,39and superiorcurrent–voltage characteristics are observed for the resultantrectifying contacts.40As with other III–V compound semi-conductors,HCl and HF can significantly reduce the oxideon GaN,41,42while the bases NH4OH and NaOH can alsodissolve the oxide.43To this point,there has been no cleardemonstration of the effect of the polarity of the epilayer onthe barrier height.Mohney and Lau18have also commented on the fact thatthere can be significant spatial differences in the quality ofSchottky barrier contacts on a n-GaN,with diodes showingideality factors ranging fromϽ1.1toу1.3on the same wa-fer.While thermionic emission is clearly the dominant cur-rent transport mechanism in most diodes,tunneling and gen-eration recombination may also be present.In many casesthe high dislocation density in material used to date are prob-ably responsible for most of the spatial variations.The thermal stability of Schottky contacts on GaN iscritically important for practical device operation.The ther-mal limits of most of the metal/GaN combinations are be-tween300and600°C,specifically300°C for Pd,44400°Cfor Pt,45575°C for Au,46and600°C for Ni.18As mentionedearlier,the silicides of Pt and Ni display greater thermalstability than the pure metals.45These contacts however maynot be in thermodynamic equilibrium with the GaN,leadingto the formation of metal gallides and silicon nitride uponprolonged annealing.There is little information on barrier heights on p-GaNdue to the general difficulty in growing high quality p-typematerial and the low hole mobility.A barrier height of2.38eV was reported for Au on p-GaN.472.n-ohmic contactsThe commonly accepted ohmic contact to n-GaN is Ti/Al,which is generally annealed to produce oxide reductionon the GaN surface.Multilevel Au/Ni/Al/Ti structures ap-pear to give wider process windows,by reducing oxidationof the Ti layer.48R C values ofр10Ϫ5⍀cm2have been pro-duced on heterostructurefield-effect transistor͑HFET͒de-vices using Ti/Al annealed at900°C for20s.15Both W and WSi x on nϩepi-GaN layers(n ϳ1019cmϪ3)produce reasonable contacts(R Cϳ8ϫ10Ϫ5⍀cm2),but extremely stable behavior49—annealing at1000°C led to shallow reacted regions ofр100Å,and in junctionfield-effect transistor structures these contacts can withstand implant activation anneals at1100°C.50Reaction with the GaN is relatively limited,although␤-W2N interfa-cial phases are found after800°C anneals,and this appears to be a barrier to Ga out-diffusion.49By contrast to the n-metal systems,the standard p-ohmiccontact to GaN is Ni–Au,with R C valuesу10Ϫ2⍀cm2. Efforts tofind a superior alternative have proved fruitless to date,51even though strong efforts have been made on multi-component alloyed contacts where one attempts to extract one of the lattice elements,replace it with an acceptor dop-ant,and simultaneously reduce the‘‘balling-up’’of the met-allization during this reaction.The model system for thistype of contact is AuGeNi/n-GaAs.A promising approach isto reduce the band gap through use of p-type InGaN on thetop of the GaN.To date there have been reports of achievingp-doping(ϳ1017cmϪ3)in compositions up toϳ15%In.The III nitrides pose a problem however,in the develop-ment of low resistance ohmic contacts because of their wideband gaps.Most of the work done in the area has been fo-cused on n-type GaN.Au and Al single metal contacts tonϩGaN and nonalloyed Au/Ti and Al/Ti were found to havecontact resistances ofϳ10Ϫ3to10Ϫ4⍀cm2.46,52–58Al-containing contacts perform best when oxidation is mini-mized.Lin et al.59reported the lowest contact resistance tonϩGaN,with Ti/Al contacts after annealing at900°C for30s in a rapid thermal annealer(R Cϭ8ϫ10Ϫ6⍀cm2).They suggested the formation of a TiN interface as important in the formation of the low resistance contact.Most of the tran-sition metal elements,including Ti,V,and Sc react with GaN to form nitrides,gallides,and metal–Ga–N ternary phases.18Thermodynamic calculation indicate that the met-als themselves are not in equilibrium with GaN under normal processing conditions,with the consequent probability of in-terfacial reactions occurring.60Both Ti and TiN have been shown to produce ohmic contacts on n-GaN,61,62with Ti consuming GaN during reaction͑a few hundred angstroms at ϳ1000°C for30s͒to form TiN.63Three phase equilibria for the Ti–Ga–N͑and V–Ga–N and Cr–Ga–N͒systems at 800°C have been reported by several groups.64–66The an-nealing ambient plays a strong role,since if the N2partial pressure is greater than that in equilibrium with a metal/GaN contact during reaction,there is a driving force to incorpo-rate nitrogen from the gas phase.18Moreover,there is quite different behavior observed for deposited TiN/GaN contacts relative to those formed by reaction of Ti with GaN in a N2 ambient.In the former case,thermionic emission appears to be the dominant conduction mechanism whereas in the latter tunneling seems to be most important.18Modification of the GaN surface by high temperatureannealing67or reactive ion etching48,68to produce preferen-tial loss of N2can improve n-type ohmic contact resistanceby increasing electron concentration in the near-surface re-gion.Many other metals can be employed to form bilayer Al/metal/GaN n-ohmic contacts,including Pd,69,70Ta,71Nd,72Sc,18and Hf.18All of these form good ohmic contacts,withspecific contact resistances in the10Ϫ5⍀cm2range.A particularly attractive method for reducing R C on de-vice structures is self-aligned implantation of Siϩto heavilydope source/drain ohmic contact regions.This approach hasbeen employed to achieve high quality contacts on hetero-structurefield-effect transistor͑HFET͒structures.73W wasfound to produce low resistance ohmic contacts to nϩGaN(R Cϭ8ϫ10Ϫ6⍀cm2)with little interaction between the semiconductor and the metal up to800°C.49WSi x on nϩGaN was found to be stable to800°C as well,with a contact resistance ofϳ10Ϫ5⍀cm2.Graded contact layers to GaN have been formed with both InN74and InGaN using WSi x metallization.Nonalloyed Ti/Pt/Au on InN produced specific contact resistance R Cϭ1.8ϫ10Ϫ7⍀cm2.74GradedIn x Ga1Ϫx As/InN contacts have been employed on GaAs/ AlGaAs heterojunction bipolar transistors,with R C as low as 5ϫ10Ϫ7⍀cm2.75For high temperature electronics applications,or for high reliability,it would be preferable to employ refractory metal contacts such as W and WSi x.Moreover,the contact resis-tance could be reduced if lower band gap In-containing al-loys͑or InN͒were used as contact layers on GaN,much as in the case of InGaAs on GaAs.However,the In-based nitrides are less thermally stable than GaN,and we need to establish the trade off between contact resistance and poorer tempera-ture stability.Recent experiments on formation of W,WSi0.44and Ti/Al contacts deposited on nϩIn0.65Ga0.35N(n ϳ1020cmϪ3),nϩInN(nϳ1020cmϪ3),and nϪIn0.75Al0.25N (nϳ1018cmϪ3)have been reported.74,75The electrical, structural,and chemical stability of these contacts were ex-amined after anneals up to900°C.It was found that InGaN allows achievement of excellent contact resistance (Ͻ10Ϫ6⍀cm2),with stability up toϳ600°C for W metal-lization.The2000-Å-thick InGaN,InN,and InAlN samples weregrown using metal organic molecular beam epitaxy ͑MOMBE͒on semi-insulating,͑100͒GaAs substrates in an Intevac Gen II system as described previously.76,77The InN,In0.65Ga0.35N,and In0.75Al0.25N were highly autodoped n type(ϳ1020,ϳ1019,and8ϫ1018cmϪ3,respectively͒due to thepresence of native defects.The samples were rinsed in H2O:NH4OH͑20:1͒for1min just prior to deposition of the metal to remove nativeoxides.The metal contacts were sputter deposited to a thick-ness of1000Åin the case of W and WSi0.44͑film composi-tion͒,and then etched in SF6/Ar in a plasma-therm reactive ion etcher͑RIE͒to create transmission line method͑TLM͒patterns.For the Ti/Al contacts,200Åof Ti and then1000Åof Al was deposited and the TLM pattern formed by lift-off of the resist mask.The nitride samples were subsequently etched in Cl2/CH4/H2/Ar in an electron cyclotron resonance ͑ECR͒etcher to produce the mesas for the TLM patterns.78 The samples were annealed at temperatures from300to 900°C for15s under a nitrogen ambient in a rapid thermal annealing͑RTA͒system͑AG-410͒.The contact resistance for W,WSi x,and Ti/Al ohmiccontacts to InGaN as a function of annealing temperature isshown in Fig.3.All contacts had similar contact resistanceas deposited,2–4ϫ10Ϫ7⍀cmϪ2.Above600°C,the Ti/Al contacts degraded rapidly,and the WSi x continued to de-grade,while R C for both samples increased up to ϳ10Ϫ5⍀cm2at900°C.The error in these measurements was estimated to beϮ10%due mainly to geometrical con-tact size effects.The widths of the TLM pattern spacings varied slightly due to processing͑maximum ofϮ5%͒as de-termined by scanning electron microscopy͑SEM͒measure-ments,which were taken into account when calculating the contact resistances.SEM micrographs of W and Ti/Al contacts on InGaN asgrown and annealed showed the W was still quite smootheven after900°C anneal,while the Ti/Al had significant pit-ting at the lowest anneal of500°C even though the contact resistance did not degrade untilу600°C.Auger electron spectroscopy͑AES͒showed that the degradation was due to out-diffusion of In and N.The contact resistance for ohmic contacts of W,WSi x, and Ti/Al to InN as a function of annealing temperature is shown in Fig.4.As-deposited samples had similar contact resistances to InGaN,indicating a similar conduction mecha-nism.WSi x contacts showed the most degradation at low temperature,with the resistance rising a factor of5after 300°C annealing and then remaining constant.Ti/Al devi-ated little from initial values,although there was severe pit-ting on samples annealed at500°C while W resistance began to degrade at500°C.In Fig.5the contact resistance is shown for W,WSi x, and Ti/Al ohmic contacts to InAlN as a function of annealing temperature.As-deposited Ti/Al had the lowest contact re-sistance on this material,R Cϳ10Ϫ4⍀cm2.Tungsten hadthe FIG.3.Contact resistance for W,WSi0.44,and Ti/Al ohmic contacts to InGaN as a function of annealingtemperature.FIG.4.Contact resistance for ohmic contacts of W,WSi x,and Ti/Al to InN as a function of annealing temperature.highest initial contact resistance,R C ϳ10Ϫ2⍀cm 2.The con-tacts showed morphological stability up to 400°C ͑Ti/Al ͒and to 800°C ͑W ͒.SEM micrographs of InAlN contacted with W,WSi x ,and Ti/Al as-grown and annealed at 800,700,and 400°C,respectively,were examined.The W on InAlN remained smooth until 800°C,and then began to form hillocks,as did the WSi x contact at 700°C.The Ti/Al began pitting at 400°C.The pitting in the Ti/Al contacts was due to diffusion of the Al through the Ti into the sample.Hillocks appear to be formed from diffusion of In from the nitride sample into the contact layer.In summary,W,WSi x ,and Ti/Al were found to produce low resistance ohmic contacts on n ϩInGaN and InN.W con-tacts proved to be the most stable,and also gave the lowest resistance to InGaN and InN,R C Ͻ10Ϫ7⍀cm 2after 600°C anneal,and 1ϫ10Ϫ7⍀cm 2after 300°C anneal,respec-tively.Significant interdiffusion of In,N,and Al,as well as Ti and W,were found after annealing.The contact resistance stability varies for each material and degraded at tempera-tures Ͼ400°C on InN,у500°C on InAlN,and у600°C on InGaN.Only W contacts remained smooth at the highest anneal temperature.3.p-ohmic contactsOne of the life-limiting factors in GaN laser diodes to date has been the p -ohmic contact.79Due to the relatively poor specific contact resistance (R C )achievable,the metal-lization will heat-up as current flows across the p -n junction,leading to metal migration down threading dislocations and eventual shorting of the junction.79Removal of the disloca-tions,such as in epitaxial lateral overgrowth structures,will greatly extend the device lifetime.79There are a number of contributing factors to the high R C values for contacts on p -GaN,including:͑i ͒The absence of a metal with a sufficiently high workfunction ͑the band gap of GaN is 3.4eV,and the electron affinity is 4.1eV,but metal work functions are typically р5eV ͒.͑ii ͒The relatively low hole concentrations in p -GaN due to the deep ionization level of the Mg acceptor ͑ϳ170meV ͒.͑iii ͒The tendency for the preferential loss of nitrogen from the GaN surface during processing,which may produce surface conversion to n -type conductivity.In the search for improved contact characteristics,a wide variety of metallizations have been investigated on p -GaN besides the standard Ni/Au,80–85including Ni,82,83,86Au,82,85,87,88Pd,82Pd/Au,88,89Pt/Au,84Au/Mg/Au,61,87Au/C/Ni,90Ni/Cr/Au,88,91and Pd/Pt/Au.84Typically Ni,Pd,or Pt is the metal in direct contact with the GaN,and the structure is annealed at 400–750°C.This produces contact resistances in the 10Ϫ1–10Ϫ3⍀cm 2range.For higher tem-peratures severe degradation in contact morphology is ob-served,usually resulting from the formation of the metal gallides.To examine thermal stability of contacts,p -type (N A ϭ1018cm Ϫ3),Mg-doped GaN layers 1␮m thick were grown on Al 2O 3substrates by molecular beam epitaxy ͑MBE ͒using solid Ga and radio-frequency ͑rf ͒plasma-activated N 2.92Strong cathodoluminescent was observed at ϳ385nm,with very little deep level emission,indicative of high quality ma-terial.Undoped GaN layers ϳ3␮m thick were grown on Al 2O 3by metal organic chemical vapor deposition,with similar cathodoluminescent properties to the MBE material.These samples were implanted with 100keV Si ϩions at a dose of 5ϫ1015cm Ϫ2,and annealed with AlN caps in place to 1400°C for 10s.93This produced a peak n -type doping concentration of ϳ5ϫ1020cm Ϫ3.W or WSi 0.45layers ϳ1000Åthick were deposited using an MRC501sputtering system.The sample position was biased at 90V with respect to the Ar discharge.Prior to sputtering,the native oxide was removed in a 201H 2O:NH 4OH solution.Transmission line patterns were defined by dry etching the exposed metal with SF 6/Ar,and forming mesas around the contact pads using BCl 3/N 2dry etching to confine the current flow.For com-parison,on the p -GaN,Au ͑1000Å͒/Ni ͑500Å͒was depos-ited by e -beam evaporation,defined by lift-off and mesas formed by dry etching.Both n -and p -type samples were annealed for 60s ͑in some experiments this was varied for 30–300s ͒at 300–1000°C under flowing N 2.From Fermi–Dirac statistics we can calculate the Fermi level position E F for p -GaN containing 1018acceptors cm Ϫ3as a function of absolute temperature T,fromN A11ϩ2exp ͓͑E a ϪE F ͒/kT ͔ϭN V exp ͓Ϫ͑E F ϪE V ͒/kT ͔,where N A is the acceptor concentration,E a ϭ171meV for Mg in GaN,and N v is the valence band density of ing this relation,we calculated the ionization efficiency for Mg as a function of sample temperature,as shown in Fig.6.Since the hole concentration in the p -GaN will increase rapidly with temperature,we would expect better ohmic con-tact properties at hightemperatures.FIG.5.Contact resistance for W,WSi x ,and Ti/Al ohmic contacts to InAlN as a function of annealing temperature.Figure 7shows annealing temperature dependence of the current–voltage (I –V )characteristics of the Ni/Au,W,and WSi on p -GaN,with the measurements made at 25°C in all cases.Note that for the optimum anneal temperatures ͑700°C for Ni/Au and W,and 800°C for WSi x ),the con-tacts are not ohmic,but are more accurately described as leaky Schottky diodes.In the case of W and WSi,we assume that annealing above the optimum temperature produces loss of N 2and poorer contact properties.The contact morphology on the W and WSi metalliza-tion remained featureless to the highest temperature we in-vestigated.This is in sharp contrast to the case of Ni/Au,as shown in Fig.8.For the latter metallization,islanding is quite severe after 700°C annealing due to reaction of the Ni with the GaN.94,95From the earlier discussion,we would expect the contact properties to improve at elevated temperatures because of the increased hole density and more efficient thermionic hole emission across the metal–GaN interface.Figure 9shows the I –V characteristics for the 700͑Ni/Au and W ͒or 800°C ͑WSi ͒annealed samples,as a function of the measurement temperature ͑25–300°C ͒.For the Ni/Au,the contacts be-come ohmic at у200°C,while for W and WSi x this occurs at 300°C.Table I shows the R C values at 300°C are 9.2ϫ10Ϫ2͑Ni/Au ͒, 6.8ϫ10Ϫ2͑W ͒,and 2.6ϫ10Ϫ2⍀cm 2͑WSi ͒.The TLM measurements showed that the substrate sheet resistance is reduced from 1.39ϫ104⍀/ᮀat 200°C,to 8470⍀/ᮀat 250°C,and 4600⍀/ᮀat 300°C,indicating that the increased hole concentration plays a major role in decreasing R C .There was not a strong dependence of the room tempera-ture I –V characteristics on annealing time.An example is shown in Fig.10for W/p -GaN,annealed at 700°C.There is little change in the characteristics for 30–120s,but the con-tacts become more rectifying for longer times,probably due to the onset of metal-semiconductor reactions.As a comparison to n -type GaN,Fig.11shows the an-nealing temperature dependence of R C for W contacts on Si-implanted ͑n type ͒GaN.The specific contact resistance improves with annealing up to ϳ950°C,which appears to correspond to the region where the ␤-W 2N interfacial phase is formed.Cole et al.49reported that W and WSi contacts on GaN annealed in the range of 750–850°C showed the mini-mum degree of metal protrusion in the interfacial regions containing the ␤-W 2N phase,whereas at lower annealing temperatures the horizontal spatial extent of this phase was smaller and allowed more protrusions to develop.Excellent structural stability of the W on GaN was shown in SEM micrographs,where a sharp interface was retained after 750°C annealing.In summary,one of the emerging applications for GaN is in ultrahigh power electronic switches,where thermal stabil-ity of the contact metallization will be of paramount impor-tance.Tungsten-based contacts on both n -and p -type GaN offer superior thermal stability to the standard metallization used in photonic devices,TiAl and Ni/Au,respectively.B.Rapid thermal processing 1.Surface protectionThe usual environment for high-temperature annealing of III nitrides is NH 3,79but this is inconvenient for processes such as rapid thermal annealing for implant activation,con-tact annealing for implant isolation.In those situationsweFIG.7.Annealing temperature dependence of I –V characteristics of WSi,W,and Ni/Au contacts on p -GaN ͑60s anneal times ͒.FIG.6.Ionization efficiency of Mg acceptors in GaN and Fermi level po-sition for GaN containing 1018cm Ϫ3Mg acceptors,as a function of tempera-ture.would like to provide some form of N 2overpressure to mini-mize loss of nitrogen from the semiconductor surface at high temperature.96With conventional III–V materials such as GaAs and InP this is achieved in several ways,97–110namely by two methods:͑i ͒placing the sample of interest face down on a substrate of the same type,101,107so that the onset of preferential As or P loss quickly suppresses further loss.The disadvantages of this method include the fact that some group V atoms are lost from the near surface.There is al-ways a possibility of mechanical abrasion of the face of the sample of interest,and contamination can easily be trans-ferred from the dummy wafer to the one of interest.The second method involves:͑ii ͒placing the wafer in a SiC-coated graphite susceptor,109,110which either has had its in-ternal surfaces coated with As or P by heating a sacrificial wafer within it,or in which granulated or powdered GaAs or InP is placed in reservoirs connected to the region in which the wafer is contained.In both cases subsequent heating of the susceptor produces an As or P vapor pressure above the surface of the process wafer,suppressing loss of the group V element.The former approach is widely used in III–V research and is known as the proximity geometry.The latter approach is widely used in industry for anneal processes for GaAs and to a lesser extent InP.It would be convenient for GaN device processing if development of a similar process for rapid thermal process-ing of III nitrides occurred,in which an overpressure of N 2is supplied to a susceptor.In this section we compare use of powdered AlN or InN as materials for use in the susceptor reservoirs,and compare the results with those obtained by simple proximity annealing.The GaN,AlN,InN,InGaN,and InAlN samples were grown using metal organic molecular beam epitaxy on semi-insulating,͑100͒GaAs substrates or Al 2O 3c -plane substrates in an Intevac Gen II system as described previously.76,77The group-III sources were triethylgallium,trimethylamine alane,and trimethylindium,respectively,and the atomic nitrogen was derived from an electron cyclotron resonance Wavemat source operating at 200W forward power.The layers were single crystal with a high density (1011Ϫ1012cm Ϫ2)of stacking faults and microtwins.The GaN and AlN were re-sistive as-grown,and the InN was highly autodoped n type (Ͼ1020cm Ϫ3)due to the presence of native defects.InAlN and InGaN were found to contain both hexagonal and cubic forms.The In 0.75Al 0.25N and In 0.5Ga 0.5N were conducting n type as grown (ϳ1020cm Ϫ3)due to residual autodoping by native defects.The samples were annealed either ͑i ͒face down on samples of the same type,i.e.,GaN when annealingGaN,FIG.8.SEM micrographs of Ni/Au contacts on p -GaN after 60s anneals at either 400͑top left ͒or 700°C ͑top right ͒,or W contacts after similar annealing at 400͑bottom left ͒or 900°C ͑bottom right ͒.。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

Journal of Financial Economics67(2003)217–248Testing the pecking order theory of capitalstructure$Murray Z.Frank a,*,Vidhan K.Goyal ba Faculty of Commerce,University of British Columbia,Vancouver,BC Canada,V6T1Z2b Department of Finance,Hong Kong University of Science and Technology,Kowloon,Hong KongReceived22March2001;accepted9January2002AbstractWe test the pecking order theory of corporate leverage on a broad cross-section of publicly traded Americanfirms for1971to1998.Contrary to the pecking order theory,net equity issues trackthefinancing deficit more closely than do net debt issues.While largefirms exhibit some aspects of pecking order behavior,the evidence is not robust to the inclusion of conventional leverage factors,nor to the analysis of evidence from the1990s.Financing deficit is less important in explaining net debt issues over time forfirms of all sizes.r2002Elsevier Science B.V.All rights reserved.JEL classification:G32Keywords:Pecking order theory;Capital structure;Financing deficit$We would like to thank Mike Barclay(the referee),Ken Bechmann,Robert Chirinko,Sudipto Dasgupta,Charles Hadlock,Keith Head,Vojislav Maksimovic,Sheridan Titman,and Karen Wruck,for helpful comments.Feedbackfrom the seminar participants at the2000European Finance Association meetings,the2000Financial Management Association meetings,the11th Annual Financial Economics and Accounting Conference(2000)at the University of Michigan,2001American Finance Association meetings,the2001Rutgers Conference on Corporate Finance,the University of Hong Kong,the University of Victoria,and the Hong Kong University of Science and Technology are appreciated.Murray Frankthank s the B.I.Ghert Family Foundation forfinancial support.We alone are responsible for the contents and any errors.*Corresponding author.Tel.:+1-604-822-8480;fax:+1-604-822-8477.E-mail addresses:murray.frank@commerce.ubc.ca(M.Z.Frank),goyal@ust.hk(V.K.Goyal). 0304-405X/02/$-see front matter r2002Elsevier Science B.V.All rights reserved.PII:S0304-405X(02)00252-01.IntroductionThe pecking order theory of capital structure is among the most influential theories of corporate leverage.According to Myers (1984),due to adverse selection,firms prefer internal to external finance.When outside funds are necessary,firms prefer debt to equity because of lower information costs associated with debt issues.Equity is rarely issued.These ideas were refined into a key testable prediction by Shyam-Sunder and Myers (1999).The financing deficit should normally be matched dollar-for-dollar by a change in corporate debt.As a result,if firms follow the pecking order,then in a regression of net debt issues on the financing deficit,a slope coefficient of one is observed.Shyam-Sunder and Myers (1999)find strong support for this prediction in a sample of 157firms that had traded continuously over the period 1971to 1989.This is an attractive and influential result.The pecking order is offered as a highly parsimonious empirical model of corporate leverage that is descriptively reasonable.Of course,157firms is a relatively small sample from the set of all publicly traded American firms.It is therefore important to understand whether the pecking order theory is broadly applicable.In this paper,we study the extent to which the pecking order theory of capital structure provides a satisfactory account of the financing behavior of publicly traded American firms over the 1971to 1998period.Our analysis has three elements.First,we provide evidence about the broad patterns of financing activity.This provides the empirical context for the more formal regression tests.It also serves as a checkon the significance of external finance and equity issues.Second,we examine a number of implications of the pecking order in the context of Shyam-Sunder and Myers’(1999)regression tests.Finally,we checkto see whether the peck ing order theory receives greater support among firms that face particularly severe adverse selection problems.The pecking order theory derives much of its influence from a view that it fits naturally with a number of facts about how companies use external finance.1Myers (2001)reports that external finance covers only a small proportion of capital formation and that equity issues are minor,with the bulkof external finance being debt.These key claims do not match the evidence for publicly traded American firms,particularly during the 1980s and 1990s.External finance is much more significant than is usually recognized in that it often exceeds investments.Equity finance is a significant component of external finance.On average,net equity issues commonly exceed net debt issues.Particularly striking is the fact that net equity issues trackthe financing deficit much more closely than do net debt issues.Shyam-Sunder and Myers (1999)focus on a regression test of the pecking order.In this test one needs to construct the financing deficit from information in the corporate accounts.The financing deficit is constructed from an aggregation of1The pecking order theory also derives support from indirect sources of evidence.Eckbo (1986)and Asquith and Mullins (1986)provide event study evidence that adverse selection is more significant for equity issues than for debt issues.Cadsby et al.(1990)provide experimental evidence of adverse selection in company financing.M.Z.Frank,V.K.Goyal /Journal of Financial Economics 67(2003)217–248218M.Z.Frank,V.K.Goyal/Journal of Financial Economics67(2003)217–248219 dividends,investment,change in working capital and internal cashflows.If the pecking order theory is correct,then the construction of thefinancing deficit variable is a justified aggregation.Under the pecking order,each component offinancing deficit should have the predicted dollar-for-dollar impact on corporate debt.The evidence does not support this hypothesis.Even if a theory is not strictly correct,when compared to other theories it might still do a better job of organizing the available evidence.The pecking order is a competitor to other mainstream empirical models of corporate leverage.Major empirical alternatives such as the model tested by Rajan and Zingales(1995)use a different information set to account for corporate leverage.It is therefore of interest to see how thefinancing deficit performs in a nested model that also includes conventional factors.The pecking order theory implies that thefinancing deficit ought to wipe out the effects of other variables.If thefinancing deficit is simply one factor among many thatfirms tradeoff,then what is left is a generalized version of the tradeoff theory.Wefind that thefinancing deficit does not wipe out the effects of conventional variables.The information in thefinancing deficit appears to be factored in along with many other things thatfirms take into account.This is true acrossfirm sizes and across time periods.Since the pecking order does not explain broad patterns of corporatefinance,it is natural to examine narrower sets offirms.According to the pecking order theory,financing behavior is driven by adverse selection costs.The theory should perform best amongfirms that face particularly severe adverse selection problems.Small high-growthfirms are often thought of asfirms with large information asymmetries. Contrary to this hypothesis,small high-growthfirms do not behave according to the pecking order theory.In fact,the pecking order works best in samples of large firms that continuously existed during the1970s and rgefirms with long uninterrupted trading records are not usually considered to befirms that suffer the most acute adverse selection problems.To understand the evidence it is important to recognize the changing population of publicfipared to the1970s and1980s,many more small and unprofitablefirms became publicly traded during the1990s.Since smallfirms generally do not behave according to the pecking order,this accounts for part of the reason that the pecking order theory is rejected.But the time period has a stronger effect than just this.Forfirms of all sizes,thefinancing deficit plays a declining role over time.Previous literature provides other evidence pertinent to a general assessment of the pecking order theory.The pecking order theory predicts that high-growthfirms, typically with largefinancing needs,will end up with high debt ratios because of a manager’s reluctance to issue equity.Smith and Watts(1992)and Barclay et al. (2001)suggest precisely the opposite.High-growthfirms consistently use less debt in their capital structure.The pecking order theory makes predictions about the maturity and priority structure of debt.Securities with the lowest information costs should be issuedfirst, before thefirm issues securities with higher information costs.This suggests thatshort-term debt should be exhausted before the firm issues long-term debt.Capitalized leases and secured debt should be issued before any unsecured debt is issued.Barclay and Smith (1995a,b)find that 50%of their firm-year observations have no debt issued with less than one-year maturity,23%have no secured debt,and 54%have no capital leases.It seems difficult to understand this evidence within a pure pecking order point of view.Chirinko and Singha (2000)question the interpretation of the Shyam-Sunder and Myers (1999)regression test.Chirinko and Singha show that equity issues can create a degree of negative bias in the Shyam-Sunder and Myers test.Suppose that firms actually follow the pecking order theory,but that these firms issue an empirically observed amount of equity.In that case,they show that the predicted regression coefficient is actually 0.74rather than one.This amount of bias is not trivial,but it still leaves the coefficient very far from the magnitudes of slope coefficients that are observed.Chirinko and Singha also point out that if,contrary to the pecking order,firms follow a policy of using debt and equity in fixed proportions,then the Shyam-Sunder and Myers regression will identify this ratio.As a result,finding a coefficient near one would not disprove the tradeoff theory.Chirinko and Singha’s cautionary note reinforces an important methodological point.Most empirical tests have various weaknesses.It is therefore important to examine the predictions of a theory from a number of points of view rather than relying solely on a single test.The structure of the rest of this paper is as follows.Section 2presents the pecking order theory and the associated empirical hypotheses.The data are described in Section 3.Section 4presents the empirical results.Conclusions are presented in Section 5.2.TheoryThe pecking order theory is from Myers (1984)and Myers and Majluf (1984).Since it is well known,we can be brief.Suppose that there are three sources of funding available to firms:retained earnings,debt,and equity.Retained earnings have no adverse selection problem.Equity is subject to serious adverse selection problems while debt has only a minor adverse selection problem.From the point of view of an outside investor,equity is strictly riskier than debt.Both have an adverse selection riskpremium,but that premium is large on equity.Therefore,an outside investor will demand a higher rate of return on equity than on debt.From the perspective of those inside the firm,retained earnings are a better source of funds than is debt,and debt is a better deal than equity financing.Accordingly,the firm will fund all projects using retained earnings if possible.If there is an inadequate amount of retained earnings,then debt financing will be used.Thus,for a firm in normal operations,equity will not be used and the financing deficit will match the net debt issues.In reality,company operations and the associated accounting structures are more complex than the standard pecking order representation.This implies that in order to test the pecking order,some form of aggregation must be used.M.Z.Frank,V.K.Goyal /Journal of Financial Economics 67(2003)217–248220We define notation as follows:DIV t cash dividends in year t;I t net investment in year t(i.e.,I t¼capital expenditures+increase in invest-ments+acquisitions+other use of fundsÀsale of PPEÀsale of investment);D W t change in working capital in year t(i.e.,D W t¼change in operating workingcapital+change in cash and cash equivalents+change in current debt);C t cashflow after interest and taxes(i.e.,C t¼income before extraordinaryitems+depreciation and amortization+extraordinary items and discontinued operations+deferred taxes+equity in net lossÀearnings+other funds from operations+gain(loss)from sales of PPE and other investments);R t current portion of the long-term debt in year t;D D t net debt issued in year t;(i.e.,D D t=long-term debt issuanceÀlong-term debtreduction);D E t Net equity issued in year t(i.e.,D E t¼sale of common stockminus stockrepurchases).Using this notation,we can use theflow of funds data to provide a partially aggregated form of the accounting cashflow identity as,DEF t¼DIV tþI tþD W tÀC t¼D D tþD E t:ð1ÞShyam-Sunder and Myers(1999)argue that under the pecking order hypothesis, after an Initial Public Offering(IPO),equity issues are only used in extreme circumstances.The empirical specification is thus given asD D it¼aþb PO DEF itþe it;ð2Þwhere e it is a well-behaved error term.In Eq.(2),the pecking order hypothesis is that a¼0and b PO¼1:Shyam-Sunder and Myers(1999)find that the pecking order model is statistically rejected.However it does provide a goodfirst-order approximation of their data.In contrast to the accounting definition,Shyam-Sunder and Myers(1999)include the current portion of long-term debt as part of thefinancing deficit beyond its role in the change in working capital.Following their argument,the relevantflow of funds deficitðDEF SSMtÞis defined asDEF SSMt¼DIV tþI tþD W tþR tÀC t:ð3ÞIf their alternative version of thefinancing deficit is to be used,then replace DEF itwith DEF SSMt in Eq.(2).We try both approaches andfind that empirically thecurrent portion of long-term debt does not appear to belong in the definition of DEF it:With the exception of column(7)of Table5,we report only the results for which the current portion of long-term debt is not included as a separate component of thefinancing deficit.This choice favors the pecking order,but it does not affect our conclusions.How is cash to be treated in thefinancing deficit?Changes in cash and cash equivalents are included with changes in working capital.Cash could be correlated with the amount of debt issued.This could arise in the presence of lumpy debt and M.Z.Frank,V.K.Goyal/Journal of Financial Economics67(2003)217–248221equity issues,with excess proceeds parked for some period of time in excess cash balances.If this takes place over a number of years,a more complex dynamic theory of leverage is needed.We report results in which the change in cash and cash equivalents are included.This choice favors the pecking order,but the conclusions are not affected.In a panel regression,one can treat all year-firm combinations as equally important independent observations.If that is done,then a simple regression can be run.If one is willing to accept the classical error term assumptions,then standard fixed-effects or random-effects panel estimators may be used.Yet another possibility is to downplay the differences across time and focus on the cross-sectional differences.One could follow Fama and MacBeth (1973)and use the average of a series of annual cross-sectional regressions as the point estimate and use the time series of these estimates to construct standard errors.This is the approach taken by Fama and French (2002).We have tried these alternatives and our conclusions are not sensitive to the choice of approach.According to theory,the specification in Eq.(2)is defined in levels.When actually estimating Eq.(2),it is conventional to scale the variables by assets or by sales.The pecking order theory does not require such scaling.Of course,in an algebraic equality if the right-hand side and the left-hand side are divided by the same value,the equality remains intact.However,in a regression the estimated coefficient can be seriously affected if the scaling is by a variable that is correlated with the variables in the equation.Scaling is most often justified as a method of controlling for differences in firm size.The reported results are based on variables scaled by net assets (total assets minus current liabilities).We replicate all the tests by scaling variables by total bookassets,by the sum of bookdebt plus mark et equity,and by sales.The results are very similar and do not affect our conclusions.There is an important econometric issue that needs to be addressed.The pecking order theory treats the financing deficit as exogenous.The financing deficit includes investment and dividends.Yet,much financial theory is devoted to attempting to understand the determinants of these factors.As a result,it is not entirely obvious that the components of the financing deficit should be properly regarded as exogenous.If they are truly endogenous,then the regression in Eq.(2)is misspecified.If a model is misspecified,then small changes to the specification may lead to large changes in the coefficient estimates.The model is also likely to be unstable across time periods and its performance would likely not generalize to other samples of firms.Such instability would itself be indicative of a failure of the model.In order to deal with these concerns,two steps are taken.First,all tests are subjected to a large number of robustness checks.In most cases the findings are robust.However,the findings are not robust on one crucial dimension.Requiring firms to have complete trading records over the period 1971–1989makes a big difference to the coefficient estimates.Second,the ability of the estimated models to predict debt issues by a holdout sample of firms is directly examined.This is a simple way to address concerns about model misspecification.A model may fit well within sample but its performance may not generalize.Such a model is,of course,much less interesting than an empirical model that also performs well out of sample.M.Z.Frank,V.K.Goyal /Journal of Financial Economics 67(2003)217–248222M.Z.Frank,V.K.Goyal/Journal of Financial Economics67(2003)217–248223The ability of eachfitted model to predict is tested on data from thefive years subsequent to(or prior to)the time period over which the model isfit.For eachfirm year in the holdout sample,we plug the actual values of the exogenous variables into thefitted equation.This provides a predicted value of the endogenous variable (usually net debt issue).In this manner,we obtainfive years offirm-specific predictions from eachfitted model.To assess the quality of these predictions,the predicted debt issues are regressed against a constant and the actual debt issues.A goodfit will be reflected in an intercept of zero,a slope of one,and a high R2:In order to save space we only report the R2that is obtained on the hold out sample.ing the same information:disaggregation of thefinancing deficitTo test the pecking order theory we need to aggregate the accounting data.Is the aggregation step justified?It seems plausible that there could be information in DEF it that helps to account for D D it;but not in the manner hypothesized by the pecking order theory.An easy way to check whether the aggregation step is justified is to run the equation on a disaggregated basis and then checkwhether the data satisfies the aggregation step.Consider the following specification,D D it¼aþb DIV DIV tþb I I tþb W D W tÀb C C tþe it:ð4ÞUnder the pecking order theory,it is DEF it itself that matters.A unit increase in any of the components of DEF it must have the same unit impact on D D it:The pecking order hypothesis is thus b DIV¼b I¼b W¼b C¼1:If that hypothesis is correct,then the aggregation in Eq.(1)is justified.If however,the significance is actually only driven by some of the individual components,then alternative coefficient patterns are possible.ing other information to account for leverageThe pecking order test implicitly makes different exogeneity assumptions and uses a different information set than is conventional in empirical research on leverage and leverage-adjusting behavior.Harris and Raviv(1991)explain the conventional set of variables and then Rajan and Zingales(1995)distill these variables into a simple cross-sectional model.The conventional set of explanatory factors for leverage is the conventional set for a reason.The variables have survived many tests.As explained below,these variables also have conventional interpretations.Excluding such variables from consideration is therefore potentially a significant omission.It is also true that including such variables potentially poses a tough test for the pecking order theory.The conventional leverage regression is intended to explain the level of leverage, while the pecking order regression is intended to explain the change rather than the level.As long as the shocks are uncorrelated across years,we can equally well run the conventional specification infirst differences.Of course,a lower R2will be obtained. The assumption of uncorrelated shocks is unlikely literally correct.When we run theconventional regression in first differences,we expect to lose some accuracy.Running the conventional regression in first differences may also bias the variable coefficients towards zero.It turns out that this bias is not large enough to alter our conclusions about the relative empirical validity of the two approaches.The benefit is that we then have an appropriate specification in which to nest the financing deficit variable.Alternatively,one could run a regression that explains the level of leverage,then use a cumulated past financing deficit variable to represent the pecking order.If that is done there is an issue about when to start the cumulating.We try such a procedure and obtain results that are very similar to those reported in Table 7.At the heart of the conventional empirical analysis is a regression of leverage on four factors:tangibility of assets (denoted T ),market-to-book ratio (denoted MTB),log sales (denoted LS),and profitability (denoted P ).Let D denote the first differences between years.Our version of the basic regression equation is thereforeD D i ¼a þb T D T i þb MTB D MTB i þb LS D LS i þb P D P i þb DEF DEF i þe i :ð5ÞEq.(5)is simply a conventional regression run in first differences but with financing deficit as an added factor.In the conventional regression,this term is not present.From the perspective of testing the pecking order,the most important of the conventional variables is tangibility.According to Harris and Raviv (1991),under the pecking order theory,one might expect that firms with few tangible assets would have greater asymmetric information problems.Thus,firms with few tangible assets will tend to accumulate more debt over time and become more highly levered.Hence,Harris and Raviv argue that the pecking order predicts that b T o 0:This is not the conventional prediction regarding the role of tangibility.A more common idea is based on the hypothesis that collateral supports debt.It is often suggested that tangible assets naturally serve as collateral.Hence,collateral is associated with increased leverage.The usual prediction is that b T >0:Firms with high market-to-book ratios are often thought to have more future growth opportunities.As in Myers (1977),there may be a concern that debt could limit a firm’s ability to seize such opportunities when they appear.Goyal et al.(2002)find that when growth opportunities of defense firms decline,these firms increase their use of debt financing.Barclay et al.(2001)present a model showing that the debt capacity of growth options can be negative.The common prediction is that b MTB o 0:Large firms are usually more diversified,have better reputations in debt markets,and face lower information costs when borrowing.Therefore,large firms are predicted to have more debt in their capital structures.The prediction is that b LS >0:The predictions on profitability are ambiguous.The tradeoff theory predicts that profitable firms should be more highly levered to offset corporate taxes.Also,in many asymmetric information models,such as Ross (1977),profitable firms are predicted to have higher leverage.But Titman and Wessels (1988)and Fama and French (2002)show that this is not a common finding.Instead,the literature finds profits and leverage to be negatively correlated.While MacKay and Phillips (2001)challenge this common finding,we expect to find that b P o 0:M.Z.Frank,V.K.Goyal /Journal of Financial Economics 67(2003)217–248224M.Z.Frank,V.K.Goyal/Journal of Financial Economics67(2003)217–248225 Fama and French(2002)note that the negative relationship between profits and leverage is consistent with the pecking order theory.But the pecking order is not the only possible interpretation of the relationship.There are at least two issues.First, current profitability can also serve as a signal of investment opportunities.There is a large macro-finance literature,including studies by Gilchrist and Himmelberg(1995) and Kaplan and Zingales(1997),in which this interpretation issue plays a key role.It is well known that it is difficult to construct a convincing proxy for investment opportunities.If Tobin’s q or the simpler measure,market-to-book assets,is measured with error,then it may not adequately control for the information content in afirm’s profitability.For an analysis of measurement error in this context,see Erickson and Whited(2000).The second issue is thatfirms may facefixed costs of adjustment.Fischer et al. (1989)analyze the effect of havingfixed costs associated with actively adjusting leverage.When afirm earns profits,debt gets paid off and leverage falls automatically.Only periodically will large readjustments be made in order to capture the tax benefits of leverage.Empirically,most of the data reflects the process of paying off the debt by using profits.Thus,profitablefirms will be less levered even if the tradeoff theory is at workand the adjustment costs are tak en into account.3.DataWe need the data from fundsflow statements to test the pecking order theory.This restricts the beginning of the sample period to1971since that is when American firms started reporting fundsflow statements.The data ends with1998.Variables are deflated to constant1992dollars.Following standard practice,financialfirms(6000–6999),regulated utilities(4900–4999),andfirms involved in major mergers(Compustat footnote code AB)are excluded.2Also excluded arefirms with missing bookvalue of assets and a small number offirms that reported format codes4,5,pustat does not define format codes4and6.Format code5is for the Canadianfile.The balance sheet and cashflow statement variables as a percentage of assets are trimmed to remove the most extreme0.50%in either tail of the distribution.This serves to remove outliers and the most extremely misrecorded data.3The balance sheet presentation of corporate assets and liabilities is reasonably consistent over time.Average common-size balance sheets for a number of years between1971and1998are presented in Table1.The asset side shows significant changes over time.Cash increases dramatically over the period,going from5%to 2Leaving in the data onfirms involved in major mergers had no material effect on our conclusions.We, therefore,do not report these results separately.3Prior to trimming,several balance sheet and cashflow statement items are recoded as zero if they were reported missing or combined with other data items in Compustat.The data is often coded as missing when afirm does not report a particular item or combines it with other data items.After examining accounting identities,we determine that recoding missing values on these items as zero respects the reported accounting identities.。

相关文档
最新文档