SCUBA sub-millimeter observations of gamma-ray bursters III. GRB 030329 The brightest sub-m

合集下载

High power pulsed magnetron sputtering A review on scientific and engineering state of the art

High power pulsed magnetron sputtering A review on scientific and engineering state of the art

High power pulsed magnetron sputtering:A review on scienti fic and engineering state of the artK.Sarakinos a ,⁎,1,J.Alami b ,1,S.Konstantinidis c ,1a Materials Chemistry,RWTH Aachen University,52074Aachen,Germanyb Sulzer Metaplas GmbH,Am Böttcherberg 30-38,51427,Bergisch Gladbach,GermanycLaboratoire de Chimie Inorganique et Analytique,Universitéde Mons,Avenue Copernic 1,7000Mons,Belgiuma b s t r a c ta r t i c l e i n f o Article history:Received 31March 2009Accepted in revised form 6November 2009Available online 19January 2010Keywords:HPPMS HiPIMSIonized PVDHigh power pulsed magnetron sputtering (HPPMS)is an emerging technology that has gained substantial interest among academics and industrials alike.HPPMS,also known as HIPIMS (high power impulse magnetron sputtering),is a physical vapor deposition technique in which the power is applied to the target in pulses of low duty cycle (b 10%)and frequency (b 10kHz)leading to pulse target power densities of several kW cm −2.This mode of operation results in generation of ultra-dense plasmas with unique properties,such as a high degree of ionization of the sputtered atoms and an off-normal transport of ionized species,with respect to the target.These features make possible the deposition of dense and smooth coatings on complex-shaped substrates,and provide new and added parameters to control the deposition process,tailor the properties and optimize the performance of elemental and compound films.©2009Elsevier B.V.All rights reserved.Contents 1.Introduction ..............................................................16622.HPPMS:power supplies and processes .................................................16623.Plasma dynamics in HPPMS ......................................................16643.1.On the HPPMS plasma .....................................................16643.2.Determination of plasma properties:a theoretical overview ....................................16653.3.The HPPMS plasma properties ..................................................16663.3.1.The effect of the pulse on/off time con figuration on the plasma properties .........................16683.3.2.Rarefaction in the HPPMS plasma ............................................16693.4.Instabilities in the HPPMS plasma ................................................16713.4.1.Plasma instabilities:an overview ............................................16713.4.2.Plasma initiation and the relationship between plasma instabilities and the chargetransport in the HPPMS plasma .....16714.Plasma-surface interactions and deposition rate .............................................16734.1.Non-reactive HPPMS ......................................................16734.1.1.Target-plasma interactions and target erosion rate ....................................16734.1.2.Transport of ionized sputtered species ..........................................16754.2.Reactive HPPMS ........................................................16755.Growth of thin films by HPPMS ....................................................16765.1.HPPMS use for deposition on complex-shaped substrates .....................................16765.2.Interface engineering by HPPMS .................................................16775.3.The effect of HPPMS on the film microstructure .........................................16785.4.Phase composition tailoring of films deposited using HPPMS ...................................16796.Towards the industrialization of HPPMS ................................................16816.1.Examples of industrially relevant coatings by HPPMS .................................. (1681)Surface &Coatings Technology 204(2010)1661–1684⁎Corresponding author.Tel.:+492418025974,+492204299390,+3265554956.E-mail address:sarakinos@mch.rwth-aachen.de (K.Sarakinos).1All authors have contributed equally to the preparation of themanuscript.0257-8972/$–see front matter ©2009Elsevier B.V.All rights reserved.doi:10.1016/j.surfcoat.2009.11.013Contents lists available at ScienceDirectSurface &Coatings Technologyj o u r n a l h om e p a g e :w w w.e l s ev i e r.c o m /l o c a t e /s u r fc o a t6.2.Up-scaling of the HPPMS process (1682)6.3.Positioning HPPMS in the coating and components market (1682)Acknowledgements (1682)References (1682)1.IntroductionThinfilms are used in diverse technological applications,such as surface protection and decoration,data storage,and optical and microelectronic devices.The increasing demand for new functional films has been a strong incentive for research towards not only understanding the fundamentals and technical aspects of thinfilm growth,but also developing new deposition techniques which allow for a better control of the deposition process.Among the various methods employed forfilm growth,Physical Vapor Deposition(PVD)techniques, such as magnetron sputtering,are widely used[1].Magnetron sputtering is a plasma-based technique,in which inert gas(commonly Ar)atoms are ionized and accelerated as a result of the potential difference between the negatively biased target(cathode)and anode.The interactions of the ions with the target surface cause ejection(sputtering)of atoms which condensate on a substrate and form afilm[1,2].The condensation and film growth processes frequently occur far from the thermodynamic equilibrium,due to kinetic restrictions[2,3].Thus,a good control of both the thermodynamic and the kinetic conditions has implications on the growth dynamics and enables the tailoring of the structural,optical, electrical,and mechanical properties of thefilms[4].One way to control the growth dynamics is by heating the substrate during the deposition [2].The magnitude of the deposition temperature affects the energy transferred to thefilm forming species(adatoms)[2].This energy is,for instance,decisive for the activation of surface and bulk diffusion processes[5,6]and enables control over thefilm morphology[7–10]. Another source of energy are the plasma particles that impinge onto the growingfilm transferring energy and momentum to the adatoms[4,5]. The bombardingflux consists both of neutral and charged gas particles,as well as sputtered species.Numerous studies have shown that the energy, theflux,the angle of incidence,and the nature of the bombarding species are of importance for the properties of the depositedfilms[4,5,11,12].In general,the plasma particles exhibit a relatively broad energy distribu-tion with a mean value of several eV[13].The particle energy,andflux are affected by the target-to-substrate distance and the working pressure[1]. When the particles are charged the control of their energy can additionally be achieved by the use of electricfields,e.g.by applying a bias voltage to the substrate[4,14],while theirflux depends on a number of factors including the plasma density,the target power,and the magneticfield configuration at the target[4,14].It is therefore evident, that a high ion fraction in a depositingflux facilitates a more efficient and accurate control of the bombardment conditions providing,thus,added means for the tuning of thefilm properties.In magnetron sputtering processes the degree of ionization of the plasma particles is relatively low[13],resulting in a low total ionflux towards the growingfilm[13].As a consequence,in many cases bias voltage values of several tenths or even hundreds of V are required in order to increase the average energy provided to the deposited atoms and significantly affect thefilm properties[5].Moreover,the degree of ionization of sputtered species is typically less than1%[13,15,16].As a result,the majority of the charged bombarding particles is made of Ar+ions[13].This fact in combination with the relatively high bias voltages may cause subplantation of the Ar atoms in thefilm[17,18], leading to a generation of lattice defects,[18,19]high residual stresses [20–23],a deterioration of the quality of thefilm/substrate interface [21],and a poorfilm adhesion[24].Thus,the increase of the fraction of the ionized sputtered species has been an objective of many research works during the recent decades.Some of the approaches that have been demonstrated include the use of an inductively coupled plasma (ICP)superimposed on a magnetron plasma[25],the use of a hollow cathode magnetron[25],and the use of an external ion source[26]. Parallel to the above mentioned approaches Mozgrin et al.[27], Bugaev et al.[28]and Fetisov et al.[29]demonstrated in the mid90s that the operation of a conventional sputtering source in a pulsed mode,with a pulse duration ranging from1µs to1s and a frequency less than1kHz,allowed for pulsed target currents two orders of magnitude higher than the average target current in a conventional sputtering technique,such as direct current magnetron sputtering (dcMS)[27–29].These high pulse currents resulted,in turn,in ultra-dense plasmas with electron densities in the order of1018m−3[27–29],which are much higher than the values of1014–1016m−3 commonly obtained for dcMS[13,30].A few years later,Kouznetsov et al.[31]demonstrated in a publication,that received much acclaim,the high power pulsed operation of the sputtered target showing that in the case of Cu the high plasma densities obtained using this pulsed technique resulted in a total ionflux two orders of magnitude higher than that of a dcMS plasma,and a sputtered material ionization of ∼70%.The new deposition technique was called High Power Pulsed Magnetron Sputtering(HPPMS).In the years after Kouzsetov's report, HPPMS received extensive interest from researchers and led to a substantial increase of the HPPMS-related publications(Fig.1).Later, several research groups[25]adopted the alternative name High Power Impulse Magnetron Sputtering(HiPIMS)for this technique.The aim of this review paper is to describe the scientific and engineering state of the art of HPPMS.For this purpose,the principle and the basic structure of the existing HPPMS power supplies are described and the various existing approaches to produce high power pulses are demonstrated.The spatial and the temporal evolution of the HPPMS plasma and its interactions with the target and the growingfilm are discussed both from a theoretical and an experi-mental perspective.Furthermore,the effect of the HPPMS operation on the deposition rate and thefilm properties is presented.Finally,the use of HPPMS in industry and the challenges that the latter faces regarding the use of HPPMS are briefly addressed.2.HPPMS:power supplies and processesThe experimental realization of HPPMS requires power supplies different than those used in conventional magnetron sputtering processes. Those power supplies must be able to provide the target with pulses of high power density(typically in the range of a few kW cm−2),while maintaining the time-averaged target power density in values similarto Fig.1.Number of HPPMS-related publications in the period1999–2008.1662K.Sarakinos et al./Surface&Coatings Technology204(2010)1661–1684those during dcMS (i.e.a few W cm −2).The low average target power density is necessary to prevent overheating of the cathode and damage of the magnets and the target.A number of research groups and companies have developed HPPMS arrangements for both laboratorial and industrial use.These devices exhibit similarities in their basic structure which is schematically depicted in Fig.2.A dc generator is used to load the capacitor bank of a pulsing unit,which is connected to the magnetron.The charging voltage of the capacitor bank ranges typically from several hundreds of V up to several kV.The stored energy is released in pulses of de fined width and frequency using transistors with a switching capability in the µs range,located between the capacitors and the cathode.The pulse width (also referred to as pulse on-time)ranges,typically,from 5to 5000µs,while the pulse repetition frequency spans from 10Hz to 10kHz.Under these conditions,the peak target current density may reach values of up to several Acm −2,which are up to 3orders of magnitude higher than the current densities in dcMS [25].The high target current densities during the pulse on-time are accompanied by differences in the electrical characteristics of the discharge,as manifested by the current –voltage (I –V )curves of a magnetron operating in a dcMS and an HPPMS mode (Fig.3).According to Thornton [32]the target and the voltage of the magnetron are linked by the power law I ∝Vnð1ÞIn dcMS processes the exponent n in Eq.(1)ranges from 5up to 15(Fig.3).In HPPMS the exponent n changes from values similar to dcMS at low target voltages to values close to the unity when high voltages are applied (Fig.3)[33,34].Despite their common basic architecture,the available HPPMS power supplies exhibit differences primarily in the width and the shape of the pulses they can deliver,and in whether they can sustain a constant voltage during the pulse on-time.Arrangements able to provide high power pulses were already developed in the mid 90s,in order to explore the feasibility of producing high-current quasi-stationary magnetron glow discharges for high rate deposition [27–29].However,it was in 1999when Kouznetsov et al.[31]emphasizedthat the high power pulsing could allow for a high ionization degree of copper vapor (70%)as well as a homogeneous filling of 1:2aspect ratio trenches (see Fig.4).The time-dependent target current and target voltage curves of the HPPMS power supply developed by Kouznetsov et al.[31]are plotted in Fig.5where it is shown that the ignition of the plasma is accompanied by a drop of the voltage from ∼1000V to ∼800V,as the current increases to its peak value ranging from 200to 230A at the highest applied working pressure [35,36].A characteristic of this power supply is that the loading voltage is not maintained constant during the pulse on-time,and its temporal evolution is rather determined by the size of the capacitor bank and the time-dependent plasma impedance.Other power generators were designed to deliver voltage pulses with a rectangular shape (see Fig.6),i.e.a constant voltage value during the pulse on-time [34,37–39].It is to be noted that pulses with widths larger than about 50µs allow for a saturation of the discharge current and establishment of a steady-state plasma [40]provided that enough energy is stored in the capacitor bank.A common problem encountered in sputtering processes is the occurrence of arcs.This phenomenon is particularly pronounced during the deposition from conducting targets covered by insulating layers and can be detrimental for the quality of the films,due to the ejection of µm sized droplets from the target [41,42].The high peak target currents during the HPPMS operation may enhance the frequency of the arc events [43].Therefore,sophisticated electronics can be used in conjunction with the HPPMS power supplies to limit their effect if formed [43].An alternative way to alleviate this problem is to operate HPPMS using short pulses with widths ranging from 5to 20µs [44].A waveform during this mode of operation is presented in Fig.7.Within this short period of time,the glow discharge remains in a transient regime [40](as manifested by the triangular form of the target current)so that an eventual glow-to-arc transition is prevented.These short pulses have been,for instance,shown to allow for a stable and an arc-free reactive deposition of metal oxide films [45,46].However,in order to eliminate the time lag for plasma ignition which would allow for the production of high-current pulses within this short period of time,a pre-ionization stage has to be implemented [44,47].The pre-ionization of the discharge ensures that a low density plasma is maintained between the pulses as the charge carriers are already present before the voltage is turned on.The plasma pre-ionization step may be achieved by super-imposing a secondary plasma (such as a dc,a microwave,or an inductively coupled plasma)on the HPPMS plasma or by running the discharge at relatively high frequencies [44,47–49].The latter pre-ionization method is extensively implemented in mid-frequency pulsed plasmas [50–53].In contrast to the short pulses described above,long pulses with a width of several thousands of µs can also be used [54](Fig.8).In this case the pulse is composed of two (or more)stages.TheFig.2.Basic architecture of an HPPMS power supply.The dc generator charges the capacitor bank of a pulsing unit.The energy stored in the capacitors is dissipated into the plasma in pulses of well-de fined width and frequency using ultra fastswitches.Fig.3.Current –voltage curves of a magnetron during operation in dcMS and HPPMS mode.The change of the slope from 7to 1at a voltage value of 650V indicates loss of the electron con finement (data taken from [33]).Fig. 4.Cross-section SEM image of two via holes with an aspect ratio 1:2homogeneously filled by Cu using HPPMS (reprinted from [31]after permission,1999©Elsevier).1663K.Sarakinos et al./Surface &Coatings Technology 204(2010)1661–1684first step of the discharge is made of a weakly ionized regime,while the second part of the pulse brings the glow discharge to the high ionization state.The weakly ionized regime allows for the formation of a stable discharge prior to entering the highly ionized regime which helps to suppress the arc formation during the high power mode of operation [54].3.Plasma dynamics in HPPMS 3.1.On the HPPMS plasmaPlasmas are partially ionized gases characterized by,for example,the charge (electron)density or the electron energy distribution function [55].A classi fication of different plasmas could be based on the charge density as shown in Fig.9,where the magnetron sputtering plasmas are shown to exhibit densities ranging from ∼1014up to 1020m −3.It is known that high density plasmas such as arc plasmas exhibit high ionization fractions [25,31].This is not the case for conventional magnetron plasmas where the electron density does not exceed the value of 1016m −3,penning ionization is the main ionization mechanism [30]and therefore the ionization of the sputtered species is at the best of few percent [13].In order to achieve a high ion fraction in a discharge it is necessary to promote the electron impact ionization [55]which is best realized by increasing the electron density and temperature.Different approaches have been made to achieve this,for example by using an electron cyclotron resonance (ECR)plasma as a plasma electron source [56],or by employing an inductively coupled plasma (ICP)[57,58].A higher frequency of electron impact ionization collisions can alsobeFig.5.Target current and voltage waveforms produced by the power supply developed by Kouznetsov et al.at working pressures of (a)0.06Pa,(b)0.26Pa,(c)1.33Pa and (d)2.66Pa.The charging voltage of 1000V drops down to 800V and the plasma is ignited.Peak target currents of up to 230A are achieved.The time lag for the ignition of the plasma decreases when the pressure is increased (reprinted from [35]after permission,2002©Elsevier).Fig.6.Rectangular shape voltage waveforms generated by the HPPMS power supply developed by MELEC GmbH.The data were recorded from an Ar –Al HPPMS discharge operating at pulse on/off time con figuration 25/475µs (K.Sarakinos,unpublisheddata).Fig.7.Temporal evolution of target voltage and current during operation using a 10µs long pulse.The target current is interrupted before the plasma enters into the steady-state regime.The data were recorded from an Ar –Ti HPPMS discharge (S.Konstantinidis unpublished data).1664K.Sarakinos et al./Surface &Coatings Technology 204(2010)1661–1684achieved by increasing the power to the sputter source and therewith the electron density.This approach is,however,limited since a power excess can result in an overheating of the target or an exceeding of the Curie temperature of the magnetron's magnets.Thus,the work of Kouznetsov et al.[31]opened a new era of magnetron sputtering deposition and plasma studies,since it allowed for the generation of highly ionized plasmas using a conventional magnetron source.The early analyses of the HPPMS plasma showed that its temporal characteristics are complex and unique [35,59,60].Understanding the discharge evolution and the mechanisms that take place in the plasma is therefore of utmost importance as this sheds light on the growth conditions during the thin film deposition.To achieve this,modeling and a number of analytical studies have been carried out,using Langmuir probes as well as optical emission,mass,and absorption spectroscopy techniques.These studies are reviewed in the following chapters and the basic properties of the HPPMS plasma are thus outlined.3.2.Determination of plasma properties:a theoretical overview Plasmas are commonly modeled as fluids,and each species is described with its fluid equation [61].The fluid approximation is suf ficiently accurate to describe the majority of observed plasma phenomena.When this is not viable,the accumulated behavior of anensemble of charged particles is described using a statistical approach,in the form of the velocity distribution function which is given as f (r,u,t )d 3r d 3u ,i.e.the number of particles inside a six-dimensional phase space volume d 3r d 3u at (r ,u )and time t [55].Here,r is the position vector and u is the speed vector.Knowing the distribution function,fundamental features of a large ensemble can be quanti fied,i.e.by integrating over the velocity,the average (macroscopic)quantities associated with the plasma are de fined as:n ðr ;t Þ=∫f ðr ;u ;t Þd u plasmanumber ðion ;electron ;orneutral Þdensityð2Þu Àðr ;t Þ=1n∫u f ðr ;u ;t Þd u plasmabulkfluid velocity ð3ÞεÀðr ;t Þ=m 2n∫ðu À−u Þ2f ðr ;u ;t Þd u plasmameankineticenergy ð4Þwhere m is the particle (ion,electron,or neutral)mass.This procedure of integrating over velocity space is referred to as taking velocity moments,with each one yielding a physically signi ficant quantity.Knowledge of the electron energy (velocity)distribution function (EEDF)permits us to determine important parameters.The EEDF can be measured using a Langmuir probe,as was demonstrated by Druyvesteyn [55]:g e ðV Þ=2m e A pr 2eV m12d 2I edV ð5Þwhere A pr is the probe area,I e is the electron current,m is the electronmass,and V =V pl −V b is the potential difference between the plasma potential and the probe potential.With the change of variables ε=eV ,the electron density n e is determined as:n e =∫∞0g e ðεÞd εð6ÞPlasma parameters,such as the electron density,are easily obtained if the mathematical description of the EEDF is Maxwellian.The Maxwellian distribution is characterized by its single temperature and is obtained when plasma electrons undergo enough collisions and are in equilibrium with other plasma components and with the electron ensemble itself.At a low pressure,the EEDF is generally non-Maxwellian,and the electron temperature is thought of as aneffectiveFig.8.Target voltage and current waveforms during the long-pulse operation.In the first 500µs a weakly ionized plasma is generated which is followed by a highly ionized steady-state plasma (reprinted from [54]after permission,Society of Vacuum Coaters ©2007).Fig.9.Plasma density is used in order to classify plasmas.The HPPMS plasma spans over a large density range depending on the pulse on/off time con figuration used for a process and presents therefore a tool for better choosing the plasma dynamics needed for the deposition process.The abbreviation ECR stands for Electron Cyclotron Resonance (taken from [36]).1665K.Sarakinos et al./Surface &Coatings Technology 204(2010)1661–1684electron temperature,T eff∼23〈ε〉,representing the mean electron energy determined from the EEDF according to[55]:〈ε〉=1n e ∫∞εg eðεÞdεð7Þ3.3.The HPPMS plasma propertiesOne of the earliest works related to the study of the properties of the HPPMS plasma was carried out by Gudmundsson et al.[35,59]. Using a Langmuir probe[62]in an Ar–Ta plasma,they measured the temporal behavior of the EEDF by recording the time-dependent probe current curves.It was shown that the EEDF evolved from a Druyvesteyn-like distribution with a broad energy distribution(high average energy)during the pulse to a double Maxwellian distribution (i.e.a distribution that is composed of two-distinct energy-distribu-tions)toward the end of the pulse,andfinally a Maxwellian-like distribution hundreds ofµs after the pulse had been switched off.This is shown in Fig.10,where the temporal evolution of the EEDF is also seen to be affected by the working pressure[35].Similar results were obtained by Pajdarova et al.[63],although the double Maxwellian distribution was observed from the beginning of the pulse,which could be explained by the nature of the target material(Cu in this case)and the HPPMS parameters used in the experiments.Further-more,it was found that the effective electron temperature peaked a fewµs after the pulse start,and at the same time,independent of measurement position[35].This is indicative of the existence of high energy electrons accelerated by the target voltage,and thus escaping the confinement region in the vicinity of the cathode[35].InaFig.10.The EEDF(electron energy distribution function)for an HPPMS plasma with a100µs on-time and a50Hz frequency at three chamber pressures(a)0.25,(b)1.33,and(c) 2.66Pa.The distribution function changes during the pulse from a Druvesteyn distribution to a double-Maxwellian distribution(reprinted from[35]after permission,Elsevier©2002).1666K.Sarakinos et al./Surface&Coatings Technology204(2010)1661–1684different study by Seo et al.[64],Langmuir probes were utilized to determine the properties of a HPPMS plasma at a higher pulsing frequency of 10kHz and a duty cycle of 10%.By using such a high frequency,the electron density reached a value of 1.5×1017m −3,which could be seen as a relatively low density HPPMS plasma.The temporal evolution of the energy distribution in this plasma showed ef ficient electron heating,which took place within the first few µs from the pulse start.This work revealed that the plasma dynamics presented in the studies by Gudmundsson et al.[35,59]are to a large extent general features of the low duty cycle pulsed plasmas regardless of the frequency.The early studies on the HPPMS plasma characteristics also showed the potentially high ionization fraction of the sputtered material that could be achieved by using this technique [16,25,33,60,65,66].Spectroscopic analyses,such as optical emission spectroscopy [67–69],mass spectroscopy [69],and absorption spectroscopy [68,69]were,therefore,employed in order to better quantify and understand the average and time-dependent evolution of the ion population in the HPPMS plasma.The optical emission studies [33,65,70–72],were performed both for the HPPMS and the dcMS discharges and showed the much higher metal ion emission intensity in HPPMS (Fig.11).However,in order to better quantify these observations,absorption spectroscopy [70,71,73,74]and mass spectroscopy [75,76]studies were performed,con firming the much higher ionization of the sputtered material and of the sputtering gas achieved in HPPMS.One well studied material in this regard is Ti.Mass spectroscopy measurements of the peak ion compositions in an Ar –Ti discharge for the pulse on-time of 100µs and a frequency of 50Hz [75]showed that the ionized part of the plasma contained 50%of Ti 1+,24%of Ti 2+,23%of Ar 1+,and 3%of Ar 2+ions.When longer pulse on-times of several hundreds of µs were used even Ti 3+and Ti 4+could also be detected [77].It is deduced that for the formation of multiply charged ions certain conditions should be ful filled [77];(i)the discharge should contain a high enough number of metal ions,(ii)the pulse on-time should be long enough (and the capacitance of the used power supply should be large enough to store enough energy for the whole pulse),(iii)the self-sputtering yield (i.e.the sputtering yield of the target material from ionized target species)should be low,and (iv)the ionization energy of each ionization step should not be too high.These two examples show clearly that the HPPMS plasma properties depend not only on the plasma density and energy distribution but also,to a large degree,on the pulse on/off time con figuration.Since the introduction of the HPPMS technique,most of the discharge studies have been performed for pulses with pulses on-time longer than 50µs.Exceptions to this can be found in the works by Konstantinidis et al.[70,73],Ganciu et al.[44],and Vasina et al.[48],where shorter pulses ranging between 5and 30µs have been used,an example of which is shown in ing time-resolved optical emission spectroscopy in a Ti –Ar discharge,an increase in the line intensity ratio of the ion-to-neutral Ti and Ar species with increasing the pulse duration was observed [70].This means that the ionization degree increases with increasing pulse length,indicating the necessity for the sputtered atoms to reside a suf ficient time in the target vicinity for the vapor to become ionized.The increase of the AremissionFig.11.Optical emission spectroscopy from a typical HPPMS plasma.The emission from the Ti and Ar species increases greatly in HPPMS indicating the high ionization of the metal and gas species (data taken from [72]).Fig.10(continued ).1667K.Sarakinos et al./Surface &Coatings Technology 204(2010)1661–1684。

一些物理里面名词的中文对照(大部分是材料领域的)

一些物理里面名词的中文对照(大部分是材料领域的)

Abelian group 阿贝尔群,又称Abel群ablation 烧蚀abnormal dispersion 反常色散Abrikosov vortex lattice 阿布里科索夫涡旋线格子Abrikocov vortex state 阿布里科索夫涡旋态absorber 吸收体absorption spectroscopy 吸收光谱abundance 丰度acceptor doping 受主掺杂acceptor impurity 受主杂质accumulation layer 累积层achromatic phase matching 消色差相位匹配achromatic wave plate 消色差波片achromatism 消色差[性]ac Josephson effect 交流约瑟夫森效应,又称交流Josephson效应acoustic compliance 声顺acoustic ohm 声欧[姆]acoustic stiffness 声劲[度]acoustic-optic tensor 声光系数张量acousto-optic effect 声光效应acousto-optic Q-switch 声光Q-开关acousto-optic signal processor 声光信号处理器acousto-optical tunable filter 声光可调滤波器actinide element 锕系元素activated tunneling 激活隧穿active device 有源器件active region 激活区addressing electrode 寻址电极adiabatic theorem,绝热定理adiabatic transformation 绝热变换adiabatic transport,绝热输运adiabaton 浸渐子,绝热子advection 平流aerodynamic sound 空气动力声aersol 气溶胶affinity potential 亲和势aggregate 聚集体aggregation 聚集Aharonov-Bohm (AB) effect AB效应,又称Aharonov-Bohm (AB) 效应Aharonov-Bohm (AB) flux AB磁通,又称Aharonov-Bohm (AB)磁通allowed state 容许态alpha decay ( -decay) 衰变alpha particle ( -particle) 粒子Altshular-Aronov-Spivak (AAS) effect AAS效应,又称Altshular-Aronov-Spivak效应amplification without inversion 无反转放大amplitude limiting 限幅amplitude transformer 变幅杆Andreev reflection 安德列也夫反射,又称Andreev反射Andreev mirror 安德列也夫镜[子],又称Andreev镜[子] Andreev scattering 安德列也夫散射,又称Andreev散射angular resolved photoemission spectroscopy 角分辨光电子谱[学] anisotropic confinement 各向异性限域anisotropic scatterer, 各向异性散射体anisotropy energy 各向异性能anomalon 反常子anomalous power laws 反常幂[次]率anomalous proximity effect,反常临近效应anomaly 反常antidot 同quantum antidot 反量子点antidodal point 腹点antigravity 反引力antihyperon 反超子anti-localization, 反局域化antimeson 反介子anti-exclusive principle 反不相容原理antiferromagnetic interaction 反铁磁相互作用antiferromagnetic semiconductor 反铁磁半导体anti-Stokes scattering 反斯托克斯散射anti-time ordered function, 反时序函数anyon 任意子aphelion 近日点, 远核点areal density 面密度armchair nanotube 扶手椅型纳米管arrayed waveguide gratings 阵列波导光栅artificial atom,人[工]构[造]原子artificial barrier 人工势垒artificial elment 人造元素atom laser 原子凝射器atom optics 原子光学atom trapping 原子陷俘,原子捕获atom waveguide 原子波导atomic clock原子钟atomic diffraction 原子衍射atomic fountain 原子喷泉atomic form factor 原子形状因子atomic time 原子时attenuation 衰减attosecond X-ray pulse 阿秒X射线脉冲Auger process 俄歇过程,又称Auger过程avalanche counter 雪崩计数器avalanche effect 雪崩效应avalanche photodiodes,apd 雪崩光电二极管azimuth 方位角back-action evasion 非干扰[测量]background radiation 本底辐射,背景辐射background temperature 本底温度, 背景温度balanced homodyne detection平衡零拍探测ballistic aggregate 弹道聚集体ballistic aggregation 弹道聚集ballistic electron injection 弹道电子注入ballistic transport弹道输运ballistics 弹道学band bending 带弯曲band index 带指标band of rotation-vibration 振转[谱]带band offset 带阶band repulsion 带排斥band theory 能带论bar 巴(压强单位),杆Barkhausen noise 巴克豪森噪声,又称Barkhausen噪声barn 巴恩(截面单位,10-24厘米2)barrier 势垒barrier curvature 势垒曲率barrier height 势垒高度barrier state 势垒态barrier tunneling 势垒隧穿base-centered orthorhombic lattice 底心正交格[子] base line 基线base material 基质base metal 碱金属basis vector 基矢beam 束,梁beam dump 束流捕集器beam focusing 束流聚焦behaviour 行为,性能Bell inequality贝尔不等式,又称Bell不等式bend resistance,弯曲电阻bent crystal 弯晶Berry phase 贝里相位,又称Berry相位βdecay β衰变βradioactivity β放射性βray β射线βspectum β谱βstability line β稳定线bevatron 吉伏质子加速器(高能质子同步稳相加速器)bicritical point 双临界点bicrystal junction 双晶结big bang model 大爆炸模型binary diffractive optical element 二元衍射光学元件bioastrophysics 天体生物物理学biochip 生物芯片bipolar junction transistor 双极[结]晶体管bit rate 比特率blackness 黑度blaze line 闪耀角bleaching effect 漂白效应blob 团迹,链滴Bloch electron 布洛赫电子,又称Bloch电子Bloch frequency,布洛赫频率,又称Bloch频率Bloch oscillation,布洛赫振荡,又称Bloch振荡Bloch theorem 布洛赫定理,又称Bloch定理blockade 阻塞Blonder-Tinkham-Klapwijk [BTK] model BTK模型body-centered cubic lattice 体心立方格[子]body-centered orthorhombic lattice 体心正交格[子]Bogoliubov [-de Gennes] equations 博戈留波夫[-得简斯]方程,又称Bogoliubov [-de Gennes]方程Boltzmann distribution 玻尔兹曼分布Boltzmann transport equation,玻尔兹曼输运方程bond-angle order 键角有序bond-orientational order 键取向有序bond polarizability 键极化性bond valence 键价boojum 布经(超流氦3中的取向织构)bosonization of field operators 场算符的波色化Bragg peak 布拉格峰,又称Bragg峰Bragg plane 布拉格平面,又称Bragg平面Bragg reflection 布拉格反射,又称Bragg 反射Bragg reflectors 布拉格反射器,又称Bragg 反射器Bragg waveguide 布拉格波导,又称Bragg 波导break junction 断裂结breathing mode呼吸模breeder 增殖反应堆breakup reaction 崩裂反应bright state 亮态brittleness 脆性buffer amplifier 缓冲放大器buffer gas 缓冲气体buffer layer, 缓冲层burn-up 燃耗Büttiker formula, 比特克公式,又称Büttiker公式buzzer 蜂鸣器C-15 structure C-15结构C[a]esium clock 铯钟calorie 卡【洛里】candle 烛光candescence 白热,又称白炽canonical commutation relation 正则对易关系canonical variable 正则变量cantact angle 接触角canted spin order倾斜自旋有序cantilever 悬臂(原子力显微镜中的)canthotaxie眼角[式]排列(另文说明)carbon cycle 碳循环(恒星内部的)carbon nanotube 碳纳米管carrier 载流子carrier concentration 载流子浓度carrier diffuse 载流子扩散carrier reservoir 载流子库Cartesian coordinates 笛卡儿坐标Cauchu-Schwarz inequality Cauchu-Schwarz不等式cavity dark state 腔暗态cavity dumping 腔倒空cavity quantum electrodynamics 腔量子电动力学cavity resonator [谐振]腔共振器14C dating 碳14测年celestial X-ray source 宇宙X 射线源center of inversion 反演中心center of moment 矩心central collision中心碰撞center-of-mass energy 质心系能量centrifuge 离心机centrifugal separation 离心分离ceramic 陶瓷chain folding 链折叠chain statistics 链统计学chalcogenide 硫属化物channel waveguide 沟道波导chaos synchronization 混沌同步chaotic communication 混沌通讯chaotic noise 混沌噪声characteristic impedance 特性阻抗characteristic curve 特征曲线charge-separated plasma 电荷分离等离子体(正负电荷在空间不同区域的等离子体) charge imbalance 电荷不平衡charge ordering 电荷有序charge parity effect,电荷宇称效应charge qubit 电荷量子比特(超导量子比特的一种)charge-phase qubit 电荷-相位量子比特(超导量子比特的一种)charge reservoir 电荷库charge stiffness 电荷劲度(衡量外场作用下电荷被自由加速的难易程度)charge-spin coupling电荷自旋耦合(用于自旋电子学)charge stripe phase 电荷条纹相charge-to-mass ratio 荷质比charge transfer insulator 电荷转移绝缘体charge transfer salt 电荷转移盐charge velocity 电荷速度(见于电荷-自旋分离现象)charging energy,充电能chemical shift 化学位移chiral liquid crystal 手征液晶chiral molecule手征分子,又称手性分子chiral symmetry broken 手征对称[性]破缺chirp啁啾chirped Gaussian pulse 啁啾高斯脉冲chirp filter 啁啾滤波器,又称线性调频滤波器,或色散延迟线chopper 斩波器circumlunar orbit 环月轨道circumsolar orbit 环日轨道circumterrestrial orbit 环地轨道cis-lunar space 月地空间clad 覆盖clamping 箝位classical fluid 经典液体clean limit [干]净极限cleaved coupling cavity 解理耦合腔cloning fidelity克隆保真度closed shell 满壳层,又称闭壳层,英文又称closure shellcluster state簇态CNO cycle 碳氮氧循环coalescence 聚合, 并合code 1,[代]码;2,密码;3,符号coding 编码codirectional coupling 同向耦合coefficient of correlation 关联系数coefficient of elasticity 弹性系数coexistence line 共存线(相图中的)coexisting phase 共存相coherence factor 相干因子coherence length,相干长度coherent atomic recoil 相干原子反冲coherent electron tunneling 相干电子隧道coherent peak 相干峰coherent photoassociation 相干光缔合coherent population oscillation相干布居振荡coherent population trapping相干布居囚禁coherent population transfer相干布居迁移coherent structure 拟序结构coherent terahertz waves相干太赫波coherent transient effects 相干暂态效应coherent trap 相干捕获cold finger 冷头cold fusion 冷聚变collective coordinate 集体坐标collective mode 集体模collective motion 集体运动collective pinning model 集体钉扎模型collinear phase matching 共线相位匹配colloid 胶体,胶质colloidal metal 胶体金属colored noise 色噪声colossal magnetoresistance [CMR] 庞磁电阻commensurate lattice 公度格子compact star 致密星compensated impurity 补偿杂质complementary metal oxide semiconductor [CMOS] 互补金属氧化物半导体complex 1,复合体;2,络合物complex analytical signal theory 复解析信号理论complex-conjugate pulses 复共轭脉冲compliance 1,柔度;2,顺度composite Fermion 复合费米子compression of ultrashort pulses 超短脉冲压缩compressor 压缩器,压机concurrence并发纠缠,又称量子并发condensate 凝聚体condensation energy 凝聚能condenser 冷凝器conductance fluctuation, 电导涨落conductance quantization 电导量子化conduction electron 传导电子confinement 1,约束(等离子);2,限域(凝聚态);3,禁闭(高能)congregating effect 聚集效应conjugate variable 共轭变量conservation of angular momentum 角动量守恒conservation of crystal momentum 晶体动量守恒conservative dislocation motion 保守位错运动(位错沿滑移面平行于Burgers矢量运动无净质量流)conservation of energy 能量守恒conservation law of flux 磁通守恒律conservation of momentum 动量守恒conservation of particle number粒子数守恒contact angle 接触角contact potential 接触势contact resistance 接触电阻continuation 延拓continuous group 连续群contour line 等值线contour map 等值线图contradirectional coupling反向耦合conventional unit cell 惯用单胞,简称单胞convergence factor 收敛因子conversion electron 内转换电子coolant moderator 载热减速剂cooperative diffusion 合作扩散Cooperon, 库珀子Cooper pair box 库珀对盒子coplanar waveguide 共面波导copolymer 共聚物core energy 芯能core nucleus 核芯[核]correlated spontaneous emission 关联自发发射correlation exponent 关联指数cosmic aerodynamics 宇宙气体动力学cosmic age 宇宙年龄cosmic constant 宇宙常量cosmic [microwave] background radiation [CMBR] 宇宙[微波]背景辐射cosmic microwave background 宇宙微波背景cosmic string 宇宙弦cosmochemistry 宇宙化学,天体化学cosmological nucleosynthesis 宇宙核合成cosmos 宇宙co-tunneling 共隧穿Couette flow 库埃特流Coulomb blockade 库仑阻塞Coulomb gap 库仑隙Coulomb interaction 库仑[相互]作用Coulomb island 库仑岛,又称单电子岛(single electron island)Coulomb potential 库仑势Coulomb repulsion 库仑斥力Coulomb staircase 库仑台阶counter telescope 计数器望远镜coupled-channels model 耦合道模型coupled mode theory 耦合模理论coupled waveguides,耦合波导coupled wells耦合阱coupling energy 耦合能coupling strength 耦合强度covalent bond 共价键creep wave蠕波,又称爬波critical assembly [核反应堆]临界装置critical density 临界密度critical dimension 临界维度cross-phase-modulation 交叉相位调制cross field 交叉场cross junction, 十字结crosstalk attenuation 串扰衰减crystal-field splitting 晶[体]场劈裂crystalline anisotropy晶态各向异性crystal symmetry class 晶体对称类cubic lattice 立方格子cuprate 铜氧化物curie 居里(非国际制放射性活度单位)current bias 电流偏置current operator 电流算符cutoff energy,截止能量cyclone 气旋cyclotron effective mass 回旋有效质量D/A converter 等于digital to analog converter 数模转换器damping radiation 阻尼辐射dark current 暗电流dark energy 暗能量dark state 暗态dark-state polariton 暗态光极化子date line 日界线dc Josephson effect 直流约瑟夫森效应,直流Josephson效应dc SQUID (superconducting quantum interference device) 直流超导量子干涉器Debye wave vector 德拜波矢decay heat 衰变热decay time,衰减时间deceleration 减速度decibel 分贝decoherence 退相干,又称消相干decoherence-free 无退相干,又称无消相干decontamination factor 去污因子decoupling epoch 退耦期decoy state 诱骗态deformation potential,形变势degeneracy collapse 简并塌缩degenerate pressure 简并压degenerate star 简并星de Gennes-Taupin length de Gennes-Taupin长度degree of order 有序度de Haas-Shubnikov effect de Haas-Shubnikov效应delay time,延迟时间demultiplexer 解复用器dendrite 1,枝晶;2,枝蔓;3,枝蔓体dense coding 密集编码dense wavelength division multiplexing 密集波分复用density correlation function,密度关联函数density distribution 密度分布density wave 密度波depairing 拆对dephasing length,退相位长度depinning 脱钉[扎]depleted Uranium 贫化铀deplation force 排空力(胶体物理用语)depletion layer 耗尽层descreening 去屏蔽deterministic equation 确定(论)的方程deuterium 氘, 即重氢deuterium oxide 重水dextrorotation 右旋diabatic approach 非绝热近似diagnostics 诊断学diagonal element 对角元diagonal matrix 对角矩阵diagonalization 对角化diamond structure 金刚石结构diblock copolymer 双嵌段共聚物dielectric response function 介电响应函数dielectric function,介电函数dielectric microcavity 介电[质]微腔dielectric reflector 介[电]质反射器differential conductance 微分电导differential input 差分输入differential rotation 较差自转(天文学用语)differential scanning calorimetry 差分扫描量热术diffraction-free beam 消衍射光束diffractive binary optics 衍射二元光学diffuseness [parameter] 弥散参数diffusion constant,扩散常数diffusion current 扩散电流diffusion region 扩散区diffusive transport,扩散输运digit 数字digital circuit 数字电路digital cross connect 数字交叉连接digit[al] to analog converter (DAC) 数模转换器digital micromirror device 数字微镜器件dilation 膨胀dilute phase 稀相dilation symmetry 伸缩对称dimensionless conductance 无量纲电导dimer 二聚体dimerization 二聚化dipole interaction 偶极相互作用dipole giant resonance 偶极巨共振Dirac braket 狄拉克括号Dirac picture 狄拉克绘景, 即相互作用绘景directed diffusion 定向扩散directional bond 定向键directional coupler 定向耦合器directional ordering 取向有序directional quantization 方向量子化direction of magnetization 磁化方向direct lattice 正格子,又称正点阵direct transition 直接跃迁dirty limit 脏极限dirty-metal regime,脏金属区discontinuity 1,不连续[性];2,突变[性] dislocation network 位错网络disordered alloy 无序合金disordered system 无序系统dispersion compensation 色散补偿dispersion-managed solitons 调控色散孤子dissipationless flow 无耗散流dissociation energy 离解能distillable entanglement 可萃取纠缠distinguishable states可区分态distributed Bragg reflector 分布布拉格反射器domain 1,畴;2,[定义]域;3,区域donor level 施主能级dopant 掺杂物doping 掺杂dosimetry 剂量学double-barrier tunneling,双势垒隧穿double exchange interaction 双交换相互作用double heterostructure DH 双异质结doublet state 双重态dressed atom 着衣原子,又称缀饰原子droplet model 小液滴模型Drude model,德鲁德模型duty ratio 占空比d-wave pairing d波配对dyad 并矢dynamical mass 动力学质量(08.02dynamic random access memory [DRAM] 动态随机存储器dynamic screening,动态屏蔽dynamically induced coherence 动态诱导相干dynamo theory 发动机理论dyne 达因early universe 早期宇宙eccentricity 偏心率eclipse 1,食;2,交食edge channel,边缘通道edge dislocation 刃[型]位错edge state,边缘态effective field theory 有效场理论effective Hamiltonian 有效哈密顿量effective mass approximation,有效质量近似Einstein-Podolsky-Rosen thought experiment EPR思想实验Einstein-Podolsky-Rosen effect EPR效应Einstein-Podolsky-Rosen pair EPR对Einstein-Podolsky-Rosen paradox EPR佯谬elastic compliance 弹性顺度elastic deformation 弹性形变electrical isolation 电绝缘electric breakdown 电击穿electric capacity 电容electric resistance 电阻electrical quadrupole moment 电四极矩electrochemical potential 电化学势electromagnetic absorption 电磁吸收electromagnetically induced absorption 电磁感生吸收electromagnetically induced transparency 电磁感生透明electromagnetic-environment effect,电磁环境效应electron backscattering pattern 电子背散射图样electron-beam lithography 电子束刻蚀electron configuration 电子组态electron density 电子密度electron-doped high temperature superconductor 电子掺杂的高温超导体electronegativity 电负性electron-electron interaction,电子-电子相互作用electron-hole pair 电子空穴对electron-hole recombination 电子-空穴复合electron hologram 电子全息术electron transition 电子跃迁electron pair 电子对electron pair tunneling 电子对隧穿electron-phonon coupling 电子声子耦合electron temperature,电子温度electron tunneling 电子隧穿electron waveguide,电子波导electron volt (eV) 电子伏electrorheological effect 电流变效应electrorheological fluid 电流变液Eliashberg equations Eliashberg方程Eliashberg theory of strong coupling Eliashberg强耦合理论elliptical orbit 椭圆轨道elliptic flow 椭圆流emittance 发射度empirical pseudopotential method 经验赝势方法empty lattice approximation 空晶格近似endohedral fullerene 内嵌原子富勒烯end-butt coupling 端面对接耦合energy relaxation length,能量弛豫长度energy transport velocity 能量传输速度ensemble average,系综平均entangled state 纠缠态entanglement 1,纠缠;2,纠缠度entanglement concentration 纠缠浓缩entanglement measure 纠缠度量entanglement monotone 单调纠缠量entanglement of formation 生成纠缠entanglement purification 纠缠纯化entanglement witness 纠缠见证entropy force 熵力envelope function,包络函数epithermal neutron 超热中子epoxy 环氧树脂erbium-doped fiber amplifier 掺饵光纤放大器error correction 纠错Esaki diode 江崎二极管evanescent state,衰逝态even-odd nucleus 偶奇核even parity 偶宇称evolution of inflation 暴涨演化Ewald construction Ewald作图法Ewald sphere Ewald球excess current 过剩电流excess neutron 过剩中子exchange-correlation hole 交换关联空穴exchange-correlation functional 交换关联泛函exchange hole 交换空穴exchange integral 交换积分excitation spectrum 激发谱excluded volume 排除体积exclusion of flux 磁通排斥exclusion principle 不相容原理exotic nucleus 奇特核expanding universe 膨胀宇宙extended [Brillouin] zone scheme 扩展[布里渊]区图式extraterrestrial life 地外生命extravehicular activity(EV A) [太空]舱外活动f-sum rule f求和规则face-centered orthorhombic lattice 面心正交格[子] face-on 正向facsimile 传真,英文简写为faxfacula 光斑Fahrenheit thermometer 华氏温度计faint object 暗天体fan diagram 扇形图F-center F中心Feno lineshape Feno线型Feno resonance Feno共振fan spin order 扇状自旋有序farad (F) 法拉(电容单位)Faraday depolarization 法拉第退偏振Faraday law of electrolysis 法拉第电解定律far-from-equilibrium system,远离平衡态系统far-side 背面(far-side of the moon, 月球背面)far-ultraviolet (FUV) 远紫外fast fission 快裂变fatigue crack 疲劳裂纹fatigue fracture 疲劳断裂fatigue strength 疲劳强度feed [source] 馈源feeder 馈线femto (f) 飞(=10-15)(01)femtosecond pulse shaping 飞秒脉冲成形Fermi age 费米[中子]年龄Fermi age-diffusion equation 费米年龄扩散方程Fermi arc 费米弧Fermi coupling constant 费米耦合常数Fermi energy 费米能量Fermi gas 费米气体Fermi golden rule 费米黄金定则Fermi liquid 费米液体Fermi liquid parameter 费米液体参数Fermi loop 费米环Fermi point 费米点Fermi transition费米跃迁Fermi vacuum 费米真空Fermi velocity 费米速度Fermi wavelength 费米波长Fermi wave vector,费米波矢Fermi’s golden rule费米黄金规则ferrielectric crystal 亚铁电晶体ferrimagnet 亚铁磁体ferroelectric 铁电体ferroelectric crystal 铁电晶体ferromagnet 铁磁体few-cycle pulse少周[期]脉冲few nucleon transfer 少[数]核子转移Feynman path,费曼路径Feynman path integral,费曼路径积分fiber cross connect 光纤交叉连接fiber grating 光纤光栅Fibonacci sequence 斐波那契序列fiducial confidence bar 置信棒fiducial point 基准点field intensity 场强field quantization 场量子化field quantum 场量子field strength 场强figure of merit,又称qualityfactor 品质因数filament 1,丝;2,丝极finite-amplitude wave 有限振幅波,又称大振幅波finite-difference method 有限差分方法finite element method 有限元法finite size effect 有限尺寸效应finite-size scaling 有限尺寸标度first approximation 一级近似first Brillouin zone 第一布里渊区first point of Aries 春分点,英文又称:vernal equinoxfirst point of Cancer 夏至点,英文又称:summer solsticefirst point of Capricornus 冬至点,英文又称:winter solsticefirst point of Libra 秋分点,英文又称:autumnal equinoxFiske steps 费斯克台阶,又称自感应台阶fissility 易裂变性fission 1,裂变;2,分裂fission isomer 裂变同质异能素fission nuclide 裂变核素fission reactor 裂变反应堆fission-spectrum neutron 裂变谱中子fission track dating 裂变径迹年代测定fitting curve 拟合曲线five-fold symmetry 5重对称fixed-range hopping 定程跳跃flash memory 闪速存储器,简称闪存flat spectrum 平谱flattening factor 扁率floating probe 浮置电极,又称浮置探针floating phase 浮置相Floquest theorem 弗洛开定理flow resistance 流阻fluctuating wall 涨落壁fluctuation 涨落(统计物理〕,又称起伏(声学〕fluence 注量fluorescence probe 荧光探剂flux,通量flux 1通量,又称流量;2,注量率;3,焊料;4 助熔剂flux bundle 磁通束flux flow amplifier (FFA) 磁通流放大器flux flow oscillator (FFO) 磁通流振荡器flux flow transistor(FFT)磁通流三极管,又称涡旋流三极管(vortex flow transistor) flux-line lattice 磁通线格子flux line 磁通线flux tube 磁流管flux quantum 磁通量子flux quantization 磁通量子化foam 泡沫focal point 焦点focal ratio 焦比focus 1,焦点;2,震源folding Brillouin zone 折叠布里渊区forbidden beta decay 禁戒b衰变forecast 预报forward bias 正向偏压four-Josephson junction logic (4JL) 四约瑟夫森结逻辑门Fourier analysis 傅里叶分析Fourier transform 傅里叶变换Fourier [transform] nuclear magnetic resonance 傅里叶[变换]核磁共振Fourier [transform] Raman spectroscopy 傅立叶[变换]拉曼谱学four probe method 四探针法four-terminal resistance,4端电阻fractional chain yield 相对链产额fractional cumulative yield 分积累产额fractional distillation 分馏fractional independent yield 分独立产额fractional statistics 分数统计法fragment 1,碎片;2,片段Franck-Condon principle弗兰克-康登原理free electron approximation 自由电子近似free electron gas 自由电子气体free energy 自由能free –free transition 自由-自由跃迁,又称自由态间跃迁freely falling body 自由落体free radical 自由基free spectral range 自由光谱范围freezing point 凝固点Frenkel exciton 弗仑克尔激子frequency conversion 频率转换Frequency division multiplexing 频分复用frequency jitter 频率抖动frequency multiplication 倍频friction 摩擦Friedel oscillation,Friedel振荡Friedel sum rule Friedel求和规则Frohlich interaction Frohlich相互作用front velocity波前速度frustrated magnet 窘组磁体fuel cell 燃料电池Fulde-Ferrell state Fulde-Ferrell态fullerene 富勒烯full moon 满月function 函数functional (1)泛函(2)功能(的)fundamental interaction 基本相互作用fundamental space-filling mode 基本空间填充模fuse (1)熔解(2)保险丝fused silica熔融石英fusion reactor 聚变[核反应]堆fuzzy information 模糊信息fussy mathematics 模糊数学gain-clamping 增益箝位gain efficiency 增益效率Galton plate 伽尔顿板-陈gamma(γ)伽马(地磁场强单位γ=nT)gamma rayγ射线gap 1,隙;2,能隙gap anisotropy 能隙各向异性gap parameter 能隙参数gaseous state 气态gate1,门;2,栅(极)gate voltage 门电压gauge symmetry 规范对称性gauss (G) 高斯(磁感应强度单位G=10-4T)Gaussian fluctuation 高斯涨落Gauss law 高斯定理Gauss surface 高斯面generalized Balmer formula 广义巴尔末公式generalized work 广义功general refractive index 广义折射率(量子信息)geomagnetic declination 地磁偏角geomagnetic inclination 地磁倾角geometrical structure factor 几何结构因子geometrization of gravitation 引力几何化German silver 德银g-factor g因子g-factor of electrons 电子的g因子g shift g移位ghost imaging 鬼成像giant magnetoresistance (GMR) 巨磁电阻Giaever tunneling 盖沃尔隧穿(单电(粒)子隧穿)Gibbs ensemble 吉布斯系综gilbert 吉尔(磁通势单位)Ginzburg-Landau coherence length 金兹堡-朗道(GL)相干长度Ginzburg-Landau equation 金兹堡-朗道(GL)方程Ginzburg-Landau-Abrikosov Go’rkov theory(GLAG)金兹堡-朗道-阿布里科索夫-高里科夫理论Glan-Thompson prism 格兰-汤普森棱镜Glan-Taylor prism 格兰-泰勒棱镜glass phase 玻璃相glassy ceramics 微晶玻璃glassy metal 玻璃态金属Glauber state Glauber态glide axis 滑移轴glide line 滑移线global phase 整体相位(量子信息)goniometer 测角器graded bandgap layer 缓变带隙层Gorter-Casimir two-fluid model 高特-卡西米尔二流体模型Graded index lens (GRIN) 梯度折射率透镜gradient of electric potential 电势梯度gram-molecule 克分子,摩尔(mole)grand free energy 巨自由能granular matter 颗粒物质granular superconductor 颗粒超导体granule 颗粒granularity 颗粒性granular metal 颗粒金属graphite 石墨graphite structure 石墨结构graph [线]图graph state 图态(量子信息)gravitational deflection of light 光线的引力偏折gravity acceleration 重力加速度Gray code 格雷码grazing angle 1,掠射角;2,擦边角greenhouse effect 温室效应group index of refraction 群折射率group theory 群论group velocity dispersion 群速度色散growth 生长growth model 生长模型guest host liquid crystal 宾主型液晶guided wave optics 导波光学gyroscopic effect 回转效应half metal 半金属half metallic magnet 半金属磁体half wave filter 半波滤波器half wave oscillator 半波振子half- wave zone method 半波带法half-wave voltage 半波电压Hall angle 霍尔角Hall coefficient 霍尔系数Hall field 霍尔电场[强度]Hall plateau 霍尔平台Hall resistance 霍尔电阻Hall voltage 霍尔电压halo nucleus 晕核halogen 卤素Hamiltonian matrix哈密顿[量]矩阵hard sphere 硬球hard sphere approximation 硬球近似harmonic generation 谐波产生Hartree-Fock electron 哈特里-福克电子H-center H心health physics 保健物理heat conductivity 1,导热性;2,热导率heat flow vector 热流矢量heat flux 热通量heat switch 热开关heavy electron 重电子heavy element 重元素heavy fermion superconductor 重费米子超导体heavy [fission] fragment 重【裂变】碎片heavy hole 重空穴heavy wall 重壁heavy water 重水hedgehog 猬缺陷height of potential barrier 势垒高度Heisenberg Hamiltonian 海森伯哈密顿量Heisenberg operators 海森伯算符Heisenberg uncertainty principle 海森伯不确定【性】原理Heitler-London theory 海特勒-伦敦理论Helfrich spontaneous curvature model 黑弗里希自发曲率模型helical spin order螺旋自旋有序helium liquefier 氦液化器heptahgedron 七面体Hermite polynomial 厄米多项式Hermitian matrix 厄米矩阵hertz (Hz) 赫兹, 频率单位heterotic superstring theory 杂化超弦理论Heusler alloy 霍伊斯勒合金hexadecapole 十六极hexahedron 六面体hexatic phase 六角相high coherence model 高相干模型high electron mobility transistor 高电子迁移率晶体管(简写:HEMT)high energy particle 高能粒子high-field domain 强场畴high-order dispersion 高阶色散high-order harmonic generation 高阶谐波产生high pass filter 高通滤波器high temperature reservoir 高温热源high temperature superconductor(HTS)高温超导体high vacuum 高真空high voltage electron microscopy 高压电子显微术Hohenberg-Kohn energy functional 霍恩伯格-科恩能量泛函hole-electron recombination 空穴-电子复合hole surface 空穴面(k空间中未占据态区的表面)hole-type high temperature superconductor 空穴型高温超导体holey fiber 多孔光纤hollow core optical fibers 空心光纤holon 空穴子homodyne零拍homodyne detection 零拍探测homolog[ue] 同系物homopolymer 单聚合物honeycomb photonic band gap fiber 蜂窝型光子带隙光纤hopping conductance 跳跃电导hopping energy,跳跃能hopping probability 跳跃概率hopping transport 跳跃输运host 基质host crystal 基质晶体,又称主晶hot carrier 热载流子h/e oscillation h/e振荡h/2e oscillation h/2e振荡Huang equations 黄[昆]方程组Huang-Rhys factor 黄昆-里斯因子Hubbard Hamiltonian 哈勃德哈密顿量Hubbard model 哈勃德模型Hubble time 哈勃时间hybrid bond 杂化键hybrid field effect 混合场效应hydrodynamics 流体[动]力学hydrodynamic mode 流体[动]力学模hydromagnetic disturbance 磁流体扰动hydromagnetic instability 磁流体不稳定性hydrophilic force 亲水力hydrophobic association 疏水缔合hydrophobic force 疏水力hyperbolic point 双曲点hypernucleus 超核hyper-Rayleigh scattering 超瑞利散射hyperspherical coordinate 超球座标hysteresis loop 1,滞后回线;2,磁滞回线hysteresis loss 1,滞后损失;2,磁滞损耗。

英文翻译 太阳能定日镜 2012

英文翻译 太阳能定日镜 2012

连续追踪的定日镜亚伯拉罕.克劳布斯陈晓东译摘要:即使太阳的运动是连续的,但是在当前的定日镜跟踪运动通常为不连续的步骤。

瞄准误差离散步骤通常大约1 毫弧度或更多。

为了有效降低跟踪误差提出了光滑连续跟踪。

就是在一个现有的定日镜中使用两个利用电子速度控制单元调整转速的交流电机。

持续跟踪系统在以色列魏兹曼科学中心定日镜领域成功的实现了。

定日镜运动的测量表明目标跟踪误差间隔几乎是消除的。

一种在定日镜运动和不稳定目标的不连续跟踪和连续跟踪的对比的报道。

关键词:跟踪误差;连续跟踪系统;交流电机;定日镜运动1 简介:许多现代的太阳能转化工艺要求的温度升高,温度一般都在1000℃以上,为了高效的操作这些高温过程需要高强度的阳光照射。

经常在相当大的程度上高于1000个太阳。

目前的技术领域的定日镜通常是用来生产平均浓度比值已达到1000。

虽然从理论上讲有这种可能性更高的浓度【1】。

由于几何误差在定日镜结构面误差和操作系统跟踪误差精度中存在,导致定日镜领域的表现是有限的。

比起原来的太阳能磁盘,他们能有效分散反射辐射在更大的角范围内,使得定日镜中集中而又潜在的错误减少了。

当在定日镜领域里不能提供必要的水平时,二级集中器就常常被用于增加聚光浓度【2-4】。

这包含额外的损失和附加的功率损失。

另外,因为定日镜的误差,二级集中器不能弥补失去的潜在热能。

所以这就有利于提高定日镜的性能,减少对二级集中器的依赖。

定日镜误差其中之一的跟踪误差存在于运转的定日镜中。

在定日镜跟踪控制中一般误差在1-1.5毫弧度,在某些特别的情况下会更高点【5】。

这些误差是油很多因素造成的。

比如:驱动公差以及齿侧间隙、太阳的位置精度模型、定日镜运动之间的间隔(由于编码器分辨率)、由重力和风造成的结构变形。

这些因素的影响可能有十分不同的时间段内。

类似于太阳模拟不正确以及重心转移的偏移的误差发生在相对较长的时间段内。

一个可能的解决方法来检测和消除这些误差是一个闭环控制系统【8】。

雅思考试阅读常见题材-11自然 detection of a meteorite lake

雅思考试阅读常见题材-11自然 detection of a meteorite lake

Detection of a meteorite LakeA A s the sun rose over picturesque Lake Bosumtwi, a team of SyracuseUniversity researchers prepared for another day of using state-of-the-art equipment to help unlock the mysteries hidden below the lake bottom.Nestled in the heart of Ghana (加纳),the lake holds an untapped reservoir of information that could help scientists predict future climate changes by looking at evidence from the past. This information will also improve the scientists’understanding of the changes that occur in a region struck by a massive meteorite (陨石).B T he project, led by earth sciences professor Christopher Scholz of the Collegeof Arts and Sciences and funded by the National Science Foundation (NSF), is the fi rst large-scale effort to study Lake Bosumtwi, which formed 1.1 million years ago when a giant meteor crashed into the Earth’s surface. The resulting crater is one of the largest and most well-preserved geologically young craters (火山口)in the world, says Scholz, who is collaborating on the project with researchers from the University of Arizona, the University of South Carolina, the University of Rhode Island, and several Ghanaian institutions. “Our data should provide information about what happens when an impact hits hard, pre-Cambrian (前寒武纪),crystalline rocks (结晶岩)that are a billion years old, he says.C E qually important is the fact that the lake, which is about 8 kilometers indiameter, has no natural outlet. The rim of the crater rises about 250 meters above the water’s surface. Streams flow into the lake, Scholz says, but the water leaves only by evaporation, or by seeping through the lake sediments.For the past million years, the lake has acted as a tropical rain gauge (测量器),fi lling and drying with changes in precipitation and the tropical climate.The record of those changes is hidden in sediment below the lake bottom. “The lake is one of the best sites in the world for the study of tropical climate (热带气候)changes,” Scholz says. “The tropics are the heat engine for the Earth’s climate. To understand global climate, we need to have records of climatechanges from many sites around the world, including the tropics.”D B efore the researchers could explore the lake’s subsurface, they needed aboat with a large, working deck area that could carry eight tons of scienti fi c equipment. The boat—dubbed R/V Kilindi—was built in Florida last year. It was constructed in modules that were dismantled, packed inside a shipping container, and reassembled over a 10-day period in late November and early December 1999 in the rural village of Abono, Ghana. The research team then spent the next two weeks testing the boat and equipment before returning to the United States for the holidays.E I n mid-January, fi ve members of the team—Keely Brooks, an earth sciencesgraduate student; Peter Cattaneo, a research analyst; and Kir am Lezzar, a postdoctoral scholar, all from SU; James McGill, a geophysical fi eld engineer;and Nick Peters, a Ph.D. student in geophysics from the University of Miami—returned to Abono to begin collecting data about the lake’s subsurface using a technique called seismic re fl ection pro fi ling. In this process, a high-pressure air gun is used to create small, pneumatic explosions in the water. The sound energy penetrates about 1,000 to 2,000meters into the lake’s subsurface beforebouncing back to the surface of the water.F T he reflected sound energy is detectedby underwater microphones—calledhydrophones—embedded in a 50-meter-long cable that is towed behind the boat asit crosses the lake in a carefully designedgrid pattern. On-board computers recordthe signals, and the resulting data are thenprocessed and analyzed in the laboratory.“The results will give us a good idea ofthe shape of the sediment are, and whenand where there were major changes insediment accumulation,” Scholz says. “Weare now developing three-dimensional perspective of the lake’s subsurface and the layers of sediment that have been laid down.”G T eam members spent about four weeks in Ghana collecting the data. Theyworked seven days a week, arriving at the lake just after sunrise. On a good day, when everything went as planned, the team could collect data and be back at the dock by early afternoon. Except for a few relatively minor adjustments, the equipment and the boat worked well. Problems that arose were primarily non-scientific— tree stumps, fishing nets, cultural barriers, and occasional misunderstandings with local villagers.H L ake Bosumtwi, the largest natural freshwater lake in the country, is sacred tothe Ashanti people, who believe their souls come to the lake to bid farewell to their god. The lake is also the primary source of fi sh for the 26 surrounding villages. Conventional canoes and boats are forbidden. Fishermen travel on the lake by fl oating on traditional planks (木板)they propel with small paddles (船桨).Before the research project could begin, Scholz and his Ghanaian counterparts had to secure special permission from tribal chiefs to put the R/V Kilindi on the lake.I W hen the team began gathering data, rumors (谣言)fl ew around the lake asto why the researchers were there. “Some thought we were dredging the lake for gold, others thought we were going to drain the lake or that we had bought the lake,” Cattaneo says. “But once the local people understood why we were there, they were very helpful.”Questions 14-18 .............................................................................Do the following statements agree with the information given in Reading Passage 1?In boxes 14-18 on your answer sheet, writeTRUE if the sataement agrees with the informationFALSE if the statement contradicts the informationNOT GIVEN if there is no information on this14W ith the analysis of the bottom of the lake, scientist will predict the climate changes in the future.15T he water stored in lake Bosumtwi was gone only by seeping through the lake sediments.16T he crater resulted from a meteorite impact is the largest and most preserved one in the world.17H istorical climate changes can be detected by the analysis of the sediment of the lake.18R esearch of scientist and co-workers had been interfered by the locals due to their indigenous believes.Questions 19-22 .............................................................................There are three steps of collecting data from the lake as followings, please filling the blanks in the Flow Chart below:Step 1a 50-meter 19 ,with many 20embededStep 2a 21 isneeded to create theexplosion into thewaterStep 3the 22enters deep intothe water andreturn backQuestions 23-27 ............................................................................. SummaryComplete the following summary of the paragraphs of Reading Passage, using no more than three words from the Reading Passage for each answer. Write your answers in boxes 23-27 on your answer sheet.The boat-double R/V Kilindi crossed the lake was dismantled and stored in a 23. The technology they used called 24; They created sound energy in to 1000-2000 metres in to the bottom of the lake, and used separate equipment to collect the returned waves. Then the data had been analyzed and processed in the 25. Scholz also added that they were now building 26view of the sediment or sub-image in the bottom of the lake. Whole set of equipment works well yet the ship should avoid tree stumps or 27fl oating on the surface of the Bosumtwi lake.。

层序地层学的名词英语解释

层序地层学的名词英语解释

层序地层学的名词英语解释GlossaryAccommodation. Another term for relative sea-level. Can be thought of as the space in which sediments can fill, defined at its base by the top of the lithosphere and at its top by the ocean surface.Basinward Shift in Facies. When viewed in cross-section, a shifting of all facies towards the center of a basin. Note that this is a lateral shift in facies, such that in vertical succession, a basinward shift in facies is characterized by a shift to shallow facies (and not a vertical shift to more basinward or deeper-water facies).Bed. Layer of sedimentary rocks or sediments bounded above and below by bedding surfaces. Bedding surfaces are produced during periods of nondeposition or abrupt changes in depositional conditions, including erosion. Bedding surfaces are synchronous when traced laterally; therefore, beds are time-stratigraphic units. See Campbell, 1967 (Sedimentology 8:7-26) for more information.Bedset. Two or more superposed beds characterized by the same composition, texture, and sedimentary structures. Thus, a bedset forms the record of deposition in an environment characterized by a certain set of depositional processes. In this way, bedsets are what define sedimentary facies. Equivalent to McKee and Weir's coset, as applied to cross-stratification. See Campbell, 1967 (Sedimentology 8:7-26) for more information.Condensation. Slow net rates of sediment accumulation. Stratigraphic condensation can occur not only through a cessation in the supply of sediment at the site of accumulation, but also in cases where the supply of sediment to a site is balanced by the rate of re moval of sediment from that site. Where net sediment accumulation rates are slow, a variety of unusual sedimentologic features may form, including burrowed horizons, accumulations of shells, authigenic minerals (such as phosphate, pyrite, siderite, glauconite, etc.), early cementation and hardgrounds, and enrichment in normally rare sedimentary components, such as volcanic ash and micrometeorites.Conformity. Bedding surface separating younger from older strata, along which there is no evidence of subaerial or submarine erosion or nondeposition and along which there is no evidence of a significant hiatus. Unconformities (sequence boundaries) and flooding surfaces (parasequence boundary) will pass laterally into correlative conformities, most commonly in deeper marine sediments.Eustatic Sea Level. Global sea level, which changes in response to changes in the volume of ocean water and the volume of ocean basins.Flooding Surface. Shortened term for a marine flooding surface.Highstand Systems Tract. Systems tract overlying a maximum flooding surface, overlain by a sequence boundary, and characterized by an aggradational to progradational parasequence set.High-Frequency Cycle. A term applied to a cycle of fourth order or higher, that is, having a period of less than 1 million years. Parasequences and sequences can each be considered high-frequency cycles when their period is less than 1 million years. Isostatic Subsidence. Vertical movements of the lithosphere as a result of increased weight on the lithosphere from sediments, water, or ice. Isostatic subsidence is a fraction of the thickness of accumulated material. For example, 100 meters of sediment will drive about 33 meters of subsidence (or less, depending on the rigidity of the lithosphere).Lowstand Systems Tract. Systems tract overlying a type 1 sequence boundary, overlain by a transgressive surface, and characterized by a progradational to aggradational parasequence set.Marine Flooding Surface. Surface separating younger from older strata, across which there is evidence of an abrupt increase in water depth. Surface may also display evidence of minor submarine erosion. Forms in response to an increase in water depth.Maximum Flooding Surface. Marine flooding surface separating the underlying transgressive systems tract from the overlying highstand systems tract. This surface also marks the deepest water facies within a sequence. This flooding surface lies at the turnaround from retrogradational to progradational parasequence stacking, although this turnaround may be gradational and characterized by aggradational stacking. In this case, a single surface defining the point of maximum flooding may not be identifiable, and a maximum flooding zone is recognized instead. The maximum flooding surface commonly, but not always, displays evidence of condensation or slow deposition, such as burrowing, hardgrounds, mineralization, and fossil accumulations. Because other flooding surfaces can have evidence of condensation (in some cases, more than the maximum flooding surface), condensation alone should not be used to define the maximum flooding surface. Meter-Scale Cycle. A term applied to a cycle with a thickness of a couple of meters or less. Parasequences and sequences can each be considered meter-scale cycles when they are thinner than a couple of meters.Parasequence. Relatively conformable (that is, containing no major unconformities), genetically related succession of beds or bedsets bounded by marine-flooding surfaces or their correlative surfaces. Parasequences are typically shallowing-upward cycles.Parasequence Boundary. A marine flooding surface.Parasequence Set. Succession of genetically related parasequences that form a distinctive stacking pattern, and typically bounded by major marine flooding surfaces and their correlative surfaces. Parasequence set boundaries may coincide with sequence boundaries in some cases. See progradational, aggradational and retrogradational parasequence sets.Peritidal. All of those depositional environments associated with tidal flats, including those ranging from the highest spring tides to somewhat below the lowest tides.Relative Sea Level. The local sum of global sea level and tectonic subsidence. Locally, a rise in eustatic sea level and an increase in subsidence rates will have the same effect on accommodation. Likewise, a fall in eustatic sea level and tectonic uplift will have the same effect on accommodation. Because of the extreme difficulty in teasing apart the effects of tectonic subsidence and eustatic sea level in regional or local studies, sequence stratigraphy now generally emphasizes relative changes in sea level, as opposed to its earlier e mphasis on eustatic sea level.Sequence. Relatively conformable (that is, containing no major unconformities), genetically related succession of strata bounded by unconformities or their correlative conformities.Sequence Boundary. Form in response to relative falls in sea level.Sequence Stratigraphy. The study of genetically related facies within a frameworkof chronostratigraphically significant surfaces.Shelf Margin Systems Tract. Systems tract overlying a type 2 sequence boundary, overlain by a transgressive surface, and characterized by a progradational to aggradational parasequence set. Without regional seismic control, most shelf margin systems tracts may be unrecognizable as such and may be inadvertently lumped with the underlying highstand systems tract as part of one uninterrupted progradational parasequence set. If this occurs, the overlying transgressive surface may be erroneously inferred to also be a type 1 sequence boundary.Systems Tract. Linkage of contemporaneous depositional systems, which are three-dimensional assemblages of lithofacies. For example, a systems tract might consist of fluvial, deltaic, and hemipelagic depositional systems. Systems tracts are defined by their position within sequences and by the stacking pattern of successive parasequences. Each sequence consists of three systems tract in a particular order. For a type 1 sequence, these are the lowstand, transgressive, and highstand systems tracts. For a type 2 sequence, these are the shelf margin, transgressive, and highstand systems tracts.Tectonic Subsidence. Vertical movements of the lithosphere, in the absence of any effects from changes in the weight of overlying sediments or water. Also called driving subsidence. Tectonic subsidence is generated primarily by cooling, stretching, loading (by thrust sheets, for example), and lateral compression of the lithosphere.Transgressive Surface. Marine flooding surface separating the underlying lowstand systems tract from the overlying transgressive systems tract. Typically, this is the first major flooding surface following the lowstand systems tract. In depositionally updip areas, the transgressive surface is commonly merged with the sequence boundary, with all of the time represented by the missing lowstand systems tract contained within the unconformity. The transgressive surface, like all of the major flooding surfaces within the transgressive systems tract, may display evidence of stratigraphic condensation or slow net deposition, such as burrowed surfaces, hardgrounds, mineralization, and fossil accumulations.Transgressive Systems Tract. Systems tract overlying a transgressive surface, overlain by a maximum flooding surface, and characterized by a retrogradational parasequence set.Type 1 Sequence Boundary. Characterized by subaerial exposure and associated erosion from downcutting streams, a basinward shift in facies, a downward shift in coastal onlap, and onlap of overlying strata. Forms when the rate of sea-level fall exceeds the rate of subsidence at the depositional shoreline break (usually at base level or at sea level). Note that this means that if such changes can be observed in outcrop and the underlying strata are marine, then the boundary is a type 1 sequence boundary.Type 2 Sequence Boundary. Characterized by subaerial exposure and a downward shift in onlap landward of the depositional shoreline break (usually at base level or at sea level). Overlying strata onlap this surface. Type 2 sequence boundaries lack subaerial erosion associated with the downcutting of streams and lack a basinward shift in facies. Forms when the rate of sea-level fall is less than the rate of subsidence at the depositional shoreline break. Note that the lack of a basinward shift in facies and the lack of a relative fall in sea level at the depositional shoreline break means that there are essentially no criteria by which to recognize a type 2 sequence boundary in outcrop.Unconformity. Surface separating younger from older strata, along which there is evidence of subaerial erosional truncation or subaerial exposure or correlative submarine erosion in some areas, indicating a significant hiatus. Forms in response to a relative fall in sea level. Note that this is a much more restrictive definition of unconformity than is commonly used or used in earlier works on sequence stratigraphy (e.g., Mitchum, 1977).Walther's Law states that "...only those facies and facies areas can be superimposed, without a break, that can be observed beside each other at the present time" (Middleton translation from German). At a Waltherian contact, one facies passes gradationally into an overlying facies, and those two facies represent sedimentary environments that were originally adjacent to one another.Water Depth. The distance between the sediment surface and the ocean surface. Water depth is reflected in sedimentary facies. A very large number of studies that purport to describe sea-level changes (both eustatic and relative) are actually only describing changes in water depth. The effects of isostatic subsidence and compaction must be removed from water depth to calculate relative sea level. This is typically done through backstripping. To calculate eustatic sea level, the rate of tectonic subsidence must then be subtracted from the relative sea-level term.。

The Detection and Characterization of cm Radio Continuum Emission from the Low-mass Protost

The Detection and Characterization of cm Radio Continuum Emission from the Low-mass Protost

a r X i v :0705.1747v 1 [a s t r o -p h ] 12 M a y 2007Draft version February 1,2008Preprint typeset using L A T E X style emulateapj v.03/07/07THE DETECTION AND CHARACTERIZATION OF CM RADIO CONTINUUM EMISSION FROM THELOW-MASS PROTOSTAR L1014-IRSYancy L.Shirley 1,Mark J.Claussen 2,Tyler M.Bourke 3,Chadwick H.Young 4,Geoffrey A.Blake 5Draft version February 1,2008ABSTRACTObservations by the Cores to Disk Legacy Team with the Spitzer Space Telescope have identified a low luminosity,mid-infrared source within the dense core,Lynds 1014,which was previously thought to harbor no internal source.Followup near-infrared and submillimeter interferometric observations have confirmed the protostellar nature of this source by detecting scattered light from an outflow cavity and a weak molecular outflow.In this paper,we report the detection of cm continuum emission with the VLA.The emission is characterized by a quiescent,unresolved 90µJy 6cm source within 0.′′2of the Spitzer source.The spectral index of the quiescent component is α=0.37±0.34between 6cm and 3.6cm.A factor of two increase in 6cm emission was detected during one epoch and circular polarization was marginally detected at the 5σlevel with Stokes V/I =48±16%.We have searched for 22GHz H 2O maser emission toward L1014-IRS,but no masers were detected during 7epochs of observations between June 2004and December 2006.L1014-IRS appears to be a low-mass,accreting protostar which exhibits cm emission from a thermal jet or a wind,with a variable non-thermal emission component.The quiescent cm radio emission is noticeably above the correlation of 3.6cm and 6cm luminosity versus bolometric luminosity,indicating more radio emission than expected.In this paper,we characterize the cm continuum emission in terms of observations of other low-mass protostars,including updated correlations of centimeter continuum emission with bolometric luminosity and outflow force,and discuss the implications of recent larger distance estimates on the physical attributes of the protostar and dense molecular core.Subject headings:radiation mechanisms:thermal,non-thermal —radio continuum:stars —stars:for-mation1.INTRODUCTIONIt is extremely difficult to identify the incipient stages of low-mass (≈1M sun )star formation be-cause dense molecular cloud cores obscure nascent pro-tostars.Submillimeter dust continuum surveys (e.g.,Ward-Thompson et al.1994,Shirley et al.2000,Visser et al.2002,Kirk et al.2005)have identified several dense cores with no apparent internal sources,based on the lack of an IRAS point source and the diffuse na-ture of submillimeter dust emission.Observations by the Cores to Disk Legacy Team (c2d)with the Spitzer Space Telescope have identified a few mid-infrared sources that are embedded near the submillimeter continuum peaks of previously classified starless cores (e.g.Young et al.2004,Bourke et al.2006).These new objects are of low-luminosity (L int ≤0.1L ⊙)and presumably low-mass since they were not previously detected by IRAS.Some of these objects may be in the earliest stages of accre-tion.These newly identified low-mass,low-luminosity protostars warrant detailed follow-up studies to deter-mine their evolutionary status.The first newly identified object detected in the c2d survey,L1014-IRS (Young,et al.2004),was mod-eled as a very low-luminosity (L int <0.1L ⊙),low-mass1Bart J.Bok Fellow,Steward Observatory,University of Ari-zona,933Cherry Ave.,Tucson,AZ 857212NRAO,P.O.Box 0,1003Lopezville Road,Socorro,NM 878013Harvard-Smithsonian Center for Astrophysics,60Garden St.MS42,Cambridge,MA 021384Nicholls State University,Thibodaux,LA 703105Division of Geological and Planetary Sciences 150-21,Califor-nia Institute of Technology,Pasadena,CA 91125(M <0.1M ⊙)object embedded within the Lynds 1014dark cloud (Lynds 1962)at a distance of approximately 200pc.This object has been classified as a VeLLO (Very Low-Luminosity Object)by the c2d team:an object with an internal protostellar luminosity ≤0.1L ⊙that is di-rectly associated with a dense molecular core.The recent study of Morita et al.(2006)has suggested a revised dis-tance estimate of 400to 900pc based on the possible age ranges of nearby T-Tauri stars that are spatially within 2◦of the L1014dense core.However,it is not clear that these T-Tauri stars are directly associated with L1014.Determining the evolutionary state of L1014-IRS has been the subject of several follow-up studies.The large scale molecular distribution in the dense core was deter-mined by the single-dish mapping survey of Crapsi et al.(2005).No evidence for a large scale CO outflow was detected;however,a molecular outflow was detected on small scales with the SMA (Bourke et al.2005).Near-infrared observations detect scattered light,presumably from the outflow cone,at 1.6µm and 2.2µm (Huard et al.2006).The SMA-detected CO outflow is aligned with the direction of the near-infrared scattered light cavity.High resolution molecular observations with BIMA indi-cate that the protostar is not at the peak of the molec-ular and dust column density in the core,but offset by about 8′′in the plane of the sky (Lai et al.,in prepa-ration).This offset is also seen in (sub)millimeter con-tinuum maps (Young et al.2004)and the near-infrared extinction map (Huard et al.2006).Despite these significant observational efforts,a single,consistent picture of the evolutionary state of L1014-IRS has not emerged.In order to better characterize the2physical nature of L1014-IRS,we have conducted cen-timeter radio continuum observations using five array configurations of the Very Large Array 6(VLA).Cen-timeter radio continuum emission is well correlated with the luminosity of protostellar sources (Anglada 1995)and is thought to arise from shock ionization from protostel-lar winds (Ghavamian &Hartigan 1998),from interac-tion of the protostellar jets with dense gas in the interface of the outflow cavity (Curiel et al.1987,1989,Shang et al.2004),or from accretion shock-driven photoionization (Neufeld &Hollenbach 1996).In this paper we report the detection and characteriza-tion of the cm radio continuum emission toward L1014-IRS (§3.1).We compare the detected cm emission with observations of other low-mass protostars (§4.1).We dis-cuss the non-detections of 22GHz H 2O masers and com-pare our upper limits to the recent maser monitoring surveys of low-mass protostars (§4.2).We also compare the derived source properties from the diverse studies of the protostar and dense core in terms of the range of distance estimates to L1014(§4.3).2.VLA OBSERVATIONSL1014-IRS was observed during 10epochs in 5ar-ray configurations (D,A,BnA,B,and C)with the Very Large Array (Table 1).All observations were centered on the published Spitzer mid-infrared source (α=21h 24m 07s .51,δ=+49◦59′09.′′0,J2000.0).Con-tinuum observations were made at 3.6cm,and 6.0cm,with two polarization pairs at adjacent frequencies,pro-viding a total equivalent bandwidth of 172MHz.We also attempted to detect H 2O masers by observing the J K a K c =616→523transition at 22.23508GHz with,typically,24.4kHz spectral resolution (0.3km/s)span-ning ±20km/s velocity coverage.For the 3.6and 6.0cm data,the data were reduced independently using the standard routines in AIPS++and AIPS .Complex gain calibration was performed by switching to the nearby quasar 2137+510,2.4◦from L1014-IRS,on time-scales of 15to 30minutes (Hamaker,Bregman,&Sault 1996a,b).The absolute flux density and bandpass calibration were determined from observa-tions of the quasars 3C48and 3C283.The Stokes I and V images were deconvolved using the Cotton-Schwab algorithm (e.g.,Schwab 1984)and Clark-H¨o gbom algorithm (H¨o gbom 1974,Clark 1980)with a few thousand iterations and interactive CLEAN regions.Imaging the L1014-IRS field was difficult due to the presence of several bright sources within the VLA pri-mary beam (Figure 1a).Special care had to be taken in the CLEANing process (e.g.,Cornwell,Braun,&Briggs 1999)and multiple reductions with variations in the CLEAN parameters were performed.We have checked the consistency of our images by also reducing the data with the standard AIPS routines,and the fluxes agree within the statistical errorbars.Generally,the images are made with natural weighting of the visibilities (Briggs 1995);however,uniform weighting was used to obtain better angular resolution for the full track (9hour)ob-servations on the days of July 1,2004and August 21,6The VLA is operated by NRAO.The National Radio Astron-omy Observatory is a facility of the National Science Foundation operated under a cooperative agreement by Associated Universi-ties,Inc.2004.3.RESULTS3.1.Radio Continuum DetectionsWe detected cm radio continuum emission from a source within 0.′′2of the Spitzer mid-IR source using the VLA at 3.6cm and 6cm (Figure 1).The initial detec-tions were made during 9hour tracks at 3.6cm and 6cm with the VLA in the D configuration.Subsequent observations detected the source at 6cm in the other three VLA configurations,and again at 3.6cm in D con-figuration.All of the 6cm observations,except for the initial 6cm detection on August 21,2004,indicate a constant flux density source with an average 6cm flux of 88±11µJy (8σ).The source is unresolved in all array configurations.The two detections at 3.6cm are also in agreement with each other despite being separated by 17months.The average 3.6cm flux is 111±8µJy (14σ).The spectral index between two wavelengths (λ2>λ1)is defined asα=ln(S λ1/S λ2)3to detect that variation due to the higher noise level in progressively shorter time blocks.A Stokes V source (circularly polarized)was marginally detected at the 5σlevel at 6cm on August 21,2004with a flux of 84±17µJy at the position of the 6cm Stokes I source (Figure 1d).The fraction of circular polarization,f c =Stokes V/I =48±16%,is quite high if the detection is significant.Since the Stokes V result was obtained during a full single track of VLA observa-tions,it will be difficult to confirm until L1014-IRS is observed with wider bandwidth (e.g.,with the eVLA).If the circular polarization signal is real,then this observa-tion indicates that the radio emission observed on August 21,2004must originate from a non-thermal mechanism (§4.1.2).3.2.Water Maser SearchWe searched for H 2O masers during seven epochs span-ning 22months by observing the J K a K c =616→523transition at 22.23508GHz.No H 2O masers were de-tected at any epoch.The combined 1σrms of the non-detection is 3.1mJy/beam with a channel spacing of 24.4kHz and a total bandwidth of 3.125MHz (∼40km/s).The individual observations are summarized in Table 1.4.DISCUSSION4.1.Centimeter Radio Continuum Emission4.1.1.Steady ComponentCentimeter continuum emission has been detected to-ward many but not all high-mass and low-mass proto-stars.The emission is most commonly thought to origi-nate from bremsstrahlung (free-free)emission from ion-ized gas,although some protostellar objects also display non-thermal emission.For high-mass protostars,the ion-ization mechanism is photoionization usually in the form of an embedded HII region (Churchwell 1990).For pro-tostars that are later in spectral type than B,the ioniz-ing radiation from the star is not enough to significantly photoionize the surrounding envelope and an alternative mechanism is needed to explain the observed emission (Rodr´ıguez et al.1989,Anglada 1995).Since nearly all low-mass,embedded protostars are known to have molecular outflows (e.g.Wu et al.2004),the ioniza-tion is postulated to arise from shocks generated from a jet (e.g.,Cohen,Bieging,&Schwartz 1982,Bieging &Cohen 1989,Curiel et al.1987,1989,Rodr´ıguez &Reipurth 1996,Shang et al.2004).Direct evidence for this hypothesis comes from observations using interfer-ometers at high angular resolution.Elongated centime-ter continuum emission regions are observed with the same orientation as the large-scale molecular outflow to-ward a few protostars (e.g.,Anglada 1995;Bontemps,Ward-Thompson,&Andr´e 1996).In addition,the ob-served outflow force (M ⊙km/s/yr)theoretically provides enough energy in the shock to explain the observed cen-timeter fluxes toward most low-mass protostars (Cabrit &Bertout 1992,Skinner et al.1993,Anglada 1995).Finally,the radio spectral index of many protostellar sources is consistent with optically thin (α≈−0.1)to partially optically thick free-free emission between 3.6cm and 6.0cm (2.0>α>−0.1,e.g.,Anglada et al.1998,Beltr´a n et al.2001).The centimeter continuum luminosity (e.g.,L 3.6=S 3.6D 2mJy kpc 2)of low-mass and intermediate-mass protostellar sources was first cataloged from the litera-ture in the review by Anglada (1995).Anglada plotted the 3.6cm luminosity against the bolometric luminosity of protostars with L bol <103L ⊙and found a well corre-lated relationship (r =0.79),L 3.6=10−2.1(L bol /1L ⊙)0.7mJy kpc 2.This relationship has formed the basis for predictions of the amount of centimeter emission that is expected in searches for new low-mass protostars (e.g.,Harvey et al.2002,Stamatellos et al.2007).The 3.6cm luminosity correlation is directly related to the well established correlation of outflow force vs.bolometric lu-minosity (Bontemps et al.1996,Wu et al.2004);higher luminosity sources drive more powerful outflows that re-sult in a larger degree of shock ionization and therefore a larger 3.6cm luminosity (Anglada 1995).Since 1995,many more centimeter observations have been made and the spectral coverage of the photometry of protostars has increased.We have used the detailed summary tables of Furuya et al.(2003;Table 4)and Anglada (1995)supplemented by the surveys of Eiroa et al.(2005)and Anglada et al.(1998)to catalog the 3.6cm,6.0cm,and bolometric luminosities of detected protostellar sources.We updated the bolometric lumi-nosity of sources observed in the submillimeter surveys of Shirley et al.(2000),Mueller et al.(2003),and Young et al (2003).The resulting sample of 58sources at 3.6cm and 40sources at 6.0cm are plotted in Figure 2.We find updated correlations oflog(L 3.6/1mJy kpc 2)=−(2.24±0.03)+(0.71±0.01)log(L bol /1L ⊙)(2)log(L 6.0/1mJy kpc 2)=−(2.51±0.03)+(0.87±0.02)log(L bol /1L ⊙),(3)with correlation coefficients of r =0.66and r =0.74respectively.This sample is not complete as there are many more protostellar centimeter detections for which no L bol has been published.Nevertheless,we have up-dated the correlation of Anglada (1995)with twice as many points at 3.6cm and plotted the correlation at 6.0cm for the first time.For comparison,we have plotted the 3.6cm and 6.0cm luminosities of L1014-IRS in Figure 2at the distances of 200,400,and 900pc.L1014-IRS is above the correlation at all distances indicating that we have detected more centimeter continuum flux than expected.The spectral index of sources between 3.6cm and 6.0cm is used to argue for the interpretation that the emission mechanism is consistent with partially optically thick free-free emission.Optically thin free-free emission is expected to have a spectral index α=−0.1at cen-timeter wavelengths.In the optically thick limit,αap-proaches 2.0with intermediate values indicative of par-tially optically thick plasmas.We have also plotted the spectral index of sources from the literature that have been detected at both wavelengths and have a published L bol determination in Figure 2.No correlation of αis observed with bolometric luminosity,probably indicat-ing that the optical depth associated with the jet’s shock ionization is not dependent on the total protostellar lu-minosity or the strength of the molecular outflow.The median spectral index is α=0.5,with most protostellar sources having flat or positive spectral indices.This me-dian value is close to the result expected for an ionized wind or jet with a 1/r 2density gradient (e.g.,Panagia &Felli 1975,Wright &Barlow 1975,Reynolds 1986).Unfortunately,most of the sources in Figure 2were not4observed at both wavelengths on the same day and vari-ability may result in significant scatter in the plot.The spectral index of the quiescent emission of L1014-IRS(0.37±0.34)is consistent with ionized free-free emis-sion with a density gradient.The spectral index agrees well with the median of the ensemble of protostellar sources measured.We have also updated the correlation between outflow force,F out(M⊙km/s/yr),and centrimetric luminosity of Anglada(1995)using the molecular outflow compilations of Bontempts et al.(1996),Furuya et al.(2003),and Wu et al.(2004).The updated correlation of44sources is weak(r=0.55),log(F out/1M⊙km/s/yr)=−(3.15±0.07)+(0.67±0.03)log(L3.6/1mJy kpc2).(4) This correlation is interpreted as evidence that jets from the molecular outflow provide enough shock ionization to account to the observed centimeter continuum emission (e.g.,Anglada1995).Assuming maximum ionization ef-ficiency,the minimum outflow force needed to create the observed level of3.6cmflux was estimated by Curiel et al.(1987,1989)to be F out=10−3.5(L3.6/1mJy kpc2). This level is shown as a dashed line in Figure2d.The observed outflow force toward L1014-IRS ranges from 0.04−2.9×10−6M⊙km/s/yr for distances of200to 900pc(see Bourke et al.2005).At a distance of200pc, the upper limit on the outflow force is a factor of2lower than the Curiel theoretical minimum outflow force.The upper limit in the observed outflow force includes esti-mates of the missingflux due to interferometric spatial filtering as well as the average opacity corrections for low-mass protostellar outflows(see Bourke et al.2005). Given the range of uncertainty in these estimates,the observed outflow force is below,but not necessarily in-consistent with the theoretical minimum outflow force needed to produce the observed3.6cm luminosity.The disagreement between the observed outflow force is more pronounced at larger distances,increasing to an order of magnitude below the theoretical curve for a distance of 900pc.The small observed molecular outflow force and the large observed centimeter continuum luminosities may indicate that another ionization mechanism is operating in L1014-IRS.There are a few other possibilities that have been discussed in the literature.We shall analyze the viability of two popular possibilities.Radiative transfer modeling of emission at24and70µm indicate aflux excess from a disk around L1014-IRS (Young et al.2004).Neufeld&Hollenbach(1996)pos-tulated that the supersonic infall of material onto a pro-tostellar disk will create an accretion shock with enough ionization to generate∼1mJy of continuum emission at centimeter wavelengths at distances of≈200pc.How-ever,this mechanism does not appear to be able to pro-vide enough ionization to explain the emission observed toward L1014-IRS.In order to produce90µJy emission at6cm at a distance of200pc,the protostellar mass has to be>2M⊙and the accretion rate onto the disk must be>10−4M⊙/yr(see Figure2of Neufeld&Hollenbach 1996).Estimates of the protostellar mass are very un-certain and highly distance dependent,however,2M⊙is likely larger than the L1014-IRS protostar and disk mass, even for the far distance estimate of900pc(see Young et al.2004).Furthermore,this accretion rate is an order of magnitude larger than the range of inferred accretion rates from the observed outflow momentumflux and the modeled internal luminosity of the source(≤3×10−5 M⊙/yr;Bourke et al.2006,Young et al.2004).Ion-ization from an accretion shock does not appear to be a likely explanation.A more plausible possibility is that there is a spheri-cal wind component.The expected continuum emission from self-shocked spherical winds have been modeled by numerous authors(e.g.,Panagia&Felli1975,Wright& Barlow1975,Reynolds1986,Gonz´a lez&Cant´o2002). The recent study of Gonz´a lez&Cant´o model a time variable wind that generates internal shocks(Raga et al. 1990)which produce ionization and centimeter contin-uum emission(see Ghavamian&Hartigan1998).Their models produce centimeter continuum emission of sev-eral hundreds ofµ-Janskys and spectral indices that are positive for mass loss rates of10−6M⊙/yr at a distance of150pc.Accounting for the larger distance estimates of200to900pc and potentially lower mass loss rates for L1014-IRS,then this type of emission may still account for theflux observed at3.6cm and6cm(∼100µJy). However,the spherical wind models is usually applied to evolved protostellar sources that are Class II(classi-cal T-Tauri stars)or later(e.g.,Evans et al.1987).If L1014-IRS is an older,more evolved protostar,then this may not be a problem(§4.3).4.1.2.Variable ComponentWhile the quiescent component of L1014-IRS is constant over5epochs with a positive spectral in-dex,the emission properties of L1014-IRS were signif-icantly different during the single epoch of August21, 2004.The6cmflux was larger by a factor of two, S6.0(21AUG2004)=173±16µJy.Unfortunately,the spectral index of elevated emission was not determined since L1014-IRS was observed at a single wavelength. Circular polarization was detected at the5σlevel indi-cating non-thermal emission.Thus,L1014-IRS has vari-able centimeter emission,although the timescale of the variability is not constrained since it was seen to vary during only a single epoch.In general variability of centimeter continuum sources has not been properly addressed since observations of sources are limited to a few epochs.There has not been a systematic,monthly monitoring campaign of deeply embedded sources to characterize their centimeter vari-ability;however,there is observational evidence for vari-ability among deeply embedded protostellar sources.For instance,the Class0source,B335,is known to vary be-tween an upper limit of<80µJy(1994December)and 390µJy(2001January)at3.6cm(Avila et al.2001, Reipurth et al.2002).Another example is the variability and purported evidence for jet precession of the centime-ter continuum sources toward IRAS16293(Chandler et al.2005).Variability contributes scatter(de-correlation) of the luminosity correlations shown in Figure2.In ad-ditional to the known variable thermal sources,several non-thermal protostellar sources are known to be highly variable(e.g.T Tauri stellarflares,see White1996); but,most of those objects are more evolved than deeply embedded protostars.The observed increase in emission on August21,2004 and5σStokes V detection is indicative of variable non-5thermal emission toward L1014-IRS.While there have been a few high-mass protostars with observed nega-tive spectral indices(e.g.,Reid et al.1995;Garay et al.1996),most embedded(≤Class I)low-mass proto-stars with cm radio emission have positive spectral in-dices(see Figure2).There are only a few embedded low-mass protostars toward which negative spectral in-dex,non-thermal emission is detected(e.g.,Shepherd &Kurtz1999,Girart et al.2002).One case,R CrA IRS5,was detected with significant circular polarization (Feigelson et al.1998).The authors postulate that the emission is due to gyrosynchrotron emission and may originate from magnetic reconnection events associated withflares(Feigelson et al.1998).The physical mecha-nism for generating the radioflare is still not well under-stood(e.g.,Basri2004)and it is questionable whether it applies to the embedded phase of low-mass protostars (see below).A second case,IRAS19243+2350,is a steep spectrum non-thermal source(α=−0.82±0.04)that is elongated in the direction of its CO outflow(Girart et al.2002).Girart et al.postulate that the emission may originate from a bi-conical synchrotron source tracing the protostellar jet,similar to observations of the high-mass source W3(OH)(Reid et al.1995,Wilner,Reid,& Menten1999).Non-thermal emission may be present in low-mass protostellar jets,but the level of emission may be dominated by the thermal,shock-ionized component of the jet.In the case of L1014-IRS,this non-thermal jet component would have to be variable.We shall dis-cuss several possibilities for the origin of the observed non-thermal emission toward L1014-IRS.In order to better understand the non-thermal emis-sion mechanism,we estimate the brightness tempera-ture of the emission to be T b=Sνλ2d2/2kR2emit= 0.14−2.8×106K for distances of200to900pc and an emitting region that is R emit=1AU in size.The brightness temperature is very sensitive to the assumed size of the emitting region.A R emit of1AU is appro-priate for a smallflare;but,smaller when compared to the solar coronal emitting region which is typically less than≈5AU(e.g.Leto et al.2000).Since the emission observed toward L1014-IRS was unresolved,even in the A-array configuration,then we can only limit the size of the emission region to<90(D/200pc)AU.For instance, if we assumed that the emitting region was45AU(half of our A-array resolution),then the brightness temper-ature drops to<100K.This is too low even for the steady,thermal free-free component(T∼104K),unless the emission was very optically thick.This cannot be the case sinceτ>>1would implyαapproaching2.0 which is not observed.While the size of the emitting re-gion is severely unconstrained,the detection of circular polarization indicates non-thermal emission probably on small(<few AU)size scales,most likely due to gyrosyn-chrotron emission(e.g.,Ramaty1969;Dulk&Marsh 1982).Radioflares are observed toward very low-mass ob-jects including brown dwarfs(Berger et al.2002,Os-ten et al.2006),late M dwarfs(Berger2006),and T-Tauri stars(Bieging&Cohen1989;White,Pallavicini, &Kundu1992).The typical level of quiescent emis-sion toward low-mass stars and brown dwarfs is100µJy and theflares are1mJy(for nearby distances of≈30pc)with significant circular polarization(f c≈50%) detected during theflaring events.The spectral lumi-nosity at6cm of L1014-IRS during the event is Lν= 4πD2Sν=4.7×1015(D/200pc)2erg s−1Hz−1.This is about2400times larger than the most luminous ob-served solarflares(Bastian2004)and about200times larger than theflares observed by Berger(2006)toward late M dwarfs.If the radio emission is due to aflare,it must be a powerfulflare since L1014-IRS is at least20 times farther away than the average distance of sources detected by Berger( d =10.6±5.1pc).However,it is not larger than the typical non-thermalflaring emission observed toward young T Tauri stars of Lν=1015−1018 erg s−1Hz−1(G¨u del2002).The ratio of spectralflare luminosity to bolometric luminosity,Lν/L bol=4×10−18 Hz−1,is also similar to the ratio observed toward classi-cal T-Tauri stars(see G¨u del2002).The timescale over which radioflares toward low-mass stars are observed tends to occur over minutes to hours. For instance,the low-mass brown dwarf,LP944-2,dis-covered by the2001NRAO summer students(Berger et al.2002),displaysflaring activity with an average timescale of10to15minutes.This is very different from the activity observed toward L1014-IRS on August21, 2004.We detected no evidence for short-term variabil-ity within our detected emission.The source appears to have a nearly constantflux that is twice as high as the steady component for at least an8hour period.This elevated emission then appears to be longer in duration than that observed duringflaring events toward low-mass (proto)stars(G¨u del2002).An alternative possibility is that the elevated emission is not due to aflare,but due to rotational modulation of a non-thermal component associated with the magnetic connection between disk and accretion onto the star(i.e., Bieging&Cohen1989).Such a mechanism has been postulated for the T-Tauri star,V410Tauri,with a rota-tional period of1.9days.The observed emission toward V410is1mJy with a negative spectral index.If the accretion spot is blocked from view for a fraction of the stellar rotation period,then it is possibly that we could have observed the source with the accretion spot is in view on August21,2004and with the accretion spot blocked from view during the other epochs.The rota-tional period of L1014-IRS must be longer than8hours since elevated emission was observed during the entire8 hour track.A negative spectral index was observed to-ward V410Tauri,while a negative spectral has not been observed toward L1014-IRS.This hypothesis is highly speculative and would require a regular monitoring cam-paign to test.Unfortunately,it is currently not possible to strongly constrain the origin of the non-thermal component to-ward L1014-IRS.The expanded bandwidth of the eVLA is needed to permit a systematic monitoring campaign with high enough signal-to-noise in only a few hour ob-servations to routinely check for a Stokes V detection and to determine an instantaneous spectral index.4.2.Water Maser Non-detectionsA compact,weak molecular outflow has been detected toward L1014-IRS(Bourke et al.2005);therefore it may be possible to detect water masers if the jet is impinging on dense knots of material near the protostar.Previous。

Photosynthetic suspended-growth systems in aquaculture

Photosynthetic suspended-growth systems in aquaculture

Photosynthetic suspended-growth systems in aquacultureJohn A.Hargreaves *Louisiana State University Agricultural Center,Aquaculture Research Station,2410Ben Hur Road,Baton Rouge,LA 70808,USAAbstractStandardized evaluation and rating of biofilters for aquaculture should be assessed in the context of the economic efficiency of ecological services (waste assimilation,nutrient recycling,and internal food production)provided by earthen ponds,and the availability and cost of land,water,and electrical energy resources required to support particular classes of production systems.In photosynthetic suspended-growth systems,water quality control is achieved by a combination of natural and mechanical processes.Natural processes include photosynthesis of oxygen,algal nutrient uptake,coupled nitrification–denitrification,and organic matter oxidation;mechanical processes include aeration and water circulation.Ammonia is controlled by a combination of phytoplankton uptake,nitrification,and immobilization by bacteria.Unlike biofilters for recirculating aquaculture systems,unit processes are combined and are an integral part of the culture unit.The important design and operational considerations for photosynthetic suspended-growth systems include temperature effects,aeration and mixing,quantity and quality of loaded organic matter,and fish water quality tolerance limits.The principle advantages of photosynthetic suspended-growth systems are lower capital costs relative to other recirculating aquaculture systems and increased control over stock management relative to conventional static ponds.The main disadvantage is the relatively low degree of control over water quality and phytoplankton density,metabolism,and community composition relative to other recirculating aquaculture systems.Examples of photo-synthetic suspended-growth systems include semi-intensive ponds,intensively aerated outdoor lined ponds,combined intensive–extensive ponds,partitioned aquaculture systems,greenwater tanks,greenwater tanks with solids removal,and greenwater recirculating aquaculture systems.#2005Elsevier B.V .All rights reserved.Keywords:Nitrification;Microbial floc;Pond aquaculture1.IntroductionBiofilters in recirculating aquaculture systems function to oxidize ammonia,1nitrite,and dissolvedorganic carbon (BOD).Additionally,one class of biofilters known as bioclarifiers capture particulate solids.All aquaculture production systems must function to prevent accumulation of these substances to levels that inhibit growth or immunocompetence of cultured organisms.Although dissolved oxygen concentration is often considered to be the primary water quality factor limiting production intensifica-tion,it is relatively easily controlled with straightfor-ward aeration or oxygenation technologies./locate/aqua-onlineAquacultural Engineering 34(2006)344–363*Tel.:+12257652848;fax:+12257652877.E-mail address:jhargreaves@.1The convention adopted in this paper will be to use the term ‘ammonia’to refer to un-ionized ammonia plus ionized ammonium (NH 3+NH 4+)expressed as nitrogen by mass unless otherwise explicitly indicated.0144-8609/$–see front matter #2005Elsevier B.V .All rights reserved.doi:10.1016/j.aquaeng.2005.08.009The general strategies for control of metabolite concentrations are removal by advection or treatment through chemical or biological transformation either within rearing units or in reactors hydraulically linked to rearing units.In open systems(e.g.,flow-through systems and net-pens),metabolites are removed by advection.In ponds and recirculating aquaculture systems well-known biogeochemical processes related to plant photosynthesis or the metabolism of heterotrophic and chemautotrophic bacteria are responsible for chemical transformations.Suspended-growth systems for aquaculture can be defined as those for which the maintenance of water quality depends on an active mass of phytoplankton, free and attached bacteria,aggregates of living and dead particulate organic matter,and microbial grazers that is maintained in suspension.This definition encompasses a broad range of aquaculture production systems,including conventional ponds,where pro-cesses related to the metabolism of phytoplankton have the primary influence on water quality,to various intensive pond,tank,and raceway systems,where a combination of autotrophic and heterotrophic micro-bial processes regulate water quality.These systems have been described with a wide variety of terms (Table1),most emphasizing the role of bacterial processes and giving the impression that processes related to phytoplankton metabolism are not impor-tant.It is more accurate to describe these systems as ‘‘mixed’’suspended-growth systems,insofar as water quality is maintained by a mixture of photosynthetic and bacterial processes.However,in most systems,the role of phytoplankton is paramount and so use of the term‘‘photosynthetic suspended-growth’’(PSG) system is appropriate and will be used here.In suspended-growth systems,substrates(e.g., dissolved organic carbon,ammonia,nitrite)are typically mixed with suspended microbes in rearing units.In attached-growth systems,substrates are transported from rearing units to specialized reactors performing a specific unit operation in a treatment train.In suspended-growth systems,the major ecological services associated with ponds—waste assimilation,nutrient recycling,and internal food production—are generally integral,as opposed to discrete unit processes that characterize systems with attached-growth biofilters.Earthen ponds are used to produce the over-whelming majority offinfish and crustacean produc-tion in aquaculture.These ponds use a suspension of phytoplankton and bacteria to maintain water quality. Thus,most aquaculture production occurs in PSG systems.There are very few examples of aquaculture industries based on attached-growth biofiltration technology.The eel culture industry in the Nether-lands is one notable example.In the1980s,many commercial shrimp ponds were operated with daily water exchange rates around10–15%to maintain water quality.Research at the Waddell Mariculture Center(Sandifer et al.,1991; Hopkins et al.,1993,1995a,b;Sandifer and Hopkins, 1996)and elsewhere indicated that water exchange was ineffective as a means to control water quality. Additionally,water exchange was identified as an important factor contributing to several disease epizootics in shrimp growing areas.As a result, shrimp producers began to reduce water exchange and to become more conservative and intensive in water use.Concerns related to biosecurity,the environ-mental impact of pond effluents,and the ineffective-J.A.Hargreaves/Aquacultural Engineering34(2006)344–363345 Table1Descriptive terms used for photosynthetic suspended-growth aquatic systemsTerm ReferenceAlgae,bacteria,zooplankton,and detritus(ALBAZOD)Soeder(1984)Organic detrital algae soup(ODAS)Serfling(2000)Zero water exchange/recycle strategy McIntosh(2000b)Zero exchange,aerobic,and heterotrophic(ZEAH)culture systems McNeil(2000)Aerated microbial reuse systems Chamberlain et al.(2001) High-rate bacterial systems McIntosh(2001)Active suspension ponds Milstein et al.(2001) Activated sludge ponds Boyd and Clay(2002) Suspended,bacterial-based treatment process Rakocy et al.(2004)ness of water exchange as a water quality management practice contributed to interest in managing ponds with long hydraulic residence times and the develop-ment of so-called‘‘zero-exchange’’systems. Although most recirculating aquaculture systems operate with very low water exchange rates(<5% dayÀ1)are thus conservative of water,many PSG systems are also intensive and conservative of water (Table2).This paper will discuss attributes of PSG systems for biofiltration in comparison tofixed-film biofiltra-tion.The paper will review ammonia dynamics in PSG systems,design,and operational considerations,the advantages and disadvantages of such systems,and conclude with descriptions of commercial and experimental applications of suspended-growth tech-nology in aquaculture.The focus of the discussion will be on ammonia and dissolved organic carbon dynamics,the two major water quality variables that are transformed in biofilters.2.Ammonia dynamics in photosynthetic suspended-growth systemsWater quality in PSG systems is controlled by a combination of natural and mechanical processes.The ultimate source of the nearly all the nitrogen in PSG systems is derived from feed loaded from outside the system boundary.Within the system,ammonia is produced by excretion by cultured organisms and mineralization of organic matter.Principal natural process involved in ammonia transformations include phytoplankton uptake,immobilization by hetero-trophic bacteria,and nitrification by chemautotrophic bacteria.Each process assumes variable importance depending on system configuration and management intensity.Although phytoplankton and bacteria play prominent roles,the activities of grazers and other components of the microbial community(e.g.,fungi and actinobacteria)interact with phytoplankton and bacteria as consumers or sources of substrates. Grazers range in size from small heterotrophic microflagellates to microzooplankton(e.g.,ciliated andflagellated protozoans)to zooplankton.Mechanical processes influence ammonia dynamics by creating the conditions that affect the rate of ammonia transformations.These include processes that provide oxygen(e.g.,aeration and oxygenation devices) and those that minimize spatial concentration gradients of substrates and maintaining particles in suspension by generating water movement(e.g.,aerators,low-speed paddlewheels,destratifiers,and water blenders).2.1.Phytoplankton uptakeIn most system configurations,processes related to photosynthesis and respiration of phytoplankton dominate oxygen dynamics and nutrient cycling. Phytoplankton uptake of ammonia varies seasonally because the rate is affected by temperature,light,and substrate concentration.Seasonal variation in perfor-mance must be considered in the design process.In most PSG settings,light limits photosynthesis because of self-shading associated with high phytoplankton densities(Smith,1988).J.A.Hargreaves/Aquacultural Engineering34(2006)344–363346Table2Estimated water requirements for the culture offish or shrimp in various photosynthetic suspended-growth systemsSystem Species Water use(m3kgÀ1)ReferenceIntensive greenhouse tank Tilapia0.03–0.04Serfling(2000) Intensive greenwater tank Tilapia0.08–0.12Rakocy et al.(2004) Greenwater recirculating raceway Shrimp0.37–0.48Moss et al.(2001) Intensive greenwater tank Tilapia0.8Milstein et al.(2001)Hybrid striped bass 1.3Combined intensive–extensive Tilapia0.85Mires et al.(1990) Semi-intensive pond Channel catfish 1.25–1.75Yoo and Boyd(1994) Zero-exchange,lined pond Penaeid shrimp1–2.26Boyd and Clay(2002) Conventional(w/water exchange)Penaeid shrimp6–64Hopkins et al.(1993)11–55Phillips et al.(1991)40–80Boyd and Clay(2002)Phytoplankton uptake of ammonia can be esti-mated on the basis of measured carbonfixation rates and Redfield stoichiometric ratios(atomic C:N:P=106:16:1).Conventional semi-intensive pondsfix between1and3g C mÀ2dayÀ1(Tucker, 1996),equivalent to ammonia uptake of176–528mg N mÀ2dayÀ1.In a model of ammonia dynamics in semi-intensive catfish ponds,phytoplank-ton uptake of ammonia during mid-summer of about 450mg N mÀ2dayÀ1was estimated(Hargreaves, 1997).Ammonia uptake rates using15N-labeled ammonia in intensively mixed lined shrimp ponds ranged from about600–1500mg N mÀ2dayÀ1(Bur-ford et al.,2003).Phytoplankton in some PSG systems canfix from10to12g C mÀ2dayÀ1(Brune et al., 2003),equivalent to ammonia uptake rates of1760–2113mg N mÀ2dayÀ1.Although phytoplankton uptake represents a powerful mechanism of ammonia removal from culture water,it represents only a temporary packa-ging of nitrogen in the form of cellular protein. Eventually,phytoplankton cells die and the organic nitrogen(and other nutrients)contained in the cells are mineralized by heterotrophic bacteria and nitrogen is recycled to the water as ammonia.2.2.Heterotrophic immobilizationAs system intensity increases,the importance of bacteria in nutrient cycling increases,although phytoplankton processes remain dominant.For exam-ple,as organic loading to a high-rate algae pond increased,the proportion of nitrogen contained in a ‘‘filamentousfloc’’fraction also increased(Cromar et al.,1992).Densities of heterotrophic bacteria in sediment(from107to109gÀ1)are typically two to three orders of magnitude greater than in the water column(Ram et al.,1982).However,bacteria densities,particularly,attached bacteria density,in the water column can be increased by maintaining organic particles in suspension.For example,total bacteria densities from about3Â107to5Â107mlÀ1 have been measured in the water column of intensively mixed lined shrimp ponds(Burford et al.,2003). Suspended detritus,most derived from settled phytoplankton,serves as a substrate for the coloniza-tion of bacteria.About40%of bacteria in the water column of intensively mixed lined shrimp ponds were associated with suspended particles(Burford et al., 2003).In many aquatic systems,the metabolic activity of attached bacteria are much greater than those of free bacteria and the proportion of attached bacteria increases with particle concentration(Kirchman and Ducklow,1987).The requirements of heterotrophic bacteria for carbon and nitrogen can be described by a stoichio-metric relationship between carbon and nitrogen of about5:1by atoms(Goldman et al.,1987).Bacteria C:N ratio varies little(between5:1and6:1)despite a wide variation in substrate C:N ratio.The regeneration or immobilization of nutrients depends on the total carbon taken up by bacteria,the gross conversion efficiency by bacteria(i.e.,the proportion of total carbon used for biosynthesis),the C:N ratio of bacteria,and the C:N ratio of substrate(Fenchel and Blackburn,1979).In empirical studies using substrates with a range of C:N ratio,regeneration of ammonia increased as C:N ratio decreased(Goldman et al.,1987).In a model of nitrogen dynamics in intensive shrimp ponds using the data of Goldman et al.(1987),regeneration rate was described as a function of C:N ratio in a polynomial equation: regeneration rate(%)=88.0À7.0(C:N ratio)+0.12 (C:N ratio)2(Montoya et al.,2002).In empirical studies,almost no(0.1%)ammonia was regenerated at a substrate C:N ratio of10:1,but70%of ammonia was regenerated at a substrate C:N of1.5:1.Thus,nitrogen limitation of bacteria growth can be promoted by providing substrates with a C:N ratio of10:1or greater.Most particulate detritus in eutrophic aquaculture ponds is derived from settled phytoplankton.The C:N ratio of phytoplankton is greater than that of bacteria and can be described by Redfield stoichiometry (atomic C:N:P ratio of106:16:1).The C:N ratio of phytoplankton is more variable than that of bacteria and ranges from about6:1to7:1.Based on the C:N ratio of phytoplankton and the carbon and nitrogen requirements of bacteria,about35%of the nitrogen in detritus derived from phytoplankton will be regener-ated as ammonia.Bacterial immobilization of ammonia in PSG systems can be promoted by manipulating input C:N ratio,a factor that has been emphasized as a driving force in PSG system performance by numerous investigators(Avnimelech et al.,1989;J.A.Hargreaves/Aquacultural Engineering34(2006)344–363347Avnimelech et al.,1994;Avnimelech,1999;Burford et al.,2004;Hari et al.,2004).The accumulation of inorganic nitrogen in PSG systems can be minimized by two approaches:(1)the addition of organic carbon sources with a wide C:N ratio and(2)reduction of feed protein content.For example,a32%crude protein feed hasaC:Nratioofabout10.6,buta22%proteinfeed hasa C:aniccarbon sourcesaddedto promote immobilization of inorganic nitrogen include glucose,cassava meal(Avnimelech andMokady,1988), cellulose,sorghum meal(Avnimelech et al.,1989),corn flour(Milstein et al.,2001),molasses(Burford et al., 2004),and tapiocaflour(Hari et al.,2004).Immobilization of ammonia by heterotrophic bacteria occurs much more rapidly than nitrification because the growth rate and microbial biomass yield per unit substrate of heterotrophic bacteria is much greater than that of nitrifying bacteria(Metcalf& Eddy Inc.and Tchobanoglous,1991).Therefore, promoting nitrogen limitation of bacteria growth by adding carbon may reduce ammonia concentration more rapidly than nitrification.One of the negative aspects of this approach to ammonia control is that the addition of organic matter as an ammonia control method increases the demand for oxygen.However, oxidation of the additional BOD can be managed through relatively simple and straightforward aeration technology.2.3.NitrificationPhilosophically,maximizing nitrogen recycling and retention within the system or maximizing nitrogen removal are opposing management goals of aquaculture production systems.Nitrogen retention within a production system can be maximized by promoting processes that encourage nitrogen recy-cling and increasing nitrogen recovery infish,such as the harvest of bacterial or algal protein by culture species able to access these potential food resources. Alternatively,nitrogen accumulation during culture can be limited by promoting processes that enhance nitrogen removal from the system.Coupled nitrifica-tion–denitrification is the major nitrogen removal process in PSG systems in either case.Ultimate removal of inorganic nitrogen from a PSG system during culture can only occur through coupled nitrification–denitrification.(V olatilization of ammo-nia from the water surface to the atmosphere is minor in relation to the potential for nitrification–denitrifica-tion as an inorganic nitrogen removal mechanism.) Ammonia mustfirst be oxidized to nitrate,which diffuses to anaerobic microzones for denitrification. Although denitrification proceeds only in the absence of oxygen,the overall process rate is ultimately oxygen-dependent and so nitrification is the rate-limiting step(Rysgaard et al.,1994).Thus,nitrogen can cycle repeatedly between inorganic and organic compartments(e.g.,fish,phytoplankton,bacteria biomass),before it is ultimately lost from the system as nitrogen gas.The factors that control nitrification rate in PSG systems are similar to those that control this reaction in most biofilters.Maintenance of dissolved oxygen concentration above a critical threshold(>2mg lÀ1)is necessary for good nitrification.In biofilms in laboratory-scale reactors,nitrification shifted from ammonia limitation to oxygen limitation when the ratio between dissolved oxygen and ammonia was less than 1.5–2.0mg O2mg TANÀ1(Nogueira et al., 1998).In suspended-growth systems,oxygen diffu-sion limitations on nitrification may not be as severe as infixed-film systems because the biofloc is suspended in a turbulent environment with dissolved oxygen concentrations above the critical threshold forfish production,which is usually greater than the threshold for nitrification.As in biofilters in recirculating aquaculture systems,nitrification is often limited by the practical consideration of maintaining low ammonia concen-tration to promote goodfish or shrimp growth.Thus, substrate concentrations often affect reaction kinetics and nitrification is likely zero-or half-order with respect to concentration(Harremo¨es,1978).Nitrifica-tion in PSG systems is also controlled by attachment sites.In conventional ponds,the thin layer at the sediment–water interface is the location for the most active nitrification.Suspension of organic and mineral particles increases nitrification by increasing the number and surface area for nitrifying bacteria attachment.Nitrification in some PSG systems may be affected by high rates of organic matter loading and long hydraulic residence times.These conditions promote the accumulation of particulate and dissolved organic matter,which inhibit nitrification infixed-filmJ.A.Hargreaves/Aquacultural Engineering34(2006)344–363 348biofilters(Bovendeur et al.,1990;Figueroa and Silverstein,1992;Ling and Chen,2005).Infixed-film biofilters,the greater growth rate of heterotrophic bacteria compared to nitrifying bacteria may result in inhibition of nitrification by competition for attach-ment sites and physical overgrowth of nitrifying bacteria.In wastewater treatment systems,the proportion of total microbial biomass that is nitrifying bacteria is inversely related to the BOD5:TKN ratio in the wastewater(Metcalf&Eddy Inc.and Tchoba-noglous,1991).At a BOD5:TKN ratio of0.5,the nitrifier fraction is35%,but reduces to<1%at a BOD5:TKN ratio of3.In suspended-growth systems, the nitrification rate(R;g N g VSSÀ1dayÀ1)declines exponentially as a function of the influent COD:N ratio according to the following equation: R=0.0323+0.334e[À1.660(COD:N)](Carrera et al., 2004).At an influent COD:N ratio of3,the nitrification rate is0.035g N g VSSÀ1dayÀ1,but at a COD:N ratio of0.5,the nitrification rate is0.178g N g VSSÀ1dayÀ1.In an intensive greenwater tank system,average nitrification was 3.2mg lÀ1dayÀ1 (Rakocy et al.,2004).Assuming a TSS concentration of500mg lÀ1and a VSS:TSS of0.25,then the nitrification rate is about0.025g N g VSSÀ1dayÀ1, suggesting that nitrification is restricted by accumu-lated dissolved organic carbon in this PSG system. Thus,nitrification in PSG systems may be inhibited below potential by dissolved organic carbon,despite otherwise favorable conditions with respect to mixing, oxygen concentration,and the presence of abundant attachment sites on suspended particles.Nitrification rates in suspended-growth and attached-growth systems are not directly comparable. In attached-growth systems,biofilters have media with a very high specific surface area in comparison with most suspended-growth systems.Furthermore,nitri-fication rates in attached-growth biofilters at substrate-limiting ammonia concentrations can be enhanced by increasing mass loading to the biofilter by reducing hydraulic retention time in the reactor.In volumetric terms,nitrification in attached-growth biofilters is orders of magnitude greater than in suspended-growth systems.However,there are two alternative approaches to achieve the same overall effect for a given loading rate:(1)a low volume but high specific surface area biofilter in attached-growth systems or(2) a relatively large volume combined rearing unit-nitrification reactor with a relatively low volumetric nitrification rate in suspended-growth systems.Nitrification rates in PSG systems are variable.In intensively aerated,lined shrimp ponds,nitrification ranged from0.1to0.8mg N lÀ1dayÀ1(Burford et al., 2003).In a greenwater recirculating tank system,in situ nitrification was about1mg N lÀ1dayÀ1(Mia, 1996).In an intensive and well-mixed greenwater tank system with solids removal,nitrification of 3.2mg N lÀ1dayÀ1was estimated from nitrate accu-mulation rate(Rakocy et al.,2004).Assuming an ammonia conversion infixed-film biofilters of200mg N mÀ2dayÀ1and a volumetric nitrification rate in PSG systems of1mg N lÀ1dayÀ1,then the nitrification rate of5m2of biofilter media surface area would be equivalent to that in1m3of water in a PSG system.2.4.Algal–bacterial interactionsInteractions between bacteria and algae contribute to the complexity of water quality dynamics in PSG systems.The relative importance of each of these groups varies,depending upon the quality and rate of organic matter loading,season(i.e.,temperature),and intensity and extent of turbulent mixing.In an intensively loaded high-rate algal pond,three fractions were separated that corresponded to algae only, bacteria only,andfloccular aggregates of algae and bacteria(Cromar and Fallowfield,2003).Algae and bacteria have a range of stimulatory or inhibitory effects on each other(Cole,1982).In general, bacteria productivity increases as a function of primary productivity.Senescent algae or algal detritus has a greater stimulatory effect on bacteria growth than living cells.However,living cells release large quantities of dissolved organic carbon in the form of simple to complex polysaccharides that are readily utilized by heterotrophic bacteria.Bacteria can stimulate phyto-plankton growth by regenerating nutrients,producing vitamins and other stimulatory products.Algae and bacteria also have inhibitory effects on each other.Both groups produce substances that are antagonistic to growth.Algae can produce antibiotics and both groups can produce allelopathic substances (Cole,1982).Some bacteria can lyse phytoplankton cells.Each group can modify the chemical micro-environment of the other group by metabolic activity, thereby affecting nutrient availability or metabolism.J.A.Hargreaves/Aquacultural Engineering34(2006)344–363349Phytoplankton compete with bacteria for sub-strates.Notably,phytoplankton will compete with chemoautotrophic nitrifying bacteria for ammonia. Phytoplankton may be more effective competitors for ammonia than nitrifying bacteria when concentrations are low during the summer,but the reverse may be the case when concentrations are higher during the winter (Hargreaves,1997).There are clear trade-offs associated with dom-inance by phytoplankton or bacteria.Phytoplankton are important oxygenators and provide a powerful mechanism for ammonia removal,but water quality is inherently more dynamic andfluctuations more extreme than in bacteria-dominated ponds.However, ponds dominated by bacteria require an external oxygen source that must be provided by a level of mechanical aeration that far exceeds the natural supply of oxygen in algae-dominated ponds.Although phytoplankton densities are lower in ponds dominated by non-algal(mineral)turbidity,cyanobacteria, including those that may cause off-flavor,are favored in the dimly lit waters of such ponds.The stability of microbial suspensions has not been well-documented, but some threshold cell density,as indexed by TSS concentration,may be necessary to maintain good water quality in bacteria-dominated ponds.3.The role of sediment in photosynthetic suspended-growth systemsMost aquaculture ponds are constructed in soil and are shallow,suggesting that sediment–water interac-tions are important regulators of water quality.The research conducted with reduced water exchange shrimp ponds at Waddell Mariculture Center used lined ponds.This research stimulated interest in the evaluation of lined ponds in a commercial setting. Therefore,sediment may or may not be included as a system component and there are advantages and disadvantages of including sediment as a component of PSG systems.In simple,low-intensity PSG systems,the sediment is the location for the most active microbial transformations involving organic matter and nutrient cycling.The sediment–water interface is character-ized by very steep chemical and physical gradients. Concentrations of dissolved oxygen,ammonia,and reactive phosphate and pH and oxidation–reduction potential can vary widely over several centimeters. The sediment serves as a sink for dissolved oxygen and nitrate and a source of ammonia and reactive phosphate,particularly under anaerobic conditions. Additionally,sediment can be a source of sulfide or other potential toxins in brackish or seawater systems. Thus,reactions occurring across the sediment–water interface can cause growth retardation of cultured organisms(Avnimelech and Zohar,1986),particularly for species with a benthic orientation(e.g.,shrimp). Furthermore,typical aquaculture pond sediments are poorly consolidated,thereby presenting a difficult working surface for harvesting.One of the greatest benefits of including sediment as a component of PSG systems is the enormous denitrification potential associated with the large volume of anaerobic sediment.Nitrate produced by active nitrification in the thin oxidized layer at the sediment–water interface diffuses into the anaerobic sediment,where it is denitrified.Nitrogen gas is then lost from the pond by gas ebulition.In PSG systems, where the culture unit is lined with an impermeable barrier and where well-mixed and well-oxygenated conditions promote nitrification,nitrate may accumu-late because anaerobic conditions are comparatively sparse.Photosynthetic suspended-growth systems that include sediment as a component have enormous phosphorus sorption capacity,especially for those constructed in clay soils.Thus,coupled nitrification–denitrification and phosphorus sorption in sediment are responsible for substantial nutrient removal in PSGs that include sediment.In lined systems,nitrogen and phosphorus may accumulate,thereby increasing the pollution potential of any effluent,but also increasing the opportunities for nutrient recovery. 4.Considerations for design and operation of photosynthetic suspended-growth systems4.1.Temperature effectsThe biofiltration capacity of PSG systems is profoundly affected by water temperature.Tempera-ture affects the kinetics of chemical reactions,nutrient uptake by phytoplankton,and microbial growth.J.A.Hargreaves/Aquacultural Engineering34(2006)344–363 350。

托福讲座单词

托福讲座单词

chaotic utilize canoe constellation zenith star latitude
混乱的 利用 独木舟 星座 天顶星 纬度
evidence breakthrough existence geomagnetic magnetic field fluctuate
indistinguishable 难以区别的 orbit 轨道
TPO30-L3 microorganism invention exploration jarosite yellowish mineral crystalline potassium hydroxide acidic planetary essential prerequisite biological incorporate detect remnant micro fabricated analyzer pure eliminate contaminate amino acide
引人注目的 准确指出 撞击 大坑 航天器 轨道 地形的 海拔 直径 太阳系 行星的 重新分配 地壳 地幔 渗透 构成 形成 盆地 精确的 流星雨 粗略的 估计 调查 自我维持的 反射的
asteroid belt a decimal point tantalizingly skip sense predicting gravitation mathematically attempt pattern telescope faint
scatter geological crucial dynamic lubricate sulfur dioxide fluctuation eruption burst intriguing Venus
分散 地质的 重要的 动力的 润滑 二氧化硫 起伏 喷发 爆发 有趣的 金星

NRL_FORMULARY_11

NRL_FORMULARY_11

3
Euler-Mascheroni constant1 γ = 0.57722 Gamma Function Γ(x + 1) = xΓ(x): Γ(1/6) Γ(1/5) Γ(1/4) Γ(1/3) Γ(2/5) Γ(1/2) = = = = = = 5.5663 4.5908 3.6256 2.6789 2.2182 √ 1.7725 = π Γ(3/5) Γ(2/3) Γ(3/4) Γ(4/5) Γ(5/6) Γ(1) = = = = = = 1.4892 1.3541 1.2254 1.1642 1.1288 1.0
C
dl · A dl × A f dg = − gdf
C
(35)
S
(d S × ∇) × A = d S · (∇f × ∇g ) =
C
(36)
S
C
DIFFERENTIAL OPERATORS IN CURVILINEAR COORDINATES5 Cylindrical Coordinates Divergence ∇·A= Gradient (∇f )r = Curl (∇ × A)r = ∂Aφ 1 ∂Az − r ∂φ ∂z ∂Ar ∂Az − ∂z ∂r 1 ∂ 1 ∂Ar (rAφ ) − r ∂r r ∂φ ∂f ; ∂r (∇f )φ = 1 ∂f ; r ∂φ (∇f )z = ∂f ∂z 1 ∂Aφ ∂Az 1 ∂ (rAr ) + + r ∂r r ∂φ ∂z
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . 26 . . . . . . . . . . . . . . . . . . . . 30 . . . . . . . . . . . . . . . . . . . . . 31 . . . . . . . . . . 40 . . . . . . . . . . . . . . . . . . . . . . 42 . . . . . . . . . . . . . . . . . . . . . 43 . . . . . . . . . . . . . . . . . . . . . . 44 . . . . . . . . . . . . . . . . . . . . 46

青藏高原及邻区大地构造单元初步划分

青藏高原及邻区大地构造单元初步划分

记录了晚古生代 — 中生代弧后扩张、多岛弧盆系发 育、 弧# 弧碰撞、 弧#陆碰撞的地质演化历史。 碰撞之 后该区的大部分地区于晚三叠世转化为陆地,并形 成碰撞后地壳伸展背景下的裂陷或裂谷盆地。 冈瓦纳北缘晚古生代 — 中生代冈底斯 # 喜马拉 雅构造区:班公湖 # 丁青 # 碧土 # 昌宁 # 孟连对接带 是冈瓦纳与劳亚 # 泛华夏大陆的分界线,亦即冈瓦 纳大陆的北界。伯舒拉岭 # 高黎贡山属于冈瓦纳晚 古生代 — 中生代前锋弧, 聂荣隆起、 嘉玉桥变质地体 等是前锋弧的残块。在前锋弧的后面 (南侧) 是晚古 生代 — 中生代冈底斯 # 喜马拉雅弧后扩张、多岛弧 盆系发育、 弧 #弧碰撞、 弧#陆碰撞的演化结果。 该区 三叠纪和侏罗纪 — 早白垩世的雅鲁藏布江蛇绿岩是 目前青藏高原乃至中国大陆内, 保存最好、 最完整的 蛇绿岩 “三位一体” 组合, 代表了特提斯洋向南俯冲 诱导出的一系列藕断丝连的弧后扩张盆地。
(或晚古 !#*&#奥依塔格 &库地 &苏巴什 &东昆中结合带 生代裂谷带) (含前寒武纪岩块、 海山和洋岛残块) !#*&!南昆仑残余弧
#+&!高喜马拉雅结晶岩带或基底逆冲带 #+&)低喜马拉雅褶冲带 #+&+锡伐利克后造山前陆盆地带
!#*&)麻扎&康西瓦&木孜塔格&西大滩晚古生代结合带 ! 泛华夏大陆晚古生代羌塘 &三江构造区 "#扬子陆块:
青藏高原具有复杂而独特的巨厚地壳和岩石圈 结构, 是一个由北部劳亚大陆、 泛华夏陆块西缘和南 部冈瓦纳大陆北缘不断弧后扩张、裂离,又互相对 接、 镶嵌构成的复杂地区, 经历了漫长的构造变动历 史, 特别是古生代以来的多岛弧盆系的形成演化, 最 终由!"多条规模不等的弧#弧、弧#陆碰撞结合带和 其间的岛弧或陆块拼贴而成。由于后期印度板块向 北强烈顶撞,在其左右犄角处分别形成帕米尔和横 断山构造结及相应的弧形弯折,在东西两端改变了 原来东西向展布的构造面貌,加之华北和扬子刚性 陆块的阻抗和陆内俯冲对原有构造,特别是深部地 幔构造的改造, 造成了青藏高原独特的构造、 地貌景 观, 形成了统一的深部幔拗和地表的隆升。

Computer simulations of interferometric imaging with the VLT interferometer and the AMBER i

Computer simulations of interferometric imaging with the VLT interferometer and the AMBER i

a r X i v :a s t r o -p h /0003474v 1 31 M a r 2000Computer simulations of interferometric imaging with the VLTinterferometer and the AMBER instrumentT.Bl¨o cker,K.-H.Hofmann,F.Przygodda and G.WeigeltMax-Planck-Institut f¨u r Radioastronomie,Auf dem H¨u gel 69,53340Bonn,GermanyABSTRACTWe present computer simulations of interferometric imaging with the VLT interferometer and the AMBER instru-ment.These simulations include both the astrophysical modelling of a stellar object by radiative transfer calculationsand the simulation of light propagation from the object to the detector (through atmosphere,telescopes,and the AMBER instrument),simulation of photon noise and detector read-out noise,and finally data processing of the interferograms.The results show the dependence of the visibility error bars on the following observational pa-rameters:different seeing during the observation of object and reference star (Fried parameters r 0,object =2.4m,r 0,ref .=2.5m),different residual tip-tilt error (δtt ,object =2%of the Airy disk diameter,δtt ,ref .=0.1%),and object brightness (K object =3.5mag and 11mag,K ref .=3.5mag).Exemplarily,we focus on stars in late stages of stellar evolution and study one of its key objects,the dusty supergiant IRC +10420that is rapidly evolving on human timescales.We show computer simulations of VLTI interferometry of IRC +10420with two ATs (wide-field mode,i.e.without fiber optics spatial filters)and discuss whether the visibility accuracy is sufficient to distinguish between different theoretical model predictions.Keywords:VLTI,AMBER,interferometric imaging,IRC +10420,dust shells,circumstellar matter,supergiants,stellar evolution1.INTRODUCTIONThe Very Large Telescope Interferometer 1(VLTI)with its four 8.2m unit telescopes (UTs)and three 1.8m auxiliary telecopes (ATs)will certainly establish a new era of optical and infrared interferometric imaging within the next few years.With a maximum baseline of up to more than 200m,the VLTI will allow the study of astrophysical key objects with unprecendented resolution opening up new vistas to a better understanding of their physics.The near-infrared focal plane instrument of the VLTI,the Astronomical MultiBEam Recombiner 2,3(AMBER),will operate between 1and 2.5µm and for up to three beams allowing the measurement of closure phases.In a second phase its wavelength coverage is planned to be extended to 0.6µm.Objects as faint as K =20mag are expected to be observable with AMBER when a bright reference star is available,and as faint as K =14mag otherwise.Among the astrophysical key issues 4are,for instance,young stellar objects,active galactic nuclei and stars in late phases of stellar evolution.The simulation of a stellar object consists in principle of two components:(i)the calculation of an astrophysical model of the object,typically based on radiative transfer calculationspredicting,e.g.,its intensity distribution.To obtain a robust and non-ambiguous model,it is of particular importance to take diverse observational constraints into account,for instance the spectral energy distribution and visibilities.(ii)the determination of the interferometer’s response to this intensity signal,i.e.the simulation of light propagationin the atmosphere and the interferometer.Often,only one of the above parts is considered in full detail.The aim of this study is to combine both efforts and to present a computer simulation of the VLTI performance for observations of one object class,the dusty supergiants.For this purpose,we calculated a detailed radiative transfer model for one of its most outstanding representatives,the supergiant IRC +10420,and carried out computer simulations of VLTI visibility measurements.The goal is to estimate how accurate visibilities can be measured with the VLTI in this particular but not untypical case,to discuss if the accuracy is sufficient to distinguish between different theoretical model predictions and to study on what the accuracy is dependent.2.THE SUPERGIANT IRC+10420:EVOLUTION ON HUMAN TIMESCALES The star IRC+10420(=V1302Aql=IRAS19244+1115)is an outstanding object for the study of stellar evolution. Its spectral type changed from F8I+a in19735to mid-A today6,7corresponding to an effective temperature increase of1000-2000K within only25yr.It is one of the brightest IRAS objects due to its very strong infrared excess by circumstellar dust and one of the warmest stellar OH maser sources known.8,9,10Large mass-loss rates,typically of the order of several10−4M⊙/yr were determined by CO observations.11,6IRC+10420is believed to be a massive luminous star currently being observed in its rapid transition from the red supergiant stage to the Wolf-Rayet phase.9,10,12,6,7Its massive nature(initially∼20to40M⊙)can be concluded from its distance(d=3–5kpc)and large wind velocity(40km/s)ruling out alternative scenarios.IRC+10420is the only object observed until now in its transition to the Wolf-Rayet phase.The structure of the circumstellar environment of IRC+10420appears to be very complex,13and scenarios proposed to explain the observed spectral features of IRC+10420include a rotating equatorial disk,12bipolar outflows,14and the infall of circumstellar material onto the star’s photosphere.15 Several infrared speckle and coronographic observations16,17,18,19,20were conducted to study the dust shell of IRC+10420.The most recent study21reports thefirst diffraction-limited73mas bispectrum speckle interferometry of IRC+10420and presents thefirst radiative transfer calculations that model both the spectral energy distribution and the visibility of this key object.In the following we will briefly describe the main results and conclusions of this study which will serve as astrophysical input for the VLTI computer simulations presented in the next section.2.1.Visibility measurements and bispectrum speckle interferometrySpeckle interferograms of IRC+10420were obtained with the Russian6m telescope at the Special Astrophysical Observatory on June13and14,1998.21The speckle data were recorded with our NICMOS-3speckle camera through an interferencefilter with a centre wavelength of2.11µm and a bandwidth of0.19µm.The observational parameters were as follows:exposure time/frame50ms;number of frames8400;2.11µm seeing(FWHM)∼1.′′0.A diffraction-limited image of IRC+10420with73mas resolution was reconstructed from the speckle interferograms using the bispectrum speckle interferometry method.22,23,24The modulus of the object Fourier transform(visibility) was determined with the speckle interferometry method.25Figure1(left panel)shows the reconstructed2.11µm0.00.20.40.60.81.0024681012NormalizedVisibilityCycles per arcsec-2.5-2.0-1.5-1.0-0.50.0050100150200Log(NormalizedIntensity)Image Space (mas)Figure1.Left:Azimuthally averaged2.11µm visibility of IRC+10420with error bars for selected frequencies. This visibility function consists of a constant plateau(visibility∼0.6)caused by the unresolved central object and a triangle-shaped low-frequency function caused by the faint extended nebula.Right:Azimuthally averaged radial plots of the reconstructed diffraction-limited2.11µm-images of IRC+10420(solid line)and HIP95447(dashed line). visibility function of IRC+10420.There is only marginal evidence for an elliptical visibility shape(position angle of the long axis∼130◦±20◦,axis ratio∼1.0to1.1).The visibility0.6at frequencies>4cycles/arcsec shows that the stellar contribution to the totalflux is∼60%and the dust shell contribution is∼40%.The Gaußfit FWHM diameter of the dust shell was determined to be d FWHM=(219±30)mas.By comparison,previous3.8m telescopeK-band observations 19found a dust-shell flux contribution of ∼50%and d FWHM =216mas.However,as will be shown later,a ring-like intensity distribution appears to be much better suited than the assumption of a Gaussian distribution whose corresponding FWHM diameter fit may give misleading sizes (see Sect.2.2).The right panel of Fig.1displays the azimuthally averaged diffraction-limited images of IRC +10420and the unresolved star HIP 95447.The K -band visibility will strongly constrain radiative transfer calculations as shown in the next section.2.2.Dust-shell modelsThe spectral energy distrubion (SED)of IRC +10420with 9.7and 18µm silicate emission features is shown in Fig.2.It corresponds to the ’1992’data set used by Oudmaijer et al.6)and combines VRI (October 1991),near-infrared (March and April 1992)and Kuiper Airborne Observatory (June 1991)photometry 12with the IRAS measurements and 1.3mm data.26Additionally,we included data for λ<0.55µm.27In contrast to the near-infrared,the optical magnitudes have remained constant during the last twenty years within a tolerance of ≈0.m 1.IRC +10420is highly reddened due to an extinction of A total V≈7m by the interstellar medium and the circum-stellar shell.From polarization studies an interstellar extinction of A V ≈6m to 7m was estimated.27,12Based on the strength of the diffuse interstellar bands,E (B −V )=1.m 4±0.m 5can be inferred for the interstellar contributioncompared to a total of E (B −V )=2.m 4.15We will use an interstellar A V of 5m adopting the method of Savage &Mathis 28with A V =3.1E (B −V ).6-16-15-14-13-12-11-10-90.11101001000l o g (λ*F λ) [W m -2]λ [µm]A=10 Y=4.5A=20 Y=4.5A=40 Y=4.5A=80 Y=4.500.20.40.60.81246810V i s i b i l i t y (λ=2.1µm )q [arcsec -1]A=10 Y=4.5A=20 Y=4.5A=40 Y=4.5A=80 Y=4.5Figure 2.SED (top)and 2.11µm visibility (bottom)for a superwind model with Y =r/r 1=4.5and different amplitudes.The inner shell obeys a 1/r 2density distribution,the outer shell a 1/r 1.7density distribution.Model parameters are:black body,T eff=7000K,T 1=1000K,τ0.55µm =7.0,silicates,32standard grain size distribution 31with a max =0.45µm,and Y out =104.The symbols refer to the observations (see text)corrected for interstellar extinction of A v =5m .0.20.40.60.81110100f l u x f r a c t i o nλ [µm]stellar radiation scattered radiationdust emission1e-081e-071e-061e-050.00010.0010.010.11-300-200-100010*******I n t e n s i t y (λ=2.11µm )ϑ [mas]Figure 3.Left:Fractional contributions of the emerging stellar radiation as well as of the scattered radiation and of the dust emission to the total flux as a function of the wavelength for a superwind model with Y =r/r 1=4.5,A =40and a 1/r 1.7density distribution in the outer shell.Model parameters are:black body,T eff=7000K,T 1=1000K,τ0.55µm =7.0,silicates,32standard grain size distribution 31with a max =0.45µm,and Y out =104.Right:Corresponding normalized intensity vs.angular displacement ϑ.The (unresolved)central peak belongs to the central star.The inner hot rim of the circumstellar shell has a radius of 35mas,and the cool component is located at 155mas.Both loci correspond to local intensity maxima.In order to model both the observed SED and 2.11µm visibility,we performed radiative transfer calculations for dust shells assuming spherical symmetry.We used the code DUSTY 29which solves the spherical radiative transfer problem utilizing the self-similarity and scaling behaviour of IR emission from radiatively heated dust.30To solve the radiative transfer problem including absorption,emission and scattering several properties of the central source and its surrounding envelope are required,viz.(i)the spectral shape of the central source’s radiation;(ii)the dust properties,i.e.the envelope’s chemical composition and grain size distribution as well as the dust temperature at the inner boundary;(iii)the relative thickness of the envelope,i.e.the ratio of outer to inner shell radius,and the density distribution;and (iv)the total optical depth at a given reference wavelength.The code has been expanded for the calculation of synthetic visibilities.We calculated various models of the dust shell of IRC +10420considering black bodies and model atmospheres as central sources of radiation,different silicates,grain-size and density distributions.We refer to Bl¨o cker et al.21for a full description of the model grid.The remaining fit parameters are the dust temperature,T 1,which determines the radius of the shell’s inner boundary,r 1,and the optical depth,τ,at a given reference wavelength,λref .We refer to λref =0.55µm.Models were calculated for dust temperatures between 400and 1000K and optical depths between 1and 12.Significantly larger values for τlead to silicate features in absorption.The near-infrared visibility strongly constrains dust shell models since it is,e.g.,a sensitive indicator of the grain size.Accordingly,high-resolution interferometry results provide essential ingredients for models of circumstellar dust-shells.In the instance of IRC +10420(assuming the central star to be a black body of T eff=7000K),the silicate 32grain sizes,a ,were found to be in accordance with a standard distribution function,31n (a )∼a −3.5,with a ranging between a min =0.005µm and a max =0.45µm.However,the observed dust shell properties cannot be fitted by single-shell models but seem to require multiple components.At a certain distance we considered an enhancement over the assumed 1/r x density distribution.The best model for both SED and visibility was found for a dust shell with a dust temperature of 1000K at its inner radius of r 1=69R ∗.At a distance of r =308R ∗(Y =r/r 1=4.5)the density was enhanced by a factor of A =40and and its density exponent was changed from x =2to x =1.7.These fits for SED and 2.11µm visibility are shown in Fig.2.The various flux contributions at 2.11µm are 62.2%stellar light,26.1%scattered radiation and 10.7%dust emission (see Fig.3),i.e.the radiation emitted by the circumstellar shell itself consists of 71%scattered radiation and 29%direct dust emission.The shell’s model intensity distribution is shown in Fig.3and was found to be ring-like.Table1.Parameters of simulationswavelength opt.throughput quantum efficiency photon number K=3.5mag15e−50msFigure4.Schematic view of interferometric imaging with the VLTI AMBER camera.Figure5.Flow chart of the simulation of VLTI AMBER interferograms(wide-field mode,i.e.withoutfiber optics spatialfiltering).star(Fried parameters r0,object=2.4m,r0,ref.=2.5m),different residual tip-tilt error(δtt,object=2%of the Airy disk diameter,δtt,ref.=0.1%),and object brightness(K object=3.5mag and11mag,K ref.=3.5mag).In the computer experiments(a)-(c),object and reference star were assumed to have the same brightness(K=3.5mag,see Table1), in experiment(d)fainter objects(K=11mag)were simulated as well.Fig.6shows the results of the simulations (a)to(c)based on the A=40intensity profile(see Fig.3)for baselines of50and100m together with the model predictions for IRC+10420.To obtain error bars each simulation was repeated three times.The insetted table lists the parameters of the simulations.Simulation(a)based on2400interferograms represents the ideal case of excellent seeing conditions(r0,object=2.5m, r0,ref.=2.5m)and an almost perfect tip-tilt correction(residual tip-tilt errorδtt=0.1%).The visibility error∆V amounts to±0.0029.Simulation(b)illustrates the influence of a larger residual object tip-tilt error(δtt=2%)leading to∆V±0.0036.Simulation(c)shows the impact of different seeing conditions for object and reference star.In the last simulation(d)of this series(shown only in the insetted table of Fig.6)an IRC+10420-like intensity distribution is assumed but the simulated K-magnitude is11mag(instead of3.5mag).Although the photon number decreases to N=4.41·103(i.e.by a factor of1000,see Table1),the visibility accuracy is still high,i.e.∆V±0.0036.4.CONCLUSIONSWe have presented computer simulations of interferometric imaging with the VLT interferometer and the AMBER instrument in the wide-field mode.These simulations include both the astrophysical modelling of a stellar object by radiative transfer calculations and the simulation of light propagation from the object to the detector and sim-ulation of photon noise and detector read-out noise.We focussed on stars in late stages of stellar evolution and examplarily studied one of its most outstanding representatives,the dusty supergiant IRC+10420.The model in-tensity distribution of this key object,obtained by radiative transfer calculations,served as astrophyiscal input for the VLTI/AMBER simulations.The results of these simulations show the dependence of the visibility error bar on various observational parame-ters.With these simulations at hand one can immediately see under which conditions the visibility data quality would allow us to discriminate between different model assumptions(e.g.the size of the superwind amplitude).Inspection of Fig.6shows that in all studied cases the observations will give clear preference to one particular model.Therefore, observations with VLTI will certainly be well suited to gain deeper insight into the physics of dusty supergiants.ACKNOWLEDGMENTSThe bispectrum speckle observations were made with the SAO6m telescope operated by the Special Astrophysical Observatory,Russia.The radiative-transfer calculations are based on the code DUSTY developed byˇZ.Ivezi´c,M. Nenkova and M.Elitzur.REFERENCES1.Glindemann,A.,et al.,2000,Proc.SPIE Conf.4006-012.Petrov,R.,Malbet,F.,Richichi,A.,Hofmann,K.-H.:1998,ESO Messenger92,113.Petrov,R.,et al.,2000,Proc.SPIE Conf.4006-074.Richichi,A.,et al.,2000,Proc.SPIE Conf.4006-085.Humphreys R.M.,Strecker D.W.,Murdock T.L.,Low,F.J.,1973,ApJ179,L496.Oudmaijer R.D.,Groenewegen M.A.T.,Matthews H.E.,Blommaert J.A.D.L,Sahu K.C.,1996,MNRAS280,10627.Klochkova V.G.,Chentsov E.L.,Panchuk,V.E.,1997,MNRAS292,198.Giguere P.T.,Woolf N.J.,Webber J.C.,1976,ApJ207,L1959.Mutel R.L.,Fix J.D.,Benson J.M.,Webber J.C.,1979,ApJ228,77110.Nedoluha G.E.,Bowers P.F.,1992,ApJ392,24911.Knapp G.R.,Morris M.,1985,ApJ292,64012.Jones T.J.,Humphreys R.M,Gehrz,R.D.et al.,1993,ApJ411,32313.Humphreys R.M.,Smith N.,Davidson K.et al.,1997,AJ114,277814.Oudmaijer R.D.,Geballe T.R.,Waters,L.B.F.M,Sahu K.C.,1994,A&A281,L33V i s i b i l i t y (λ=2.1µm )q [arcsec -1]Figure 6.Dependence of the error bars of simulated AT-VLTI/AMBER (wide-field mode)observations of IRC +10420at 2.11µm on (i)seeing differences of object and reference star observations (Fried parameter,r 0,differ-ences);(ii)different residual tip-tilt errors,δtt ,of object and reference star observations;and (iii)object brightness.Lines refer to radiative transfer models of different superwind amplitudes A .The symbols refer to AT-VLTI/AMBER simulations (each with N =2400interferograms)of the A =40intensity profile (see Fig.3)for baselines of 50and 100m (100m:q =227arc sec −1).To better distinguish between the simulations,the data points belonging to one baseline are somewhat shifted with respect to the spatial frequency.The error bars are based on a triple iteration of each simulation.The table gives the parameters of the simulations (a)to (d).Note that though simulation (d)corresponds to an IRC +10420-like intensity distribution,it refers to a much fainter object of K =11mag (instead of 3.5mag).15.Oudmaijer R.D.,1998,A&AS129,54116.Dyck H.,Zuckerman B.,Leinert C.,Beckwith S.,1984,ApJ287,80117.Ridgway S.T.,Joyce R.R.,Connors D.,Pipher J.L.,Dainty C.,1986,ApJ302,66218.Cobb M.L.,Fix J.D.,1987,ApJ315,32519.Christou J.C.,Ridgway S.T.,Buscher D.F.,HaniffC.A.,McCarthy Jr.D.W.,1990,Astrophysics with infraredarrays,R.Elston(ed.),ASP conf.series14,p.13320.Kastner J.,Weintraub D.A.,1995,ApJ452,83321.Bl¨o cker T.,Balega Y.,Hofmann K.-H.,Lichtenth¨a ler J.,Osterbart R.,Weigelt G.,1999,A&A348,80522.Weigelt G.,1977,Optics Commun.21,5523.Lohmann A.W.,Weigelt G.,Wirnitzer B.,1983,Appl.Opt.22,402824.Hofmann K.-H.,Weigelt G.,1986,A&A167,L15beyrie A.,1970,A&A6,8526.Walmsley C.M.,Chini R.,Kreysa E.et al.,1991,A&A248,55527.Craine E.R.,Schuster W.J.,Tapia S.,Vrba F.J.,1976,ApJ205,80228.Savage B.D.,Mathis J.S.,1979,ARA&A17,7329.Ivezi´cˇZ.,Nenkova M.Elitzur M.,1997,User Manual for DUSTY,University of Kentucky(/˜moshe/dusty)30.Ivezi´cˇZ.,Elitzur M.,1997,MNRAS287,79931.Mathis J.S.,Rumpl W.,Nordsieck K.H.,1977,ApJ217,42532.Draine B.T.,Lee H.M.,1984,ApJ285,8933.Roddier F.,1981,Progress in Optics XIX,28134.Tallon M.,Tallon-Bosc I.,1992,A&A253,64135.Malbet F.,et al.,2000,VLTI AMBER Instrument Analysis Review.。

斯仑贝谢所有测井曲线英文名称

斯仑贝谢所有测井曲线英文名称

斯仑贝谢所有测井曲线英文名称斯仑贝谢所有测井曲线英文名称OCEAN DRILLING PROGRAM ACRONYMS USED FOR WIRELINE SCHLUMBERGER TOOLSACT Aluminum Clay ToolAMS Auxiliary Measurement SondeAPS Accelerator Porosity SondeARI Azimuthal Resistivity ImagerASI Array Sonic ImagerBGKT Vertical Seismic Profile ToolBHC Borehole Compensated Sonic ToolBHTV Borehole TeleviewerCBL Casing Bond LogCNT Compensated Neutron ToolDIT Dual Induction ToolDLL Dual LaterologDSI Dipole Sonic ImagerFMS Formation MicroScannerGHMT Geologic High Resolution Magnetic Tool GPIT General Purpose Inclinometer ToolGR Natural Gamma RayGST Induced Gamma Ray Spectrometry Tool HLDS Hostile Environment Lithodensity Sonde HLDT Hostile Environment Lithodensity Tool HNGS Hostile Environment Gamma Ray Sonde LDT Lithodensity ToolLSS Long Spacing Sonic ToolMCD Mechanical Caliper DeviceNGT Natural Gamma Ray Spectrometry Tool NMRT Nuclear Resonance Magnetic Tool QSST Inline Checkshot ToolSDT Digital Sonic ToolSGT Scintillation Gamma Ray ToolSUMT Susceptibility Magnetic ToolUBI Ultrasonic Borehole ImagerVSI Vertical Seismic ImagerWST Well Seismic ToolWST-3 3-Components Well Seismic Tool OCEAN DRILLING PROGRAM ACRONYMS USED FOR LWD SCHLUMBERGER TOOLSADN Azimuthal Density-NeutronCDN Compensated Density-NeutronCDR Compensated Dual ResistivityISONIC Ideal Sonic-While-DrillingNMR Nuclear Magnetic ResonanceRAB Resistivity-at-the-BitOCEAN DRILLING PROGRAM ACRONYMS USED FORNON-SCHLUMBERGER SPECIALTY TOOLSMCS Multichannel Sonic ToolMGT Multisensor Gamma ToolSST Shear Sonic ToolTAP Temperature-Acceleration-Pressure Tool TLT Temperature Logging ToolOCEAN DRILLING PROGRAM ACRONYMS AND UNITS USED FOR WIRELINE SCHLUMBERGER LOGS AFEC APS Far Detector Counts (cps)ANEC APS Near Detector Counts (cps)AX Acceleration X Axis (ft/s2)AY Acceleration Y Axis (ft/s2)AZ Acceleration Z Axis (ft/s2)AZIM Constant Azimuth for DeviationCorrection (deg)APLC APS Near/Array Limestone Porosity Corrected (%)C1 FMS Caliper 1 (in)C2 FMS Caliper 2 (in)CALI Caliper (in)CFEC Corrected Far Epithermal Counts (cps) CFTC Corrected Far Thermal Counts (cps) CGR Computed (Th+K) Gamma Ray (API units)CHR2 Peak Coherence, Receiver Array, Upper DipoleCHRP Compressional Peak Coherence, Receiver Array, P&SCHRS Shear Peak Coherence, Receiver Array, P&SCHTP Compressional Peak Coherence, Transmitter Array, P&SCHTS Shear Peak Coherence, Transmitter Array, P&SCNEC Corrected Near Epithermal Counts (cps) CNTC Corrected Near Thermal Counts (cps)CS Cable Speed (m/hr)CVEL Compressional Velocity (km/s)DATN Discriminated Attenuation (db/m)DBI Discriminated Bond IndexDEVI Hole Deviation (degrees)DF Drilling Force (lbf)DIFF Difference Between MEAN and MEDIAN in Delta-Time Proc. (microsec/ft)DRH HLDS Bulk Density Correction (g/cm3) DRHO Bulk Density Correction (g/cm3)DT Short Spacing Delta-Time (10'-8' spacing; microsec/ft)DT1 Delta-Time Shear, Lower Dipole (microsec/ft)DT2 Delta-Time Shear, Upper Dipole (microsec/ft)DT4P Delta- Time Compressional, P&S (microsec/ft)DT4S Delta- Time Shear, P&S (microsec/ft)) DT1R Delta- Time Shear, Receiver Array, Lower Dipole (microsec/ft)DT2R Delta- Time Shear, Receiver Array,Upper Dipole (microsec/ft)DT1T Delta-Time Shear, Transmitter Array, Lower Dipole (microsec/ft)DT2T Delta-Time Shear, Transmitter Array, Upper Dipole (microsec/ft)DTCO Delta- Time Compressional (microsec/ft) DTL Long Spacing Delta-Time (12'-10' spacing; microsec/ft)DTLF Long Spacing Delta-Time (12'-10' spacing; microsec/ft)DTLN Short Spacing Delta-Time (10'-8' spacing; microsec/ftDTRP Delta-Time Compressional, Receiver Array, P&S (microsec/ft)DTRS Delta-Time Shear, Receiver Array, P&S (microsec/ft)DTSM Delta-Time Shear (microsec/ft)DTST Delta-Time Stoneley (microsec/ft)DTTP Delta-Time Compressional, Transmitter Array, P&S (microsec/ft)DTTS Delta-Time Shear, Transmitter Array,P&S (microsec/ft)ECGR Environmentally Corrected Gamma Ray (API units)EHGR Environmentally Corrected High Resolution Gamma Ray (API units)ENPH Epithermal Neutron Porosity (%) ENRA Epithermal Neutron RatioETIM Elapsed Time (sec)FINC Magnetic Field Inclination (degrees) FNOR Magnetic Field Total Moment (oersted) FX Magnetic Field on X Axis (oersted)FY Magnetic Field on Y Axis (oersted)FZ Magnetic Field on Z Axis (oersted)GR Natural Gamma Ray (API units)HALC High Res. Near/Array Limestone Porosity Corrected (%)HAZI Hole Azimuth (degrees)HBDC High Res. Bulk Density Correction(g/cm3)HBHK HNGS Borehole Potassium (%)HCFT High Resolution Corrected Far Thermal Counts (cps)HCGR HNGS Computed Gamma Ray (APIunits)HCNT High Resolution Corrected Near Thermal Counts (cps)HDEB High Res. Enhanced Bulk Density (g/cm3) HDRH High Resolution Density Correction(g/cm3)HFEC High Res. Far Detector Counts (cps)HFK HNGS Formation Potassium (%)HFLC High Res. Near/Far Limestone Porosity Corrected (%)HEGR Environmentally Corrected High Resolution Natural Gamma Ray (API units) HGR High Resolution Natural Gamma Ray (API units)HLCA High Res. Caliper (inHLEF High Res. Long-spaced Photoelectric Effect (barns/e-)HNEC High Res. Near Detector Counts (cps) HNPO High Resolution Enhanced Thermal Nutron Porosity (%)HNRH High Resolution Bulk Density (g/cm3) HPEF High Resolution Photoelectric Effect(barns/e-)HRHO High Resolution Bulk Density (g/cm3) HROM High Res. Corrected Bulk Density(g/cm3)HSGR HNGS Standard (total) Gamma Ray (API units)HSIG High Res. Formation Capture Cross Section (capture units)HSTO High Res. Computed Standoff (in) HTHO HNGS Thorium (ppm)HTNP High Resolution Thermal Neutron Porosity (%)HURA HNGS Uranium (ppm)IDPH Phasor Deep Induction (ohmm)IIR Iron Indicator Ratio [CFE/(CCA+CSI)]ILD Deep Resistivity (ohmm)ILM Medium Resistivity (ohmm)IMPH Phasor Medium Induction (ohmm)ITT Integrated Transit Time (s)LCAL HLDS Caliper (in)LIR Lithology Indicator Ratio [CSI/(CCA+CSI)] LLD Laterolog Deep (ohmm)LLS Laterolog Shallow (ohmm)LTT1 Transit Time (10'; microsec)LTT2 Transit Time (8'; microsec)LTT3 Transit Time (12'; microsec)LTT4 Transit Time (10'; microsec)MAGB Earth's Magnetic Field (nTes)MAGC Earth Conductivity (ppm)MAGS Magnetic Susceptibility (ppm)MEDIAN Median Delta-T Recomputed (microsec/ft)MEAN Mean Delta-T Recomputed (microsec/ft)NATN Near Pseudo-Attenuation (db/m)NMST Magnetometer Temperature (degC)NMSV Magnetometer Signal Level (V)NPHI Neutron Porosity (%)NRHB LDS Bulk Density (g/cm3)P1AZ Pad 1 Azimuth (degrees)PEF Photoelectric Effect (barns/e-)PEFL LDS Long-spaced Photoelectric Effect (barns/e-)PIR Porosity Indicator Ratio [CHY/(CCA+CSI)]POTA Potassium (%)RB Pad 1 Relative Bearing (degrees)RHL LDS Long-spaced Bulk Density (g/cm3)RHOB Bulk Density (g/cm3)RHOM HLDS Corrected Bulk Density (g/cm3)RMGS Low Resolution Susceptibility (ppm)SFLU Spherically Focused Log (ohmm)SGR Total Gamma Ray (API units)SIGF APS Formation Capture Cross Section (capture units)SP Spontaneous Potential (mV)STOF APS Computed Standoff (in)SURT Receiver Coil Temperature (degC)SVEL Shear Velocity (km/s)SXRT NMRS differential Temperature (degC) TENS Tension (lb)THOR Thorium (ppm)TNRA Thermal Neutron RatioTT1 Transit Time (10' spacing; microsec)TT2 Transit Time (8' spacing; microsec)TT3 Transit Time (12' spacing; microsec)TT4 Transit Time (10' spacing; microsec) URAN Uranium (ppm)V4P Compressional Velocity, from DT4P (P&S; km/s)V4S Shear Velocity, from DT4S (P&S; km/s)VELP Compressional Velocity (processed from waveforms; km/s)VELS Shear Velocity (processed from waveforms; km/s)VP1 Compressional Velocity, from DT, DTLN, or MEAN (km/s)VP2 Compressional Velocity, from DTL, DTLF, or MEDIAN (km/s)VCO Compressional Velocity, from DTCO (km/s)VS Shear Velocity, from DTSM (km/s)VST Stonely Velocity, from DTST km/s)VS1 Shear Velocity, from DT1 (Lower Dipole; km/s)VS2 Shear Velocity, from DT2 (Upper Dipole; km/s)VRP Compressional Velocity, from DTRP (Receiver Array, P&S; km/s)VRS Shear Velocity, from DTRS (Receiver Array, P&S; km/s)VS1R Shear Velocity, from DT1R (Receiver Array, Lower Dipole; km/s)VS2R Shear Velocity, from DT2R (Receiver Array, Upper Dipole; km/s)VS1T Shear Velocity, from DT1T (Transmitter Array, Lower Dipole; km/s)VS2T Shear Velocity, from DT2T (TransmitterArray, Upper Dipole; km/s)VTP Compressional Velocity, from DTTP (Transmitter Array, P&S; km/s)VTS Shear Velocity, from DTTS (Transmitter Array, P&S; km/s)#POINTS Number of Transmitter-Receiver Pairs Used in Sonic ProcessingW1NG NGT Window 1 counts (cps)W2NG NGT Window 2 counts (cps)W3NG NGT Window 3 counts (cps)W4NG NGT Window 4 counts (cps)W5NG NGT Window 5 counts (cps)OCEAN DRILLING PROGRAMACRONYMS AND UNITS USED FOR LWD SCHLUMBERGER LOGSAT1F Attenuation Resistivity (1 ft resolution; ohmm)AT2F Attenuation Resistivity (2 ft resolution; ohmm)AT3F Attenuation Resistivity (3 ft resolution; ohmm)AT4F Attenuation Resistivity (4 ft resolution; ohmm)AT5F Attenuation Resistivity (5 ft resolution; ohmm)ATR Attenuation Resistivity (deep; ohmm)BFV Bound Fluid Volume (%)B1TM RAB Shallow Resistivity Time after Bit (s)B2TM RAB Medium Resistivity Time after Bit (s)B3TM RAB Deep Resistivity Time after Bit (s)BDAV Deep Resistivity Average (ohmm)BMAV Medium Resistivity Average (ohmm)BSAV Shallow Resistivity Average (ohmm)CGR Computed (Th+K) Gamma Ray (API units)DCAL Differential Caliper (in)DROR Correction for CDN rotational density (g/cm3).DRRT Correction for ADN rotational density (g/cm3).DTAB AND or CDN Density Time after Bit (hr)FFV Free Fluid Volume (%)GR Gamma Ray (API Units)GR7 Sum Gamma Ray WindowsGRW7+GRW8+GRW9-Equivalent to Wireline NGT window 5 (cps)GRW3 Gamma Ray Window 3 counts(cps)-Equivalent to Wireline NGT window 1GRW4 Gamma Ray Window 4 counts(cps)-Equivalent to Wireline NGT window 2GRW5 Gamma Ray Window 5 counts(cps)-Equivalent to Wireline NGT window 3GRW6 Gamma Ray Window 6 counts(cps)-Equivalent to Wireline NGT window 4GRW7 Gamma Ray Window 7 counts (cps) GRW8 Gamma Ray Window 8 counts (cps) GRW9 Gamma Ray Window 9 counts (cps) GTIM CDR Gamma Ray Time after Bit (s) GRTK RAB Gamma Ray Time after Bit (s) HEF1 Far He Bank 1 counts (cps)HEF2 Far He Bank 2 counts (cps)HEF3 Far He Bank 3 counts (cps)HEF4 Far He Bank 4 counts (cps)HEN1 Near He Bank 1 counts (cps)HEN2 Near He Bank 2 counts (cps)HEN3 Near He Bank 3 counts (cps)HEN4 Near He Bank 4 counts (cps)MRP Magnetic Resonance PorosityNTAB ADN or CDN Neutron Time after Bit (hr)PEF Photoelectric Effect (barns/e-)POTA Potassium (%) ROPE Rate of Penetration (ft/hr)PS1F Phase Shift Resistivity (1 ft resolution; ohmm)PS2F Phase Shift Resistivity (2 ft resolution; ohmm)PS3F Phase Shift Resistivity (3 ft resolution; ohmm)PS4F Phase Shift Resistivity (4 ft resolution; ohmm)PS5F Phase Shift Resistivity (5 ft resolution; ohmm)PSR Phase Shift Resistivity (shallow; ohmm)RBIT Bit Resistivity (ohmm)RBTM RAB Resistivity Time After Bit (s)RING Ring Resistivity (ohmm)ROMT Max. Density Total (g/cm3) from rotational processingROP Rate of Penetration (m/hr)ROP1 Rate of Penetration, average over last 1 ft (m/hr).ROP5 Rate of Penetration, average over last 5 ft (m/hr)ROPE Rate of Penetration, averaged over last 5 ft (ft/hr)RPM RAB Tool Rotation Speed (rpm)RTIM CDR or RAB Resistivity Time after Bit (hr)SGR Total Gamma Ray (API units)T2 T2 Distribution (%)T2LM T2 Logarithmic Mean (ms)THOR Thorium (ppm)TNPH Thermal Neutron Porosity (%)TNRA Thermal RatioURAN Uranium (ppm)OCEAN DRILLING PROGRAM ADDITIONAL ACRONYMS AND UNITS(PROCESSED LOGS FROM GEOCHEMICAL TOOL STRING)AL2O3 Computed Al2O3 (dry weight %)AL2O3MIN Computed Al2O3 Standard Deviation (dry weight %)AL2O3MAX Computed Al2O3 Standard Deviation (dry weight %)CAO Computed CaO (dry weight %)CAOMIN Computed CaO Standard Deviation (dry weight %)CAOMAX Computed CaO Standard Deviation (dry weight %)CACO3 Computed CaCO3 (dry weight %)CACO3MIN Computed CaCO3 Standard Deviation (dry weight %)CACO3MAX Computed CaCO3 Standard Deviation (dry weight %)CCA Calcium Yield (decimal fraction)CCHL Chlorine Yield (decimal fraction)CFE Iron Yield (decimal fraction)CGD Gadolinium Yield (decimal fraction)CHY Hydrogen Yield (decimal fraction)CK Potassium Yield (decimal fraction)CSI Silicon Yield (decimal fraction)CSIG Capture Cross Section (capture units) CSUL Sulfur Yield (decimal fraction)CTB Background Yield (decimal fraction) CTI Titanium Yield (decimal fraction)FACT Quality Control CurveFEO Computed FeO (dry weight %)FEOMIN Computed FeO Standard Deviation (dry weight %)FEOMAX Computed FeO Standard Deviation(dry weight %)FEO* Computed FeO* (dry weight %)FEO*MIN Computed FeO* Standard Deviation (dry weight %)FEO*MAX Computed FeO* Standard Deviation (dry weight %)FE2O3 Computed Fe2O3 (dry weight %)FE2O3MIN Computed Fe2O3 Standard Deviation (dry weight %)FE2O3MAX Computed Fe2O3 Standard Deviation (dry weight %)GD Computed Gadolinium (dry weight %)GDMIN Computed Gadolinium Standard Deviation (dry weight %)GDMAX Computed Gadolinium Standard Deviation (dry weight %)K2O Computed K2O (dry weight %)K2OMIN Computed K2O Standard Deviation (dry weight %)K2OMAX Computed K2O Standard Deviation (dry weight %)MGO Computed MgO (dry weight %)MGOMIN Computed MgO Standard Deviation (dry weight %)MGOMAX Computed MgO Standard Deviation (dry weight %)S Computed Sulfur (dry weight %)SMIN Computed Sulfur Standard Deviation (dry weight %)SMAX Computed Sulfur Standard Deviation (dry weight %)SIO2 Computed SiO2 (dry weight %)SIO2MIN Computed SiO2 Standard Deviation (dry weight %)SIO2MAX Computed SiO2 Standard Deviation (dry weight %)THORMIN Computed Thorium Standard Deviation (ppm)THORMAX Computed Thorium Standard Deviation (ppm)TIO2 Computed TiO2 (dry weight %)TIO2MIN Computed TiO2 Standard Deviation (dry weight %)TIO2MAX Computed TiO2 Standard Deviation (dry weight %)URANMIN Computed Uranium Standard Deviation (ppm)URANMAX Computed Uranium Standard Deviation (ppm)VARCA Variable CaCO3/CaO calcium carbonate/oxide factor。

A COMPARISON OF VERTICAL MOTIONS OB

A COMPARISON OF VERTICAL MOTIONS OB

5. TERM
NEGLECTS DEFORMATION
5. EXTENDED
6. WIND
INCLUDES AGEOSTROPHIC
• Surface/Upper-Air Data: 0000 UTC Nov. 2nd to 0000 UTC Nov. 3rd
• 2-pass Barnes Analysis Scheme
• SYNOPTIC PATTERNS OMEGA vs. SYNOPTIC FEATURES
• PRECIPITATION PATTERNS % PRECIP IN UPWARD MOTION
12 UTC 2 Nov. 1999
(TOP): SL Pressure & 6hr. Precip.
(MIDDLE): 500-mb Height & Abs. Vort.
• Case study: Strahl and Smith (MWR:2001) – Explosive cyclone development accompanied by upper-air wave merger over North America (Nov. 2-3 1999)
INTRODUCTION
4. Isentropic – Pressure Advection along Potential Temperature Surface
5. Omega Equation – Combine Vorticity Equation and 1st Law of Thermodynamics
QG THEORY
• QG APPROXIMATION WORKS QUITE WELL EVEN FOR A CASE OF EXPLOSIVE CYCLONE DEVELOPMENT.

柴达木盆地咸湖相烃源岩特征——以英西地区下干柴沟组上段为例

柴达木盆地咸湖相烃源岩特征——以英西地区下干柴沟组上段为例

断块油气田FAULT-BLOCK OIL&GAS FIELD第28卷第2期2021年3月doi:10.6056/dkyqt202102007柴达木盆地咸湖相H源岩特征——以英西地区下干柴沟组上段为例舒豫川,胡广,庞谦,胡朝伟,夏青松,谭秀成(西南石油大学地球科学与技术学院,四川成都610500/摘要盐湖环境对油气生成具有重要的地质意义!国内外有诸多含油气盆地发育咸湖相?源岩!自新生代以来,在青藏高原的隆升作用下,柴达木盆地的海拔升高、湖盆逐渐封闭、持续的干寒气候及盐源的充分供应,使得柴达木盆地逐渐变成了一个典型的高原咸化湖盆,其中,下干柴沟组为柴达木盆地西部英西地区的主力?源岩层系!文中通过有机岩石学、有机地球化学及生物标志物地球化学研究,对英西地区下干柴沟组上段进行了?源岩评价,并探讨了其沉积环境及生?母质类型!研究结果表明:英西地区下干柴沟组上段?源岩总有机碳质量分数在0.31%〜1.49%,平均值为0.82%;氯仿沥青“A”质量分数在562.16x10-6~2999.65X10"6,平均值为1619.59x10^o?源岩生?母质为藻类、细菌和高等植物,以低等水生生物为主!有机质类型主要为!型,其次为"型,?源岩成熟度介于低成熟一成熟!英西地区下干柴沟组?源岩综合评价为好的?源岩!生物标志化合物特征表明,下干柴沟组上段?源岩沉积于水体盐度较高、还原性较强的湖相环境之中!关键词?源岩;生物标志化合物;下干柴沟组;英西地区;柴达木盆地中图分类号:TE13513;P618.13文献标志码:ACharacteristics of source rocks of salt lake facies in Qaidam Basin:taking XiaganchaigouFormation in Yingxi region as an exampleSHU Yuchuan?HU Guang,PANG Qian?HU Chaowei,XIA Qingsong?TAN Xiucheng(School of Geoscience and Technology?Southwest Petroleum University?Chengdu61O5OO?China) Abstract:Salt lake environment has important geological significance for oil and gas generation.There are many oil-bearing basins containing salt lake facies source rocks at home and abroad.Since Cenozoic era,with the uplift of the Tibetan Plateau,the Qaidam Basin has gradually become a typical plateau salt lake basin because of the altitude rise,the lake basin gradually closed,the persistent dry and cold climate,and the adequate supply of salt.In the basin,the Xiaganchaigou Formation is the main source rock series in the Yingxi region which located in the western part of Qaidam Basin.Based on study of organic petrology,organic geochemistry and biomarker geochemistry,the source rocks of the Xiaganchaigou Formation in the Yingxi region were evaluated, and their sedimentary environment and the types of hydrocarbon generating parent material were discussed.The results show that the total organic carbon content of source rocks in the Xiaganchaigou Formation in the Yingxi region is0.31%to1.49%,with an average value of0.82%.The value range of chloroform bitumen"A H is562.16X10"6to2,999.65X10"6,with an average of1,619.59x10"6.The source rocks are mainly composed of algae,bacteria and higher plants,mainly low aquatic organisms.The organic matter type is mainly type!,followed by type".The maturity of source rock is between low maturity and prehensive evaluation shows that the source rocks of the Xiaganchaigou Formation in the Yingxi region are good source rocks.The characteristics of biomarkers indicate that the source rocks of Xiaganchaigou Formation are deposited in the lacustrine environment which has high salinity and strong reducibility.Key words:source rocks;biomarker compounds;Xiaganchaigou Formation;Yingxi region;Qaidam Basin近年来,在咸湖相碳酸盐油气地质理论和三维地震攻关技术的助推下,柴达木盆地英西地区古近系盐下油气勘探获得重大突破,多口钻井在下干柴沟组上段盐间和盐下湖相碳酸盐岩中获得了日产超千吨的高产油气流,揭示了柴达木盆地古近系沉积体系所具有收稿日期:2020-07-30"改回日期:2021-01-07o第一作者:舒豫川,1993年生,男,在读硕士研究生,从事沉积岩石学方面的研究工作。

天文词汇中英文对照

天文词汇中英文对照

天文词汇中英文对照(源自网络,仅供参考)天文词汇检索(按英文字母序)(A)absolute bolometric magnitude绝对热星等:量测恒星的绝对星等时,如果我们能同时侦测到所有的波长,所量测到就是绝对热星等。

absolute visual magnitude (Mv)绝对视星等:恒星的真正亮度。

定义为恒星在距离我们 10 秒差距 (32.6光年) 时的视星等。

absolute zero绝对零度,绝对零点:温度的最低点。

在绝对零度时,不可能再抽取物质、原子或分子里的粒子的动能。

absorption line吸收谱线:光谱里的暗线。

来自天体的光,被原子或分子选择性的吸收,导致那部分的光从星光中被消去,留下一条条的暗线。

absorption spectrum (dark line spectrum)吸收光谱:含有吸收谱线的光谱。

acceleration 加速度:速度对时间的变率。

速度在方向或速率的改变,就代表有加速度的存在。

参考速度(velocity) 。

acceleration of gravity重力加速度:量测天体附近重力强度的物理量。

accretion吸积:粘合较小的固体粒子,形成更大粒子的过程。

accretion disk吸积盘:指白矮星、中子星或黑洞等致密天体周围,由气体所形成的盘状天体。

Achernar水委一:波江座α星。

achondrites无球粒陨石:不含球粒(chondrules) 或挥发性物质的石陨石。

achromatic lens消色差透镜:由两种不同材质的透镜组合而成,消色差透镜的用途是把两种不同颜色的光聚焦到同一点,或称为修正色像差(chromatic aberration)。

active galactic nucleus (AGN)活跃星系核:活跃星系中央的能量源。

active galaxy活跃星系:一种会发生极大量无线电波、X射线、珈玛射线辐射的星系。

active optics活动光学:由计算机控制,来改变光学组件的位置和形状,以补偿大气扰动所产生的散焦效应之光学系统。

地质学专业术语翻译

地质学专业术语翻译

凹坑pool,chute 白垩chalk白云化灰岩dolomitized limestone白云母muscovite白云石dolomite白云岩dolostone, dolomitite斑点构造mottled structure斑礁patch reef斑晶porphyrotope斑脱岩bentonite搬运作用transportation板状晶体laths半深海helmipelagic半自形晶subhedral棒形蛀木虫迹teredolites clavatus包卷层理convolute lamination包卷构造comvolute structure堡礁(障壁礁) barrier reef抱球虫软泥globigerina ooze比拟试验模式scaled experimental model碧玉jasper边滩point bar辫状河braided channel标志层key bed标准偏差standard deviation表面结构surface texture表皮鲕superficial Oolits表栖动物epifauna表生成岩作用epidiagenesis表生作用epigenesis表土regolith滨湖亚相lake-shore滨面shoreface滨线shoreline冰雹痕hail impression冰川glacier冰-海沉积物glacier-marine sediment 1/21冰川相glacial facies冰碛物till冰碛岩tillite冰水沉积outwash波长wave length波高wave height波痕ripple波状层理wavy bedding剥离线理parting lineation不等粒inequigranular不整合unconformity槽痕flute槽模flute casts侧向沉积作用lateral sedimentation侧向加积lateral accretion测井沉积学well-logging sedimentology 测井相electrofacies层bed层间流interflow层孔虫stromatoporoid层理stratification,bedding,lamination 层流laminar flow层内溶解作用intrastratal solution层系set层系组coset层序地层学sequence stratigraphy次长石净砂岩subarkose次生孔隙secondary porosity次岩屑净砂岩sublitharenite长石feldspar长石砂岩feldspar sandstone长石净砂岩arkose长石净砂岩类arkosic arenite长石杂砂岩feldspathic greywacke长石杂砂岩类arkosic greywacke潮道tidal channel潮间带intertidal潮控三角洲tidal-dominated delta2/21潮坪tidal flat潮上带supratidal潮汐层理tidal- bedding潮汐流ebb current潮汐三角洲tidal delta潮下带subtidal zone沉积分异作用sedimentary differentiation沉积构造sedimentary structure沉积后作用postdeposition沉积环境depositional environment沉积模式sedimentary model沉积盆地sedimentary basin沉积铁矿sedimentary iron ores沉积物sediment沉积物重力流sediment gravity flow沉积相sedimentary(depositional)facies 沉积学sedimentology沉积岩sedimentary rocks沉积岩石学sedimentary petrology沉积作用sedimentation,deposition成分成熟度compositional matuarity成岩作用diagenesis赤铁矿hematite冲积平原alluvial plain冲积扇alluvial fan冲裂avulsion冲流swash冲流回流带zone of swash and backwash冲刷痕scour mark冲洗交错层理swash bedding 冲溢扇washover fan储层reservoir储油窗oil–windows床沙bed床沙形体bed form床沙载荷bed load吹程fetch垂向沉积作用vertical accretion 3/21丛藻迹chondrites粗砾cobble粗枝藻门dasyclads大潮spring tide大陆架continental shelf大陆隆continental rise大陆坡continental slope大洋盆地basin带穴迹taenidium带状迹teanidium袋状钻孔迹rogerella单杯迹monocraterion单成份砾岩oligomictic conglomerate 岛(弧)-海(沟)体系arc-trench system 岛弧island arc等深流contour current等深积岩contourite等趾迹isopodichnus等足迹isopodichnus低镁方解石low-magnesian calcite低弯度河道low-sinuosity channel低位体系域low-stand systems tract 堤岸bank底痕sole mark底积层bottomset底模solecast底栖生物benthos底载河道bedload channel地层strata地层倾斜录井dipmeter log地层水formation water地球化学geochemistry地球物理geophysics地震地层学seismic stratigraphy地震勘探seismic exploration递变层graded layer递变悬浮作用graded suspension 4/21典型浊积岩clastic turbidite电测井electric log叠层石stromatolite叠瓦状排列imbricate arrangement 顶积层topsets动态模式dynamic model动藻迹zoophycus豆粒pisolite段member对称波痕symmetrical ripple鲕粒oolite二叶石迹cruziana反递变层理reverse grading反沙丘层理antidune bedding反向递变的inversely graded泛滥盆地flood basin泛滥平原flood plain方沸石Analcime方解石calcite放射虫Radiolarians放射性测井radioactivity log非补偿(饥饿)沉积starved basin 非显晶质的aphanitic沸石zeolite分流河道distributary channel分流河口砂坝distributary-mouth bar 分流间湾interdistributary分选系数So分选作用sorting粉砂silt粉砂岩siltite,siltstone粉砂页岩silt shale球粒pellet风暴大潮storm surge风成eolian风成波痕eolian ripple,wind ripple风成交错层理aeolian cross bedding 风成岩eolianite5/21风化作用weathering风生流wind-drift current风蚀deflation峰度curtosis缝合线stylolite弗劳德数Froude number浮游生物plankton腐泥sapropel负载(荷)构造load structure复成分砾岩polymictic conglomerate 复理石flysch钙结层caliche干酪根kerogen干盐湖playa高岭石kaolinite高弯度河道high-sinuosity channel 高位体系域high-stand systems tract根珊瑚迹rhizocoralliun耕作迹agrichnia沟模groove casts沟蠕形迹taphrhelminthopsis古拉迹paleohelcura古生态学paleoecology古生物paleontological古水流标志paleocurrent indicator古土壤paleosols古网迹palaeodictyon古斜坡paleoslope骨架岩framestone固着底栖生物sessile benthos光合作用photosynthesis硅镁层sima硅土,二氧化硅silica硅质silic-硅质碎屑siliciclastic国际沉积学家学会IAS(International Association of Sedimentologists) 文档来自于网络搜索海(湖)泛面flooding surface6/21海(湖)滩beach海岸带coastal zone海底扇submarine fan海底峡谷submarine canyon海沟trench海进体系域transgressive systems tract 海浪sea wave海绿石glauconite海平面变化sea-level changes海平面升降eustasy海侵transgression海侵层序transgressive sequence海生迹thalassinoides海滩marine beach海滩脊ridge of a beach海相三角洲marine delta海洋ocean旱谷arroyo, wadi河床(道)channel河床滞留lag河底滞留沉积channel floor lag河控三角洲fluvial delta河口湾estuary河流相fluvial河漫overbank核形石oncilite褐铁矿limonite黑色页岩black shale黑云母biotite红藻植物门rhodophyta后滨backshore后成作用katagenesis,anadiagenesis 弧间盆地interarc basin弧内盆地intra-arc basin弧前盆地forearc basin湖lake湖泊环境lacustrine environment湖泊相lacustrine7/21湖面波动seiche湖沼学limnology滑[动]痕slide mark滑塌构造slump structure滑塌岩slumps还原硫酸盐细菌desulfouibrio,desulfomaculmi 还原作用reduction环境沉积学environmental sedimentology环线迹spirorhaphe黄铁矿pyrite黄土loess灰泥micrite灰泥岩mudstone灰质白云岩calcareous dolostone回流backflow,backwash混杂沉积岩diamicton混杂堆积物mélange混载荷流mixed-load channel火山弹bomb火山灰ash火山灰流ash flow火山角砾岩volcanic breccia火山块block火山泥石流lahar火山碎屑pyroclastic debris,tephra火山碎屑流pyroclastic flow火山碎屑岩pyroclastic rock火焰构造flame structure鸡笼状的chickenwire集合粒aggregate集块岩agglomerate加积accretion加积作用aggradation假亮晶pseudospar尖峰度的jeptokurtic剪切shearing交错层cross beds交错层理cross bedding; cross stratification 8/21交错纹crossed-lamellar交错纹层cross laminae胶结物cement胶结作用cementation胶磷矿collophane礁reef礁核reef core礁后相backreef facies礁麓堆积reeftalus礁前相forereef facies礁相reef facies角度不整合angular unconformity 角砾岩breccias节藻迹phycodes结晶碳酸盐岩crystalline carbonate 界面bounding surfaces金藻门chrysophyta津(波)浪tsunami近滨nearshore进潮口tidal inlet进积作用progradation进食构造feeding structures进食迹fodinichnia晶面crystal face晶体crystal井形穴shafts净砂岩arenite静态模式static model居住构造dwelling structure居住迹domichnia居住潜穴dwelling burrows巨画迹megagrapton巨砾boulder决口扇crevasse,crevasse splay均匀悬浮uniform suspension颗粒grain,particle颗粒流grain flow颗粒岩grainstone9/21颗粒支撑的particle -supported,grain-supported 可劈性fissility可容空间accommodation space克拉通craton孔隙(度)pore,porosity块状层理massive bedding块状砂岩massive sandstone浪成波痕wave-generated ripple浪成沙纹层理wave- ripple bedding浪基面wave base浪控三角洲wave-dominated delta雷诺数Reynolds number-Re类沙蚕迹nereites累积频率曲线cumulative-frequency curve 离岸(裂)流rip current犁沟迹aulichnites丽线迹cosmorhaphe砾gravel砾漠reg砾屑灰岩calcirudite砾屑岩rudite砾岩conglomerate粒度grain size粒度概率图probability plot粒度中值Md粒间孔隙度intergranular porosity粒屑灰岩grainstone粒序层理graded bedding亮晶sparite亮晶方解石sparry calcite亮晶灰岩sparite亮煤clarain裂隙fracture临滨shoreface磷灰岩phosphorite磷酸盐岩phosphate rock菱铁矿siderite流动体制flow regimes10/21流痕rill impression流态flow regime流体化作用fluidization硫酸盐sulphate漏斗形funnel shape陆表海epicontinental sea陆架边缘体系域shelf margin systems tract 陆相terrestrial陆缘海marginal sea陆源碎屑岩terrigenous clastic rocks陆表海epeiric sea绿泥岩chlorite绿藻门chlorophyta,卵石质砂岩pebbly sandstone螺管迹gyrolithes洛伦茨迹corenzinia脉状层理flaser bedding漫流sheet flow漫游迹planoliies煤coal煤化作用coalification锰质岩manganess rock觅食迹pascichnia觅食拖迹grazing trails密度底流density underflow密度流density current密度跃层pycnocline磨拉石molasses母岩source rock, parent rock母岩区provenance,source area 内成岩endogenetic内栖动物infauna内碎屑intraclast钠长石albite钠长石化albitization能量带energy belt泥mud泥板岩argillite11/21泥晶micrite泥晶套micrite envelope泥粒灰岩(泥质颗粒岩)packstone 泥裂mud crack泥流mud flow泥石流debris flow泥质岩mudrock泥页岩argillutite泥支撑的mud supported泥质的argil拟蠕形迹helminthopsis逆行沙波regressive sand wave逆行沙丘antidune年代地层单位chronostratigraphic unit鸟眼构造birds-eyes structure鸟足状三角洲bird-foot delta凝灰岩tuff凝聚颗粒aggregate grains凝块灰岩thrombolite牛顿流体Newtonian fluid牛轭湖ox-bow lake爬迹repichnia爬升沙纹层理climbing ripple bedding 爬行crawling拍岸浪surf盘旋虫科spirillinidae螃蟹迹psilonichnus跑动running泡沫痕foam impression盆内岩intrabasinal rock盆外岩石extrabasinal rock片流sheetwash偏度Sk1,skewness漂砾erratic平均粒径Mz平均值mean平行层理parallel bedding平直河straight channel12/21破浪breaker破浪带breaker zone葡萄石grapestones气候带climate belts牵引[作用] traction牵引流tractive current牵引毯状沉积traction carpet前滨foreshore前积层foresets前三角洲沉积prodelta deposit浅海区neritic province浅滩shoal墙形迹teichichus丘状交错层理hummocky cross bedding 球度sphericity球粒pellet,pellet peloids球状粒peloid曲带迹curuolithus曲流meander曲流裁直neck cutoff曲流河meandering river取心钻井core drilling全球海平面变化eustatic change of sea level 缺氧事件anoxic event群group群体生物colony(organisms)热云nuee ardente热云岩ignimbrite溶解载荷solutional load熔结凝灰岩ignimbrite蠕形迹helminthoida软流圈asthenosphere软泥ooze弱光带disphotic zone萨布哈(盐沼)sabkha三角洲delta三角洲平原delta plain三角洲前缘delta front13/21三联点triple junction扫毛虫sabellarians沙坝bar沙波sand wave沙流sand flow沙漠desert沙漠相desert facies沙丘dune沙蜀迹arenicolites沙体Sand body沙纹ripples沙嘴spit砂sand砂火山sand volcanoes砂流sand flow砂屑detritus砂屑灰岩calcarenite砂岩sandstone砂枕构造pillow structure 筛状沉积sieve deposition山鹿冲积平原bajada山麓piedmont珊瑚coral扇砾岩fanglomerate扇三角洲fan delta扇状三角洲fan-like delta 舌菌迹glossiofunhites蛇形迹ophiomorpha深海pelagic, deep sea深海带abyssal zone深海沟deep-sea trench/2114深海平原abyssal深海扇deep-sea fan深泓线thalweg深水bathymetric深水盆地bathymetric basin 渗透率permeability生长构造growth structure 生物层biostrome生物带biozone生物地层学biostratigraphy生物灰岩biolithite生物礁reef生物颗粒bioclast生物亮晶灰岩biosparite生物内碎屑bioclastic生物泥晶灰岩biomicrite生物丘bioherm生物扰动作用bioturbation生物遗迹化石trace fossil石膏gypsum石化作用lithification石灰岩limestone石漠hamada石英quarts石英砂岩quartz sandstone,silicarenita 石英净砂岩quartzarenite石英杂砂岩quartzwack石针迹skolithos石枝藻litho-thammium始成岩作用eodiagenesis示顶底标志geopetal criterion事实模式actual model数学模式mathematical model刷模brush casts双杯迹diplocraterion双形迹dimorphichus双趾迹diplichifes水道channel15/21水流current,aqueous flow水流痕current mark水流上涌upwelling水流线理current lineation水平层理horizontal bedding水体分层stratification of the water column 水下泥石流subaqueous debis flow水下岩崩subaqueous rockfall水悬浮作用aqueous suspension斯科阳迹scoyenia似海生迹thalassinoides似海星迹asteriacites似蠕形迹helminthopsis似藻迹phycodes碎屑颗粒detrital particle燧石flint,chert他生沉积物allochthonous sediments它形晶anhedral塌积角砾岩collapse braccia塔斯曼迹tasmanadia台地platform苔藓虫fenestrate bryozoans滩bank滩脊beach ridge碳酸盐carbonate碳酸盐岩carbonate rock, carbonatite碳酸盐岩补偿深度carbonate compensation depth 碳酸盐岩隆(建隆)carbonate buildup逃逸构造escape structure逃逸迹fugichnia体系域systems tract天然堤levee,natural levee填集作用packing跳模bounce casts跳跃作用saltation铁岩ironstone铁质岩ferruginous rock停息迹cubchnia16/21同沉积滑塌构造syndepositinal slump structure 同生白云岩syngenetic dolostone同生作用syngenesis透光带photic zone透镜状层理lenticular bedding湍流turbulent flow团块lumps退覆层序offlap sequence洼状交错层理swaley cross bedding外成岩exogenetic rock外来岩块exotic block外生沉积物exogentic sediments弯曲锚形迹ancorichnus coronus晚成岩作用telodiagenesis网状(交织)河流anastomosing stream,anastomosed river 微亮晶microspar微斜长石microcline胃形钻孔迹gastrochaenolites温跃层thermoclime纹层laminae,lamination纹层状的laminated纹理lamination纹泥,季节泥varve雯石(霰石)aragonite紊流turbulent flow无光带aphotic zone无粘结性颗粒cohesionless particle物源区provenance矽藻diatom细砾granule线管迹tigllites相facies相变facies change相标志facies mark相对海平面relative sea-level相序facies sequence(succession)向上变粗coarsening-upward向上变浅shoaling-upward,shallowing--upward 17/21小潮neap tide小间断diastem斜坡地形clinoform泄水构造water-escape structure泻湖lagoon新生变形作用neomorphism新月形沙丘barchan dune信风trade winds星瓣迹asterosoma休息迹resting traces悬浮(移)负载suspended load选择交代作用selective replacement压扁层理flaser bedding压刻痕toll mark压溶作用pressure solution压实作用compaction盐水楔saltwater wedge岩石地层单位lithostratigraphic unit岩石地层学lithostratigraphy岩石圈lithosphere岩石学petrology岩相lithofacies岩相古地理lithofacies paleogeography岩屑rock fragment,lithic岩屑砂岩lithic sandstone,rock fragment sandstone 岩屑净砂岩litharenite岩屑杂砂岩lithic greywacke岩屑质长石净砂岩lithic arkose岩性lithological沿岸流longshore current沿岸沙坝longshore bar盐度salinity盐湖salt lakes盐碱滩salina氧化作用oxidation页岩shale液化liguefaction液化沉积物流fluidized sediment flow18/21液化作用liquefaction伊利石illite遗迹化石ichnofossil,trace fossil遗迹相ichofacies异化颗粒allochem粘结灰岩bindstones硬砂岩greywacke硬石膏anhydrite涌浪swells尤尔斯特龙图解Hiulstrom diagram油气储层地质建模reservoir geological Modeling 油页岩oil shale游泳生物nekton organisms原生孔隙primary porosity有效孔隙度effective porosity羽状交错层理herringbone cross bedding雨痕rain impression原始古网迹protopalaeodictyon圆度roudness越岸沉积overbank deposit运动迹moving traces韵律层理rhythmic bedding杂基matrix杂基支撑matrix supported杂砂岩graywacke,wacke载荷load再作用面(复活面)reactivation surface 藻管迹phycosiphon藻类algae藻团块algae lumps藻团粒algal poliods藻席Algal mats造礁珊瑚hermatypic coral粘度viscosity粘结岩boundstone粘土级clay size粘土矿物clay mineral粘土岩claystone19/21涨潮流flood current帐篷构造tepees障壁barrier障壁岛barrier island障壁沙咀barrier spit障积岩bufflestone沼泽相swamp振荡波痕oscillation ripples 蒸发岩evaporites正常鲕normal Oolits正递变的normally graded正化颗粒orthochems正态分布normal distribution 直观模式visual model指状沙坝bar-finger sand指状钻孔迹trypanites中成岩作用mesodiagenesis 中峰度的mesokertic中砾pebble中值median钟形bell shape众数mode重晶石barite重流层heavy-fluid layer皱饰迹rusophycus蛀木虫迹teredolites铸模mold铸模孔隙moldic porosity铸型cast锥模prod casts坠积物colluvium准层序parasequence准层序组parasequence set浊积石灰岩allodapic limestone浊积岩turbidite浊流turbidity current自生沉积物autochthonous sediments 自生的authigentic20/21自形晶euhedral自悬浮autosuspension组formation组构fabric钻孔迹trypanites21/21。

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

a r X i v :a s t r o -p h /0503573v 1 26 M a r 2005Astronomy &Astrophysics manuscript no.astroph˙scuba˙grb3˙1February 2,2008(DOI:will be inserted by hand later)SCUBA sub-millimeter observations of gamma-ray burstersIII.GRB 030329:The brightest sub-millimeter afterglow to dateI.A.Smith 1,R.P.J.Tilanus 2,N.Tanvir 3,R.A.M.J.Wijers 4,P.Vreeswijk 5,E.Rol 6,and C.Kouveliotou 7,81Department of Physics and Astronomy,Rice University,6100South Main,MS-108,Houston,TX 77005-1892USA e-mail:iansmith@2Joint Astronomy Centre,660N.Aohoku Place,Hilo,HI 96720USA3Centre for Astrophysics Research,University of Hertfordshire,College Lane,Hatfield,Herts AL109AB,UK 4Astronomical Institute ‘Anton Pannekoek’,University of Amsterdam and Center for High-Energy Astrophysics,Kruislaan 403,1098SJ Amsterdam,The Netherlands 5European Southern Observatory,Alonso de C´o rdova 3107,Casilla 19001,Santiago 19,Chile 6Department of Physics and Astronomy,University of Leicester,Leicester LE17RH,UK7NASA Marshall Space Flight Center,SD-50,NSSTC,320Sparkman Drive,Huntsville,AL 35805USA 8Universities Space Research AssociationReceived ;acceptedAbstract.We discuss our ongoing program of Target of Opportunity (ToO)sub-millimeter observations of gamma-ray bursts(GRBs)using the Sub-millimetre Common-User Bolometer Array (SCUBA)on the James Clerk Maxwell Telescope (JCMT).In this paper,we present all the ToO observations of GRB 030329.This was by far the brightest sub-millimeter afterglow seen to date.The flux density at 850µm was approximately constant up to a break that took place ∼7days after the burst.This was consistent with being a jet break.The 850µm results agree with those at longer wavelengths that show a brighter flux ∼7days after the burst,right at the time of the break.No short-lived large-scale brightenings were detected in the sub-millimeter light curve.However,the 850µm light curve may have had a drop <∼16days after the burst.The peak of the afterglow emission was at ∼90GHz in the days before the break in the light curve.A simple modeling is consistent with the spectral indices remaining the same as the afterglow evolved,with the breaks in the spectrum moving to longer wavelengths at later times and the flux at the peak falling.No significant sub-millimeter emission was detected from the host galaxy.Key words.gamma rays:bursts –submillimeter1.IntroductionThe discovery of localized transients in the error boxes of gamma-ray burst (GRB)sources has led to intense multi-wavelength campaigns that have revolutionized our under-standing of these sources.For reviews see Van Paradijs et al.(2000)and M´e sz´a ros (2002).Sub-millimeter observations form a key element of the multiwavelength observations of the afterglow.They provide “clean”measures of the source intensity,una ffected by scin-tillation and extinction.We have therefore been performing Target of Opportunity (ToO)sub-millimeter observations of GRB counterparts using the Sub-millimetre Common-User Bolometer Array (SCUBA)on the James Clerk Maxwell Telescope (JCMT)on Mauna Kea,Hawaii.Including the results shown here,SCUBA has so far per-formed ToO observations of 23bursts (Smith et al.1999,2000,2001,2005;Bloom et al.1998;Galama et al.1999;Kulkarni2I.A.Smith et al.:SCUBA sub-millimeter observations of gamma-ray burstersnearby made it one of the brightest bursts and afterglows ever recorded.GRB030329was an X-ray rich burst,and a distinct, bright,soft X-ray component may have been present during the burst.The X-ray afterglow seen by the Rossi X-Ray Timing Explorer and XMM-Newton was bright(Marshall&Swank 2003;Marshall et al.2003;Tiengo et al.2003,2004).However, the X-ray afterglow light curve was sparsely sampled.It might extrapolate back to match the soft X-ray tail seen during the burst.There appeared to be a break in the decay of the X-ray light curve∼0.5days after the burst,and there may have been aflattening of the X-ray light curve∼40days after the burst.Optical observations simultaneous with the burst did not detect the source,placing a limit of V∼5.5on the reverse shock emission(Torii et al.2003).However,a bright optical afterglow was found,with R∼12.4at67minutes after the burst(e.g.Price et al.2003;Sato et al.2003;Uemura et al. 2003;Urata et al.2004).The optical afterglow was polarized at the0.3−2.5%level, and the polarization was variable on time scales down to hours (Greiner et al.2003b).The radio afterglow was not signifi-cantly polarized(Finkelstein et al.2004;Taylor et al.2005).The optical afterglow decay can broadly be described using a broken power law with a break∼3−8days after the burst (Lipkin et al.2004).However,extensive observations of the op-tical afterglow light curve found a complex evolution of bumps and wiggles superimposed on the overall fading(e.g.Greiner et al.2003b;Smith et al.2003;Matheson et al.2003;Bloom et al. 2004;Lipkin et al.2004).These could have been due to struc-tures in the forward and/or reverse shocks,or repeated energy injection from the central engine(e.g.Sari&M´e sz´a ros2000; Granot et al.2003),or inhomogeneities in the ambient medium that thefireball expanded into(e.g.Berger et al.2000).The in-terpretation of the various features in the light curve is therefore complicated.The radio and millimeter afterglows were the brightest recorded for any burst to date(Berger et al.2003b;Sheth et al. 2003;Kuno et al.2004;Taylor et al.2004;Kohno et al.2005). The radio and millimeter light curves were initially approxi-matelyflat or slowly rising.After this,they broke to a power law Fν∝t−δwith a steep decay ofδ∼2.This is consistent with a jet break(e.g.Sari et al.1999).The light curves for the observations from43to250GHz all broke at∼8days after the burst.In the single broken power law model,this would suggest that the jet break took place at∼8days after the burst rather than at∼3days which is allowed from the optical observa-tions.The peak occurred at increasingly later times for longer wavelengths,for example,it was at∼30days at4.86GHz. These longer wavelengths were below the synchrotron self-absorption frequency(νa)and/or below the typical synchrotron frequency for the lowest energy electron in the power law(νm). Interstellar scintillation was an increasingly important factor below15GHz for the early observations when the source was compact.At late times,the radio light curve showed an achro-maticflattening which is consistent with the blast wave becom-ing transrelativistic(Frail et al.2005).A more complex model for the multiwavelength emission used breaks in the early optical afterglow decay at∼0.25and ∼0.5days(e.g.Sato et al.2003;Burenin et al.2003).Thefirst break could be due to the cooling break frequency(νc)passing through the optical band.The second break would be similar to the one seen in X-rays;since it is achromatic,this suggests it is a jet break.However,an additional source of optical emis-sion would then be required∼1day after the burst.A pos-sible model for the multiwavelength afterglow then uses two jets with different opening angles(Berger et al.2003b;Sheth et al.2003;Lipkin et al.2004;Tiengo et al.2004).The ini-tial part of the optical and X-ray afterglow would come from a narrow ultra-relativistic jet.The late X-ray afterglow and op-tical brightening would come from a wide,mildly relativistic jet,with a jet break∼10days after the burst.The plateau seen at longer wavelengths before the break could then be a con-volution of the fallingflux density(as t−1/3)from the narrow angle jet and a risingflux density(as t1/2)from the wide an-gle jet.Although this model describes the overall trends in the light curves,the multiple bumps and wiggles in the optical light curve still remain to be explained.Beginning∼6days after the burst it became evident that there was an underlying component from the Type Ic su-pernova SN2003dh.The broad spectral features indicated a large expansion velocity similar to the Type Ic hypernovae SN 1998bw and SN1997ef(Stanek et al.2003;Hjorth et al.2003; Kawabata et al.2003;Matheson et al.2003;Mazzali et al. 2003;Kosugi et al.2004).This was direct spectroscopic evi-dence that at least some classical GRBs originate from core-collapse supernovae.The SN2003dh observations also indi-cated that the explosion was aspherical.The host is believed to be a low to moderate metallicity starburst dwarf galaxy with V∼22.7and R>22.5(Fruchter et al.2003;Matheson et al.2003).3.SCUBA observing detailsSCUBA is the sub-millimeter continuum instrument for the JCMT(for a review see Holland et al.1999).The observing, calibration,and reduction techniques used for GRB030329 were the same as has been used on other bursts and are de-scribed in detail in Smith et al.(1999,2001,2005).SCUBA uses two arrays of bolometers to simultaneously observe the same region of sky.The arrays are optimized for operations at450and850µm.Photometry observations are made using a single pixel of the arrays.The other bolometers in the arrays are used to perform a good sky noise subtraction (Archibald et al.2002).During an observation,the secondary is chopped between the source and sky.The term“integration time”always refers to the“on+off”time,including the amount of time spent off-source.An18sec integration thus amounts to a9sec on-source observation time.A typical measurement consists of50integrations of18seconds;we refer to this as a“run.”Each observation of a source in general consists of several such runs,with focus,pointing,and calibration obser-vations in between.Version1.6of SURF1(Jenness&Lightfoot2000)was used to combine the integrations,remove anomalous spikes,flatfieldI.A.Smith et al.:SCUBA sub-millimeter observations of gamma-ray bursters3 Table1.SCUBA850µm(350GHz)afterglow observations ofGRB030329.Burst SCUBA observing times Time since Integrationτ850Afterglow850µmStart Stop burst(days)time(sec)flux density(mJy)a Not corrected for small pointing errorthe array,and subtract the signal from the sky bolometers.The zenith sky opacity was determined using“skydips”in which the sky brightness temperature was measured at a range of ele-vations.This was used to apply atmospheric extinction correc-tions to the observed targetfluxes.At least one standard cali-bration target was observed each night to determine the abso-luteflux of the GRB.A typical integration time of2hours gives an rms∼1.5mJy at850µm.However,the sensitivity depends significantly on the weather and the elevation of the source;since our ToO ob-servations are done on short notice,sometimes these factors are less than ideal.Based on observed variations of the gain factor and signal levels we estimate typical systematic uncertainties in the absoluteflux calibrations of10%at850µm.In general the rms errors of the observations presented here are larger than this uncertainty.The pointing of the JCMT is checked several times during the night to ensure that it is reliable.The pointing accuracy is usually a few arcsec.However,an error in the track model used between2000August25and2003April25resulted in pointing errors that were non-negligible(shifts larger than4′′) for targets with elevations above60◦and over small ranges of azimuth.Since the850µm bolometric pixel has a diffraction limited resolution of14′′,the target remains well within the bolometric pixel.However,there will be an error in theflux that is measured.For GRB030329,the only problem was with a pointing calibration observation of3C273on2003April5. However,this did not significantly affect the result for that day.4.SCUBA observations of GRB030329 Unfortunately,SCUBA was not on the telescope at the time of the burst and it was several days before thefirst SCUBA obser-vation could be made.The weather was not good for most of the observations leading to larger rms uncertainties than nor-mal.Despite these difficulties,GRB030329was still easily de-tected at850µm in the initial observations.Table1summarizes all the850µm SCUBA ToO obser-vations for GRB030329.These expand on and supercede the results in Hoge et al.(2003).The start times are the times when the GRB wasfirst observed on each day,and the stop times are when the last observation of the GRB was completed.The time since the burst uses the mid-point between these start and stop times(although it should be noted that focus,pointing,and calibration observations may have also taken place in be-tween the start and stop times);the time is the elapsed time in the Earth frame,not the rest frame of the host galaxy.The “integration time”for each observation is the“on+off”time; only half of this is spent on-source.The zenith optical depth at 850µm(τ850)is given for the time of the observation.The∼30mJy source initially detected for GRB030329 was by far the brightest SCUBA afterglow seen to date.For all the previous ToO observations performed by SCUBA,the850µmflux density from the afterglow and/or reverse shock was <10mJy.4.1.SCUBA850µm light curveThe850µm light curve is shown in Figure1.This is similar to the observations at longer wavelengths.Theflux density is ap-proximately constant initially,before breaking with a rapid de-cay.There is no evidence for an underlying dusty host galaxy.The weighted mean for thefirst three SCUBA observations was30.7±1.8mJy at850µm.This gives the constant dashed curve shown in Figure1for thefirst week.Unfortunately,we do not have high quality observations during the decay phase,due to the terrible weather at that time. However,the lack of detection17days after the burst allows us to use simple qualitative models to constrain the break time and/or the decay slope.The short dashed curve in Figure1shows a break at8days after the burst and a decay ofδ=2.This would predict aflux density of7.0mJy at16.8days after the burst.This is therefore inconsistent with our non-detection at the3.7σlevel.To maintainδ=2for the whole decay,it would be neces-sary for the break to occur earlier.The solid curve in Figure1 places the break at7days.This is just consistent with our3σlimit at16.8days.To have the jet break occur at later times,it would be neces-sary to have a faster decay.For example,the long dashed curve in Figure1uses a break at8days andδ=3.Such a steep decay(on average)has not been seen in the other longer wave-length observations of this burst.However,as shown in§4.2,it is possible that a break moved through the sub-millimeter band between7and17days after the burst,and this might account for a faster decay.Alternatively,Kohno et al.(2005)found that theflux density at93GHz was approximately constant from 12to17days after the burst before a rapid drop took place be-4I.A.Smith et al.:SCUBA sub-millimeter observations of gamma-rayburstersFig.1.The 850µm light curve for GRB 030329.The same data and models are plotted in the linear-linear (top)and log-log (bottom)versions,except that the last two data points are plot-ted as 3σupper limits in the latter.The initial dashed curve is a constant with F ν=30.7mJy.This breaks to F ν∝t −δ,with the break at 8days and δ=2(short dashed),the break at 8days and δ=3(long dashed),and the break at 7days and δ=2(solid).tween 17and 18days after the burst.It is therefore possible that a similar rapid drop took place <∼16days after the burst at 850µm.The SCUBA flux density at 850µm before thebreak in the light curve is consistent with being constant.Sheth et al.(2003)found a small (<20%)drop in the flux density at 100GHz starting April 2followed by a rise to its mean level on April 5.This is similar to our observation on April 4being lower than on April 3and 5.On April 5,the flux density at 43.3GHz was the highest recorded,and the flux density at the four other longer wavelengths measured was also temporarily higher on that day (Berger et al.2003b).Thus this variability may have been real.If so,it complicates the determination of the actual time for the break.The observations during each SCUBA run show no signif-icant variability.One possible explanation for the fluctuations in the optical light curve is that the central engine had mul-tiple ejections,rather than a single episode (e.g.Granot et al.2003).This would lead to refreshed shocks when the faster moving ejecta caught up with slower moving material.A pre-diction of this model is that there should be short-lived reverse shocks when the ejecta collide producing significant longer wavelength flares.The lack of large-scale variability at 850µm agrees with the results from other observations (Sheth et al.2003),and argues against this model.However,we caution that our data only sparsely sample the light curve.Fig.2.The longer wavelength spectrum for GRB 030329be-fore the break in the light curve.The data points show typi-cal averaged values for ∼6days after the burst.Added to the 850µm (350GHz)and 450µm (670GHz)data from this paper are the results at 4.86,8.46and 15GHz (Berger et al.2003b),23and 43GHz (Berger et al.2003b;Kuno et al.2004),90GHz (Kuno et al.2004),100and 250GHz (Sheth et al.2003).The upper limit at 450µm is 3σ.The solid curve is F ν=0.15ν2up to 15GHz,F ν=13.7ν1/3up to a peak of 61mJy at 90GHz,and F ν=582ν−1/2after the peak.4.2.Evolution of the longer wavelength spectrumDue to the poor observing conditions,the source was not de-tected at 450µm in any of the individual observations.The source was also not detected at 450µm when the results from separate days are combined.The weighted mean for the first three days was −69±33mJy at 450µm.Figure 2illustrates the longer wavelength spectrum before the break in the light curve using the average fluxes for the observations ∼6days after the burst.Since data has been combined from di fferent observatories,and the source may be variable,the points will have uncertainties of up to ∼20%.However,this does give a good qualitative picture for the longer wavelength spectrum at the time before the break.The solid curve in Figure 2shows one qualitative fit to the data.The rise uses F ν∝ν2up to 15GHz and F ν∝ν1/3up to a peak of 61mJy at 90GHz.The fall is F ν∝ν−1/2.The short dashed curve extrapolates the ν1/3curve to highlight the break at longer wavelengths.Spectral indices such as these are commonly found in syn-chrotron fireball models (e.g.Sari et al.1998;Piran 1999;Wijers &Galama 1999;Chevalier &Li 2000;Granot et al.2000;Granot &Sari 2002;Panaitescu &Kumar 2004).For ex-ample,the breaks could be νa =15GHz and νm =90GHz.For the R band (4.2×105GHz)this fit would extrapolate to give 0.9mJy or R ∼16.The source was that bright ∼1dayI.A.Smith et al.:SCUBA sub-millimeter observations of gamma-ray bursters5Fig.3.The longer wavelength spectrum for GRB 030329afterthe break in the light curve.The data points show typical av-eraged values for ∼17days after the burst.Added to the 3σ850µm (350GHz)and 450µm (670GHz)limits from this paper are the results at 4.86,8.46,15,and 43GHz (Berger et al.2003b),23GHz (Kuno et al.2004),93GHz (Kohno et al.2005),100and 250GHz (Sheth et al.2003).The solid curve is F ν=0.294ν2up to 8GHz,F ν=9.41ν1/3up to a peak of 27mJy at 24GHz,and F ν=133ν−1/2after the peak.The short dashed curve shows the fit from Figure 2.after the burst,and was at R ∼18at 6days after the burst.In this case,it would then be necessary to have another break in the spectrum between the sub-millimeter and optical bands.For example,this could be the cooling break νc .Figure 3illustrates the longer wavelength spectrum after the break in the light curve using the average fluxes for the observations ∼17days after the burst.Again,the uncertainties in the points may be up to ∼20%.The solid curve in Figure 3shows one qualitative fit to the data,with F ν∝ν2up to 8GHz,F ν∝ν1/3up to a peak of 27mJy at 24GHz,and a fall with F ν∝ν−1/2.The short dashed curve in Figure 3shows the fit from Figure 2.Qualitatively,the spectral indices have stayed the same,with the two longer wavelength breaks (νa and νm in the example above)moving to longer wavelengths and the flux at the peak falling.However,the model predicts a flux density of 7.1mJy at 850µm at ∼17days after the burst.This is inconsistent with our lack of detection at the 3.7σlevel if we use the limit from April 15,and at the 4.4σlevel if we combine our last two SCUBA observations.One way to reduce the disagreement between the model and the SCUBA observation at ∼17days after the burst is if the third break (νc in the example above)passed through the sub-millimeter band shortly before our last observations.ThisTable 2.Host galaxy fluxes for GRB 030329determined by combining the SCUBA 850µm and 450µm observations taken on the last two days.BurstRedshiftHost flux density (mJy)850µm 450µmcould also help to explain why the decay of the 850µm light curve was faster than δ=2.However,the spectral break would need to be sharp to fully explain the SCUBA limit.Also,in the simple fireball models,νc should remain constant after the jet break.Thus this would require a more complex model,such as having the magnetic energy change with time (e.g.Yost et al.2003).An alternative explanation for the lack of detection at 850µm at ∼17days after the burst is that there could have been a rapid drop in the light curve,as discussed in §4.1.The drop must have taken place earlier than 16.8days after the burst at 850µm.Figure 3shows that the 250GHz point was still on the ν−1/2curve at 16.3days after the burst.The drop took place between 17and 18days after the burst at 93GHz.This in-dicates that if the drop is real,it did not occur at exactly the same time for all wavelengths.Instead,it swept from shorter to longer wavelengths.4.3.SCUBA host galaxy limitsThe results in Table 2show the weighted means for the last two days of SCUBA observations.Although a small amount of af-terglow emission may be present,this contamination should be small.There is no indication for any sub-millimeter emission from the host galaxy.It is not surprising that GRB 030329is not a bright qui-escent sub-millimeter source given that (1)it is at a low red-shift,(2)the host is probably a dwarf galaxy,and (3)the X-ray and optical observations show relatively little excess absorp-tion from dust above the Galactic value (Tiengo et al.2003;Matheson et al.2003).However,given the short observing times and less than ideal observing conditions,the limits on the host galaxy fluxes are less restrictive than for some of the other bursts whose host galaxies have been observed by SCUBA (Barnard et al.2003;Berger et al.2003a;Tanvir et al.2004;Smith et al.2005).The SCUBA host galaxy compilation of Tanvir et al.(2004)focused on hosts with 850µm rms flux densities <1.4mJy.Since our rms is slightly higher than this,and the redshift of GRB 030329is <2,the lack of detection of a sub-millimeter host galaxy does not shed any new light on the question of whether GRBs are closely linked to the most luminous dusty star-forming galaxies.5.SummaryGRB 030329was by far the brightest sub-millimeter afterglow seen to date.Despite our limited observations,we were still able to constrain the break and /or decay light curve.Assuming6I.A.Smith et al.:SCUBA sub-millimeter observations of gamma-ray burstersthat the decay is not steeper thanδ=2,the break must have taken place∼7days after the burst at850µm.This jet break could have been at a slightly later time if another break(per-haps the cooling break)passed through the sub-millimeter band between7and17days after the burst,or if there was a drop in the850µmflux at<∼16days after the burst.No short-lived large-scale brightenings were detected in the sub-millimeter light curve.Theflux density at850µm was con-sistent with being constant up to the jet break.However,the 850µm results also agree with those at longer wavelengths that show a slightly brighterflux∼7days after the burst,right at the time of the break.The peak of the afterglow emission was at∼90GHz in the days before the break in the light curve.The spectral indices were qualitatively given by+2and+1/3up to the peak and −1/2beyond it.A simple modeling is consistent with the spec-tral indices remaining the same as the afterglow evolved,with the spectral breaks moving to longer wavelengths at later times and theflux at the peak falling.GRB030329did not have any significant sub-millimeter contribution from the host galaxy.The observations of GRB030329show that for bright af-terglows it is important to monitor theflux evolution carefully to look forfluctuations in the light curve and to determine the breaks.Observations of new bursts are continuing to produce surprises,and there is much left to learn about GRB afterglows and host galaxies.To obtain a complete picture of their nature will require the careful study of many bursts to expand our sam-ple.Sub-millimeter observations with a∼mJy sensitivity are a key component to the multi-wavelength coverage.To this end, our program of ToO observations using SCUBA is ongoing. Acknowledgements.The James Clerk Maxwell Telescope is operated by The Joint Astronomy Centre on behalf of the Particle Physics and Astronomy Research Council of the United Kingdom,the Netherlands Organisation for Scientific Research,and the National Research Council of Canada.We thank the JCMT Director Gary Davis for authorizing the ToO observations.We are indebted to all the observers whose time was dis-placed by these observations,and acknowledge the dedicated efforts of the JCMT telescope operators Jim Hoge and Jonathan Kemp for their valuable assistance with the observations.We are grateful to Scott Barthelmy and Paul Butterworth for maintaining the GRB Coordinates Network(GCN),and to the other ground-based observers for the rapid dissemination of their burst re-sults.We acknowledge the data analysis facilities provided by the Starlink Project which is run by CCLRC on behalf of PPARC.The work at Rice University was supported in part by AFOSR/NSF grant number NSF AST-0123487.ReferencesArchibald,E.N.,Jenness,T.,Holland,W.S.,et al.2002,MNRAS, 336,1Barnard,V.E.,Blain,A.W.,Tanvir,N.R.,et al.2003,MNRAS,338, 1Berger,E.,Sari,R.,&Frail,D.A.2000,ApJ,545,56Berger,E.,Cowie,L.L.,Kulkarni,S.R.,et al.2003a,ApJ,588,99 Berger,E.,Kulkarni,S.R.,Pooley,G.,et al.2003b,Nature,426,154Bloom,J.S.,Frail,D.A.,Kulkarni,S.R.,et al.,1998,ApJ,508,L21 Bloom,J.S.,van Dokkum,P.G.,Bailyn,C.D.,et al.2004,ApJ,127, 252Burenin,R. A.,Sunyaev,R. A.,Pavlinksy,M.N.,et al.2003, Astronomy Letters,29,573Chevalier,R.A.,&Li,Z.-Y.2000,ApJ,536,195Finkelstein,A.M.,Ipatov,A.V.,Gnedin,Y.N.,et al.2004,Astronomy Letters,30,368Frail,D.A.,Berger,E.,Galama,T.,et al.2000,ApJ,538,L129 Frail,D.A.,Bertoldi,F.,Moriarty-Schieven,G.H.,et al.2002,ApJ, 565,829Frail,D.A.,Yost,S.A.,Berger,E.,et al.2003,ApJ,590,992 Frail,D.A.,Soderberg,A.M.,Kulkarni,S.R.,et al.2005,ApJ,619, 994Fruchter,A.,Levan,A.,Hook,R.,et al.2003,GCN2243 Galama,T.J.,Briggs,M.S.,Wijers,R.A.M.J.,et al.1999,Nature, 398,394Granot,J.,&Sari,R.2002,ApJ,568,820Granot,J.,Piran,T.,&Sari,R.2000,ApJ,534,L163Granot,J.,Nakar,E.,&Piran,T.2003,Nature,426,138Greiner,J.,Peimbert,M.,Estaban,C.,et al.2003a,GCN2020 Greiner,J.,Klose,S.,Reinsch,K.,et al.2003b,Nature,426,157 Hjorth,J.,Sollerman,J.,Møller,P.,et al.2003,Nature,423,847 Hoge,J.C.,Meijerink,R.,Tilanus,R.P.J.,&Smith,I.A.2003,GCN 2088Holland,W.S.,Robson,E.I.,Gear,W.K.,et al.1999,MNRAS,303, 659Jenness,T.,&Lightfoot,J.F.2000,Starlink User Note216,Starlink Project,CLRCKawabata,K.S.,Deng,J.,Wang,L.,et al.2003,ApJ,593,L19 Kohno,K.,Tosaki,T.,Okuda,T.,et al.2005,PASJ,57,147 Kosugi,G.,Mizumoto,Y.,Kawai,N.,et al.2004,PASJ,56,61 Kulkarni,S.R.,Frail,D.A.,Sari,R.,et al.1999,ApJ,522,L97 Kuno,N.,Sato,N.,Nakanishi,H.,et al.2004,PASJ,56,L1 Lipkin,Y.M.,Ofek,E.O.,Gal-Yam,A.,et al.2004,ApJ,606,381 Marshall,F.E.&Swank,J.H.2003,GCN1996Marshall,F.E.,Markwardt,C.,&Swank,J.H.2003,GCN2052 Matheson,T.,Garnavich,P.M.,Stanek,K.Z.,et al.2003,ApJ,599, 394Mazzali,P.A.,Deng,J.,Tominaga,N.,et al.2003,ApJ,599,L95M´e sz´a ros,P.2002,ARA&A,40,137Panaitescu,A.,&Kumar,P.2004,MNRAS,350,213Piran,T.1999,Physics Reports,314,575Price,P.A.,Fox,D.W.,Kulkarni,S.R.,et al.2003,Nature,423,844 Sari,R.,&M´e sz´a ros,P.2000,ApJ,535,L33Sari,R.,Piran,T.,&Narayan,R.1998,ApJ,497,L17Sari,R.,Piran,T.,&Halpern,J.P.1999,ApJ,519,L17Sato,R.,Kawai,N.,Suzuki,M.,et al.2003,ApJ,599,L9 Schnoor,P.W.,Welch,D.L.,Fishman,G.J.,Price,A.2003,GCN 2176Sheth,K.,Frail,D.A.,White,S.,et al.2003,ApJ,595,L33 Smith,D.A.,Rykoff,E.S.,Akerlof,C.W.,et al.2003,ApJ,596, L151Smith,I.A.,Tilanus,R.P.J.,Van Paradijs,J.,et al.1999,A&A,347, 92Smith,I.A.,Van Paradijs,J.,Tilanus,R.P.J.,et al.2000,in Gamma-Ray Bursts:5th Huntsville Symposium,ed.R.M.Kippen,R.S.Mallozzi,&G.J.Fishman(New York:AIP),326Smith,I.A.,Tilanus,R.P.J.,Wijers,R.A.M.J.,et al.2001,A&A, 380,81Smith,I.A.,Tilanus,R.P.J.,Tanvir,N.,et al.2005,A&A,submitted Stanek,K.Z.,Matheson,T.,Garnavich,P.M.,et al.2003,ApJ,591, L17。

相关文档
最新文档