The magnetic field of pulsars and the gravito-magnetic theory
- 格式:pdf
- 大小:634.79 KB
- 文档页数:20
宇宙知识英语作文Title: Exploring the Vast Expanse of the Universe;Introduction:The universe, with its infinite expanse and countless wonders, has fascinated humanity for centuries. Exploring the mysteries of space has been a continuous endeavor, and our understanding of the cosmos has grown exponentially over time. This essay aims to provide a glimpse into the vastness of the universe and the remarkable knowledge we have acquired about it.The Expanding Universe:One of the most fundamental discoveries in cosmology is the fact that the universe is expanding. Scientists have observed that galaxies are moving away from each other, indicating that the universe is constantly growing. This led to the formulation of the Big Bang theory, suggesting that the universe originated from a single, incredibly dense point billions of years ago.Galaxies and Stars:The universe is comprised of billions of galaxies, each containing billions of stars. Galaxies come in various shapes and sizes, from spiral and elliptical to irregular formations. They are held together by gravity and play a crucial role in shaping the cosmic landscape.Stars, the building blocks of galaxies, are fascinating celestial bodies. They range in size, temperature, and luminosity. Our own star, the Sun, is just one of billions in the Milky Way galaxy. Stars are born from vast clouds of gas and dust, and their lifecycle involves fusion reactions that release tremendous amounts of energy.Black Holes and Neutron Stars:Among the most enigmatic entities in the universe are black holes and neutron stars. Black holes are regions of spacetime with such intense gravitational forces that nothing, not even light, can escape their grasp. They are formed from the remnants of massive stars that have collapsed under their own gravity.Neutron stars, on the other hand, are incredibly dense stellar objects composed mainly of neutrons. They are the remnants of massive stars that have undergone supernova explosions. Neutron stars possess powerful magnetic fields and emit beams of radiation, which can be observed as pulsars.Cosmic Microwave Background Radiation:The discovery of cosmic microwave background radiation (CMB) has provided strong evidence for the Big Bang theory. CMB is faint radiation that permeates the entire universe and is a remnant of the intense heat generated during the early stages of the universe's formation. By studying the characteristics of CMB, scientists have gained valuable insights into the evolution of the universe. Dark Matter and Dark Energy:Despite our significant progress, there are still some mysteries surrounding the universe. Dark matter and dark energy, for instance, remain elusive. Dark matter is believed to make up a significant portion of the universe, exerting gravitational forces but not interacting with light. Dark energy, on the other hand, is thought to be responsible for the accelerated expansion of the universe.Conclusion:The universe is an awe-inspiring realm, filled with endless wonders waiting to be explored. Our understanding of the cosmos has come a long way, thanks to the relentless pursuit of knowledgeand advancements in technology. However, there is still much more to learn and discover. As we delve deeper into the mysteries of the universe, we gain a greater appreciation for its grandeur and our place within it.contribute to its preservation. Join us on and embark on an exciting journey to explore and appreciate the beauty and importance of our natural environment.。
你听说过脉冲星吗英语作文In the vast universe, there exists a mysterious and fascinating celestial object - the pulsar. Have you ever heard of it?A pulsar is a highly distinctive and intriguing object. It is a type of star that exhibits unique properties and behaviors.Pulsars emit regular pulses of radiation, which make them stand out in the universe. These pulses can be detected and studied using specialized instruments.They are incredibly compact and have extremely high densities. Their gravitational pull is immense, making them some of the most powerful objects in the cosmos.The study of pulsars provides valuable insights into various aspects of astronomy. It helps us understand the nature of matter, the behavior of stars, and the properties of extreme environments.Pulsars also play a crucial role in our understanding of the universe's evolution. They offer clues about the formation and evolution of stars, as well as the processes occurring in deep space.Discovering and studying pulsars requires advanced technological tools and expertise. Astronomers use radio telescopes and other sophisticated instruments to detect and analyze their signals.。
麦婷婷,高梦祥,李利,等. 磁场调控丝状真菌生长发育及代谢产物合成的研究进展[J]. 食品工业科技,2023,44(13):450−457.doi: 10.13386/j.issn1002-0306.2022090173MAI Tingting, GAO Mengxiang, LI Li, et al. Research Progress on the Regulation of Magnetic Field on the Growth and Development of Filamentous Fungi and the Synthesis of Metabolites[J]. Science and Technology of Food Industry, 2023, 44(13): 450−457. (in Chinese with English abstract). doi: 10.13386/j.issn1002-0306.2022090173· 专题综述 ·磁场调控丝状真菌生长发育及代谢产物合成的研究进展麦婷婷1,高梦祥1,李 利1,张佳兰2, *,王劲松3(1.长江大学生命科学学院,湖北荆州 434025;2.长江大学动物科学学院,湖北荆州 434025;3.荆楚理工学院生物工程学院,湖北荆门 444800)摘 要:磁场作为一种普遍存在的环境因子,影响着微生物的生长及代谢。
丝状真菌是一类重要的异养型真核生物,在食品工业和生物医药领域有广泛应用。
目前,采用磁场处理丝状真菌已成为工业重要研究目标,而红曲霉、黑曲霉和黄曲霉是丝状真菌常用的典型菌种。
本文分别介绍了磁场调控红曲霉、黑曲霉、黄曲霉等丝状真菌生长及代谢的研究进展,通过不同磁场类型、作用时间、磁场强度等多种磁场参数分析三种丝状真菌磁场磁效应,阐述了磁场对三种典型丝状真菌的生长发育及其代谢物影响规律和代谢关系。
a r X i v :a s t r o -p h /0601246v 2 26 O c t 2006Invited talk,47th APS DPP Meeting,Denver,Oct.24–28,2005;Phys.Plasmas 13,056501(2006)[astro-ph/0601246]Turbulence,magnetic fields and plasma physics in clusters of galaxiesA.A.Schekochihin 1,∗and S.C.Cowley 2,31DAMTP,University of Cambridge,Cambridge CB30WA,UK2Department of Physics and Astronomy,UCLA,Los Angeles,California 90095-15473Plasma Physics Group,Imperial College,Blackett Laboratory,Prince Consort Road,London SW72BW,UK(Dated:February 5,2008)Observations of galaxy clusters show that the intracluster medium (ICM)is likely to be turbulent and is certainly magnetized.The properties of this magnetized turbulence are determined both by fundamental nonlinear magnetohydrodynamic interactions and by the plasma physics of the ICM,which has very low collisionality.Cluster plasma threaded by weak magnetic fields is subject to firehose and mirror instabilities.These saturate and produce fluctuations at the ion gyroscale,which can scatter particles,increasing the effective collision rate and,therefore,the effective Reynolds number of the ICM.A simple way to model this effect is proposed.The model yields a self-accelerating fluctuation dynamo whereby the field grows explosively fast,reaching the observed,dynamically important,field strength in a fraction of the cluster lifetime independent of the exact strength of the seed field.It is suggested that the saturated state of the cluster turbulence is a combination of the conventional isotropic magnetohydrodynamic turbulence,characterized by folded,direction-reversing magnetic fields and an Alfv´e n-wave cascade at collisionless scales.An argument is proposed to constrain the reversal scale of the folded field.The picture that emerges appears to be in qualitative agreement with observations of magnetic fields in clusters.I.INTRODUCTIONClusters of galaxies are vast and varied objects that have long attracted the attention of both observers (in recent decades spectacularly aided by X-ray and radio telescopes)and theoreticians.The observed properties of clusters have proved far from easy to explain as new data has confounded many old theories.The overall budget of the cluster constituents is roughly as follows:∼75%of cluster mass is dark mat-ter,whose sole function is assumed to be to provide the gravitational well,∼20%of cluster mass is the diffuse X-ray-emitting plasma (the intracluster medium,or ICM),while the galaxies have an all but negligible mass.The plasma that makes up the ICM is made of hot and ten-uous ionized hydrogen:temperatures are in the range of 1−10keV,number densities ∼10−1−10−3cm −3,with colder,denser material found in the cool cores and hotter,more diffuse one in the outlying regions.70Necessarily,the first observations and the first the-oretical models of clusters concerned what one might call large-scale features,such as the overall profiles of mass and temperature,the structure formation and the role of the central objects.As observations increased in accuracy and resolution,the ICM was revealed to be much richer than simply a dull cloud of X-ray glow smoothly petering out with distance from the center.A great panoply of features has been detected:bub-bles,filaments,ripples,edges,shocks,sound waves,etc.,as well as very chaotic density,temperature and abun-dance distributions.1,2,3,4,5It is particularly the presence of chaotic fields and the evidence that this chaos exists in a range of scales 4,6,7that makes one expect that ICM,like so many other astrophysical plasmas,is in a turbulent state.It is essential to know the properties of this turbu-lence in order to predict the current and future statisticalmeasurements of the plasma and magnetic fields (spectra,correlation and distribution functions,etc.)and to model correctly the transport processes in the ICM that deter-mine,for example,the overall temperature profiles.8,9We shall assume that the physics of small scales is,at least to some degree,independent of large-scale cir-cumstances and that we can,therefore,gain some use-ful understanding of the turbulence in clusters by ignor-ing large-scale features and considering a homogeneous subvolume of the ICM.In what follows,after reviewing briefly what is known about fluid motions (Sec.II)and magnetic fields (Sec.III)in clusters,we describe,mostly in qualitative terms,what we consider to be the essential aspects of the small-scale physics of the ICM (for plasma physics,see Sec.IV).This will lead us to a tentative overall picture of the structure of the ICM turbulence (Sec.VI)as well as of the origin of its magnetic compo-nent (Sec.V).We must emphasize that the current state of the debates on the nature of turbulence in clusters (and,indeed,on its very existence 10)is such that even a fundamental,conceptual view of the problem has not yet been agreed,and,therefore,a quantitative theory —even an incomplete one such as exists for hydrodynamic turbulence —remains a matter of future work.II.TURBULENCEThe failure of one of the instruments on the ASTRO-E2satellite has set offthe planned direct detection of cluster turbulence 11,12into the (probably not very dis-tant)future.However,indirect evidence of turbulent gas motions does exist:three recent examples are the broad spectrum of pressure fluctuations measured in the Coma cluster 4,detection of subsonic gas motions in the core of the Perseus cluster,13and,again for the Perseus cluster,2TABLE I:Cluster ParametersParameter Expression Cool cores a Hot ICMT observed3×107K108Kn observed6×10−2cm−310−3cm−3v th,i(2T/m i)1/2700km/s1300km/sνii1.5nT−3/2b5×10−13s−12×10−15s−1λmfp v th,i/νii0.05kpc30kpcµ v th,iλmfp1028cm2/s1031cm2/sη3×1013T−3/2b200cm2/s30cm2/sU inferred250km/s300km/sL inferred10kpc200kpcL/U inferred4×107yr7×108yrRe UL/µ 702Rm UL/η4×10276×1029t visc(L/U)Re−1/25×106yr5×108yrl visc L Re−3/40.4kpc100kpcl res L Rm−1/25000km8000kmΩi,eq eB eq/cm i0.3s−10.04s−1ρi,eq v th,i/Ωi,eq3000km30,000kmB0B eqρi,eq/λmfp5×10−17G2×10−19G B1Eq.(13)c3×10−14G2×10−17G B2Eq.(15)c8×10−7G2×10−7GB visc B eq Re−1/49×10−6G4×10−6GB eq(8πm i nU2/2)1/23×10−5G4×10−6Gβeq8πnT/B2eq820l⊥(B2/B eq)L0.2kpc7kpcl B observed1kpc10kpca These numbers are based on the parameters for the Hydra A cluster given in Ref.20.b In these expressions,n is in cm−3,T in Kelvin.c We usedα=3/2to get specific numbers,but the outcome is not very sensitive to the value ofα.the broadening(assumed to be caused by turbulent dif-fusion)of abundance peaks associated with the brightestcluster galaxies.14These and other studies and models based on observational data appear to converge in ex-pecting turbulentflows with rms velocities in the range U∼102−103km/s at the outer scales L∼102kpc. The energy sources for this turbulence are probably the cluster and subcluster merger events and/or,especiallyfor the turbulence in cool cores,the active galactic nu-clei(AGN).The aforementioned observational estimates of the strength and scale of the turbulence are in order-of-magnitude agreement with the outcomes of numerical simulations of cluster formation11,15,16and of the buoy-ant rise of radio bubbles generated by the AGN.17,18Fur-ther discussion and references on the stirring mechanisms for cluster turbulence can be found in Refs.19,20,21. While estimates of the turbulence parameters appear robust roughly to within an order of magnitude,a more quantitative set of numbers is elusive,partly because of the indirect and difficult nature of the observations,partly because the conditions vary both in different clus-ters and within each individual cluster.Here we shalladopt twofiducial sets of parameters:one for cool cores and one for the bulk of hot cluster plasma.These aregiven in Table I(along with some theoretical quantities that will arise in Sec.V and Sec.VI).They will allowus to make estimates that will have the virtue of being consistent and systematic but must not be interpreted asprecise quantitative predictions.They are representative of the range of conditions that can be present in clus-ters.The turbulence is assumed to be stirred at the outer scale L(with rms velocity U at this scale)and to havea Kolmogrov-type cascade below this scale.The small-scale cutoffis determined by the microphysical proper-ties of the cluster plasma.In Table I,we give the value of the particle mean free pathλmfp,which can be usedin a naive estimate of viscosity:µICM∼v th,iλmfp,where v th,i=(2T/m i)1/2is the ion thermal speed.We seethat this gives fairly low values for the Reynolds number, Re∼UL/µICM.It is this feature of the ICM that con-tinues to fuel doubts about its ability to support turbu-lence,at least in the strict,hydrodynamic high-Reynolds-number sense.71However,one should be cognizant of the fact that whatever type of turbulence might exist in the cluster plasma,it is certainly not hydrodynamic,because this plasma is highly electrically conducting72and mag-netized.The presence of the magneticfields not only has a dynamical effect on the turbulence(due to the action of the Lorentz force,the medium acquires a certain elas-tic quality),but also changes the transport properties of the plasma itself:the viscosity,in particular,becomes strongly anisotropic.22These issues will constitute the main subject of this paper,butfirst let us briefly de-scribe what is known about magneticfields in clusters.III.MAGNETIC FIELDSThefirst observed signature of cluster magneticfields was the diffuse synchrotron radio emission in the Coma cluster detected in1970.23Starting from early1990’s, increasingly detailed measurements of the Faraday Ro-tation in the emission from intracluster radio sources have made possible quantitative estimates of the mag-neticfield strength and scales in a large number of clusters.24,25,26Randomly tangled magneticfields with rms strength of order B rms∼1−10µG are consistently found,with thefields in the cool cores of the cooling-flow clusters somewhat stronger than elsewhere.This is fairly close to the value B eq that corresponds to magnetic energy equal to the energy of turbulent motions(see Ta-ble I).Thus,the magneticfield must be dynamically important.The estimates for the tangling scale l B of the field are usually arrived at by assuming that direction reversals along the line of sight(probed by the Faraday Rotation measure)can be described as a random walk with a single step size equal to l B(the estimate of B rms is obtained in conjunction with this model).This gives3 l B∼1−10kpc.The single-scale model is almost certainly not a correctdescription on any but a very rough level.Fortunately,much more detailed information on the spatial structureof the clusterfields is accessible.First,using certain sta-tistical assumptions(most importantly,isotropy),it ispossible to compute magnetic-energy spectra from themaps of the Faraday Rotation measure associated withextended radio sources(the radio lobes of the jets emerg-ing from the AGN—these can be as large as∼102kpcacross).6,27This has been done most thoroughly for a ra-dio lobe located in the cool core of the Hydra A cluster.7The spectrum has a peak at k≃2kpc−1followed bywhat appears to be a power tail consistent with k−5/3down to the resolution limit of k≃10kpc−1.The rmsmagneticfield strength is B rms=7±2µG.The second source of information on the clusterfieldstructure is the polarized synchrotron emission,whichprobes the magneticfield in the plane perpendicular tothe line of sight.28Such data,while widely used for Galac-tic magneticfield studies,29has until recently not beenavailable for clusters.This is now changing:thefirstanalysis of polarized emission from a radio relic in thecluster A2256reveals the presence of magneticfilamentswithfield reversals probably on∼20kpc scale,which,however,is dangerously close to the resolution scale.30 This data is representative of the situation in the bulk of the ICM,rather than in the cores.Statistical analysis of such data will make possible quantitative diagnosis of thefield structure and its dynamical role.31Thus,our knowledge of the magneticfields in clusters, while far from perfect,is more direct and more detailed than that of the turbulent motions of the ICM.It is also due to improve dramatically with the arrival of new radio telescopes such as LOFAR and SKA.73IV.PLASMA PHYSICSThe key property of the ICM as plasma is that it is only weakly collisional and magnetized:given the ob-served values of the magneticfield,the ion gyroradius isρi∼104km,which is much smaller than the mean free path.Asρi≪λmfp already for dynamically veryweakfields(B≫B,see Table I),this is true both inthe observed present state of the ICM and during mostof its hypothetical past,when the magneticfield was be-ing amplified from some weak seed value.In a plasmawithρi≪λmfp,the equations for theflow velocity u and for the magneticfield B may be written in the followingform,valid at time scales≫Ω−1i(Ωi=eB/m i c is the ion cyclotron frequency)and spatial scales≫ρi=v th,i/Ωi,ρd u2 +∇· ˆbˆb p⊥−p +B2 ,(1)d B4πhas been absorbed intoB,and the resistive term has been omitted in Eq.(2)inview of the tiny value of the resistivity.The turbulentmotions in clusters are subsonic(U<v th,i),so we maytake∇·u=0and setρ=1.The magneticfield is,thus,in units of velocity,pressure in units of velocity squared.The proper way to compute p⊥and p is by a kineticcalculation.In the collisional limit,this was done in Bra-ginskii’s classic paper.22It is instructive to obtain hisresult in the following heuristic way that highlights thephysics behind the formalism.32Charged particles mov-ing in a magneticfield conserve theirfirst adiabatic in-variantµ=m i v2⊥/2B.Whenλmfp≫ρi,this conserva-tion is only weakly broken by collisions.As long asµisconserved,any change in thefield strength causes a pro-portional change in p⊥:summing up thefirst adiabaticinvariants of all particles,we get p⊥/B=const.Then1dt∼1dt−νiip⊥−pBdBdt u2 2=−µ |ˆbˆb:∇u|2=−µ 1dt 2 .(5)Thus,the Braginskii viscosity only dissipates such mo-tions that change the strength of the magneticfield.Mo-tions that do not affect B are allowed at subviscous scales.In the weak-field regime,these motions take the form ofplasma instabilities.When the magneticfield is strong,acascade of Alfv´e n waves can be set up below the viscousscale.Let us elaborate.The simplest way to see that pressure anisotropies leadto instabilities is as follows.32Imagine that the large-scaleenergy sources stir up a“fluid”turbulence with u,p⊥,p ,B at time and spatial scales above viscous.Wouldsuch a solution be stable with respect to much higher-frequency and smaller-scale perturbations?Linearizing4 Eq.(1)and denoting perturbations byδ,we get−iωδu=−i k(δp⊥+BδB)+ p⊥−p +B2 δK+iˆb k δp⊥−δp − p⊥−p −B2 δBβ αΩi,(8)whereα=α1+α2>0.This changes the characteristicsof the turbulence:the effective mean free path of theparticles isλmfp,eff∼v th,i/νeff,the effective(parallel)viscosity of the ICM isµ ,eff∼v th,iλmfp,effand,therefore,the effective Reynolds number isRe eff∼UL v2νeff.(9)th,iOn the other hand,using Eq.(4)with effective viscosity,we getµ ,effand|∇u|∼(U/L)Re1/2eff|∆|∼ UαBρi,eq U√B2eqFIG.2:The mechanism of the fluctuation dynamo.Equation (11)models qualitatively the assumed out-come of an as yet inexistent proper theory of the viscos-ity of magnetized ICM.Its solution is plotted inFig.1.Weshalluse this model shortly in our discussion of the fluctuation dynamo in clusters.V.FLUCTUATION DYNAMOIn Sec.III,we reviewed the observational evidence that testified to the presence of a dynamically significant ran-domly tangled magnetic field in clusters.What is the origin of this field?There are numerous physical rea-sons to expect that a certain amount of seed magnetic energy predates structure formation and was,therefore,already present at the birth of clusters.40,41Typical val-ues given for the strength of such field are in the range of B seed ∼10−21−10−17,although this may be an underestimate.42It then falls to the random motions of the cluster plasma to amplify the field to its observed magnitude of a few µG.77This,indeed,they should be able to do by means of the fluctuation (or small-scale)dynamo mechanism:the random stretching of the field.It is a fundamental property of a succession of random (in time)linear shears that it leads on the average to expo-nential growth of the energy of the magnetic field frozen into the medium.43,44,45The rate of growth is roughly equal to the rate of strain (shear,or stretching rate)of the random flow.While the mathematical theory of this process can be nontrivial,46the physics of it is basically illustrated by Fig.2.In Kolmogorov turbulence,the rate of strain is domi-nated by the viscous scale,so |∇u |∼t −1visc∼(U/L )Re 1/2.In fact,what is relevant for the growth of the magnetic field is not the full rate-of-strain tensor but its “paral-lel”component,ˆbˆb :∇u .Since this is exactly the type of motion damped by Braginskii viscosity [see Eqs.(4–5)],we can,for the purposes of the fluctuation dynamo,ignore any subviscous-scale velocity fluctuations.Thus,the magnetic field should grow according to 781dt=ˆbˆb :∇u ∼U Lv th ,i(1−t/t c )2+α,(14)where B (0)∼the greater of B seed and B 1and t c =(2+α)(L/U )Re −1/2[B 1/B (0)]1/(2+α)is at most (for B seed <B 1)the viscous turnover time t visc associated with the collisional Re.The explosive stage contin-ues until the amplified field starts suppressing the in-stabilities,i.e.,when βdrops to values comparable to(v th ,i /U )2Re 1/2eff.This happens at B ∼B 2,whereB 2=B eqv th ,iL1/(5+2α).(15)Thus,we have a mechanism that amplifies the field by many orders of magnitude from any strength above B 1to B 2in finite,cosmologically short,time ∼t c .The value of B 2turns out to be only just over an order of magnitude below B eq (see Table I).Further growth of the field is algebraically slow (B ∼t 1/2),but it does not have to go on for a very long time because B 2is already quite close to the observed field strength.To be precise,there are two algebraic regimes.During the first,Re effis still controlled by thesecond term in Eq.(11)as B hovers just below B eq Re −1/4effwhile B is increasing and Re effis decreasing.Eventually,B ∼B visc =B eq Re −1/4and Re effis returned back to Re (plasma instabilities are suppressed).79As this is also the field strength at which the field has energy comparable to the energy of the viscous-scale motions,any further growth of the field is a nonlinear process,in which the back reaction of the field on the flow has to be taken intoFIG.3:Evolution of the magnetic-field strength for the cool-core parameters of Table I.account.This can be done by assuming that,as thefield grows,it can no longer be stretched by motions whose energy it exceeds and that,therefore,at any given time, the dominant stretching is exerted by motions whose en-ergy is equal to that of thefield.Denoting their velocity and scale by u l and l and using u2l∼B2,we have58,59 dB2∼u3ll1−A 1FIG.5:A schematic illustration of the structure of cluster turbulence proposed in Sec.VI.ratio l /l⊥,the larger is the contrast between the strong field in the straight segments of the folds and the weakfield in the corners.We assume that the rms value of thefield is determined by the straight segments becausethese are thefields that are stretched by turbulence.It istheir growth that we studied in Sec.V.In the saturated state,we expect B straight∼B eq.For such a strongfield, the plasma instabilities are suppressed.It is intuitively clear that they must be suppressed in the regions of theweakfield as well,i.e.,thefield there cannot be weaker than B2[Eq.(15)].Indeed,as we saw in Sec.V,the explosive dynamo mechanism brings anyfield up to this value nearly instantaneously(for B>B2,the growth is much slower).We may conjecture that the maximum aspect ratio of the folds is set by the maximum contrast in thefield strength l /l⊥∼B eq/B2.81Substituting the numbers,wefind(Table I)that this prediction gives the field reversal scale no smaller than a few per cent of the outer scale L(taking l ∼L).This is in passable agree-ment with the observational evidence,which is the best one can expect,given the highly imprecise nature of our argument,of most observational inferences,and of the definitions of such quantities as l B,l⊥,l and L.For comparison of our model with observational data for a number of individual clusters,see Ref.20.It is fair to acknowledge that the above argument,while providing a useful constraint,falls short of a sat-isfactory explanation of thefield structure.One might argue that if,in the course of the turbulent stretch-ing/shearing of thefield,a region offield strength below B2appears(as explained above,in the corner of a fold), Re effthere becomes very large and a localized spot of high-Reynolds-number turbulence is formed.This should have two principal effects.Thefirst is akin to that of a locally enhanced turbulent resistivity,so thefield that violates our constraint is continuously destroyed.The second is a burst of explosivefluctuation dynamo in the spot,which produces more foldedfield with B>B2 and thus shuts itself down.These folds are then further stretched,sheared,etc.,again subject to the constraint that they are destroyed and replaced by new ones wher-ever a spot of weakfield appears.We do not currently have a more detailed mechanistic scenario of how exactly the folded structure with l /l⊥∼B eq/B2is established. It may be feasible to test these ideas numerically by solv-ing MHD equations with viscosity locally determined by the magnetic-field strength according to Eq.(11).Let us assume that clusterfields do indeed have a folded structure with a direction reversal scale l⊥∼(0.01...0.1)L,possibly determined by the argument given above.82The magnetic-energy spectrum then peaks at k∼1/l⊥.What is the structure of the turbulence above and below this scale?At scales l≪l⊥,the magnetic field reversing at the scale l⊥will appear uniform and,in accordance with the old idea of Kraichnan63could sup-port a cascade of Alfv´e n waves.This cascade can rig-orously be shown to be described by the equations of Reduced MHD at collisionless scales all the way down to the ion gyroscale.37The currently accepted theory of such a cascade,primarily associated with the names of Goldreich and Sridhar,64is based on the conjecture that at each scale,the Alfv´e n frequency is equal to the turbu-lent decorrelation rate.The result is a k−5/3spectrum of Alfv´e nicfluctuations—this possibly explains what ap-pears to be a k−5/3tail in the observed spectrum of mag-netic energy for the Hydra A core.7We cannot embark on a detailed discussion of the theory of the Alfv´e n-wave cascade here,so the reader is referred to Ref.65for a re-view and to Ref.37for the theoretical basis of extending this theory to collisionless scales.Above the reversal scale,l≫l⊥,the cluster turbu-lence should resemble the saturated state of isotropic MHD turbulence:a magnetic-energy spectrum with a positive spectral index corresponding to foldedfields and a kinetic-energy spectrum populated in the inertial range by a peculiar type of Alfv´e n waves that propagate along the folds(i.e.,simultaneously perturbing the antiparallel magneticfield lines).58,59This type of turbulence is also reviewed in Ref.65.It is probably of limited relevance for clusters because the Reynolds number in the ICM is not large enough to allow a well-developed inertial range. Fig.5summarizes the—admittedly,rather specula-tive—picture of cluster turbulence proposed above.We offer this sketch in lieu of conclusions.While we believe that the set of physical arguments that has led to it is not without merit,it is clear that much analytical,numeri-cal and observational work is needed before a conclusion can truly be reached in the study of turbulence,magnetic fields and plasma physics in clusters of galaxies.AcknowledgmentsHelpful discussions with T.Enßlin,G.Hammett, R.Kulsrud,E.Quataert,and P.Sharma are gratefully acknowledged.This work was supported by a UKAFF Fellowship,a PPARC Advanced Fellowship,King’s Col-lege,Cambridge(A.A.S.)and by the DOE Center for Multiscale Plasma Dynamics. A.A.S.also thanks the Royal Society(UK)and the APS for travel support.∗Electronic address:as629@1A.C.Fabian,J.S.Sanders,S.W.Allen,C.S.Crawford, K.Iwasawa,R.M.Johnstone,R.W.Schmidt,and G.B.Taylor,Mon.Not.R.Astron.Soc.344,L43(2003).2A.C.Fabian,J.S.Sanders,G.B.Taylor,S.W.Allen,C.S.Crawford,R.M.Johnstone,and K.Iwasawa,Mon.Not.R.Astron.Soc.366,417(2006).3A.C.Fabian,J.S.Sanders,G.B.Taylor,and S.W.Allen,Mon.Not.R.Astron.Soc.(2005),in press(astro-ph/0503154).4P.Schuecker,A.Finoguenov,F.Miniati,H.B¨o hringer, and U.G.Briel,Astron.Astrophys.426,387(2004).5E.Churazov,W.Forman, C.Jones,and H.B¨o hringer, Astrophys.J.590,225(2003).6C.Vogt and T.A.Enßlin,Astron.Astrophys.412,373 (2003).7C.Vogt and T.A.Enßlin,Astron.Astrophys.434,67 (2005).8T.J.Dennis and B.D.G.Chandran,Astrophys.J.622, 205(2005).9L.M.Voigt and A.C.Fabian,Mon.Not.R.Astron.Soc.347,1130(2004).10A.C.Fabian,J.S.Sanders,C.S.Crawford,C.J.Con-selice,J.S.Gallagher,and R.F.G.Wyse,Mon.Not.R.Astron.Soc.344,L48(2003).11R.A.Sunyaev,M.L.Norman,and G.A.Bryan,Astron.Lett.29,783(2003).12N.A.Inogamov and R.A.Sunyaev,Astron.Lett.29,791 (2003).13E.Churazov,W.Forman, C.Jones,R.Sunyaev,andH.B¨o hringer,Mon.Not.R.Astron.Soc.347,29(2004). 14P.Rebusco,E.Churazov,H.B¨o hringer,and W.Forman, Mon.Not.R.Astron.Soc.359,1041(2005).15M.L.Norman and G.L.Bryan,Lect.Notes Phys.530, 106(1999).16P.M.Ricker and C.L.Sarazin,Astrophys.J.561,621 (2001).17Y.Fujita,Astrophys.J.631,L17(2005).18E.Churazov,M.Br¨u ggen,C.R.Kaiser,H.B¨o hringer,and W.Forman,Astrophys.J.554,261(2001).19K.Subramanian,A.Shukurov,and N.E.L.Haugen,Mon.Not.R.Astron.Soc.366,143(2006).20T.A.Enßlin and C.Vogt,Astron.Astrophys.(2006),in press(astro-ph/0505517).21B.D.G.Chandran,Astrophys.J.632,809(2005).22S.I.Braginskii,Rev.Plasma Phys.1,205(1965).23M.A.G.Willson,Mon.Not.R.Astron.Soc.151,1(1970). 24P.P.Kronberg,Rep.Prog.Phys.57,325(1994).25C.L.Carilli and G.B.Taylor,Ann.Rev.Astron.Astro-phys.40,319(2002).oni and L.Feretti,Int.J.Mod.Phys.D13,1549 (2004).27T.A.Enßlin and C.Vogt,Astron.Astrophys.401,835 (2003).28B.J.Burn,Mon.Not.R.Astron.Soc.133,67(1966). 29M.Haverkorn,P.Katgert,and A.G.de Bruyn,Astron.Astrophys.403,1045(2003).30T.E.Clarke and T.A.Enßlin,Astron.J.(2006),in press (astro-ph/0603166).31T.A.Enßlin,A.Waelkens,C.Vogt,and A.A.Schekochi-hin,Astron.Nachr.327,626(2006).32A.A.Schekochihin,S.C.Cowley,R.M.Kulsrud,G.W.Hammett,and P.Sharma,Astrophys.J.629,139(2005). 33M.N.Rosenbluth(1956),Los Alamos Scientific Labora-tory Report LA-2030.34S.Chandrasekhar,A.N.Kaufman,and K.M.Watson, Proc.R.Soc.London,Ser.A245,435(1958).35E.N.Parker,Phys.Rev.109,1874(1958).36A.A.Vedenov and R.Z.Sagdeev,Sov.Phys.Doklady3, 278(1958).37A.A.Schekochihin,S.C.Cowley,W.D.Dorland,G.W.Hammett,G.G.Howes,and E.Quataert,Astrophys.J.(2006),submitted.38V.D.Shapiro and V.L.Shevchenko,Sov.Phys.JETP18, 1109(1964).39K.B.Quest and V.D.Shapiro,J.Geophys.Res.101, 24457(1996).40D.Grasso and H.R.Rubinstein,Phys.Rep.348,163 (2001).41N.Y.Gnedin,A.Ferrara,and E.Zweibel,Astrophys.J.539,505(2000).42R.Banerjee and K.Jedamzik,Phys.Rev.Lett.91,251301 (2003).43G.K.Batchelor,Proc.R.Soc.London,Ser.A201,405 (1950).44Y.B.Zeldovich,A.A.Ruzmaikin,S.A.Molchanov,andD.D.Sokoloff,J.Fluid Mech.144,1(1984).45Y.B.Zeldovich,A.A.Ruzmaikin,and D.D.Sokoloff,The Almighty Chance(World Scientific,Singapore,1990).46V.I.Arnold and B.A.Khesin,Topological Methods in Hydrodynamics(Springer,Berlin,1999).47W.Jaffe,Astrophys.J.241,925(1980).48J.Roland,Astron.Astrophys.93,407(1981).49A.Ruzmaikin,D.Sokoloff,and A.Shukurov,Mon.Not.R.Astron.Soc.241,1(1989).50D.S.D.Young,Astrophys.J.386,464(1992).51O.Goldshmidt and Y.Rephaeli,Astrophys.J.411,518 (1993).52F.J.S´a nchez-Salcedo,A.Brandenburg,and A.Shukurov, Astrophys.Space Sci.263,87(1998).53K.Roettinger,J.M.Stone,and J.O.Burns,Astrophys.J.518,594(1999).54K.Dolag,M.Bartelmann,and H.Lesch,Astron.Astro-phys.348,351(1999).55K.Dolag,M.Bartelmann,and H.Lesch,Astron.Astro-phys.387,383(2002).56K.Dolag,D.Grasso,V.Springel,and achev,J.Cos-mol.Astropart.Phys.1,9(2005).57M.Br¨u ggen,M.Ruszkowski,A.Simonescu,M.Hoeft,andC.Dalla Vecchia,Astrophys.J.631,L21(2005).58A.A.Schekochihin,S.C.Cowley,G.W.Hammett,J.L.Maron,and J.C.MacWilliams,New J.Phys.4,84(2002). 59A. A.Schekochihin,S. C.Cowley,S. F.Taylor,J.L.Maron,and J.C.MacWilliams,Astrophys.J.612,276 (2004).60E.Ott,Phys.Plasmas5,1636(1998).61A.Schekochihin,S.Cowley,J.Maron,and L.Malyshkin, Phys.Rev.E65,016305(2002).62A.Brandenburg and K.Subramanian,Phys.Rep.417,1 (2005).63R.H.Kraichnan,Phys.Fluids8,1385(1965).64P.Goldreich and S.Sridhar,Astrophys.J.438,763(1995).。
武汉2024年06版小学六年级上册英语第六单元自测题考试时间:90分钟(总分:110)A卷考试人:_________题号一二三四五总分得分一、综合题(共计100题)1、填空题:I saw a ________ flying high in the sky.2、What is the name of the famous American author known for "The Adventures of Tom Sawyer"?A. Mark TwainB. Ernest HemingwayC. F. Scott FitzgeraldD. John Steinbeck答案: A3、填空题:The turtle hides in its _________. (壳)4、填空题:She is a lawyer, ______ (她是一名律师), who helps people.5、填空题:My ________ (玩具名称) is a treasure from my childhood.6、填空题:The __________ (历史的讨论) can lead to greater understanding.7、听力题:The chemical symbol for lithium is _______.8、填空题:The __________ (颜色) of a solution can indicate its chemical properties.9、听力题:The Great Red Spot is a giant storm on the planet ______.I bring a towel to dry off after swimming in the ______ (河).11、填空题:I enjoy watching ________ (比赛) at the stadium.12、What is the first month of the year?A. FebruaryB. JanuaryC. MarchD. April13、选择题:What do we call the act of taking care of someone?A. CaringB. NurturingC. SupportingD. Assisting14、听力题:The _____ (desk/table) is made of wood.15、What do you call the line separating day and night?A. EquatorB. HorizonC. MeridianD. Axis答案:B16、填空题:The elephant is known for its ______ (记忆).17、填空题:We have a ______ (愉快的) gathering for birthdays.18、How many months have 30 days?A. 4B. 5C. 6D. 7答案: A19、填空题:My pet hamster loves to explore its ______ (笼子).The chemical formula for arsenic trioxide is ______.21、听力题:We eat ______ (snacks) during recess.22、填空题:I find it ________ (有趣) to explore nature.23、What do you call a person who studies the effects of chemicals on living organisms?A. ChemistB. ToxicologistC. PharmacologistD. All of the above答案: D24、Which country is famous for kangaroos?A. CanadaB. AustraliaC. BrazilD. India答案: B25、community resource exchange) shares local assets. 填空题:The ____26、What is the capital of Mongolia?a. Ulaanbaatarb. Erdenetc. Darkhand. Choibalsan答案:a27、选择题:What do we call a person who studies the Earth?A. GeographerB. GeologistC. CartographerD. All of the above28、填空题:The ________ (杏树) in our garden gives us sweet fruit.29、What do you use to write on paper?a. Brushc. Scissorsd. Ruler答案:B30、What is the name of the famous mouse created by Walt Disney?A. TomB. JerryC. MickeyD. Donald答案:C31、听力题:She _____ (enjoy/enjoys) dancing.32、What do we call a young kangaroo?A. CalfB. JoeyC. KidD. Pup33、听力题:The chemical symbol for tungsten is ________.34、听力题:A rabbit's strong back legs help it to ______ quickly.35、听力题:I want to ___ an adventure. (have)36、听力题:The ____ is known for its agility and speed.37、填空题:A ____(refugee) flees their home due to conflict or danger.38、What do you call the process of removing the outer layer of fruit?A. SlicingB. PeelingC. ChoppingD. Dicing39、What is the weather like when it snows?A. HotB. ColdC. Humid答案:B40、填空题:The ______ (根系) structure of plants is complex and important.41、选择题:Which animal is known for its long neck?A. ElephantB. GiraffeC. ZebraD. Rhino42、What shape has four equal sides?A. TriangleB. RectangleC. SquareD. Circle答案: C43、填空题:The capital of Bangladesh is _____.44、听力题:My ______ loves to engage in discussions.45、听力题:A saturated fat is typically found in ______.46、填空题:Herbs are often grown in ______ (窗台) pots.47、填空题:The _____ (pistachio) tree produces nuts.48、填空题:The _______ (老虎) is striped.49、填空题:The three states of matter are solid, liquid, and _______. (气体)50、听力题:The nucleus of an atom contains protons and _____.51、What is the chemical symbol for water?A. H2OB. O2C. CO2D. NaCl答案:A52、Which planet is known as the Red Planet?A. EarthB. MarsC. JupiterD. Venus答案: B53、听力题:The ________ (activity) enhances learning.54、填空题:A ________ (冰川) can create valleys as it moves.55、What do you call a person who flies airplanes?A. PilotB. NavigatorC. EngineerD. Steward56、听力题:The study of Earth’s history helps us understand ______ changes.57、填空题:The __________ is a famous city known for its ancient history. (开罗)58、听力题:A mixture that does not settle over time is called a ______.59、填空题:The capital of Cuba is ________ (哈瓦那).60、填空题:We planted ________ in our garden.61、听力题:Plastics are made from long chains of _____, called polymers.62、填空题:When a solid dissolves in a liquid, it is called _______. (溶解)The ________ was a time of great change in Europe.64、What is the capital city of Hungary?A. BudapestB. DebrecenC. SzegedD. Miskolc65、How many colors are in a standard rainbow?A. 5B. 6C. 7D. 8答案:C66、听力题:They are _____ (watering) the plants.67、听力题:A solution where no more solute can dissolve is called ______.68、听力题:The main gas we breathe is ______.69、选择题:What is 6 + 4?A. 10B. 11C. 12D. 1370、What is the capital of Argentina?A. SantiagoB. MontevideoC. Buenos AiresD. Lima答案:C71、填空题:The ______ (植物的营养需求) varies greatly between species.72、听力题:We are going ________ a trip.The __________ (古埃及的象形文字) was used for writing on temple walls.74、Which animal is known for its ability to change colors?a. Chameleonb. Dogc. Catd. Rabbit答案:a75、填空题:_____ (土壤) quality is essential for plant growth.76、听力题:He is playing with his ___. (friends)77、What color is the sky?天空是什么颜色的?A. RedB. BlueC. GreenD. Yellow答案: B78、填空题:The flowers in the garden are _______ and cheerful, spreading happiness.79、What is 25 ÷ 5?A. 3B. 4C. 5D. 6答案:C80、听力题:A __________ is a natural elevation of the Earth's surface.81、What is the hardest natural substance on Earth?A. GoldB. IronC. DiamondD. Silver答案:C82、填空题:The __________ is a large archipelago in Southeast Asia. (印度尼西亚)Glacial deposits can create __________ features.84、How many days are there in a leap year?A. 364B. 365C. 366D. 367答案:C. 36685、What do you call a book with stories?A. NovelB. DictionaryC. EncyclopediaD. Journal86、听力题:The Earth's crust is divided into continental and ______ crust.87、听力题:The ____ has a unique way of moving and can hop very high.88、听力题:The Earth's surface is influenced by natural and ______ factors.89、ts can live in ______ (盐碱地). 填空题:Some pla90、选择题:What is the capital of Chile?A. SantiagoB. ValparaisoC. ConcepcionD. La Serena91、How many legs does a typical arachnid have?A. SixB. EightC. TenD. Twelve答案: B92、填空题:The dolphin jumps high out of the _________. (水)93、Which fruit is red and round?A. BananaB. AppleC. OrangeD. Grape答案:B94、填空题:The ______ (狐狸) is very clever and sly.95、听力题:The sky is _______ (clear) with no clouds.96、听力题:The ________ (lobster) is a type of seafood.97、What is the color of a typical mint leaf?A. GreenB. YellowC. BlueD. Red答案:A98、What is the capital of Somalia?a. Mogadishub. Hargeisac. Kismayod. Baidoa答案:a99、填空题:The _____ (草原) is home to many wildflowers.100、听力题:__________ change involves a change in state, not in composition.。
Magnetic D evices for a Beam Energy RecoveryTHz Free Electron LaserR. R. S. Caetano¹, G. Cernicchiaro², R. M. O. Galvão³1 2Universidade Federal do Rio de Janeiro,Macaé, R J, Bra z il, rcaetano@macae.ufrj.brCoordenação de Física Aplicada, Centro Brasileiro de Pesquisas Físicas, Rio de Janeiro, RJ, Bra z il, geraldo@cbpf.br3Instituto de Física USP, São Paulo, SP, Bra z il.Abstract:This paper presents a numerical analysis of magnetic devices, dipole, quadrupole and undulator and a THz Free Electron Laser (FEL) electron-beam recovery system. Free Electron Laser are an important source of coherent radiation being used in the study of chemical properties of substances, thus being an important tool for various fields of science such as condensed matter physics, chemistry, biology and medicine. The magnetic device of this simulation is to contribute to the proposed deployment of a national laboratory for multiuser application and development of a recovery system FEL to operate in the far infrared range between 0.3 to 1.2 THz.Keywords: Free Electron Laser THz, magnetic device, simulation in COMSOL.1.IntroductionFree electron laser can be used are in spectroscopy thus have applications in scientific fields such as medicine, chemistry, condensed matter and biology [1]. Thus the Brazilian Center for Physics Research (CBPF) proposed a construction project of a Free Electron Laser (FEL) using the components of a Free Electron Laser of the College of Optics & Photonis (CREOL). The free-electron laser operating in the far infrared range working with wavelengths in the range of 200-600 micrometers. The equipment consists of a linear accelerator that generates electrostatic energy up to 1.7 MeV, a magnetic undulator, which is designed with permanent magnets made from neodymium iron boron (NdFeB) with 185 periods with a wavelength of 8 mm length of 1486 mm and a distance between the undulator cassette (gap) of 6 mm. Has magnetic dipoles and quadrupoles working in optics from the electron beam. The vacuum system is formed by mechanical, ionic pumps and turbomolecular [2], control is being developed for the LabView platform.This paper focuses on the numerical analysis using the COMSOL software of magnetic elements that are the dipole, quadrupole and undulator. These components work changing the trajectory and the size of the electron beam.2.Theory2.1 DipoleT he magnetic dipole is an element that has the function deflect the electron beam [3]. The dipole magnetic field is generated by the electric current passing through the coils. This current in solenoid generates a magnetic flux in the iron core creating a unidirectional magnetic field in between irons given by the right hand rule. The magnitude of the magnetic field can be extracted by Ampere's law gives us the relationship between the electric current and the magnetic field.Figure 1: Magnetic dipole.The dipole magnetic field is given by equation 1 where n is the number of turns, the electrical current I, is the distance between the irons and is the permeability of air.( 1)2.2 QuadrupoleT his component consists of four poles formed by rectangular hyperbolas with alternating magnetic fields, your objective is to change the diameter of the electron beam. The geometry of the quadrupole makes the magnetic field in the core is zero and the magnetic field is generated by modulo magnitude of the field in the x-axis and y-axis [3].Figure 2: Magnetic quadrupole.The quadrupole magnetic field can be calculated from the integration by Ampere's law.( 2)The magnetic field consists of two quadrupole components in the x and y axis, there is a resulting gradient g [T/ m]. Thus, the magnetic field in the x and y axes can be given by the equations , where , and the equation of the resulting magnetic field is equal to:( 3)2.3 UndulatorThe undulator is a mechanical structure consisting of periodic magnets with alternating poles, separated by a distance called (gap). These magnets are made of magnetic material pure permanent magnetic (PPM). This structure causes synchronous radiation be concentrated in a beam, thus reducing the radiation loss.Figure 3: Description of the undulator magnetic.The magnetic undulator field is perpendicular to the x axis, and the direction of the z-axis. Because the poles are alternating magnetic field in the undulator has a sinusoidal behavior as can be seen in equation 4. [4]. Where is the wavelength of the undulator and is the initial magnetic field. In equation 5 we have the initial magnetic field which is a ratio between the gap and the wavelength.(4)(5) e of COMSOL MultiphysicsThe numerical analysis of the dipole, quadrupole and undulator were built using the module AC/DC COMSOL multiphysics. To calculate the magnetic field in the dipole and quadrupole magnetic field the tool was used with the equation 6. In the undulator, the calculation of the initial magnetic field was carried out with the magnetic field in the current tool with the equation 7.(6)(7)3.1 DipoleThe geometry of the dipole was constructed in 3D and has two types of materials are iron, which comprises all the dipole and charge air by the cylindrical shell. Infinite element domain was used for the cylindrical shell for emulating an infinite open space causing all numerical analysis considers the limited space as being infinite.Figure 4. 3D geometry of the magnetic dipole.The interface was used for the dipole Magnetic Field <Ampere Law for the iron structure and for calculating coil was used in Magnetic Field <muti-turn coil <Coil Current Calculation which calculate the magnetic field in the coil according to the current electrical and number of turns. The dipole has 877 turns and the maximum current is 2.5 A.Figure 5. Coil Current Calculation in dipole magnetic.Dipole in the fine mesh, which corresponds to 61114 domain element, boundary element 9728, 1039 edge element was used. The study used for the calculation of the dipole was the Stationary Parametric Sweep and the electric current was applied in order to have values of -2.5 A to +2.5 A.Figure 6: mesh magnetic Dipole.3.2 QuadrupoleThe quadrupole was developed in 3D for its construction it was done primarily in 2D and using the Work Plane tool, after it was extruded to 3D format.Figure7. Quadrupole magnetic geometric in 2D.Figure 8. 3D geometry of the magnetic quadrupole.The materials used in the quadrupole are air and iron and infinite element domain was appliedto the cylindrical shell. To calculate the magneticfield of the quadrupole physical Magnetic Field <Ampere Law, which calculates the magnetic field from the magnetization and uses equation was used , where N is the number of turns, I is the electric current, A is the area and V is the volume of the coil. In the case of quadrupole have 344 turns and the electrical current worth 2.5 A. The magnetization in the quadrupole is oriented in the x and z axes have that each pole has a different orientation.Quadrupole in the mesh extremely fine,which corresponds to 749885 domain element, boundary element 58352, 2772 edge element was used. The study used for the calculation of the quadrupole was the Stationary Parametric Sweep and the electric current was applied in order to have values of -2.5 A to +2.5 A.Figure 9: mesh magnetic quadrupole.3.3 UndulatorFor this work was a session of 3D undulator built with 400 mm in length and can change the wavelength and the distance of the gap. Undulator possessed these variations to resemble the original equipment. The initial magnetic field generated in the undulator does not depend on the distance of the gap but wavelength (equation 5)and thus it is possible to make a comparison between the simulation and the experiment.The materials used were a NdFeB magnetic material is (33SH) and those manufactured with the air gaps 1010 steel, their averages are: x = 10.5 mm, y = 1.3 mm and z = 25 mm and x = 30mm, y = 2 7 mm and z = 13 mm, respectively [5]. On the external surface that is a spherical vacuum volume 220 mm radius of 5mm layer was used.Figure 10. 3D geometry of the magnetic Undulator.T he magnetic field of the undulator was obtained using the Magnetic Field, No current Physical<Magnetic Flux. When applying a magnetic field in a material, the resulting field B is the sum of the applied field H and the field of the magnetized material, as in Equation 8.( 8)Magnetic materials have hysteresis magnetization curves called reduce to zero the magnetic field applied, as can be seen in equation 9 [6].( 9)The magnetization is added in Magnetic flux Conservation, as the magnetization of the magnets are oriented alternately is necessary to indicate the direction of magnetization in this case, the direction is along the x axis. The mesh used was extra fine, possessing 5842814 domain boundary elements 681 250 and 80 793 edge andthe study used was stationary.Figure 11. Mesh Undulator magnetic.6.0 Results 6.1 DipoleThe simulation dipole was constructed to compare the value of the modulus of the magnetic field with the experimental value in Figures 12 and 13 have the simulation results.Figure 12. Magnetic Flux density in an dipolemagnetic device.Figure 13. Graph magnetic field x current electric indipole.Table 1: Comparison of the magnetic field obtained from the simulation and experimental analysis in the dipole magnetic device.The magnetic field values are obtained for the maximum value of electric current of 2.5 A. Thus the table 1 it can be seen that the threevaluesarecloseto a percentage error of 3% between experiment and COMSOL. 6.2 QuadrupoleThe simulation quadrupole was constructed to compare the value of the modulus of the magnetic field with the experimental value in Figures 14 and 15 have the simulation results.Figure 14. Magnetic Flux density in an quadrupolemagnetic device.Figure 15. Graph magnetic field x current electric inquadrupole.Table 2: Comparison of the magnetic field obtained from the simulation and experimental analysis in the quadrupole magnetic device.The magnetic field values are obtained for the maximum value of electric current of 2.5 A. Thus the table 2 it can be seen that the three values are close to a percentage error of 1,9 % between CREOL and experiment.6.3 UndulatorThe simulation study of the undulator was developed for the x axis and y axis. In Figures 16, 17, 18 and 19 have a magnetic field initial described in the two axes respectively.Figure 16. Magnetic Flux density in an undulator magnetic device (front view). Figure17. Magnetic field of the undulator gap measured in the z axis direction with respect to distance in the y direction.Figure 18. Magnetic Flux density in an undulator magnetic device (side view).In Figure 19 we have the magnetic field in the undulator which is a sine function according to equation 10, so this simulation is in agreement with theory. In Table 3it canbe observedthat the values stipulated by the project, experiment and COMSOL, so it is possible that there is an error 1% compared to COMSOL and experiment.Figura 19.Undulator magnetic field measured at the gap in the z direction as a function of distance in the x direction.Table 3: Comparison of the magnetic field obtained from the simulation and experimental analysis in the undulator magnetic device.7.0 ConclusionIn this paper we present 3D simulations of the dipole, quadrupole and undulator magnetic elements that are owned by a Thz Free Electron Laser. The results have to numerical simulation results are in agreement with experiment validating the paper. To the next module using the practical tracing in COMSOL, an electron beam will be added with the aim of studying the behavior of the electron beam in magnetic elements.8 Reference[1] S rinivas Krishnagopal*, Vinit Kumar†,Sudipta Maiti, S. S. Prabhu and S. K. Sarkar, "Free-electron lasers.," CURRENT SCIENCE, VOL. 87,, pp. NO. 8, 25, OCTOBER 2004.[2] M. Tecimer, Time –Domanin analysis andtechnology of THz Free Electron Lasers, Tel –Aviv University. : Faculty of Engineering.Departamento f Electical Engineering –Physical Electronics. , 2005.[3] J. Tanabe, Iron Dominated ElectromagnetsDesign, Fabrication, Assembly and Measurements, 2006.[4] M. D. J. R. Peter Schmüser, Ultraviolet andSoft X-Ray Free-Electron Lasers - Introduction to Physical Principles, Experimental Results, Technological Challenges, Springer, 2009.[5] J. G. L. R. E. Paul P. Tesch, "Finalconstruction of the CREOL 8 mm period hybrid undulator," Nuclear Instruments and Methods in Physics Research A375, pp. 504-507, 1996. [6] R. N. Faria, L.F.C.P. Lima, Introdução aomagnetismo dos materiais, São Paulo: Livraria da Fisica, 2005.。
&KDSWHU+DOO (IIHFW 6HQVRUVIntroductionThe Hall effect was discovered by Dr. Edwin Hall in 1879 while he was a doctoral candidate at Johns Hopkins University in Baltimore. Hall was attempting to verify the theory of electron flow proposed by Kelvin some 30 years earlier. Dr. Hall found when a magnet was placed so that its field was perpendicular to one face of a thin rectangle of gold through which current was flowing, a difference in potential appeared at the opposite edges. He found that this voltage was proportional to the current flowing through the conductor, and the flux density or magnetic induction perpendicular to the conductor. A l-though Hall’s experiments were successful and well received at the time, no applications outside of the realm of theoretical physics were found for over 70 years.With the advent of semiconducting materials in the 1950s, the Hall effect found its first applications. However, these were severely limited by cost. In 1965, Everett Vorthmann and Joe Maupin, MICRO SWITCH Sensing and Control senior d e-velopment engineers, teamed up to find a practical, low-cost solid state sensor. Many different concepts were examined, but they chose the Hall effect for one basic reason: it could be entirely integrated on a single silicon chip. This breakthrough resulted in the first low-cost, high-volume application of the Hall effect, truly solid state keyboards. MICRO SWITCH Sensing and Control has produced and delivered nearly a billion Hall effect devices in keyboards and sensor products.Theory of the Hall EffectWhen a current-carrying conductor is placed into a magnetic field, a voltage will be generated perpendicular to both the current and the field. This principle is known as the Hall effect.Figure 2-1 illustrates the basic principle of the Hall effect. It shows a thin sheet of semicon-ducting material (Hall element) through which a current is passed. The output connections are perpendicular to the direction of current. When no magnetic field is present (Figure 2-1), current distribution is uniform and no potential difference is seen across the output.When a perpendicular magnetic field is present,as shown in Figure 2-2, a Lorentz force is exerted on the current. This force disturbs the currentdistribution, resulting in a potential difference (voltage) across theoutput. This voltage is the Hall voltage (V H ). The interaction of the magnetic field and the current is shown in equation form as equ a-tion 2-1.Hall effect sensors can be applied in many types of sensing devices. If the quantity (parameter) to be sensed incorporates or can incorporate a magnetic field, a Hall sensor will perform the task.V H ∝ I × BFormula (2-1)Figure 2-4 Basic Hall effect sensorA differential amplifier with these characteristics can be readily integrated with the Hall element using standard bipolar transistor technology. Temperature compensation is also easily i n tegrated.As was shown by equation 2-1, the Hall voltage is a function of the input current. The purpose of the regulator in Figure 2-4 is to hold this current constant so that the output of the sensor only reflects the intensity of the magnetic field. As many systems have a regulated supply available, some Hall effect sensors may not include an internal regul a tor.Figure 2-12 NPN (Current sinking) . . . Digital output sensor。
The magnetic field of pulsars and the gravito-magnetictheoryJacob Biemond*Vrije Universiteit, Amsterdam, Section: Nuclear magnetic resonance, 1971-1975AbstractMany authors have considered a gravitational origin of the magnetic field ofcelestial bodies. Especially, the so-called Wilson-Blackett formula has beeninvestigated, both theoretically and observationally. It appeared possible to deducethis formula from general relativity, e. g., by application of a special interpretationof gravito-magnetic theory. More consequences of the latter theory for pulsarswill be considered in this work.As an example, the standard quadrupolar charge density for pulsars can bededuced from the gravito-magnetic theory. Moreover, a new magnetic dipolemoment from electromagnetic origin is found, generated in the basic magneticfield from gravito-magnetic origin.In general, for thirteen accreting, slowly rotating, binary pulsars theagreement between the observed magnetic field and the gravito-magneticprediction is better than between the observed value and the value from thestandard magnetic dipole radiation model. At present, an analogous comparisonfor five isolated pulsars appears to be difficult.For a sample of 100 pulsars the averaged (gravito-)magnetic field, extractedfrom the magnetic dipole spin-down model, is in fair agreement with the gravito-magnetic prediction. Unfortunately, the (gravito-)magnetic field has not yetdirectly been measured.Finally, it is found that the first and second order braking indices onlydepend on the magnetic field from electromagnetic origin.1. IntroductionSince 1891 many authors have discussed a gravitational origin of the magneticfield of rotating celestial bodies. Especially, the so-called Wilson-Blackett formula wasconsidered [1–8]M = – 1/2 βc–1 G1/2 S. (1) Here M is the magnetic dipole moment of the massive body with angular momentum S, cis the velocity of light in vacuum, G is the gravitational constant and β is a dimensionlessconstant of order unity. Attempts to derive equation (1) from a more general theory havebeen given by several authors [8–14]. Equation (1) has often been considered as a conse-quence of general relativity [12–14], for example, in a special version of the gravito-magnetic theory [8]. In the latter approach the so-called “magnetic-type” gravitationalfield is identified as a common magnetic field, resulting into the magnetic dipole momentM(gm) of (1).Experiments on rotating masses in the laboratory in order to test (1) have beenperformed by Blackett [15] and others [16–18]. Available observations and theoreticalconsiderations with respect to the relation (1) and other explanations of the origin of themagnetic field of celestial bodies have been reviewed by Biemond [8]. FollowingWoodward [19], the validity of the Wilson-Blackett (or Schuster-Blackett) formula forpulsars has again been investigated in this work, but now from the gravito-magnetic pointof view._______________________________________________________________________________________*Postal address: Sansovinostraat 28, 5624 JX Eindhoven, The Netherlands.The angular momentum S for a sphere of radius R can be calculated from the relationS = I Ω, or S = I Ω = 2/5 f m R2Ω, (2)where m is the mass of the sphere, Ω = 2πP–1 its angular velocity (P is the rotationalperiod), I its moment of inertia and f is a dimensionless factor depending on the homo-geneity of the mass density in the sphere (for a homogeneous mass distribution f = 1).The value of the magnetic dipole moment M may be calculated fromM = ½ r3 B p, or M = ½ r3B p . (3)Here B p is the magnetic induction field at the pole of the sphere and r is the distance fromthe centre of the sphere to the field point where B p is measured. The magnetic moment Mhas been derived for r > R, but in this work values of B p measured at r≈R have beenintroduced into (3).Combination of (1)–(3) yields the following gravito-magnetic prediction for B pB p(gm) = – β c–1 G1/2I R–3Ω, or B p(gm) = – 5.414x1013P–1, (β = + 1) (4)where the minus sign reflects that the vectors B p(gm) and Ωpossess opposite directions forβ = + 1. Neither the sign nor the value of β follow from the gravito-magnetic theory. There-fore, absolute values of B p(gm) are given in this work. The representative values I = 1045g.cm2 and R≈r = 10 km have been inserted in (4) and in other formulas in this work. It isnoticed that (4) imposes no restriction on the maximum of B p(gm), i.e., B p(gm) may belarger than the critical field strength m e2c3/eħ = 4.414x1013 Gauss, at which electron-positron pair creation processes may become probable.As pointed out earlier [8, pp. 12, 14, 19–20 and 49], moving electric charge in themagnetic field from gravito-magnetic origin may cause an additional magnetic field fromelectromagnetic origin. The latter field may partly or completely compensate the magneticfield from gravito-magnetic origin. It is noticed that the magnetic field generated byrotating neutral mass is generally much smaller than the magnetic field generated bymoving charge. For example, for the electron one may compare the following magneticmoment to angular momentum ratios: (M/S)gravito-magnetic and (M/S)electromagnetic. Calculationshows: (M/S)gravito-magnetic/(M/S)electromagnetic = G1/2m e e–1 = 4.899x10–22. Therefore, as a rule,proposed magnetic fields from gravito-magnetic origin are extremely small and difficultto isolate from fields of electro-magnetic origin.When a magnetic induction field B(gm)due to gravito-magnetism is present, thetotal magnetic induction field B(tot) is given byB(tot) = B(gm) + B(em), (5) where B(em) is the magnetic induction field due to electromagnetism. It may be helpfulto define the dimensionless quantities βeff and β′ by the relationsB(tot) ≡βeff B(gm), and B(em) ≡β′B(gm). (6)Since the signs of βeff and β′ are unknown, absolute values of these quantities have beengiven in all tables below. Combination of (5) and (6) yieldsB(em) = (βeff – 1) B(gm), and βeff = 1 + β′. (7)When the total field B(tot) is only due to gravito-magnetic origin, B(em) = 0 and βeff and β′reduce to βeff = 1 and β′ = 0, respectively.Since charges may move in different ways in pulsars and other rotating bodies, onecan hardly expect that βeff is a constant. Indeed, different results for βeff have been foundfor about fourteen rotating bodies: metallic cylinders in the laboratory, moons, planets,stars and the Galaxy. However, for this series a mean value of βeff = 0.076 has been calcu-lated [8]. Although the agreement with the gravito-magnetic prediction βeff = 1 is onlyapproximate, the correct order of magnitude of βeff for so many, strongly different , rota-ting bodies is amazing (the values of the parameters m , R , Ω and B p occurring in (4) (see(1) and (6a)) can differ by many decades ). Such a result may reflect the validity of thegravito-magnetic hypothesis.In general, many properties of pulsars are determined by the emission model used.Therefore, we will first pay attention to some proposed emission models in section 2. Insections 3 and 6 the influence of the magnetic field on first and second order brakingindices is considered, both theoretically and observationally, respectively. Apart from theWilson-Blackett formula, other consequences of the gravito-magnetic theory for pulsarsare deduced in section 4. In section 5 observational data for different kinds of pulsars aresummarized and compared with theoretical predictions. A summary of the conclusions isgiven in section 7.2. Emission models of pulsarsThe emission model of pulsars is an important factor in the prediction of the mag-netic field of pulsars. Unfortunately, no generally accepted model is available. Therefore,we will discuss some alternative formulas for the emission of pulsars.The rotational energy E = ½ I Ω2 of the pulsar may change by magnetic dipole radi-ation, accretion or other mechanisms according toĖ = I Ω Ωֹ = – 4π2 I P –3 Pֹ , (8) where Ė, Ωֹ and Pֹ are the time-derivatives of E , Ω and P , respectively (It has been assumed that the moment of inertia I is no function of time). The contribution of the(electro)magnetic dipole radiation in vacuum to Ė of (8) is given by (see, e. g., [20, pp.188–189; 21, pp. 176–177])Ė = – ⅔ c –3 sin 2θ M (em)2 Ω4 = – (32π4/3) c –3 sin 2θ M (em)2 P –4, (9)where M (em) is the (electro)magnetic dipole moment of the pulsar and θ is the anglebetween M (em) and S . Note that the (gravito-)magnetic dipole moment M (gm) does notoccur in (9), since (gravito-)magnetic dipole radiation does not exist. Gravito-magneticquadrupole radiation (for two point masses its energy formula coincides with the familiarexpression for the gravitational quadrupole radiation (see, [8, ch. 3; 20]) may be present,however, but its influence is not considered in this work. When Ė in (8b) is taken equal to Ė in (9b), it follows that Pֹ is positive . So, magnetic dipole radiation leads to spin-down . Combination of (8b) and (9b) further yields forM (em)⎛3⎞1/2 3 c I M (em) = ⎜⎯⎯⎯⎯ ⎟ (P P ֹ )1/2. (10) ⎝ 8π2 sin 2θ ⎠Choosing θ = 90º in (10), the following magnetic induction field at the pole of the pulsarcan be calculated from (3b)3 ⎜⎯⎯⎯⎟ (P P ֹ )1/2 = 6.399x1019 (P P ֹ )1/2. (11) ⎛c 3I ⎞1/2 B p (sd) = ⎝ 2π2 R 6 ⎠It is noticed that, as a rule, θ < 90º, so that B p (sd) in (11) must be considered as the lowerlimit.Assuming deformation of magnetic field lines, or outflow of high-energy chargedparticles, approximations to Ė may be obtained by taking the product of the energydensity of the equatorial (electro)magnetic field at the light cylinder B L (em), the effectivearea of the light cylinder and the velocity of light (see, e. g., [21, pp. 188–189])B (em)2 L Ė = – ⎯⎯⎯ 4π R L 2 c , (12) 8πwhere R L = c /Ω is the radius of the light cylinder. B L (em) is taken equal toB L (em) = ½ B p (em) (R /R L )p , (13)where B p (em) is the (electro)magnetic field at the pole of the pulsar. Using (3b), combi-nation of (12) and (13) leads toĖ = – ⅛ c 3– 2p R 2p B p (em)2 Ω2p – 2 = – ½ c 3– 2p R 2p – 6 M (em)2 Ω2p – 2. (14)For example, p = 3 yieldsĖ = – ½ c – 3 M (em)2 Ω4. (15)Although the derivations of (9a) and (15) start from different assumptions, both equationsnearly coincide, apart from a different coefficient and a factor sin 2θ. Since the derivationof (9a) is more rigorous than that of (15), the former expression is usually used as thestandard formula for the determination of the (electro)magnetic field of pulsars.When p = 2, equation (14) yields for ĖĖ = – ½ c – 1 R – 2 M (em)2 Ω2 = – ⅛ c – 1 R 4 B p (em)2 Ω2. (16)Combination of (8b) and (16b) then leads to the following expression for the deformedmagnetic induction field, denoted by B p (def)⎛8 c I ⎞1/2 ⎛ P ֹ ⎞1/2 ⎛ֹ ⎞1/2 P B (def) = ⎜⎯⎯ ⎟ ⎜⎯⎟ = 1.549x1016 ⎜p ⎝ R 4 ⎠ ⎝ P ⎠ ⎝ P ⎠ ⎯⎟ . (17)The calculated values of B p (gm) of (4), B p (sd) of (11) and B p (def) of (17) will be com-pared with each other and with available observational data below.3. Magnetic field dependence of braking indicesIn order to calculate so-called braking indices, equations (8a) and (14a) are com-bined toΩֹ = – ⅛ c 3– 2p I –1 R 2p B p (em)2 Ω2p – 3. (18)Alternatively, (18) may be generalized toΩֹ = – k p c 3– 2p I –1 R 2p B p (em)2 Ω2p – 3, (19)where k p is a constant factor. In the sequel of this work it will be assumed that thequantity p is no function of time. Direct combination of (3b), (8a), (9a) and (19)shows that k 3 = 1/6 sin 2θ for p = 3.Instead of using (18) and (19), Kaspi et al . [22] and Johnston and Galloway [23]followed Manchester and Taylor [21] and considered the generalized relationΩֹ = – K Ωn 0 = – k B p (em)2 Ωn 0, (20)where the quantities k and n 0 are assumed not to depend on time. From (18)–(20) the quantity Ωֹ ֹ and the so-called first order braking index n can becalculated by differentiation (compare with the integration method discussed in section 6).From (18) and (19) one obtainsΩ Ωֹ ֹ B ֹ (em) Ω p n ≡ ⎯⎯ = 2p – 3 + 2 ⎯⎯⎯⎯(21) , Ωֹ 2 B p (em) Ωֹwhere the first order braking index n is defined in terms of the observable quantities Ω, Ωֹ and Ωֹ ֹ . In deriving (21) it has been assumed that none of the following parametersdepend on time : the moment of inertia I , the radius R , the mass M , or the angle θ of thepulsar. For convenience sake, these assumptions will be also used in the sequel of thispaper. Note that for a fixed value of p the braking index n will become smaller for B p (em) < 0 and Bֹ p (em) < 0, since Ω > 0 and usually Ωֹ < 0. From (20b) follows, analogously to (21)Ω Ωֹ ֹ B ֹ (em) Ω p n ≡ ⎯⎯ = n + 2 ⎯⎯⎯⎯ . (22) 0 Ωֹ 2 B p (em) ΩֹComparison of (21) and (22) shows that n 0 = 2p – 3. Note that, in principle, n and n 0 neednot to be integers.It is noticed that the braking indices n of (21) and (22) neither directly depend on the gravito-magnetic field B p (gm) nor on Bֹ p (gm). Therefore, the validity of the gravito-mag-netic hypothesis cannot been tested by considering the braking index n .In the evaluation of (18)–(20) it is often assumed that B p (em) is no function of time .In that case n in (21) reduces to n = 2p – 3, whereas (22) simplifies to n = n 0.In addition, one may identify the field B p (em) in (22) with the field B p (sd) from (11) that depends on the quantities Ω and Ωֹ . Differentiation of B p (em) = B p (sd) with respect to time followed by evaluation yields for the ratio Bֹ p (sd) to B p (sd) B ֹ (sd) Ωֹ 1 p ⎯⎯⎯ = ½ (n – 3) ⎯ = ¼ (3 – n ) –. (23) B p (sd) Ω τc Here the quantity τc is defined by τc ≡ – ½ Ω/ Ωֹ = ½ P /P ֹ . When n = 3, it follows from (23) that B p (sd) is independent of time, but for n ≠ 3 Bֹ p (sd) ≠ 0. When n < 3 and B p (sd) < 0, B p (sd) becomes more negative in time, whereas for n > 3 and B p (sd) < 0 B p (sd) becomesless negative in time.Combination of (22) and (23a) shows that n 0 = 3. Therefore, the quantity n 0 hasbeen called the true braking index (see, e. g., [22]). Moreover, from the relation n 0 = 2p – 3follows p = 3 in this case. Summing up, (electro)magnetic dipole radiation implies n 0 = 3and p = 3, but the value of n is determined by the time dependence of B p (sd) (or more precisely by the ratio B ֹ p (sd)/B p (sd): see equation (23)). Note that in favourable cases the braking index n can be calculated from n = Ω Ωֹ ֹ /Ωֹ 2, where n need not to be an integer.When B p (em) is identified with B p (sd), it is appears possible to deduce anexpression for the factor β′ in (6b) by combining (4), (6b) and (11)β′ = (– 6 c 5 G –1 I –1 Ω– 5 Ωֹ )1/2. (24)This result shows that β′ depends on the quantities Ω and Ωֹ , and on the particular choiceB p (em) = B p (sd).Another possibility is to identify B p (em) with B p (def) from (17). Differentiation ofB p (em) = B p (def) with respect to time yieldsB ֹ (def) Ωֹ 1 p ⎯⎯⎯ = ½ (n – 1) ⎯ = ¼ (1 – n ) –. (25) B p (def) Ω τcWhen in this case n = 1, it follows from (25) that B p (def) is independent of time. When n< 1 and B p (def) < 0, B p (def) becomes more negative in time, whereas for n > 1 andB p (def) < 0 B p (def) becomes less negative in time. Combination of (22) and (25a) showsthat now n 0 = 1. From the relation n 0 = 2p – 3 then follows that p = 2. Summing up, in thecase of deformation of magnetic field lines, or outflow of high-energy charged particles,n 0 and p may have the values n 0 = 1 and p = 2. Again the observed braking index n neednot to be an integer.In table 1 some results for n and n 0 from (21) and (22) with respect to the timedependence of B p (em) are summarized. See the text for further information.Table 1. Braking indices n and n 0 depending on the emission model (p -value) and thefield B p (em).B p (em) ≠ f(t ) B p (em) = B p (sd) B p (em) = B p (def) General casep = pn = 2p – 3 n = n 0 n 0 = 2p – 3 n 0 = 2p – 3 Magnetic dipole radiationp = 3n = n = 3 0 n 0 = 3 – Deformation of field linesp = 2 n = n = 1 0 –n 0 = 1In order to get an idea of the evolution a pulsar, the true age t of a pulsar is animportant parameter. If B p (em) is no function of time , integration of (18) or (19) yields thefollowing expression for this quantity (compare with, e.g., [21, p. 111])– Ω ⎧ ⎛ Ω ⎞n – 1 ⎫ – Ω 2 τ c t = ⎯⎯⎯⎯ ⎨1 – ⎜⎯ ⎟ ⎬ ≈ ⎯⎯⎯⎯ = ⎯⎯ , (n ≠ 1) (26) (n – 1) Ωֹ ⎩ ⎝Ωi ⎠ ⎭ (n – 1) Ωֹ n – 1 where Ωi is the initial value of Ω at time t = 0 and τc ≡ – ½ Ω/Ωֹ = ½ P /P ֹ . According to (21), n then equals to n = 2p – 3 and neither n nor p need not to be integers. From (26b) itcan be seen that the characteristic time τc may be used as an estimate for t . Therefore, τc has been added to tables 2 through 4 when P and Pֹ are known.Following Kaspi et al . [22], a second order braking index m can be calculated from(20b). Taking the time dependence of B p (em) into account, one findsΩ2 Ωֹ ֹ ֹ B ֹ (em) Ω B ֹ (em)2 Ω2 B ֹ ֹ (em) Ω2p p p m ≡ ⎯⎯ = 2n 2 – n + 6n ⎯⎯⎯⎯ 2 ⎯⎯⎯⎯⎯ 2 ⎯⎯⎯⎯⎯(27) 000 + + . Ωֹ 3 B p (em) Ωֹ B p (em)2 Ωֹ 2 B p (em) Ωֹ 2The quantity m has been defined in terms of the observable quantities Ω, Ωֹ a nd Ωֹ ֹֹand neednot to be an integer.When B p (em) is no function of time , m in (27) reduces tom 0 = 2n 02 – n 0 . (28)If B p (em) in (27) is identified with a time independent B p (sd) from (11), n = n 0 = 3 (seecomment following (23)) and m 0 obtains the valuem 0 = 15. (29)If B p (em) in (27) is identified with a time dependent B p (sd) from (11), n ≠ 3 and n 0 =3 (see again comment following (23)). Substitution of B p (sd) of (11) and its time deri-vatives B ֹ p (sd) and B ֹ ֹ p (sd) calculated from (11) into (27) yields an identity and no additional expressions for m or n . When the observable quantities Ω, Ωֹ , Ωֹ ֹa nd Ωֹ ֹֹare known, n and mcan be calculated from (21a) and (27a), respectively. When the value of n is known the ratio B ֹ p (sd)/B p (sd) can be calculated from (23). Substitution of the calculated values of m , n 0 = 3 and Bֹ p (sd)/B p (sd) into (27) yields the ratio B ֹ ֹ p (sd)/B p (sd). Since the quantity Ωֹ ֹֹ has only been measured for PSR B1509–58 (see, [22]), only results for this pulsar will be discussedin section 6.4. Magnetic field from gravito-/electromagnetic originPulsars mainly consist of electrically neutral matter, probably neutrons, whereassome charged particles may also be present. In this section the charge distribution inpulsars induced by the magnetic field from gravito-magnetic origin will be investigated.Therefore, some elements of the gravito-magnetic theory will be given here. In thestationary case the magnetic induction field B can be calculated from the simplifiedgravito-magnetic equations [8, ch. 2]∇ x B = – 4π β c – 1 G 1/2 ρ v and ∇ . B = 0, (30)where ρ is the homogeneous mass density of a sphere of radius R moving with a velocityv = Ω x r (0 ≤ r ≤ R ). A dipolar magnetic field B at distance r > R can be calculated from(30)3 M . r M B = ⎯⎯ r – ⎯ (31) ⎯ , r 5 r 3where the gravito-magnetic dipole moment M = M (gm) is given by (1) (see also (2)). It isassumed that β = + 1, so that M and S posses opposite directions. To my knowledge theorigin of the basic magnetic field B of pulsars has never been explained. The predictionof a such a field B , i. e, B (gm) of (31) from (30) may be considered as a first merit of thegravito-magnetic theory. Usually, it is thought that the magnetic field of pulsars is due tocirculating charge (see, e. g., Reisenegger [24] for a recent discussion on the origin of themagnetic fields of pulsars). In that case, the (electro)magnetic field B in (31) can be writtenas B (em) and the magnetic dipole moment M as M (em). So far, no generally acceptedmodel for the calculation of the field B (em) has, however, emerged.The components of B = B (gm) of (31) in spherical coordinates (r , θ and φ) aregiven by2 M cos θ M sin θ B = ⎯⎯⎯⎯ e , B = ⎯⎯ e and B φ = 0, (r > R ) (32) r r θ⎯θ r3 r 3where e r and e θ are unit vectors. The field B in (31) has been calculated assuming r > R ,but below we also need knowledge of the field B inside the sphere. The latter field hasnot yet been calculated from (30), but for a sphere with homogeneous mass density ρ themagnetic field at the centre, B c , can be shown to beB c = 5 M /R 3. (33)The field B inside the sphere and its components, however, may approximately be calcu-lated from (31) and (32), respectively, for all values of 0 < r ≤ R , if the magnetic dipolemoment M of the sphere with homogeneous mass density ρ shrinks to an ideal magneticdipole moment located at the centre of the sphere. As can be seen from (33), the singularityt r = 0 in (31) and (32) does not occur in the full solution for B c of the equations of (30). a If v is the velocity of charge in the pulsar, the current density j is given byj = σ (E + c – 1 v x B ′), where B ′ ≡ B + B (em). (34)Here σ is the electrical conductivity of the plasma, B = B (gm) is the magnetic inductionfield from gravito-magnetic origin, B (em) the contribution from the moving charge, B ′ =B (tot) (see (5)) and E the electric field inside the sphere. Since the conductivity inside thepulsar may be extremely high, the following well-known relation from plasma physicsmay be usedE = – c – 1 v x B ′ = – c – 1 (Ω x r) x B ′. (35)If accretion is present, the assumption j = 0 underlying (35) is probably unjustified. Thelatter case will, however, not be considered in this work. When in a first approximationB (em) = 0 is chosen in (34) and (35), combination of (32) and (35) yields the followingcomponents of E in spherical coordinates– Ω M sin 2θ 2 Ω M sin θ cos θ E = ⎯⎯⎯⎯⎯ e , E = ⎯⎯⎯⎯⎯⎯⎯ e and E φ = 0. (36) r r θθ c r 2 c r 2Then, the charge density ρe inside the pulsar can be calculated from the Maxwell equation∇ . E = 4π ρe . (37)Combination of (36) and (37) yields for the quadrupolar charge distribution ρe ρ⎯⎯⎯2θ ⎯⎯⎯0 ≤ r ≤ R ) (38) ΩM Ω.B e = (3 cos – 1) = , ( 2π c r 3 2π cwhere M still equals M (gm). A number of authors presented results related to (38). Seefor a discussion, e. g., Michel and Li [25]. Note that the right hand sides of (38) do notcontain the usual minus sign. The latter sign is obtained, however, if β = + 1 in (1) is re-placed by β = – 1. The prediction of the much applied quadrupolar charge distribution ρe of(38) may be considered as a second merit of the gravito-magnetic theory.Integration of the quadrupolar charge distribution (38) yields a positive and a nega-tive charge Q + and Q – , respectively31/2 Ω M Q + = ∫ ρe dV ⎯⎯ ln 9 c= ⎯⎯ R /R 0 = – Q – , (R 0 ≤ r ≤ R ) (39)where V is the volume of the sphere, R 0 is a chosen radius much smaller than R . So, for asmall value of R 0 the net internal charge from quadrupolar origin Q satisfies Q = Q + + Q –= 0. Note that in our case the charge-separation of available plasma is caused by the (gra-vito-)magnetic field B = B (gm) of (34). Usually (see, e.g., [25]), this charge-separation iscalculated from the (electro)magnetic field B (em) of unspecified origin in (34).Analogous to the derivation given by Michel and Li [25], an electricmonopole in the pulsar should be adopted. Combination of E r from (36) and (37)followed by inte-gration over a sphere using Gauss’ law yields for the magnitudeof the charge of this monopole2 Ω M ρe dV ∇o ⎯⎯⎯(40)Q monopole = ∫ = 1/(4π) ∫ ( . E ) dV = 1/(4π) ∫ E . dS = – . 3 cContrary to Michel and Li [25], a negative charge is obtained in our calculation. A positivecharge is obtained, however, when β = + 1 in (1) is replaced by β = – 1.In addition, integration over the whole sphere (0 ≤ r ≤ R ) shows (see, e. g. [20, § 44])that a net magnetic dipole moment M (em) generated by the quadrupolar charge distribu-tion (38) arises (in this case no singularity problem occurs)2 Ω 2 R 2 M (em)/(2 c ) ∫⎯⎯ M 15 c 2 = 1 ρe r x v dV = ⎯(gm). (41)Note that the magnetic dipole moment M (em) of (41), due to charge-separation has thesame direction as the dipole moment M (gm) from gravito-magnetic origin. When β = + 1in (1) is replaced by β = – 1, the signs on the right hand sides of (38) change, but the sign in(41) remains the same. So, M (em) in (41) always enforces M (gm). The field B (em) cor-responding to M (em) of (41) is given by an expression analogous to (31).As an aside: if the magnetic moment M in (31) would exclusively be due to circula-ting charge in the pulsar, M = M (em) and M (gm) = 0. Analogously to (41), an expressionfor the additional (electro) magnetic dipole moment M ′(em) can be deduced in that case2 Ω 2 R 2 ∫ ρe r x v dV = ⎯(em), (42) M ′(em) = 1/(2 c ) ⎯⎯ M 15 c 2where again the resulting magnetic moment M ′(em) always enforces the original M (em).It seems difficult, however, to determine the contributions M ′(em) and M (em) separately.Returning to the gravito-magnetic case, by using (41), the total magnetic momentM ′ of the pulsar can be written as2 Ω 2 R 2 (gm) + M (em) = (1 + ⎯(gm). (43) M ′ = M ⎯⎯) M 15 c 2In chapter 4 of [8] a theoretical example is discussed in which circulating charge weakensa basic magnetic field from gravito-magnetic origin. In that case the resulting magneticmoments M (em) and M (gm) in (43) are anti-parallel and the total magnetic moment M ′ isweakened by M (em).Analogous to the derivation of (38), by using (43), the following higher order qua-drupolar charge distribution ρe ′ can be deducedΩ M ′ Ω . B ′ ρe ′ = (3 cos 2π c r 3 2π c ⎯⎯⎯2θ – 1) =⎯⎯⎯ . (0 ≤ r ≤ R ) (44)Moreover, analogous to the derivation of (41), the following higher order magnetic dipolemoment M ′(em) can be deduced2 2 ) = ⎯(gm). (45) Ω2R 2 Ω2R 2 M ′(em ⎯⎯ (1 + ⎯⎯⎯) M 15 c 2 15 c 2Although other mechanisms (for example: accretion, or the mechanism mentioned in thecomment following (43)) may be the main cause of the (electro)magnetic dipole momentof (10), the magnetic dipole moment M′(em) of (45) may also contribute. In any case, thenewly deduced magnetic dipole moment M′(em) of (45) may be considered as a thirdconsequence of the gravito-magnetic theory.Note that combination of (3a), (6b) and (45) yields an additional expression for β′β′ = + 2/15 (ΩR/c)2 {1 + 2/15 (ΩR/c)2}. (46) An upper limit of β′ = 2/15 (1 + 2/15) = 0.1333 x 1.1333 = + 0.1511 is obtained from(46), when the ratio ΩR/c approaches unity value for pulsars with a very short period. It isnoticed that still higher order values for β′can be calculated, but the series convergesquickly to + 0.1538. From the comment following (41) it can be concluded that the sign ofβ′ is positive for both β = + 1 and β = – 1.It may be possible, that equation (44) remains valid for values of r somewhat largerthan R. In that case the quantity R in (45) and (46) must be replaced by r. A larger r re-sults into an increasing value of β′, but the upper limit of β′ remains 0.1538 when r equalsr =c/Ω, the classical radius of the light cylinder. A charge density according to (38)outside the pulsar radius R has been considered by a number of authors. See discussionsfrom, e.g., [21, pp. 178–181], Michel and Li [25] and Michel [26] (for recent theoreticaland observational developments).5. Observational dataIn table 2 the values for the observed magnetic field B p(tot) of a number of accre-tion powered X-ray, binary pulsars are taken from Coburn et al. [27] and others [28–30], whereas the corresponding values of P and Pֹ are taken from the same references orother ones [31–36]. Values of B p(tot) are obtained from so-called cyclotron resonancespectral features (CRSFs) that are attributed to resonant scattering of photons by elec-trons. Unfortunately, no generally accepted model is available for the shape of the conti-nuum part of the X-ray spectrum, whereas such a model is necessary to isolate the line-like spectral features. Coburn et al. [27], however, applied the same continuum model(power-law with high-energy cutoff) to ten pulsars in table 2 and therefore their resultshave been used as much as possible.More explicitly, in the presence of a gravitational field the field B p(tot) has beencalculated from the fundamental energy E cycl⎛eħ⎞⎛2Gm ⎞1/2E = ⎜⎯⎯⎟⎜1 – ⎯⎯⎟B = 11.58x10–12 0.7657 B(G) keV, (47)cycl⎝mc ⎠⎝c2r⎠ewhere m e is the mass of the electron, m = 1.4 m u is the mass of the pulsar and r = 10 km.When the scattering of the photons by the electrons takes place at the polar cap of thepulsar, the field B(G) in Gaussian units may be taken equal to B p(tot). Since the fieldB p(tot) can more or less directly be obtained, it is generally denoted as the observed field.Firstly, the observed magnetic field B p(tot) and the calculated field B p(gm) from (4)for β = 1 in table 2 may be compared. Agreement between these fields is only approximate.If the field B p(tot) would be equal to B p(gm), then the quantity βeff from (6a) would yieldunity value. It is noticed that larger values for B p(gm) and βeff are obtained, if the angularvelocity in the core, Ωcore, of the pulsar is larger than the angular velocity at the surface,Ωcrust, as a result of retrograde accretion. In that case Ωcrust in (4a) must be replaced byΩcore. The observed spin-up episodes in many of the pulsars in table 2 may be an indica-tion of the validity of such an interpretation.。