Localization transition on complex networks via spectral statistics
- 格式:pdf
- 大小:152.19 KB
- 文档页数:6
PMSF(100mM)产品简介:PMSF即Phenylmethanesulfonyl fluoride,中文名为苯甲基磺酰氟。
分子式为C7H7FO2S,分子量为174.19,纯度>99%。
常用生化试剂,用于抑制蛋白酶,进口试剂配制。
第一次使用时把一管PMSF晶体全部倒入到PMSF溶剂中,溶解混匀,即配制成10毫升100mM PMSF。
保存条件:室温保存,配成PMSF溶液后需-20ºC保存。
注意事项:有毒,注意防护。
本产品仅限于专业人员的科学研究用,不得用于临床诊断或治疗,不得用于食品或药品,不得存放于普通住宅内。
为了您的安全和健康,请穿实验服并戴一次性手套操作。
使用本产品的文献:1.Cai H, Yin D, Zhang L, Yang X, Xu X, Liu W, Zheng X, Zhang H, Wang J, Xu Y,Cheng D, Zheng M, Han Y, Wu M, Wang Y. Preparation and biological evaluation of 2-amino-6-[18F]fluoro-9-(4-hydroxy-3-hydroxy-methylbutyl) purine (6-[18F]FPCV) as a novel PET probe for imaging HSV1-tk reporter gene expression.Nucl Med Biol. 2007 Aug;34(6):717-25.2.Zhang JQ, Zhao XK, He J, Zhu L, Wang XJ. Express and role of FGFR2 in bladdertransitional cell carcinoma. China Medical Engineering. 2007 Sep;V ol.15(9):720-5.3.Sun YG, Wang XW, Yang SM, Zhou G, Wang WQ, Wang HB, Wang RQ, Fang DC.Inhibition of nucleostemin upregulates CDX2 expression in HT29 cells in response to bile acid exposure: implications in the pathogenesis of Barrett's esophagus. J Gastrointest Surg. 2009 Aug;13(8):1430-9.4.Hua W, Jiang J, Rong X, Wu R, Qiu H, Zhang Y, Chen Q. The dual role of thecystathionine gamma-lyase/hydrogen sulfide pathway in CVB3-induced myocarditis in mice. Biochem Biophys Res Commun. 2009 Oct 23;388(3):595-600.5.Li H, Zhang L, Huang Q. Differential expression of mitogen-activated proteinkinase signaling pathway in the hippocampus of rats exposed to chronic unpredictable stress. Behav Brain Res. 2009 Dec 14;205(1):32-7.6.Cai J, Yu C, Liu Y, Chen S, Guo Y, Yong J, Lu W, Ding M, Deng H. Generation ofhomogeneous PDX1(+) pancreatic progenitors from human ES cell-derived endoderm cells. J Mol Cell Biol. 2010;2(1):50-60.7.Li R, Zhang J, Zhang L, Cui Q, Liu H. Angelica injection promotes peripheralnerve structure and function recovery with increasedexpressions of nerve growth factor and brain derived neurotrophic factor in diabetic rats. Curr Neurovasc Res.2010 Aug;7(3):213-22.8.Yu LN, Yu J, Zhang FJ, Yang MJ, Ding TT, Wang JK, He W, Fang T, Chen G, YanM. Sevoflurane postconditioning reduces myocardial reperfusion injury in rat isolated hearts viaactivation of PI3K/Akt signaling and modulation of Bcl-2 family proteins. J Zhejiang Univ Sci B. 2010 Sep;11(9):661-72.9.Wang C, Xie H, Liu X, Qin Q, Wu X, Liu H, Liu C. Role of nuclear factor-κB involatile anaesthetic preconditioning with sevoflurane during myocardial ischaemia/reperfusion. Eur J Anaesthesiol. 2010 Aug;27(8):747-56.10.Wang W, Wang L, Li XX, Chen X, Zhang HY, He Y, Wang JJ, Zhao YY, Zhang BL,Xu YX. Effect of interrupted endogenous BMP/Smad signaling on growth and steroidogenesis of porcinegranulosa cells. J Zhejiang Univ Sci B. 2010 Sep;11(9):719-27.11.Hou T, Shi Y, Cheng S, Yang X, Li L, Xiao C. Nogo-A expresses on neural stemcell surface. Int J Neurosci. 2010 Mar;120(3):201-5. 12.Wang X, Yang Q, Wang P, Luo L, Chen Z, Liao B, Li G. Derp2-mutant genevaccine inhibits airway inflammation and up-regulates Toll-like receptor 9 in anallergic asthmatic mouse model. Asian Pac J Allergy Immunol. 2010 Dec;28(4):287-93.13.Feng DD, Zhang H, Zhang P, Zheng YS, Zhang XJ, Han BW, Luo XQ, Xu L, ZhouH, Qu LH, Chen YQ. Down-regulated miR-331-5p and miR-27a are associated with chemotherapy resistance andrelapse in leukaemia. J Cell Mol Med. 2011 Oct;15(10):2164-75.14.Zhang L, Jiang H, Gao X, Zou Y, Liu M, Liang Y, Yu Y, Zhu W, Chen H, Ge J. Heatshock transcription factor-1 inhibits H2O2-induced apoptosis via down-regulation of reactive oxygen species in cardiac myocytes. Mol Cell Biochem. 2011 Jan;347(1-2):21-8.15.Ding S, Peng H, Fang HS, Zhou JL, Wang Z. Pulsed electromagnetic fieldsstimulation prevents steroid-induced osteonecrosis in rats. BMC Musculoskelet Disord. 2011 Sep 29;12:215.16.Wu J, Song R, Song W, Li Y, Zhang Q, Chen Y, Fu Y, Fang W, Wang J, Zhong Z,Ling H, Zhang L, Zhang F. Chlorpromazine protects against apoptosis induced by exogenous stimuli in the developing ratbrain. PLoS One. 2011;6(7):e21966.17.Liu ML, Wen JJ, Xu XF, Zhao DM. Neurotoxic effect of the complex of the ovineprion protein (OvPrP(C)) and RNA on the cultured ratcortical neurons. Neurochem Res. 2011 Oct;36(10):1863-9.18.Guo L, Zhang J, Yan Q, Yin M. Establishment and characterization of RNA-editedserotonin 2C receptor isoform cell models andalteration of amyloid precursor protein ectodomain secretion in HEK293 APPSwe cells. Hum Cell. 2011 Jun;24(2):104-11.19.Feng S, Cong S, Zhang X, Bao X, Wang W, Li H, Wang Z, Wang G, Xu J, Du B, QuD, Xiong W, Yin M, Ren X, Wang F, He J, Zhang B. MicroRNA-192 targeting retinoblastoma 1 inhibits cell proliferation and induces cell apoptosis inlung cancer cells. Nucleic Acids Res. 2011 Aug;39(15):6669-78.20.Tian Y, Hu Y, Wang Z, Chen K, Zhang L, Wang L, Ren M, Huang A, Tang H.Hepatitis B virus regulates Raf1 expression in HepG2.2.15 cells by enhancing its promoter activity. Arch Virol. 2011 May;156(5):869-74.21.Yang L, Ping YF, Yu X, Qian F, Guo ZJ, Qian C, Cui YH, Bian XW. Gastric cancerstem-like cells possess higher capability of invasion and metastasis in associationwith a mesenchymal transition phenotype. Cancer Lett. 2011 Nov 1;310(1):46-52.22.He Y, Zhang H, Teng J, Huang L, Li Y, Sun C. Involvement of calcium-sensingreceptor in inhibition of lipolysis through intracellular cAMP andcalcium pathways in human adipocytes. Biochem Biophys Res Commun. 2011 Jan 7;404(1):393-9.23.Jiang S, Zu Y, Wang Z, Zhang Y, Fu Y. Involvement of mitochondrial permeabilitytransition pore opening in 7-xylosyl-10-deacetylpaclitaxel-induced apoptosis.Planta Med. 2011 Jul;77(10):1005-12.24.Chen X, Wu X, Zhao Y, Wang G, Feng J, Li Q, Qian G. A novel binding protein ofsingle immunoglobulin IL-1 receptor-related molecule: Paralemmin-3. Biochem Biophys Res Commun. 2011 Jan 28;404(4):1029-33.25.Zhang D, Jin T, Xu YQ, Lu Y, Wu Q, Zhang YK, Liu J. Diurnal-and sex-relateddifference of metallothionein expression in mice. J Circadian Rhythms. 2012 Jul 24;10(1):5.26.Wang L, Qin Y, Tong L, Wu S, Wang Q, Jiao Q, Guo Z, Lin L, Wang R, Zhao W,Zhong Z. MiR-342-5p suppresses coxsackievirus B3 biosynthesis by targeting the 2C-coding region. Antiviral Res. 2012 Feb;93(2):270-9.27.Liu S, Gao S, Wang XY, Wang DB. Expression of miR-126 and Crk inendometriosis: miR-126 may affect the progression ofendometriosis by regulating Crk expression. Arch Gynecol Obstet. 2012 Apr;285(4):1065-72.28.Chen X, Gao YD, Yang J. Elevated interferon regulatory factor 4 levels in patientswith allergic asthma. J Asthma. 2012 Jun;49(5):441-9.29.Chen Y, Wang Z, Xie Y, Guo X, Tang X, Wang S, Yang S, Chen K, Niu Y, Ji W.Folic acid deficiency inhibits neural rosette formation and neuronal differentiation from rhesusmonkey embryonic stem cells. J Neurosci Res. 2012 Jul;90(7):1382-91.30.Li FF, Shen J, Shen HJ, Zhang X, Cao R, Zhang Y, Qui Q, Lin XX, Xie YC, ZhangLH, Jia YL, Dong XW, Jiang JX, Bao MJ, Zhang S, Ma WJ, Wu XM, Shen H, Xie QM, Ke Y. Shp2 Plays an Important Role in Acute Cigarette Smoke-Mediated Lung Inflammation. J Immunol. 2012 Sep 15;189(6):3159-67.31.Jiao J, Hong S, Zhang J, Ma L, Sun Y, Zhang D, Shen B, Zhu C. Opsin3 sensitizeshepatocellular carcinoma cells to 5-fluorouracil treatment by regulating theapoptotic pathway. Cancer Lett. 2012 Jul 1;320(1):96-103.32.Ye XZ, Xu SL, Xin YH, Yu SC, Ping YF, Chen L, Xiao HL, Wang B, Yi L, WangQL, Jiang XF, Yang L, Zhang P, Qian C, Cui YH, Zhang X, Bian XW. Tumor-associated microglia/macrophages enhance the invasion of glioma stem-like cells via TGF-β1 signaling pathway. J Immunol. 2012 Jul 1;189(1):444-53.33.Deng J, Huang Q, Wang F, Liu Y, Wang Z, Wang Z, Zhang Q, Lei B, Cheng Y. Therole of caveolin-1 in blood-brain barrier disruption induced by focused ultrasound combinedwith microbubbles. J Mol Neurosci. 2012 Mar;46(3):677-87.34.Yu X, Tang J, Wang Q, Ye W, Tao K, Duan S, Lu C, Yang X, Dong S, Zheng X,Wang Y. The RxLR effector Avh241 from Phytophthora sojae requires plasma membrane localization toinduce plant cell death. New Phytol. 2012 Oct;196(1):247-60.35.Yu Y, Mu J, Fan Z, Lei G, Yan M, Wang S, Tang C, Wang Z, Yu J, Zhang G. Insulin-like growth factor 1 enhances the proliferation and osteogenic differentiation of humanperiodontal ligament stem cells via ERK and JNK MAPK pathways.Histochem Cell Biol. 2012 Apr;137(4):513-25.36.Xu XY, Shi Y, Zhang PP, Zhang F, Shen YY, Su X, Zhao BL. E-cadherin mediatesadhesion and endocytosis of Aspergillus fumigatus blastospores in humanepithelial cells. Chin Med J (Engl). 2012 Feb;125(4):617-21.37.Liu M, Zhou G, Song W, Li P, Liu H, Niu X, Fan Y. Effect of nano-hydroxyapatiteon the axonal guidance growth of rat cortical neurons. Nanoscale. 2012 May 21;4(10):3201-7.38.Lv XH, Chen JW, Zhao G, Feng ZZ, Yang DH, Sun WW, Fan JS, Zhu GH. N-mycdownstream-regulated gene 1/Cap43 may function as tumor suppressor in endometrialcancer. J Cancer Res Clin Oncol. 2012 Oct;138(10):1703-15.39.Hu LF, Qian AQ, WangY, Di SM, Shang P. Inhibitory Effect of SimulatedMicrogravity on Differentiating Preosteoblasts. Advances in Space Research. 2013 Jan 1;V ol.51(1):107-114.40.Li P, Huang JJ, Ni JJ, Sun FY. VEGF evokes reactive astroglia to convert intoneuronal cells by affecting the biological functionof MeCP2 in adult rat brain after cerebral ischemia. Neurochem Int. 2012 Jul 20. pii: S0197-0186(12)00231-8.41.Du YF, Liang L, Shi Y, Long QZ, Zeng J, Wang XY, He DL. Multi-target siRNAbased on DNMT3A/B homologous conserved region influences cell cycle andapoptosis of human prostate cancer cell line TSU-PR1. Genet Mol Biol. 2012 Jan;35(1):164-71.42.Lu T, Luo Y, Sun H, Qin W, Li Y. Electroacupuncture improves behavioral recoveryand increases SCF/c-kit expression in a ratmodel of focal cerebralischemia/reperfusion. Neurol Sci. 2013 Apr;34(4):487-95.43.Chen HY, Wang JM, Wang HY, Zhang YX, Liu W, Pan L, Wang WH, Chen SF, JinWG, Wang L. Effect of short hairpin RNA-induced CXCR4 silence on ovarian cancer cell. Biomed Pharmacother. 2012 Oct;66(7):549-53.44.Jia P, Xu YJ, Zhang ZL, Li K, Li B, Zhang W, Yang H. Ferric ion could facilitateosteoclast differentiation and bone resorption through the production ofreactive oxygen species. J Orthop Res. 2012 Nov;30(11):1843-52.45.Yin C, Zhou S, Jiang L, Sun X. Mechanical injured neurons stimulate astrocytes toexpress apolipoprotein E through ERKpathway. Neurosci Lett. 2012 Apr 25;515(1):77-81.46.Zhuo C, Ji Y, Chen Z, Kitazato K, Xiang Y, Zhong M, Wang Q, Pei Y, Ju H, WangY. Proteomics analysis of autophagy-deficient Atg7-/- MEFs reveals a close relationship between F-Actin and autophagy. Biochem Biophys Res Commun.2013 Aug 2;437(3):482-8.47.Wang J, Sai K, Chen FR, Chen ZP. miR-181b modulates glioma cell sensitivity totemozolomide by targeting MEK1. Cancer Chemother Pharmacol. 2013 Jul;72(1):147-58.48.Li X, Dong X, Zheng S, Xiao J. Expression and localization of TASK-1, -2 and -3channels in MG63 human osteoblast-like cells. Oncol Lett. 2013 Mar;5(3):865-869.49.Wang J, Li J, Huang Y, Song X, Niu Z, Gao Z, Wang H. Bcl-3 suppresses Tax-induced NF-κB activation through p65 nuclear translocation blockage in HTLV-1-infectedcells. Int J Oncol. 2013 Jan;42(1):269-76.50.Shen X, Tang J, Hu J, Guo L, Xing Y, Xi T. MiR-424 regulates monocyticdifferentiation of human leukemia U937 cells by directly targeting CDX2.Biotechnol Lett. 2013 Nov;35(11):1799-806.51.Wang C, Ma W, Su Y. NF-κB pathway contr ibutes to cadmium-induced apoptosisof porcine granulosa cells. Biol Trace Elem Res. 2013 Jun;153(1-3):403-10.52.Gao C, Wang X, Chen L, Wang JH, Gao ZT, Wang H. Knockdown of Bcl-3 inhibitscell growth and induces DNA damage in HTLV-1-infected cells. Asian Pac J Cancer Prev. 2013;14(1):405-8.53.Yu ZY, Bai YN, Luo LX, Wu H, Zeng Y. Expression of microRNA-150 targetingvascular endothelial growth factor-A is downregulated under hypoxiaduring liver regeneration. Mol Med Rep. 2013 Jul;8(1):287-93.54.Yan B, Zhou Y, Feng S, Lv C, Xiu L, Zhang Y, Shi J, Li Y, Wei P, Qin Z. β-Elemene-Attenuated Tumor Angiogenesis by Targeting Notch-1 in Gastric Cancer Stem-Like Cells. Evid Based Complement Alternat Med. 2013;2013:268468.55.Li WL, Chen ML, Liu SS, Li GL, Gu TY, Liang P, Qin YM, Zhan YH, Quan Y,Zhang GH. Sweet Preference Modified by Early Experience in Mice and the Related Molecular Modulations on the PeripheralPathway. J Mol Neurosci. 2013 Sep;51(1):225-36.56.Wei Q, Shi F. Cleavage of poly (ADP-ribose) polymerase-1 is involved in theprocess of porcine ovarian follicular atresia. Anim Reprod Sci. 2013 May;138(3-4):282-91.57.Liu JH, Xu TT, Liu YJ, Zhu WY, Mao SY. A high-grain diet causes massivedisruption of ruminal epithelial tight junctions in goats. Am J Physiol Regul Integr Comp Physiol. 2013 Aug 1;305(3):R232-41.58.Tang H, Ge L, Shao W, Qiu Y, Cui D. Effect of targeted silencing of hTERT mRNAby lentivirus-mediated siRNA on A549 lung cancer cells in vitro. Mol Biol Rep.2013 Jan;40(1):605-16.59.Wang P, Wang X, Yang X, Liu Z, Wu M, Li G. Budesonide suppresses pulmonaryantibacterial host defense by down-regulating cathelicidin-related antimicrobial peptide in allergic inflammation mice and in lung epithelial cells. BMC Immunol.2013 Feb 6;14:7.60.Wang J, Wang X, Gao C, Song X, Niu Z, Gao Z, Qin Z, Chang J, Wang H. Thepyrimidine analog FNC inhibits cell proliferation and viral protein synthesis in HTLV‑1‑infected cells. Mol Med Rep. 2013 May;7(5):1656-60.61.Chen Z, Yang Y, Yang X, Zhou C, Li F, Lei P, Zhong L, Jin X, Peng G. Immuneeffects of optimized DNA vaccine and protective effects in a MPTP model of Parkinson's disease. Neurol Sci. 2013 Sep;34(9):1559-70.注:更多使用本产品的文献请参考产品网页。
Physica A 374(2007)483–490Identification of overlapping community structure in complexnetworks using fuzzy c -means clusteringShihua Zhang a,Ã,Rui-Sheng Wang b ,Xiang-Sun Zhang aa Academy of Mathematics &Systems Science,Chinese Academy of Science,Beijing 100080,Chinab School of Information,Renmin University of China,Beijing 100872,ChinaReceived 28June 2006Available online 7August 2006AbstractIdentification of (overlapping)communities/clusters in a complex network is a general problem in data mining of network data sets.In this paper,we devise a novel algorithm to identify overlapping communities in complex networks by the combination of a new modularity function based on generalizing NG’s Q function,an approximation mapping of network nodes into Euclidean space and fuzzy c -means clustering.Experimental results indicate that the new algorithm is efficient at detecting both good clusterings and the appropriate number of clusters.r 2006Elsevier B.V.All rights reserved.Keywords:Overlapping community structure;Modular function;Spectral mapping;Fuzzy c -means clustering;Complex network1.IntroductionLarge complex networks representing relationships among set of entities have been one of the focuses of interest of scientists in many fields in the recent years.Various complex network examples include social network,worldwide web network,telecommunication network and biological network.One of the key problems in the field is ‘How to describe/explain its community structure’.Generally,a community in a network is a subgraph whose nodes are densely connected within itself but sparsely connected with the rest of the network.Many studies have verified the community/modularity structure of various complex networks such as protein-protein interaction network,worldwide web network and co-author network.Clearly,the ability to detect community structure in a network has important practical applications and can help us understand the network system.Although the notion of community structure is straightforward,construction of an efficient algorithm for identification of the community structure in a complex network is highly nontrivial.A number of algorithms for detecting the communities have been developed in various fields (for a recent review see Ref.[1]and a recent comparison paper see Ref.[2]).There are two main difficulties in detecting community structure.The first is that we don’t know how many communities there are in a given network.The usual drawback in many /locate/physa0378-4371/$-see front matter r 2006Elsevier B.V.All rights reserved.doi:10.1016/j.physa.2006.07.023ÃCorresponding author.E-mail addresses:zsh@ (S.Zhang),wrs@ (R.-S.Wang),zxs@ (X.-S.Zhang).algorithms is that they cannot give a valid criterion for measuring the community structure.Secondly,it is a common case that some nodes in a network can belong to more than one community.This means the overlapping community structure in complex networks.Overlapping nodes may play a special role in a complex network system.Most known algorithms such as divisive algorithm [3–5]cannot detect them.Only a few community-detecting methods [6,7]can uncover the overlapping community structure.Taking into account the first difficulty,Newman and Girvan [8]has developed a new approach.They introduced a modularity function Q for measuring community structure.In order to write the context properly,we refer to a similar formulation in Ref.[5].In detail,given an undirected graph/network G ðV ;E ;W Þconsisting of the node set V ,the edge set E and a symmetric weight matrix W ¼½w ij n Ân ,where w ij X 0and n is the size of the network,the modularity function Q is defined asQ ðP k Þ¼X k c ¼1L ðV c ;V c ÞL ðV ;V ÞÀL ðV c ;V ÞL ðV ;V Þ 2"#,(1)where P k is a partition of the nodes into k groups and L ðV 0;V 00Þ¼P i 2V 0;j 2V 00w ði ;j Þ.The Q function measuresthe quality of a given community structure organization of a network and can be used to automatically select the optimal number of communities k according to the maximum Q value [8,5].The measure has been used for developing new detection algorithms such as Refs.[5,9,4].White and Smyth [5]showed that optimizing the Q function can be reformulated as a spectral relaxation problem and proposed two spectral clustering algorithms that seek to maximize Q .In this study,we develop an algorithm for detecting overlapping community structure.The algorithm combines the idea of modularity function Q [8],spectral relaxation [5]and fuzzy c -means clustering method[10]which is inspired by the general concept of fuzzy geometric clustering.The fuzzy clustering methods don’t employ hard assignment,while only assign a membership degree u ij to every node v i with respect to the cluster C j .2.MethodSimulation across a wide variety of simulated and real world networks showed that large Q values are correlated with better network clusterings [8].Then maximizing the Q function can obtain final ‘optimal’community structure.It is noted that in many complex networks,some nodes may belong to more than one community.The divisive algorithms based on maximizing the Q function fail to detect such case.Fig.1shows an example of a simple network which visually suggests three clusters and classifying node 5(or node 9)intoFig.1.An example of network showing Q and e Qvalues for different number k of clusters using the same spectral mapping but different cluster methods,i.e.k -means and fuzzy c -means,respectively.For the latter,it shows every node’s soft assignment and membership of final clusters with l ¼0:15.S.Zhang et al./Physica A 374(2007)483–490484two clusters at the same time may be more appropriate intuitively.So we introduce the concept of fuzzy membership degree to the network clustering problem in the following subsection.2.1.A new modular functionIf there are k communities in total,we define a corresponding n Âk ‘soft assignment’matrix U k ¼½u 1;...;u k with 0p u ic p 1for each c ¼1;...;k and P kc ¼1u ic ¼1for each i ¼1;...;n .With this we define the membership of each community as ¯V c ¼f i j u ic 4l ;i 2V g ,where l is a threshold that can convert a soft assignment into final clustering.We define a new modularity function e Q as e Q ðU k Þ¼X k c ¼1A ð¯V c ;¯V c ÞA ðV ;V ÞÀA ð¯V c ;V ÞA ðV ;V Þ 2"#,(2)where U k is a fuzzy partition of the vertices into k groups and A ð¯V c ;¯V c Þ¼P i 2¯V c ;j 2¯V c ððu ic þu jc Þ=2Þw ði ;j Þ,A ð¯V c ;V Þ¼A ð¯V c ;¯V c ÞþP i 2¯V c ;j 2V n ¯V c ððu ic þð1Àu jc ÞÞ=2Þw ði ;j Þand A ðV ;V Þ¼P i 2V ;j 2V w ði ;j Þ.This of coursecan be thought as a generalization of the Newman’s Q function.Our objective is to compute a soft assignment matrix by maximizing the new Q function with appropriate k .How could we do?2.2.Spectral mappingWhite and Smyth [5]showed that the problem of maximizing the modularity function Q can be reformulated as an eigenvector problem and devised two spectral clustering algorithms.Their algorithms are similar in spirit to a class of spectral clustering methods which map data points into Euclidean space by eigendecomposing a related matrix and then grouping them by general clustering methods such as k -means and hierarchical clustering [5,9].Given a network and its adjacent matrix A ¼ða ij Þn Ân and a diagonal matrix D ¼ðd ii Þ,d ii ¼P k a ik ,two matrices D À1=2AD À1=2and D À1A are often used.A recent modification [11]uses the top K eigenvectors of the generalized eigensystem Ax ¼tDx instead of the K eigenvectors of the two matrices mentioned above to form a matrix whose rows correspond to original data points.The authors show that after normalizing the rows using Euclidean norm,their eigenvectors are mathematically identical and emphasize that this is a numerically more stable method.Although their result is designed to cluster real-valued points[11,12],it is also appropriate for network clustering.So in this study,we compute the top k À1eigenvectors of the eigensystem to form a ðk À1Þ-dimensional embedding of the graph into Euclidean space and use ‘soft-assignment’geometric clustering on this embedding to generate a clustering U k (k is the expected number of clusters).2.3.Fuzzy c-meansHere,in order to realize our ‘soft assignment’,we introduce fuzzy c -means (FCM)clustering method [10,13]to cluster these points and maximize the e Qfunction.Fuzzy c -means is a method of clustering which allows one piece of data to belong to two or more clusters.This method (developed by Dunn in 1973[10]and improved by Bezdek in 1981[13])is frequently used in pattern recognition.It is based on minimization of the following objective functionJ m ¼Xn i ¼1X k j ¼1u m ij k x i Àc j k 2,(3)over variables u ij and c with P j u ij ¼1.m 2½1;1Þis a weight exponent controlling the degree of fuzzification.u ij is the membership degree of x i in the cluster j .x i is the i th d -dimensional measured data point.c j is the d -dimensional center of the cluster j ,and k Ãk is any norm expressing the similarity between any measured data and the center.Fuzzy partitioning is carried out through an iterative optimization of the objective function shown above,with the update of membership degree u ij and the cluster centers c j .This procedure converges to a local minimum or a saddle point of J m .S.Zhang et al./Physica A 374(2007)483–4904852.4.The flow of the algorithmGiven an upper bound K of the number of clusters and the adjacent matrix A ¼ða ij Þn Ân of a network.The detailed algorithm is stated straightforward for a given l as follows:Spectral mapping:(i)Compute the diagonal matrix D ¼ðd ii Þ,where d ii ¼P k a ik .(ii)Form the eigenvector matrix E K ¼½e 1;e 2;...;e K by computing the top K eigenvectors of thegeneralized eigensystem Ax ¼tDx .Fuzzy c -means:for each value of k ,2p k p K :(i)Form the matrix E k ¼½e 2;e 3;...;e k from the matrix E K .(ii)Normalize the rows of E k to unit length using Euclidean distance norm.(iii)Cluster the row vectors of E k using fuzzy c -means or any other fuzzy clustering method to obtain a softassignment matrix U k .Maximizing the modular function:Pick the k and the corresponding fuzzy partition U k that maximizes e QðU k Þ.In the algorithm above,we initialize FCM such that the starting centroids are chosen to be as orthogonal as possible which is suggested for k -means clustering method in Ref.[12].The initialization does not change the time complexity,and also can improve the quality of the clusterings,thus at the same time reduces the need for restarting the random initialization process.The framework of our algorithm is similar to several spectral clustering methods in previous studies[5,9,12,11].We also map data points (work nodes in our study)into Euclidean space by computing the top K eigenvectors of a generalized eigen system and then cluster the embedding using a fuzzy clustering method just as others using geometric clustering algorithm or general hierarchical clustering algorithm.Here,we emphasize two key points different from those earlier studies:We introduce a generalized modular function e Q employing fuzzy concept,which is devised for evaluating the goodness of overlapping community structure. In combination with the novel e Qfunction,we introduce fuzzy clustering method into network clustering instead of general hard clustering methods.This means that our algorithm can uncover overlapping clusters,whereas general framework:‘‘Objective function such as Q function and Normalized cut function+Spectral mapping+general geometric clustering/hierarchical clustering’’cannot achieve this.3.Experimental resultsWe have implemented the proposed algorithm by Matlab.And the fuzzy clustering toolbox [14]is used for our analysis.In order to make an intuitive comparison,we also compute the hard clustering based on the original Q -function,spectral mapping (same as we used)and k -means clustering.We illustrate the fuzzy concept and the difference of our method with traditional divisive algorithms by a simple example shown in Fig.1.Just as mentioned above,the network visually suggests three clusters.But classifying node 5(or node 9)simultaneously into two clusters may be more reasonable.We can see from Fig.1that our method did uncover the overlapping communities for this simple network,while the traditional method can only make one node belong to a single cluster.We also present the analysis of two real networks,i.e.the Zachary’s karate club network and the American college football team network for better understanding the differences between our method and traditional methods.S.Zhang et al./Physica A 374(2007)483–490486S.Zhang et al./Physica A374(2007)483–490487 3.1.Zachary’s karate clubThe famous karate club network analyzed by Zachary[15]is widely used as a test example for methods of detecting communities in complex networks[1,8,16,3,4,17,9,18,19].The network consists of34members of a karate club as nodes and78edges representing friendship between members of the club which was observed over a period of two years.Due to a disagreement between the club’s administrator and the club’s instructor, the club split into two smaller ones.The question we concern is that if we can uncover the potential behavior of the network,detect the two communities or multiple groups,and particularly identify which community a node belongs to.The network is presented in Fig.2,where the squares and the circles label the members of the two groups.The results of k-means and our analysis are illustrated in Fig.3.The k-means combined with Q function divides the network into three parts(see in Fig.3A),but we can see that some nodes in one cluster are also connected densely with another cluster such as node9and31in cluster 1densely connecting with cluster2,and node1in cluster2with cluster3.Fig.3B shows the results of our method,from which we can see that node1,9,10,31belong to two clusters at the same time.These nodes in the network link evenly with two clusters.Another thing is that the two methods both uncover three communities but not two.There is a small community included in the instructor’s faction,since the set of nodes5,6,7,11,17only connects with node1in the instructor’s faction.Note that our method also classifies node1into the small community,while k-means does not.work of American college football teamsThe second network we have investigated is the college football network which represents the game schedule of the2000season of Division I of the US college football league.The nodes in the network represent the115 teams,while the links represent613games played in the course of the year.The teams are divided into conferences of8–12teams each and generally games are more frequent between members of the same conference than between teams of different conferences.The natural community structure in the network makes it a commonly used workbench for community-detecting algorithm testing[3,5,7].Fig.4shows how the modularity Q and e Q vary with k with respect to k-means and our method,respectively. The peak for k-means is at k¼12,Q¼0:5398,while for our algorithm at k¼10,e Q¼0:4673with l¼0:10. Both methods identify ten communities which contain ten conferences almost exactly.Only teams labeled as Sunbelt are not recognized as belonging to a same community for both methods.This group is classified as well in the results of Refs.[3,19].This happens because the Sunbelt teams played nearly as many games against Western Athletic teams as they played in their own conference,and they also played quite a number of gamesagainst Mid-American team.Our method identified11nodes(teams)which belong to at least twoFig.2.Zachary’s karate club network.Square nodes and circle nodes represent the instructor’s faction and the administrator’s faction, respectively.Thisfigure is from Newman and Girvan[8].communities (see Fig.5,11red nodes).These nodes generally connect evenly with more than one community,so we cannot classify them into one specific community correctly.These nodes represent ‘fuzzy’points which cannot be classified correctly by employing current link information.Maybe such points play a ‘bridge’role in two or more communities in complex network of other types.4.Conclusion and discussionIn this paper,we present a new method to identify the community structure in complex networks with a fuzzy concept.The method combines a generalized modularity function,spectral mapping,and fuzzy clustering technique.The nodes of the network are projected into d -dimensional Euclidean space which is obtained by computing the top d nontrivial eigenvectors of the generalized eigensystem Ax ¼tDx .Then the fuzzy c -means clustering method is introduced into the d -dimensional space based on general Euclidean distance to cluster the data points.By maximizing the generalized modular function e QðU d Þfor varying d ,we obtain the appropriate number of clusters.The final soft assignment matrix determines the final clusters’membership with a designated threshold l .Fig.3.The results of both k -means and our method applied to karate club network.A:The different colors represent three different communities obtained by k -means and the right table shows values of NG’Q versus different k .B:Four red nodes represent the overlap of two adjacent communities obtained by our method and the right table shows values of new Q versus different k with l ¼0:25.3.t y 0510152000.10.20.30.40.50.6k-meansK N G ' Q 0510152000.10.20.30.40.5fuzzy c-means K N e w Q Fig.4.Q and e Qvalues versus k with respect to k -means and fuzzy c -means clustering methods for the network of American college football team.S.Zhang et al./Physica A 374(2007)483–490488Although spectral mapping has been comprehensively used before to detect communities in complex networks (even in clustering the real-valued points),we believe that our method represents a step forward in this field.A fuzzy method is introduced naturally with the generalized modular function and fuzzy c -means clustering technique.As our tests have suggested,it is very natural that some nodes should belong to more than one community.These nodes may play a special role in a complex network system.For example,in a biological network such as protein interaction network,one node (protein or gene)belonging to two functional modules may act as a bridge between them which transfers biological information or acts as multiple functional units [6].One thing should be noted is that when this method is applied to large complex networks,computational complexity is a key problem.Fortunately,some fast techniques for solving eigensystem have been developed[20]and several methods of FCM acceleration can also be found in the literature [21].For instance,if we adopt the implicitly restarted Lanczos method (IRLM)[20]to compute the K À1eigenvectors and the efficient implementation of the FCM algorithm in Ref.[21],we can have the worse-case complexity of O ðmKh þnK 2h þK 3h Þand O ðnK 2Þ,respectively,where m is the number of edges in the network and h is the number of iteration required until convergence.For large sparse networks where m $n ,and K 5n ,the algorithms will scale roughly linearly as a function of the number of nodes n .Nonetheless,the eigenvector computation is still the most computationally expensive step of the method.We expect that this new method will be employed with promising results in the detection of communities in complex networks.AcknowledgmentsThis work is partly supported by Important Research Direction Project of CAS ‘‘Some Important Problem in Bioinformatics’’,National Natural Science Foundation of China under Grant No.10471141.The authors thank Professor M.E.J.Newman for providing the data of karate club network and the college football team network.Fig.5.Fuzzy communities of American college football team network (k ¼10and e Q¼0:4673)with given l ¼0:10(best viewed in color).S.Zhang et al./Physica A 374(2007)483–490489References[1]M.E.J.Newman,Detecting community structure in networks,Eur.Phys.J.B 38(2004)321–330.[2]L.Danon,J.Duch,A.Diaz-Guilera,A.Arenas,Comparing community structure identification,J.Stat.Mech.P09008(2005).[3]M.Girvan,M.E.J.Newman,Community structure in social and biological networks,A 99(12)(2002)7821–7826.[4]J.Duch,A.Arenas,Community detection in complex networks using extremal optimization,Phys.Rev.E 72(2005)027104.[5]S.White,P.Smyth,A spectral clustering approach to finding communities in graphs,SIAM International Conference on DataMining,2005.[6]G.Palla,I.Derenyi,I.Farkas,T.Vicsek,Uncovering the overlapping community structure of complex networks in nature and society,Nature 435(2005)814–818.[7]J.Reichardt,S.Bornholdt,Detecting fuzzy community structures in complex networks with a Potts model,Phys.Rev.Lett.93(2004)218701.[8]M.E.J.Newman,M.Girvan,Finding and evaluating community structure in networks,Phys.Rev.E 69(2004)026113.[9]L.Donetti,M.A.Mun oz,Detecting network communities:a new systematic and efficient algorithm,J.Stat.Mech.P10012(2004).[10]J.C.Dunn,A fuzzy relative of the ISODATA process and its use in detecting compact well-separated clusters,J.Cybernet.3(1973)32–57.[11]D.Verma,M.Meila,A comparison of spectral clustering algorithms.Technical Report,2003,UW CSE Technical Report 03-05-01.[12]A.Ng,M.Jordan,Y.Weiss,On spectral clustering:analysis and an algorithm,Adv.Neural Inf.Process.Systems 14(2002)849–856.[13]J.C.Bezdek,Pattern Recognition with Fuzzy Objective Function Algorithms,Plenum Press,New York,1981.[14]Fuzzy Clustering Toolbox-h http://www.fmt.vein.hu/softcomp/fclusttoolbox/i .[15]W.W.Zachary,An information flow model for conflict and fission in small groups,J.Anthropol.Res.33(1977)452–473.[16]M.E.J.Newman,Fast algorithm for detecting community structure in networks,Phys.Rev.E 69(2004)066133.[17]F.Radicchi,C.Castellano,F.Cecconi,V.Loreto,D.Parisi,Defining and identifying communities in networks,Proc.Natl.Acad.A 101(9)(2004)2658–2663.[18]F.Wu,B.A.Huberman,Finding communities in linear time:a physics approach,Eur.Phys.J.B 38(2004)331–338.[19]S.Fortunato,tora,M.Marchiori,A method to find community structures based on information centrality,Phys.Rev.E 70(2004)056104.[20]Z.Bai,J.Demmel,J.Dongarra,A.Ruhe,H.Vorst (Eds.),Templates for the Solution of Algebraic Eigenvalue Problems:A PracticalGuide,SIAM,Philadelphia,PA,2000.[21]J.F.Kelen,T.Hutcheson,Reducing the time complexity of the fuzzy c -means algorithm,IEEE Trans.Fuzzy Systems 10(2)(2002)263–267.S.Zhang et al./Physica A 374(2007)483–490490。
Proceedings of the New3SC International Conference: to be published in Physica CInhomogeneous Charge State in HTSC Cuprates and CMR Manganites T. EgamiDepartment of Materials Science and Engineering and Laboratory for Research on the Structure of Matter, University of Pennsylvania, Philadelphia, PA 19104(Received: Jan. 15, 2001, revised: Feb. 14, 2001)Abstract: Recent measurements of neutron elastic and inelastic scattering suggest that the charge states are spatially inhomogeneous at two lengthscales, atomic and nanometer scales, in both the high-temperature superconducting (HTSC) cuprates and colossal mag-netoresistive (CMR) manganites. We suggest that the two-phonon mechanism that con-trols the charge localization in CMR manganites is also at work in HTSC cuprates, and may hold a key to understanding the mechanism of superconductivity.Key words: neutron scattering, charge inhomogeneity, phonons, CMR manganitese-mail: egami@1. IntroductionGenerally local lattice distortion and charge inhomogeneity are considered to be harmful to superconductivity, including high-temperature superconductivity (HTSC). Thus in spite of numerous reports of local lattice distortions in cuprates [1] the majority view has been that they are unimportant to the mechanism of superconductivity, and are merely the consequences of strong interactions and short coherence length. However, the perception appears to be changing due to the observation of the spin-charge stripes in non-superconducting cuprates with 1/8 charge density and the conjecture that dynamic stripes might be present in the superconducting state [2]. In this paper we discuss recent results by neutron inelastic scattering that indicate the presence of charge inhomogemei-ties at two lengthscales, atomic and nanometer scales, which are in many ways similar to those in the manganites that show colossal magnetoresistance (CMR). We suggest that the two-phonon mechanism of charge localization in CMR manganites must be in opera-tion also in cuprates, and such phonon involvement may be an integral part of the HTSC phenomenon.2. Inelastic Neutron ScatteringWe have made a series of inelastic neutron scattering measurements on the high-energy LO phonons in cuprates, La1.85Sr0.15CuO4 [3] and YBa2Cu3O6+x [4,5]. The pho-non branch we studied was the Cu-O bond-stretching mode as shown in Fig. 1. It mainly involves oxygen motion, and at the zone-center, q = 0, it is the ferroelectric mode, while at the zone-edge, q = (π, 0, 0) in the unit of 1/a (a is the Cu-Cu distance), it is a half-breathing mode. This mode induces the hybridization between the Cu-d and O-p orbitals to change and results in charge transfer. For this reason it is expected to interact strongly with charge carriers [6,7]. Indeed this is the only mode that shows strong softening near the zone-edge upon doping [8], and the kink in the electron dispersion determined by photoemission is likely due to this phonon mode [9]. The neutron results revealed a rather anomalous nature of this branch:1. At low temperatures the branch splits into two, a high-energy branch and a low en-ergy-branch, both of which show little dispersion.2. The splitting of the phonon dispersion becomes weak at higher temperatures. InYBa2Cu3O6.95 the temperature dependence of the change in the dispersion agrees with the superconducting order parameter.3. As the doped charge density is increased in YBCO, one would expect the LO phonondispersion to soften gradually near the zone-edge. Surprisingly the energy at the zone-edge does not change much with doping and the phonon dispersion is always split into two. But the intensity is transferred from the high energy-branch to the low-energy branch with increased doping.The point 1 is best explained by doubling of the unit cell that halves the Brillouin zone, such as the charge pattern shown in Fig. 2 that has the periodicity of 2a [3]. The usual stripe pattern with the periodicity of 4a fails to account for the dispersion. The 2a perio-dicity appears to reflect the underlying Peierls instability, and may give rise to the cou-pled ladder state [10]. As point 2 shows this charge pattern is temperature dependent. The point 3 is best interpreted in terms of an inhomogeneous structure; microscopic seg-regation into charge-rich and charge-poor regions. The high-energy branch must be as-sociated with charge-poor regions, since it agrees with the dispersion in the undoped sample [8], and the low-energy branch with charge-rich regions, since its intensity grows with doping. The lengthscale for the segregation must be a nanometer at most; otherwise the q dependence of the phonon structure factors for the two branches must be similar. This micro-phase segregation was seen also by a neutron scattering PDF study [11]. Thus the results suggest that the charge state of cuprates is inhomogeneous at two lengthscales, atomic and nanometer scales.3. Polaron stability in CMR manganitesIn manganites charges can become localized forming polarons, depending on tem-perature and composition [12,13]. Nanometer-scale charge inhomogeneity is also seen, even in the metallic state just above the insulator-to-metal transition, suggesting the charge segregation tendency [14,15]. There is a striking parallelism between polaron formation in manganites (0-d) and stripe formation in cuprates (1-d) as the mechanism of charge localization, and the presence of nanometer charge inhomogeneity in both. It is well understood that, in manganites, charge localization is controlled by the balance be-tween the localization forces (electron-lattice coupling and antiferromagnetic or orbital ordering) and the delocalizaton forces (kinetic energy and elastic energy) [12,13]. These forces can be affected by the ionic size of the A-site ions in the perovskite, <r A> [16].While the standard explanation is that the ionic size modulates the bandwidth, this has been questioned by a number of results [15]. Our alternate interpretation is based upon the local structure of polarons [15,17]. In the case of simple lattice polarons the spatial extensions of the charge and the lattice distortion are about the same. However, in man-ganites and cuprates spins play a large role in confining the charge, so that a hole can be localized on a single Mn site, while the elastic field that a polaron generates is long-range. This long-range stress field increases the polaron self-energy by a factor of about 2, according to the continuum elasticity theory [18,19]. But if the value of <r A> is small enough the Mn-O-Mn bonds are buckled. Then, when the local Jahn-Teller distortion is removed within the polaron and the long Mn-O bond locally becomes short, this change can be accommodated by unbuckling the bond (Fig. 3) [17]. In other words this unbuck-ling screens the local polaronic strain, and the long-range field will not be generated. Thus the renormalization of the polaron energy due to the long-range field does not hap-pen, reducing the polaron energy by a factor of 2 compared to the previous case. In man-ganites this change in the polaron self-energy apparently is sufficient to stabilize polaron when <r A> is small, and destabilize when <r A> is large. At the crossover polarons are marginally stable, giving rise to the CMR phenomenon.4. Electron-phonon coupling in cupratesIn our view an exact parallel should be observed for cuprates [20]. An increase in the local hole density will cause the Cu-O bond to be reduced in length, or in the phonon language, excite the in-plane LO phonons. If the Cu-O-Cu bond is buckled, the in-plane LO phonons couple to the c-axis oxygen mode, as in Fig. 3. The coupling is proportional to the deviation of the Cu-O-Cu bond angle from π. This coupling will soften the in-plane mode, by the screening mechanism described above. There are two possibilities of the in-plane modes that are softened by this coupling. The first obvious choice is the stripe mode (q = π/2), since when the buckling is strong charges will be localized into the stripe state, driving the system insulating. This is consistent with the effect of the CuO6 tilting [21], and the need of Nd in (La,Nd)1.875Sr0.125CuO4 to produce the stripe structure [2], since Nd3+ is smaller than La3+ (1.27 vs. 1.36 Å [22]). The second possibility is the Peierls mode (q = π). Our neutron results strongly suggest that this mode is present, per-haps competing against the stripe mode and consequently not developing into the stable state. This could produce a local spin-singlet state, suppress the antiferromagnetic order, and promote superconductivity [10]. At the same time the two-phonon mixing could produce multi-band HTSC [23,24]. While this scenario is highly speculative at this mo-ment, it explains quite well the anomaly in the pulsed neutron PDF at T C involving the out-of-plane oxygen motion [25], and the dependence of T C on the “micro-strain” pro-posed recently by Bianconi et al. [26]. Indeed a remarkable similarity between the phase diagram in Fig. 6 of Ref. 26 and that for manganites (Fig. 4) [15,17] suggests that similar mechanisms must be at work. In Ref. 26, however, the stripe mode itself is considered to generate HTSC. In the present scenario the stripe mode merely terminates HTSC by CDW formation, while the competition between the Peierls state and the stripe state plays a crucial role in HTSC. Further experimental works are being conducted to examine these scenarios.5. ConclusionsIn CMR manganites the stability of polarons is determined by the local structure, i.e. screening of the local in-plane LO phonons by the c-axis oxygen phonons. An exact par-allel should hold in cuprates where the coupling between LO phonons and c-axis phonons could play a crucial role in determining the stability of the stripe structure, and possibly superconductivity.Acknowledgments: The author is grateful to the collaborators of the phonon and CMR projects, R. J. McQueeney, M. Yethiraj, D. Louca, Y. Petrov, M. Arai, J. F. Mitchell, H.A. Mook, G. Shirane and Y. Endoh. He is also thankful to A. Bianconi, S. J. L. Billinge, A. R. Bishop, A. Bussmann-Holder, J.B. Goodenough, L. P. Gor’kov, V. Kresin, A. Lanzara, K. A. Müller, J.C. Phillips, S. Sachdev, Z.-X. Shen, S. R. Shenoy, and M. Tachiki for useful discussions. Research at the University of Pennsylvania was supported by the National Science Foundation through DMR96-28136.References:1. T. Egami and S. J. L. Billinge, in Physical Properties of High Temperature Super-conductors V, ed. D. Ginsberg (Singapore, World Scientific, 1996) p. 265.2. J. M. Tranquada et al., Nature375, 561 (1995).3. R. J. McQueeney, Y. Petrov, T. Egami, M. Yethiraj, G. Shirane and Y. Endoh, Phys.Rev. Lett. 82, 628 (1999).4. T. Egami, R. J. McQueeney, Y. Petrov, M. Yethiraj, G. Shirane, and Y. Endoh, AIPConf. Proc. 483, 231 (1999).5. Y. Petrov, T. Egami, R. J. McQueeney, M. Yethiraj, H. A. Mook, and F. Dogan,cond-mat/0003414.6. S. Ishihara, T. Egami and M. Tachiki, Phys. Rev. B55, 3163 (1997).7. Y. Petrov and T. Egami, Phys. Rev. B58, 9485 (1998).8. L. Pintschovius, et al., Physica C185-189, 156 (1991).9. A. Lanzara, et al., cond-mat/0102227.10. S. Sachdev, Science288, 475 (2000); S. Sachdev and N. Read, Int. J. Mod. Phys. B5, 219 (1991).11. E. S. Bozin, G. H. Kwei, H. Takagi and S. J. L. Billinge, Phys. Rev. Lett. 84, 5856(2000).12. A. J. Mills, P. B. Littlewood, and B. I. Shraiman, Phys. Rev. Lett. 74, 5144 (1995).13. H. Röder, J. Zang and A. R. Bishop, Phys. Rev. Lett., 76, 1356 (1996).14. D. Louca, T. Egami, E. L. Brosha, H. Röder, and A. R. Bishop, Phys. Rev. B56,R8475 (1997).15. T. Egami and D. Louca, J. Superconductivity13, 247 (2000).16. H. Y. Hwang, S.-W. Cheong, P. G. Radaelli, M. Marezio, B. Batlogg, Phys. Rev.Lett.75, 914 (1995).17. T. Egami and D. Louca, J. Superconductivity12, 23 (1999).18. J. D. Eshelby, Proc. Roy. Soc. A241, 376 (1957).19. T. Egami and D. Louca, unpublished.20. T. Egami, AIP Conf. Proc., in press.21. B. Büchner, et al., Phys. Rev. Lett.73, 1841 (1994).22. R. Shannon, Acta Cryst. A32, 751 (1976).23. S. R. Shenoy, V. Subrahmanyam and A. R. Bishop, Phys. Rev. Lett.79, 4657 (1997).24. A. Bussmann-Holder, K. A. Müller, R. Micnas, H. Büttner, A. Simon, A. R. Bishopand T. Egami, submitted.25. B. H. Toby, T. Egami, J. D. Jorgensen and M. A. Subramanian, Phys. Rev. Lett. 64,2414 (1990).26. A. Bianconi, G. Bianconi, S. Caprara, D. Di Castro, H. Oyanagi and N. L. Saini, J.Phys.: Condens. Matter12, 10655 (2000).Figure captions:Figure 1. (a) High-energy LO phonons at the zone-center, and (b) at the zone-edge, (π, 0, 0).Figure 2. A possible charge pattern with the 2a periodicity [3]. Large circles denote oxygen, with different charge densities indicated by grades of lightness. For clarity we show the pattern with a long-range order, but in reality the pattern should be dynamic and short-range, with the correlation length about 20 × 8 Å.Figure 3. Accommodation of the reduction in the local Mn-O distance due to the presence of a hole on Mn, (a) when the Mn-O-Mn bond is straight this causes the neighboring Mn-O bonds to stretch, (b) when the Mn-O-Mn bond is buckled, unbuckling will locally accommodate the change in the Mn-O distance, without producing a long-range stress field [13,15].Figure 4. Phase diagram for the CMR phenomenon in A1-x A’x MnO3, plotted for the ionic radius (9-coordinated), <r A>, against the charge density, x [13,15]. When <r A> is greater than 1.24 Å polarons are not stable, resulting in the metallic state. When <r A> is less than 1.18 Å polarons are stable, making the system insulating. The CMR phenome-non is observed in the crossover region.(b)Fig. 1Fig. 2(a) (b) Fig. 3Fig. 4。
第 22卷第 10期2023年 10月Vol.22 No.10Oct.2023软件导刊Software Guide一种基于物体追踪的改进语义SLAM算法杜小双,施展,华云松(上海理工大学光电信息与计算机工程学院,上海 200093)摘要:在视觉同步定位与建图(SLAM)算法中,使用语义分割和目标检测以剔除异常点的方法成为主流,但使用中无法对物体语义信息进行充分追踪。
为此,提出一种基于物体追踪的改进语义SLAM算法,通过YOLACT++网络分割物体掩码,提取物体特征点后,利用帧间匹配实现物体追踪。
该方法对匹配特征点进行深度、重投影误差和极线约束三重检测后判断物体动静态,实现物体追踪并判断运动状态。
通过对TUM RGB-D数据集测试,实验表明该方法可有效追踪物体,且轨迹估计精度优于其他SLAM算法,具有较好实用价值。
关键词:视觉SLAM;语义分割;物体追踪;动态场景;几何约束DOI:10.11907/rjdk.222298开放科学(资源服务)标识码(OSID):中图分类号:TP301 文献标识码:A文章编号:1672-7800(2023)010-0205-06 An Improved Semantic SLAM Algorithm Based on Object TrackingDU Xiaoshuang, SHI Zhan, HUA Yunsong(School of Optical-Electrical and Computer Engineering,University of Shanghai for Science and Technology,Shanghai 200093,China)Abstract:In the visual SLAM (simultaneous localization and mapping), the method of using semantic segmentation and object detection to detect dynamic objects and remove outliers has become the mainstream, but its disadvantage is that it is unable to fully track the semantic in⁃formation of objects. Therefore, this paper proposes an improved semantic SLAM algorithm based on object tracking, which uses YOLACT++ network to segment object mask, extract object feature points, and use inter frame matching to achieve object tracking. The method detects the depth, reprojection error and epipolar constraint of the matched feature points, and then judges the dynamic and static state of the object to achieve object tracking and judge the motion state. After testing the TUM RGB-D dataset, the experiment shows that the method can effective⁃ly track objects, and the trajectory estimation accuracy is better than other SLAM algorithms, which has practical value.Key Words:visual SLAM; semantic segmentation; object tracking; dynamic environment; geometric constraint0 引言随着机器人技术、无人驾驶、增强现实等领域的发展与普及,视觉SLAM作为其应用的基础技术之一[1],得到了学者们的广泛关注与研究,并成为机器人定位与建图研究领域的一个热点[2]。
§5 超导体的电磁性质1本节主要内容: 1. 超导体特性之一:零电阻 2. 超导体特性之二:完全抗磁性(Meissner 迈斯纳效应) 3. 超导体的电动力学性质 4. 超导环的磁通俘获和磁通量子化现象2气体液化与低温环境的获得 1892年,发明了杜瓦瓶(中间抽真空,内胆涂有银 的双层玻璃瓶) 1899年,杜瓦(James Dewar)在伦敦皇家研究所成 立100周年庆典上,展示氢气(H2)的液化实验3水银超导体的发现Heike Kamerlingh Onnes (1853-1926) Dutch Physicist、 1882年,进入Leiden大学,研 究低温气体; 1908年,将液体的温度降低到 大约1K,成功将氦气液化; 1911年,开始研究金属在极 低温下的性质; 1912年,发现了水银的超导 电性, 1913年,获Nobel奖4The discovery of superconductivityNotebook 56, 8 April, 1911 Notebook 57, 26 October, 1911“Mercury[‘s resistance] practically zero [at 3 K] ……repeated with gold…”2014/11/5The historic plot. Superconducting transition at 4.2k in mercury5Meissner effectFritz Walther Meissner (1882-1974) 1933 Robert Ochsenfeld (1901-1993)German physicists2014/11/5Perfect diamagnetism below Tc6Londons’ theoryHeinz Fritz Wolfgang London London (1907-1970) (1900-1954) Londons’ Equation: (1935)Ampère's law:German Physicists2014/11/5 7Ginzburg-Landu theoryLev Landau (1908-1968) Vitaly Ginzburg (1916-2009) 1950 The free energy density:Complex order parameterU(1) gauge symmetry broken Soviet physicists2014/11/5其它几种超导体 元素 Al(铝) In(铟) Sn(锡) Pb(铅) Nb(铌) 1911 超导转变温度 1.2 K 3.4 K 3.7 K 7.2 K 9.2 K 198691986年,Muller和Bednorz发现:陶瓷氧化物 LaBaCuO的转变温度可达到35K。
cell theory细胞学说:一切生物从单细胞到高等动、植物都是由细胞组成;细胞是生物形态结构和功能活动的基本单位。
cytology细胞学:是研究细胞生命现象的科学。
其研究范围包括:细胞的形态结构和功能、分裂和分化、遗传和变异以及衰老和病变等。
cell biology细胞生物学:应用现代物理、化学、实验生物学及分子生物学技术和方法, 从细胞整体层次、亚细胞层次和分子层次三个层面、并将其有机地结合起来去研究探索细胞结构及其基本生命活动。
Ribozyme核酶:具有催化功能的RNA分子。
它具有高度专一内切核酸酶的活性。
miRNA小RNA:是广泛存在于真核生物中的一组短小的、不编码蛋白质的RNA家族,它们是由21-25个核苷酸组成的单链RNA。
unit membrane单位膜:电镜下,膜有三层结构组成,两层深色带,中间隔有一层疏松的浅色带,此三层结构作为一种单位,称为单位膜。
amphipathic molecule兼性分子或双亲媒性分子:即含有亲水的头又含有疏水的尾. 其头部具有亲水性,尾部疏水且寻求与其他疏水分子相聚.Fluid mosaic model液态镶嵌模型:生物膜是球形蛋白质和脂类分子呈二维状态排列的液态体,不是静止状态,而是具有流动性特点的结构. 脂类双分子层既具有液态分子的流动性、又具有固体分子排列的有序性,即流动的脂类双分子层构成的连续性膜性结构,各种球状蛋白质镶嵌于脂类双分子层中. 蛋白质分子的非极性部分嵌入脂类双分子层的疏水区;极性部分露于膜表面,似一群岛屿一样,无规则地分散在脂类的海洋中.Lipid Rafts Model 脂筏模型:即在生物膜的脂双分子层的外层有很多胆固醇和鞘磷脂富集而成的有序的、脂相的微结构域,它是一种动态结构,位于质膜的外层. 这些区域结构致密,介于无序液体与液晶之间(称为有序液体),是鞘磷脂与胆固醇的动态集合,如同“脂筏”一样,不同脂筏载有特定的蛋白,脂筏可根据胞内外不同刺激而改变自身大小和组成,以利于与另一特定蛋白相互作用,激活相应信号开关.脂筏间可彼此合并,使信号放大,脂筏就像蛋白质的停泊平台,与膜的信号转导、蛋白质分选均有密切的关系。
A fair days workA reactionA type personalityAAAAAabacusabandonment of aged abandonment reaction abandonment threat abarognosiaabaseabasementabasiaAbbe condenserabbe effectabbreviated drawing abbreviated formulaABC of reaction time ABC ABC of sense impression abderiteabdomenabdominal breath abdominal ganglion abdominal reflex abdominal wallabducent nerveabductionabductorabecedarianabepithymiaaberglaubeaberranceaberrantaberrant behavioraberrationaberration of lightaberration rateaberrometerabhorabhorrenceabhorrent to reasonabhorrerabient behaviorabient driveabient responseabilityability characteristicsability compensationability developmentability distributionability groupability groupingability levelability of cognizanceability of imaginationability of managersability of memoryability of taking criminal responsibility ability of thinkingability ratingability taskability testability trait abiogenesisabiogenistabionergy abiophysiologyabiosisabioticabiotic surround abirritationAbklingenablationablazeablepsiaablushablutomaniaabmodalityabnegateabnegationabnerval current Abney effectabnormalabnormal behavior abnormal child abnormal condition abnormal curve abnormal distribution abnormal endabnormal fertilization abnormal fixation abnormal habitude abnormal impulse workabnormal metabolism abnormal mindabnormal mind of convict abnormal personality abnormal phenomena abnormal psychology abnormal results abnormal sensation abnormal series abnormal sexual behavior abnormal sexuality abnormal spoilage abnormal time abnormalityabnormityabolition of reflex abominateabominationabominatorabortabortionabortiveabouliaaboulomaniaabove threshold abreactionabridged learning abridgment of response abridgment of table abrogateabruptionabscissaabscissa axisabscissionabsenceabsence of mindabsence of restrictionsabsent behaviorabsent responseabsenteeismabsentiaabsentia epilepticabsolute activityabsolute advantagesabsolute and relativeabsolute atmosphereabsolute beingabsolute brilliance limen absolute certitudeabsolute chanceabsolute choiceabsolute cognitive consciousness absolute conceptabsolute constructionabsolute differenceabsolute discriminationabsolute dispersionabsolute diversityabsolute egoabsolute electrostatic unit absolute equalityabsolute errorabsolute evaluation absolute existence absolute extremesabsolute fieldabsolute freedomabsolute grading standards absolute humidityabsolute ideaabsolute identityabsolute idiocyabsolute impression absolute intelligence absolute judgmentabsolute knowledge absolute lethal genes absolute light threshold absolute limenabsolute localization absolute luminosity absolute magnitude absolute mass unit absolute measurement absolute memoryabsolute mindabsolute moralityabsolute motion parallax absolute nothingabsolute notionabsolute numberabsolute pitchaposterioriAQCommunication capability index hypererosiaOAPophidiophiliaophidiophobeophidiophobiaophtalmodonesis ophthalmencephalon ophthalmocopia ophthalmodiastimeter ophthalmograph ophthalmography ophthalmologicophthalmologic examination ophthalmologyophthalmometer ophthalmophantom ophthalmoreaction ophthalmoscope ophthalmoscopy ophthalmostasisophthalmostat ophthalmotonometer ophthalmotrope ophthalmotropometer ophthalmotropometryopiateopiate receptoropinionopinion leaderopinion pollopinion pollingopinion surveyopinion testopinionableopinionaireopinionalopinionateopinionationOpinionation Scaleopinionativeopinionatoropiomaniaopiomaniacopisthencephalonopisthoporeiaopisthotonosopiumopium denopligationOppel Kundt figureOppel's illusionopponent color theoryopponent process theoryopponents theory of color vision opportune confabulationopportune moment for interrogation opportunityopportunity classopportunity factoropportunity to performopposites test opposites viewpoint oppositionopsialgia opsigenesopsinopsiometer opsomaniaopsonic action optesthesia opthalmusopticoptic angleoptic aphasiaoptic apraxiaoptic atrophyoptic axisoptic canaloptic celloptic centreoptic chiasmaoptic cupoptic discoptic glandoptic grooveoptic lobeoptic nerveoptic nerve fiber optic nystagmus optic papillaoptic radiationoptic segmentoptic stalkoptic thalamusoptic tractoptic typeoptic zoneoptica axisopticaloptical aberrationoptical activityoptical aidoptical alexiaoptical allachesthesiaoptical antimeroptical artoptical axisoptical centeroptical character recognition optical communicationoptical computeroptical correctionoptical dispersionoptical frequencyoptical illusionoptical imageoptical instrumentoptical isomerismoptical landing systemoptical lanternoptical measurementoptical mechanical systemoptical microscopeoptical parallaxoptical phenomenaoptical rangeoptical rangefinderoptical readeroptical revolving power optical rotationoptical systemoptically activeopticociliaryopticokineticopticokinetic reflex opticokineticsoptimaloptimal allocationoptimal arousal level optimal arousal level theory optimal choiceoptimal coloroptimal conditionoptimal controloptimal control model optimal decisionoptimal decodingoptimal functioningoptimal groupingoptimal illuminationoptimal movementoptimal performanceoptimal personality traits, optimal propertyoptimal proportionoptimal ratiooptimal samplingoptimal sampling plan optimal scalingoptimal seeking method optimal sequence trainer optimal state of performance optimal stimulationoptimal strategyoptimal work spaceoptimeteroptimismoptimizationoptimizeoptimizing modeloptimumoptimum cd ratiooptimum conditionoptimum controloptimum decision making optimum degreeoptimum estimateoptimum loadoptimum rangeoptimum selectionoptimum speedoptimum stateoptimum state of performance optimum strategyoptimum temperatureoption experimentoptionaloptional altruistic suicide optional routeoptional variantoptionsoptoblastoptogramoptogravit illusion optogyral illusion optokinetic response optomeninxoptometeroptometryoptomotoroptomyometeroptophoneORoraora serrataoracyoral aggressionoral anxietyoral arithmeticoral behaviororal cavityoral characteroral communicationoral complexoral contraceptive oral dependence oral driveoral erotismoral examination oral fixationoral incorporation oral languageoral literature oral methodoral optimismoral orientation oral personality oral pessimismoral phaseoral primacyoral readingoral reportoral sadismoral sensationoral sensory stage oral stageoral temperature oral testoralismorality orangoutang orangutanorbOrbeli phenomenonOrbison illusionorbitaorbitofrontal cortex orderorder codeorder disorder transition order effectorder of birthorder of dominanceorder of meritorder of merit method order of merit scale order statisticsordered alternative ordered change theory ordered metric scale ordered recallorderingordering relation orderlinessordinal dispersionordinal interaction ordinal measureordinal positionordinal scaleordinal variableordinaryordinateorectic processesOrestes complex oreximaniaorexisorganorgan inferiority organ languageorgan neurosisorgan of Cortiorgan of equilibration organ of language organ of locomotion organ of special sense organ of speechorgan of surveyorgan of tasteorgan of visionorgan sensationorgan sensuumorgan specificity organ transplantation organaorganicorganic amnesia organic anxiety organic brain damage organic brain disorder organic brain syndrome organic deafness organic dementia organic disorder organic factororganic feeling organic imitationorganic melancholiaorganic memoryorganic modelorganic psychosisorganic repressionorganic rulesorganic selectionorganic sensationorganic theory of organization organic theory of societyorganic therapyorganic traitorganic variableorganicismorganismorganismal theoryorganismic ageorganismic experienceorganismic psychologyorganismic theoryorganismic variableorganizationorganization behaviororganization chartsorganization cultureorganization developmentorganization development intervention organization effectorganization manorganization of animal society organization of dataorganization structureorganization theoryorganizational behaviorOrganizational Behavior and Human Decision Processes organizational changeorganizational climateorganizational commitmentOrganizational Commitment Scaleorganizational communicationorganizational crisisorganizational cultureorganizational departmentalizationorganizational designorganizational developmentorganizational diagnosisorganizational effectorganizational effectivenessorganizational factororganizational goalorganizational identificationorganizational maintenanceorganizational matrixorganizational mergerorganizational objectiveorganizational planningorganizational politicsorganizational processorganizational productivityorganizational psychologyorganizational reformorganizational shape organizational span of control organizational stress organizational structure organized criminal gang organized fermentorganized playorganized supplementary play organized wholeorganizerorganizingorganogenesisorganogenicorganolepticorganoleptic research organologyorganonorganon auditusorganon gustusorganon olfactusorganon visusorganonomyorganopathyorganotherapyorganphysiologyorganuleorganum gustatoriusorganum olfactoriumorganum sensuumorgasmorgasmic dysfunctionorgasmphaseorgyorientating reflexorientationorientation disturbanceorientation motion detector orientation of consuming requirement orientation of motor learning orientation of purchasing power orientation perceptionorientation reactionorientation reflexorientation responseorientation trainingorientational trainingorienting reflexorienting responseorienting response in animal orienting taskoriginorigin of consciousnessorigin of infants imaginationorigin of infants self consciousness origin of internal speechorigin of languageorigin of lifeorigin of mindorigin of neonatal attentionorigin of neonatal memoryorigin of speciesorigin of speechorigin of thinkingorigin of volitional movementoriginaloriginal civilizationoriginal dataoriginal descriptionoriginal groundoriginal learningoriginal natureoriginal qualityoriginal recordoriginal surveyoriginalityoriginality of thinkingoriginateoriginationoriginatorOrleans Hanna Algebra Prognosis Test Orleans Hanna Geometry Prognosis Test ornithineornithine decarboxylaseornithophobiaorphanorphanageorphanhoodorthergasiaortho raterortho sleeporthobiosisorthocephalyorthochoreaorthodoxorthodox psychoanalysis orthodoxsleeporthodromicorthodromic conduction orthogenesisorthogenicsorthogonal arrayorthogonal comparison orthogonal component orthogonal contrast orthogonal factororthogonal factor model orthogonal polynomial orthogonal rotation orthogonalityorthographic deep hypothesis orthographyorthomolecular psychiatry orthomutationorthophonyorthophoriaorthophreniaorthopsychiatryorthosisorthostatism orthosympatheticorthotonosorthropsiaoscillationoscillogramoscillographoscillometeroscillometryoscitancyOseretsky Tests of Motor Proficiency Osgoods surfaceosmaticosmesthesiaosmoceptorosmologyosmometerosmometric thirstosmonosologyosmoreceptorosmoscopeosmosisosmotic thirstosmylosphresiologyosphresiometerosphresisosphreticossicleossificationossifyOSTosteoanesthesiaosteogenesisosteologyOstwald color circleOstwald color pyramidOstwald color systemotaphoneother directedother perceptionother suggestionoticOtis Lennon Mental Ability TestOtis Quick Scoring Mental Ability TestOtis Self Administering Tests of Mental Ability otoconiteotoganglionotographyotohemineurastheniaotolaryngologyotolithotolithiasisotologyotomyastheniaotoneurastheniaotoneurologyotopathyotophoneotopiesisotorhinolaryngologyotorhinologyotosalpinxotosisOTUout groupout of balanceout of body experienceout of controlout of employout of stockout patientout patient treatment outbreakoutbreedingoutcomeoutcome goaloutcome researchouter causeouter directedouter directed behaviorouter earouter hair cellsouter nuclear layerouter plexiform layerouter segmentoutflow theoryoutgroupoutingoutline making in composition outlookoutlook on lifeoutnessoutputoutput of learningoutrageoutward manifestation outward senseoutward thing outwardnessoval windowovarian obesity ovariumovaryover adaptationover anxious reaction over arousalover compensationover constancyover developmentover exerciseover fatigueover learningover self confidence over the counter drugs overachievement overachieveroveractoverageoverall rating method overall reaction overcautionovercome overcompensation overconfidence overconfident overconstancy overcontrolled overcorrectionovercredulityovercredulousovercriticalovercrowovercrowdoverdeterminationoverdeterminedoverdevelopoverdoseoveremphasizeoverestimationoverexcitationoverexcitementoverextensionoverfatigueoverfeedingoverflow activityovergeneralizationovergrowthoverheard communicationoverinclusionoverinclusive thinkingoverindulgenceoverinhibitionoverjoyedoverjustificationoverjustification effectoverlapoverlappingoverlapping control of eating and drinking overlearningoverlearning reversal effect overloadoverload imitationoverload modelovermanovermanyovermasterovermaximalovermotivationoverpersuadeoverpopulationoverproductionoverproductivityoverprotectionoverrateoverreactionoverregulationoverresponseovershadowingoverstimulationoverstocking of commodities overstrainoverstrain of excitatory process overstrain of inhibitory process overstrain of nervous process overtovert behaviorovert crimeovert cultureovert homosexualityovert intrusionovert needovert responseovertimeovertoneovertrainingovertraining reversal effect overweightoverwrokoverwroughtoverzealousovigenesisovulationovumownownerownerlessownershipownguardownnessOxford indexoxidative enzymeoxyacoiaoxyblepsiaoxycodoneoxyecoiaoxyesthesiaoxygenoxygen consumptionoxygen debtoxygen deficitoxygen demandoxygen deprivation oxygen intake oxygen intake rate oxygen lack oxygen toxicity oxygeusia oxyhemograph oxylalia oxymorphone oxynervoneoxyopiaoxyopteroxyosmiaoxypathiaoxytocinPOMSponderal growth ponesiatrics pongidae ponographponspons cerebralis pons oblongata pons varolii pontibrachium ponticpontinepontine tegmentum pontobulbar pontocerellar pontopeduncularPonzo figurePonzo illusionpooled interdependencepooled observationspooling of errorpooling procedurepoorly nourisionpoppypopular childpopular responsepopularitypopulationpopulation analysispopulation by age group population by educational level population by sexpopulation composition population controlpopulation densitypopulation distribution population effectpopulation experiment population geneticspopulation growthpopulation meanpopulation of samples population parameterpopulation pressurepopulation registerpopulation regulation population revolutionpopulation situationpopulation statisticspopulation structurepopulation studypopulation theorypopulation vital statistics,population waiting for employmentPORPorch Index of Communicative Ability in Children poriomaniapornographyporphrismusporphyropsinportentPorters lawPorteus MazePorteus Maze Testporus acusticus externusporus acusticus internuspositionposition alternation problemposition constancyposition effectposition factorposition habitposition powerposition preferenceposition reactionposition receptorposition relationposition responsepositional leaderpositionalitypositionerpositioning analysispositioning movementpositive accelerationpositive adaptationpositive after imagepositive after potentialpositive and negative cases method positive attention seekingpositive attitudepositive attitude changepositive cathexispositive conceptpositive conditioned reflex positive contrastpositive contrast of reinforcement positive correlationpositive cuepositive errorpositive feedbackpositive goalpositive heuristicpositive identitypositive incentivepositive inducementpositive inductionpositive instancepositive interactionpositive interestpositive itempositive linear relationshippositive negative ambivalent quotient positive negative conflictpositive phasepositive phototaxispositive punishmentpositive reactionpositive recollectionpositive reference grouppositive regardpositive regencypositive reinforcementpositive reinforcerpositive reinforcing stimuluspositive rewardpositive sanctionpositive self feelingpositive skewnesspositive stimuluspositive symptompositive time errorpositive transferpositive transferencepositive valencepositivismpositivistic philosophypositivity biaspositron emission tomography posodynicspossessivenesspossibilitypossible errorpossible judgmentpossible outcomepost adaptationpost columnpost conventional level post enumeration survey post formation theorypost hoc fallacypost hoc testpost hornpost hypnotic amnesiapost hypnotic behaviorpost hypnotic suggestion post industrialized society post infections psychosis post partus bluespost surveypost term infantpost testpost traumatic amnesiapost traumatic delirium post traumatic dementia post voluntary attention postadolescence postadolescentpostage stamp experiment postcentral gyruspostconsummatory behavior postconventional level postconventional level of morality postconventional morality postcornupostcranialpostdecision dissonance postdormitumpostepilepticposteriorposterior association area posterior columnposterior commissureposterior comparisonsposterior hornposterior lobeposterior lobe hypophysisposterior lobe of the hypophysis posterior nuclear complexposterior parietal areaposterior pituitary glandposterior probabilityposterioritypostfebrile dementiapostfigurative culture postformation theory postganglionic fibreposthypnotic amnesiaposthypnotic effectsposthypnotic suggestion posthypophysispostinfectious depression postinfectious psychosis postjunctional potential postmaturepostmature infantspostmaturitypostmortem analysispostnatalpostnatal developmentpostnatal environment postpartum depression postpartum psychiatric reaction postpartum psychosis postpuberalpostpuberal phasepostpubertalpostpubertypostpubescencepostpubescentpostretinalpostretinal fiber postschizophrenic depression postsynapticpostsynaptic cellpostsynaptic inhibition postsynaptic membrane postsynaptic neuron postsynaptic potentialposttest designposttetanic facilitation posttetanic potentiationposttreatment follow up postulatepostulational method postural coordination postural reflex postural setpostural substrate postural swaypostureposture expression posture sense posturogaphy posturometer potamophobiapotencypotency dimension potentialpotential ability potential competition potential concept stage potential consumers potential criminal potential defect potential delinquent potential demand potential difference potential intellect potential learning potential market potential of action potential reasonpotential stimuluspotential unemploymentpotential utilitypotential variabilitypotential vocabularypotentialitypotentializationpotentiationpotentiometerpotentiometric methodpotentiometryPoter's generic competitive strategies potomaniapoundpourpouring in process in instruction poverty of movementpoverty of thoughtpowderpower coercive strategypower curvepower dependence theorypower distancepower driverpower efficiencypower factorpower functionpower function lawpower motivepower needpower of imagerypower of influencepower of intellectpower of soulpower of statistical test power of testpower relationpower spectrumpower structurepower strugglepower tacticspower testpowerlessnessPPSpractical activitypractical agopractical experiencepractical frameworkpractical intelligence practical knowledgepractical levelpractical philosophypractical problempractical psychologypractical purposepractical reasoningpractical situation practicalnesspractice curvepractice curve of action practice curve of performance practice effectpractice errorpractice limitpractice periodpractice supervisionpractice teachingpractice theorypractice theory of playpradictionpraecoxpraecox ejaculationpraecox feelingpragmatagnosiapragmatamnesiapragmatic value theorypragmatismpraiseprandial drinkingpraxeologypraxisPREpre adolescencepre after comparison methodpre competition emotional statepre competition mental preparation pre competition mental trainingpre competition psychological state pre competition statepre determinationpre drawing periodpre equatingpre established harmonypre feeding behaviorpre heatedpre inductionpre schematic stagepre school acepreadaptationpreadaptive mutationpreadaptive phasepreadolescencepreadultprearrival stagepreataxicpreattentive processprecariousprecausal thinkingprecausalityprecedenceprecedence effectprecedence methodprecentral gyrusprecisionprecision of samplingprecision of simple random sampling precociousprecocious childprecocious developmentprecocious pubertyprecocityprecodingprecognitionprecoidprecollege counselingprecollege guidanceprecompetition apathy state precompetition downcast symptom precompetition mental preparation precompetition overconfidence state precompetition overexciting state preconceptpreconceptionpreconceptual thoughtpreconditionpreconditioningpreconditioning adaptation preconsciouspreconscious process preconsciousnesspreconsummatory behavior preconventional level preconventional level of morality preconventional moralityprecoxprecursorprecursor of neocortexpredator modelpredatorypredatory behaviorpredecessorpredeterminationpredeterminepredetermined harmony predeterminismpredicate intersection model* predictabilitypredicted variablepredictionprediction of criminal mentality prediction of criminalityprediction of delinquencyprediction of early illegal behavior prediction of first offense prediction of recidivismprediction of second offense prediction to motor objectpredictivepredictive abilitypredictive attentionpredictive displaypredictive efficiencypredictive indexpredictive method of criminality predictive periodpredictive studypredictive validitypredictive valuepredictive variablepredictor variablepredilectionpredispositionpredominancepredominatepredormitiumpredreampreferencepreference methodpreference scalepreference testpreferential learning preferred colorpreferred orientation preferred orientation column preferred rewardprefoetal period preformationpreformation theory preformationismprefrontalprefrontal areaprefrontal cortexprefrontal leucotomy prefrontal lobeprefrontal lobotomy preganglionic fibre pregenital phasepregenital stage pregenitalitypregnancepregnancyprehensionprehension program prehension testpreinductionprejudgmentprejudicepreliminary analysispreliminary datapreliminary forecastingpreliminary investigation preliminary measurespreliminary mentorpreliminary observationPreliminary Scholastic Aptitude Test preliminary studypreliminary surveypreliminarytrialprelinguistic behaviorprelinguistic developmentprelogical thinking prelumirhodopisinPremack principlepremarital counselingprematurationprematurepremature childpremature ejaculationpremature infantpremature repeningprematuritypremeditated crimepremenstrual syndromepremenstrual tensionpremisepremiumpremorbidpremotor areaprenatalprenatal behaviorprenatal development prenatal development stage prenatal diagnosis prenatal environment prenatal influence prenatal periodprenatal training prenubilepreoccupationpreographpreogurative culture preoperation period preoperation thought preoperational stage preoptic areapreoptic regionpreorganic evolution preorganizerpreparationpreparation theory of play preparative period preparative stage preparatory interval preparatory measure preparatory response preparatory set preperception preponderance preponderantprepossessionprepotencyprepotentprepotent reflex prepotent response prepotent stimulus prepotential reflexes prepsychoticprepuberal stage prepuberty prepubescence prerecognition hypothesis prerecordprerelease prerepresentation prereproductive prerequisitepreretired education prerogativepresagepresbyacusispresbyatricspresbyopepresbyophrenia presbyopiaPresbyterian Presbyterianismpresbytia preschizophrenic preschoolpreschool child。
a r X i v :c o n d -m a t /0509275v 2 [c o n d -m a t .m e s -h a l l ] 13 N o v 2005Localization transition on complex networks via spectral statisticsM.Sade,T.Kalisky,S.Havlin and R.BerkovitsThe Minerva Center,Department of Physics,Bar-Ilan University,Ramat-Gan 52900,Israel(Dated:Oct.21,2005,version 3.1)The spectral statistics of complex networks are numerically studied.The features of the Anderson metal-insulator transition are found to be similar for a wide range of different networks.A metal-insulator transition as a function of the disorder can be observed for different classes of complex networks for which the average connectivity is small.The critical index of the transition corresponds to the mean field expectation.When the connectivity is higher,the amount of disorder needed to reach a certain degree of localization is proportional to the average connectivity,though a precise transition cannot be identified.The absence of a clear transition at high connectivity is probably due to the very compact structure of the highly connected networks,resulting in a small diameter even for a large number of sites.PACS numbers:89.75.Hc,72.15.RnThe Anderson transition predicts a transition from ex-tended (metallic)to localized (insulating)eigenstates as a function of the disorder or energy of a quantum sys-tem [1].This second-order phase transition turns out to be a very general property related to the transport of quantum and classical waves in disordered systems and signatures of it have been observed for electrons in met-als,microwaves in waveguides,light in liquids and gels and acoustical waves in the earth crust[2].For all these systems the clearest signature of the transition is the dif-ferent spread of a signal injected into the system.While in a metallic phase the injected wave will spread all over the system,in the insulating phase it will be localized in the vicinity of the injection point.This effect is a result of the constructive interference between time-reversed path throughout the system.Thus,the details of the tran-sition are strongly influenced by the dimensionality and topology of the system.The lower critical dimension,below which the system is localized for all values of disorder,is believed to be two [3],since the probability of returning to the origin (i.e.,constructive interference due to time reversal symmetry)is finite below d =2.The upper critical dimension [4](beyond which the critical exponents reach their mean-field values)remains uncertain although it is argued to be infinity [5].The parameters defining the transition are traditionally given as the critical disorder W c expressed in terms of the width of the distribution from which the on-sites energies in the Anderson model are drawn,and the critical exponent ν(for definitions see Sec.II).For square lattices with dimensionality d =3,4the values of W c and νare well established :for d =3,W c ∼16.5and ν∼1.5[6,7],while for d =4,W c ∼35and ν∼1[8,9].For higher dimensions the following extrapolation was of-fered [9]:ν∼0.8/(d −2)+0.5and W c ∼16.5(d −2),which was obtained by studying the transition on bifrac-tal topologies.Thus,the mean field critical index value of 1/2is obtained in the the upper critical dimension d =∞.There has been recently much interest in the proper-ties of random scale-free networks [10,11,12,13,14,15,16,17,18,19].These networks are characterized by the fact that each node is connected with some finite probability to any other node in the graph,which is very different from the usual topology of a real space lattice in which nodes are connected only to their neighbors.This leads to a very interesting behavior of the graphs when properties such as percolation,cluster structures,paths length etc.are considered [11].Interest in the influence this unusual topology has on the properties of wave inter-ference (i.e.,the Anderson transition of these graphs)is rising.Indeed,recently the Anderson transition in partic-ular networks,namely the small-world networks [20,21]and the Cayley tree [22]were studied.Older work on the localization properties of sparse random matrices [23]is directly relevant to Erd¨o s-R´e nyi graphs.Since scale free networks have an unusual topology for which anomalous classical properties have been found [24,25]it is of partic-ular interest to study the Anderson localization in these networks.Essentially,the probability to return to the origin de-fines the dimensionality of a system for the Anderson transition.Therefore,one may speculate that random graphs,which have only very long closed trajectories,correspond to systems with an infinite dimension [26].On the other hand,the critical disorder which depends roughly on the number of nearest neighbors Z is expected to follow [1]W c ∼Z ,which for a random graph with an average degree k (i.e.,the average number of connec-tions per node)corresponds to W c ∼ k .Thus,one may expect here an interesting situation in which the criti-cal index ν,which is determined by the dimensionality is close to a half,while the critical disorder is determined by k .This is very different than the situation described by the extrapolation given above,where the critical disorder for an infinite dimension should also be infinite.Beyond the general interest in investigating the Ander-son transition on scale free networks,and the fresh out-look it might provide on the localization phenomenon,the metal-insulator transition can provide insights into the functionality of complex networks.Consider for ex-ample an optical communication network.In such a net-2work the edges of the graph represent opticalfibers in which light propagates and the nodes represents a beam splitter which redistributes the incoming wave into the outgoing bonds connected to the node.Since for high quality opticalfiber there are essentially no losses or de-coherence on the bonds,the amplitude of the transmitted wave will not depend on the bond length(on the other hand phase will depend on the length).This network may be mapped on a tight binding Hamiltonian of the type described in Sec.II[27].An interesting question for such a scale free network is whether a wave injected into one of the nodes will produce a signal at all other node (which in the language of the Anderson localization is equivalent to the question is the system metallic)or only at afinite set of other nodes(i.e.,the system is insulat-ing).A metal-insulator transition in this network corre-sponds to a phase transition between those two phases as function of the properties of the nodes.Generally,for any complex network in which information propagates in a wave-like fashion and interference is possible,the An-derson transition will limit the spread of the information throughout the network.The paper is organized in the following way:In the next section(Sec.I)we describe the different networks which were considered,while in Sec II the spectral statis-tics method which was applied in order to identify the metal-insulator transition is outlined.In Sec.III the re-sults are depicted and some general characteristics of the localization on complex networks are discussed.I.CHARACTERISTICS OF THE DIFFERENTNETWORKSOur main goal is to study the Anderson transition for different complex networks.In this section we shall define the characteristics of the networks which will be consid-ered.A.Random GraphA random graph(or-random regular graph)is a graph with N nodes,each is connected to exactly k random neighbors[11].The diameter of a graph is the maximal distance between any pair of its nodes.In a random graph the diameter d is proportional to ln N.In Sec.III we shall present results of the level spacing distribution for random-regular graphs with k=3.B.Erd¨o s-R´e nyi GraphsIn their classical model from1959Erd¨o s and R´e nyi (ER)[28]describe a graph with N nodes where every pair of nodes is connected with probability p resulting in k = Np.For a large random graph the degree distribution follows the Poisson distribution:P(k)=e− kk k2−λ×K2−λ−m2−λ3[22]we checked a treein which5%of its nodes have higher degree (k =4)resulting in an average connectiv-ity 3.05and creating few closed trajectories -loops.II.METHODNow we turn to the calculation of the spectral statis-tics of these networks.First,one must construct the appropriate network structure,i.e.,to determine which node is connected to which.This is achieved using the following algorithm [14,17]:1.For each site choose a degree from the required distribution.2.Create a list in which each site is repeated as many times as its degree.3.Choose randomly two sites from the list and connect this pair of site as long as they are different sites.4.Remove the pair from the list.Return to 3.The diameter of a graph is calculated by building shells of sites [29].The inner shell contains the node with the highest degree,the next contains all of its neighbors,and so on.Of course,each node is counted only once.The diameter of the system is then determined by the number of shells.Two more options which were considered are defining the diameter by the most highly populated shell,or by averaging over the shells.The diameter obtained by the various methods are quite similar.The energy spectrum is calculated using the usual tight-binding Hamiltonian,H =iεi a †i a i − i,ja †j a i ,(3)where first term of the Hamiltonian stands for the dis-ordered on-site potential on each node i of the network.The on-site energies,εi are uniformly distributed over the range −W/2≤εi ≤W/2.The second term corresponds to the hopping matrix element which is set to 1,and i,j denotes nearest neighbor nodes which are determined ac-cording to the network structure.We diagonalize the Hamiltonian exactly,and obtain N eigenvalues E i (where N is the number of nodes in the graph)and eigenvectors ψi .Then we calculate the distribution P (s )of adjacent level spacings s ,where s =(E i +1−E i )/ E i +1−E i ,and ... denotes averaging over different realizations of dis-order and when relevant also over different realizations of node connectivities.One expects the distribution to shift as function of the on-site disorder from the Wigner surmise distribution (characteristic of extended states):P W (s )=πs4,(4)at weak disorder to a Poisson distribution (characteristic of localized states)at strong disorder:P P (s )=e −s .(5)12345s00.20.40.60.81P (s )W=2W=5W=7W=10W=15W=20W=25W=35Wigner PoissonFIG.1:The distribution P (s )for a 500sites scale-free graphwith λ=4(m =2).A clear transition from Wigner to Poisson is observed as a function of disorder.An example for such a transition is presented in Fig.1where a scale-free graph with λ=4and m =2was con-sidered.As W increases P (s )shifts toward the Poisson distribution.Additional hallmark features of the Ander-son transition,such as the fact that all curves intersects at s =2and the peak of the distribution ”climbs”along the Poisson curve for larger values of W are also appar-ent.Similar transition from Wigner to Poisson statis-tics is seen also for the other networks considered in this study.The transition point can be determined more accu-rately from calculating [30]γ= ∞2P (s )ds − ∞2P w (s )ds W c−1L 1/ν,(7)W0.60.70.80.91γFIG.2:γas function of W for different SF graphs sizes (λ=4and m =2).The typical behavior for finite size transition is seen,where a crossing in the size dependence of γbetween the metallic (small values of W )and localize (large value of W )regime is seen.01020304050|W/W c −1|L1/ν0.20.40.60.81γFIG.3:The scaling of γaccording to Eq.(7)for different SF,λ=4,m =2,networks sizes.Two branches,corresponding to the metallic and localized regimes,appear.where C is a constant.This relation enables us to extract both the critical disorder W c and the critical index ν.Scaling of the numerical data according to Eq.(7)yields two branches corresponding to the metallic and localized regimes,that are clearly seen in Fig. 3.The estimated values of νand W c (see Table I)are extracted by fitting the branches to a 4th order polynomial.Networklν2.9715.7±0.9Random-Regular (RR)11.80.66±0.08320.5±0.23Cayley-Tree100.51±0.0453.0512.4±0.1RR ”double peak”10.280.85±0.41TABLE I:Networks showing the localization transition.The value of l is for N =1000.III.RESULTS AND DISCUSSIONThe calculations for all networks mentioned above are performed for M different realizations,where M =1000,400,200,...,50for the corresponding num-ber of nodes:N =200,500,1000,...,4000.Ex-cept for the Cayley-tree networks for which M =4000,2000,1000,...,125,64for the corresponding tree sizes:N =63,127,255,...,2047,4095or L =6,7,8,...,11,12(where L is the number of ”generations”of the tree).Another exception is for Erd¨o s-R´e nyi graphs in which k is between 3and 3.5.The low connectivity of the graphs,results in one main cluster and relatively large number of not-connected nodes (about 5%).Thus,the calculations are made only for the largest cluster of each realization,since a procedure that considers all the nodes is skewed by the eigenvalues of small disconnected clusters [31].A clear localization transition is observed for a group of graphs which are all characterized by an average degree k smaller than 3.1,and an averaged last occupied shell l (for N =1000sites)larger or equal to 9.45.The results are summarized in Table I.The results for all the graphs (including those which show no clear signs of transition)can be scaled according to their average degree k .The higher the value of k is -the higher is the value of W needed in order to obtain a specific value of γ.Thus,the higher the average degree,the more metallic the system is,which makes sense.A cross section at γ=0.6of all curves is shown in Fig.4as a function of k .The k of the networks studied in Fig.4as well as the averaged last occupied shell l for N =1000sites are presented in Table II.The following observations can be gleaned out of the data for the different networks:(1)For all the networks that show a metal-insulator tran-sition,νis of order 1/2except for the Random-Regular ”double-peak”network which is the one with the highest value of connectivity that still shows a clear transition.A critical index of ν=0.5is expected for a system of infinite dimensionality.At k =3.1the value of νis significantly higher,but so is the estimate of the error bar.On the other hand for the Erd¨o s-R´e nyi graph with k =3.1no clear transition is observed.(2)All networks with connectivity above 3.1do not show clear signs of a metal-insulator transition.Nevertheless,Networkl (for N =1000)Scale-Free ,λ=4,m =212.463Random-Regular ”double peak”(p =0.95→k =3,p =0.05→k =5)10.283Cayley-Tree with loops103Scale-Free ,λ=3.5,m =29.363.1Erd¨o s-R´e nyi9.033.4Erd¨o s-R´e nyi8.333.35Erd¨o s-R´e nyi7.514Scale-Free ,λ=4,m =36.054.4Erd¨o s-R´e nyi6.317.5Erd¨o s-R´e nyi4.1TABLE II:The average connectivity k of all the networks considered in this study,as well as the averaged last occupied shell l for N =1000sites.246810<k>102030405060W (γ=0.6)FIG.4:The γcurves of the checked networks can be ordered by their k .The values of their W (γ=0.6)are presented as a function of k .One can notice the increasing W with k .one should be rather careful in interpreting this observa-tion since,as is clear from Table II,larger values of klead to smaller size,l ,of the network for the same number of nodes.Moreover,from the two networks which have the same k =3.1,only the one with the higher value of l shows clear signs of the metal insulator transition.Thus,the absence of transition may be an artifact of the small size of networks with high average connectivity.(3)The critical disorder W c fluctuates in the range of 12−20(Table I).Due to the small range of k (2.97−3.1),it is hard to determine any relation between k and W c .(4)On the other hand,there is a clear relation between the amount of disorder needed in order to reach a par-ticular value of γ(i.e.,the value of W needed to reach a certain degree of localization)and k .As can be seen in Fig.4,a linear dependence W (γ=0.6)∝ k is ob-served.Thus,the gross features of the Anderson metal-insulator transition are similar for a wide range of dif-ferent networks.The critical index for all the networks studied here are within the range expected for a system of infinite dimensionality,and the connectivity influences the degree of localization.Thus,complex networks are an example of a topology for which the critical index follows the mean field prediction,but the critical disorder can be tuned by the connectivity.On the other hand,the fact that networks with high connectivity are very compact raises the problem of identifying the transition point.It is hard to extend the usual finite-size scaling method to networks with high connectivity since the number of sites grows very rapidly with size,while for small network sizes the crossover behavior of the γcurves is very noisy.This results in an inability to clearly identify the Anderson transition,although it can not be ruled out the possibility that there is a critical connectivity for complex networks above which no transition exist.We acknowledge very useful discussions with J.W.Kantelhardt and helpful correspondence with Y.V.Fy-odorov,as well as support from the Israel Academy of Science (Grant 877/04),the ONR,and the Israeli Center for Complexity Science.This work was also supported by the European research NEST Project No.DYSONET 012911.[1]P.W.Anderson,Phys.Rev.109,1492(1958); D.J.Thouless,Phys.Rep.13,93(1974).[2]see these reviews and references within:G.Bergmann,6Phys.Rep.107,1(1984);P.A.Lee and T.V.Ramakr-ishnan,Rev.Mod.Phys.57,287(1985);R.Berkovits and S.Feng,Phys.Rep.238,135(1994);M.C.W.van Rossum and T.M.Nieuwenhuizen,Rev.Mod.Phys.71, 313(1999).[3]E.Abrahams,P.W.Anderson,D.C.Licciardello,andT.V.Ramakrishnan,Phys.Rev.Lett.42,673(1979).[4]C.Castellani,C.DiCastro,and L.Peliti,J.Phys.A19,L1099(1986);H.Kunz and B.Souillard,J.Phys.Lett.44,L503(1983);J.P.Straley,Phys.Rev.B28,5393 (1983);T.Lukes,J.Phys.C12,L797(1979);K.B.Efetov,Sov.Phys.JETP61,606(1985).[5]A.D.Mirlin and Y.V.Fyodorov,Phys.Rev.Lett.72,526(1994).[6]B.Kramer and A.MacKinnon,Rep.Prog.Phys.56,1496(1993).[7]K.Slevin,T.Ohtsuki and T.Kawarabayashi,Phys.Rev.Lett.84,3915(2000).[8]I.K.Zharekeshev and B.Kramer,Ann.Phys.(Leipzig)7,442(1998).[9]M.Schreiber and H.Grussbach,Phys.Rev.Lett.76,1687(1996).[10]R.Albert,H.Jeong,and A.L.Barab´a si,Nature,406,6794(2000).[11]For a recent review see:R.Albert and A.L.Barab´a si,Rev.Mod.Phys.74,47(2002).[12]S.N.Dorogovtsev and J.F.F.Mendes,Evolution of Net-works-From Biological Nets to the Internet and WWW (Oxford University Press,2003)[13]R.Pastor-Satorras and A.Vespignani,Evolution andStructure of the Internet:A Statistical Physics Approach (Cambridge University Press,2004)[14]R.Cohen,K.Erez, D.ben Avraham,and S.Havlin,Phys.Rev.Lett.85,4626(2000);Phys.Rev.Lett.86, 3682(2001).[15]R.Cohen and S.Havlin,Phys.Rev.Lett.90,058701(2003).[16]B.Bollobas and O.Riordan in Handboook of graphs andnetworks,eds.S.Bornholdt and H.G.Schuster(Wiley-VCH,Berlin,2002)[17]M.Molloy and B.Reed,Combinatorics,Probability andComputing7,295(1998).[18]D.S.Callaway,M.E.J.Newman,S.H,Strogatz,and D.J.Watts,Phys.Rev.Lett.85,5468(2000).[19]M.E.J.Newman,S.H,Strogatz,and D.J.Watts,Phys.Rev.E64,026118(2001).[20]C.P.Zhu and S.-J.Xiong,Phys.Rev.B62,14780(2000).[21]O.Giraud, B.Georgeot and D.L.Shepelyansky,quant-ph/0503188.[22]M.Sade and R.Berkovits,Phys.Rev.B68,193102(2003).[23]Y.V.Fyodorov and A.D.Mirlin,Phys.Rev.Lett.672049(1991);S.N.Evangelou and E.N.Economou,Phys.Rev.Lett.68,361(1992);Y.V.Fyodorov et.al.,J.Phys.I (France)2,1571(1992);Y.V.Fyodorov and A.D.Mirlin, Phys.Rev.B55,R16001(1997).[24]L.A.Braunstein,S.V.Buldyrev,R.Cohen,S.Havlin,and H.E.Stanley,Phys.Rev.Lett.91,168701(2003).[25]E.Lopez,S.V.Buldyrev,S.Havlin,and H.E.Stanley,Phys.Rev.Lett.94,248701(2005).[26]E.M.Bollt and D.ben-Avraham,cond-mat/0409465.[27]B.Shapiro,Phys.Rev.Lett.48,823(1982).[28]P.Erd¨o s and A.R´e nyi,Publ.Math.Debrecen6,290(1959).[29]T.Kalisky,R.Cohen,D.ben Avraham,and S.Havlin,in Lecture Notes in Physics:Proceedings of the23rd LANL-CNLS Conference,”Complex Networks”,Santa-Fe,2003,eds. E.Ben-Naim,H.Frauenfelder,and Z.Toroczkai(Springer,Berlin,2004).[30]B.I.Shklovskii,B.Shapiro,B.R.Sears,mbrianides,and H.B.Shore,Phys.Rev.B47,11487(1993). [31]R.Berkovits and Y.Avishai,Phys.Rev.B53,16125(1996).。