Instruction: Write an abstract for the following acadamic paper.

Input:
Paper: we have recorded the flight of a fly during take off and landing using digital high speed photography . it is shown that the dynamics of flexible wings are different for these procedures . during this observation fly flew freely in a big box and it was not tethered .
Abstract: in this fluid dynamics video , we demonstrated take off and landing of a fly . the deformation of wings is in focus in this video .

Paper: the standard model ( sm ) is a very successful effective theory . over the last few years it has emerged that the ckm mechanism is likely the dominant mechanism and new physics ( np ) beyond the sm would appear as corrections . if np were observed at the lhc its flavour sector must be studied . here we discuss the plans for an upgraded lhcb experiment . in section  [ sec : lhcb ] we present the lhcb experiment and physics programme , and motivate running at higher luminosities . the physics reach of a 100  data sample is discussed in section  [ sec : physics ] for selected key measurements : the @xmath4 mixing phase  from the tree decay @xmath7 , the @xmath8 penguin decay @xmath9 , as well as the ckm angle  @xmath6 . in section  [ sec : detector ] we present the plans for upgrading the lhcb trigger and detectors . the report ends with conclusions . the goal of the lhcb experiment is to perform precise measurements of cp violation and rare decays of @xmath10 mesons . the lhcb detector is a forward arm spectrometer which is equipped with state - of - the - art vertexing capabilities , and is designed to trigger efficiently on muons , electrons and photons as well as hadronic final states . two ring imaging cherenkov ( rich ) detectors allow excellent charged particle identification . the trigger and flavour tagging of lhcb are discussed elsewhere in these proceedings  @xcite . the lhcb experiment will start taking data in 2007 at the large hadron collider ( lhc ) at cern . more details about the current status of lhcb and its first year data - taking plans are described in  @xcite . the lhcb experiment is designed to operate at a luminosity @xmath11 . this choice maximises the number of bunch crossings with exactly one @xmath12 interaction and reduces the radiation damage to the vertex detector . this luminosity is a factor of @xmath13 below the lhc design value which is achievable by using a tunable amplitude function , @xmath14 , at the lhcb interaction point . therefore this luminosity will very likely be reached during the first lhc physics run . in a nominal year of @xmath15 , a data sample of @xmath16 will be collected by lhcb .    [ cols="^,<,^,^,^",options="header " , ]      +    the large cross section of @xmath17b for the production of @xmath18-quark pairs at the lhc will allow lhcb to collect much larger data samples of selected @xmath10 meson decays than previously available . the sensitivities of these data - sets to cp violating asymmetries and to other physics observables have been studied with simulated events . in table [ tab : sensitivities ] we present a selection of signal yields @xmath19 , background to signal ratios @xmath20 and precisions on ckm angles @xmath6 , @xmath21 , @xmath22 and @xmath5 and other observables . detailed presentations on many of the channels are presented elsewhere in these proceedings  @xcite . the plan is to operate the lhcb experiment for five years at the design luminosity and to collect a data sample of @xmath23 to @xmath24 . a major goal is to exploit the physics in the @xmath25 system . this includes the observation of cp violation in @xmath25 mesons and precision measurements of the mass difference between the @xmath4 mass eigenstates , @xmath26 , and the lifetime difference @xmath27 . other main aims are to improve the error on the ckm angle  @xmath6 by a factor of five , to probe np in rare @xmath10 meson decays with electroweak , radiative and hadronic penguin modes , and to make the first observation of the very rare decay @xmath28 . after the first five years of operation , the precision of many lhcb physics results will remain limited by the statistical error of the collected data . the following questions arise : is lhcb exploiting the full potential for @xmath10 physics at hadron colliders and is there a science case for collecting even larger data samples ? note that lhcb is the only dedicated heavy flavour experiment approved to run after 2010 . in the remainder of this report we will try to answer these questions . the lhcb experiment has commenced studying the feasibility of upgrading the detector such that it can operate at a luminosity @xmath29 , which is ten times larger than the design luminosity . this upgrade would allow lhcb to collect a data sample of about @xmath30 during five years of running . this increased luminosity is achievable by decreasing @xmath14 at the lhcb interaction point . it does not require a lhc luminosity upgrade ( super - lhc ) as the lhc design luminosity is @xmath31 although it could operate at super - lhc . thus an upgrade of lhcb could be implemented as early as 2013 . the number of interactions per beam crossing will increase to @xmath32 which will require improvements to the lhcb sub - detectors and trigger . flavour - changing neutral currents ( fcncs ) are particularly sensitive to np . new particles appear virtually in loop diagrams which lead to deviations from sm predictions . the interference in @xmath4 mixing and decay to cp eigenstates is sensitive to the cp violating weak mixing phase @xmath5 . this can be studied in the time - dependent asymmetry of flavour - tagged  decays . the sm prediction for @xmath5 is very small : @xmath34rad where @xmath35 and @xmath36 are the usual wolfenstein parameters of the ckm matrix  @xcite . hence @xmath5 is a very sensitive probe of np and will provide a stringent test for non - minimal flavour violation ( nmfv )  @xcite . currently the only direct  result is from d0 who measure @xmath37rad in the untagged decay time spectrum of  decays  @xcite . using a full simulation , the lhcb experiment expects to collect 131k  decays with a 2  data sample  @xcite . the expected precision on  is estimated with many toy monte carlo experiments based on the yields and resolutions obtained in the full simulation . we obtain @xmath38rad where the  mass difference has been set to @xmath39 . new physics in @xmath4 mixing can be parameterised as @xmath40 where the sm phase of @xmath41 is equal to @xmath42 and @xmath43 and @xmath44 are , respectively , the amplitude and the phase of the np contribution . in fig . [ fig : hsss ] we show the expected sensitivity for @xmath43 and @xmath44 for lhcb with 2  of data  @xcite .    with five years of data , lhcb expects to achieve a precision of @xmath45rad which corresponds to about a @xmath46 significance for the sm value of . this precision is expected to be statistically limited , with systematic errors likely to be much smaller . the precision of the sm prediction for  is better than 0.01rad .    a significantly larger data - set would allow lhcb to probe new physics in @xmath10 mesons at an unprecedented level of a few percent . here we present estimates for sensitivities with a 100  data sample which are based on scaling with luminosity . potential trigger efficiency improvements are not included . while these estimates have large uncertainties , these are extremely useful to motivate simulation studies with the aim to improve or at least maintain trigger and selection efficiencies at ten times higher luminosities . an upgrade of lhcb has the potential to measure the sm value of   with @xmath47 precision ( @xmath48rad ) in  decays . new physics can be probed by studying fcncs in hadronic @xmath49 transitions . one approach is to compare the time - dependent cp asymmetry in a hadronic penguin loop decay , where unknown massive particles could enter with a decay based on a tree diagram which generally is insensitive to np and which has the same weak phase . the b - factories measure the cp asymmetry @xmath50 in the penguin decay . a value for @xmath50 which is different from @xmath51 measured in  would signal physics beyond the sm . within the current available precision , all @xmath50 measurements are in reasonable agreement with the sm , but all central values are lower than expected . a naive average indicates a @xmath52 discrepancy  @xcite . the above approach can also be applied to  mesons and this will be exploited by lhcb . within the sm the weak mixing phase  is expected to be almost the same when comparing the time - dependent cp asymmetry of the hadronic penguin decay  with the tree decay . due to a cancellation of the @xmath25 mixing and decay phase , the sm prediction for  is actually very close to zero  @xcite . thus any measurement of @xmath53 is a clear signal for nmfv . lhcb expects to collect 1.2k  decays in 2  of data  @xcite . by comparing with the  yield we estimate a sensitivity of @xmath54rad for 10  of data this  precision is statistically limited . scaling the sensitivity up to a data sample of 100   we estimate a precision of @xmath55rad . this sensitivity presents a very precise and exciting np probe . a first - level detached vertex trigger is required to reconstruct this hadronic decay mode . a new lhcb simulation study on   can be found in  @xcite . lhcb will perform direct measurements of the ckm angle  @xmath6 using two interfering tree processes in neutral and charged  decays . the interference arises due to decays which are common to @xmath56 and @xmath57 mesons such as @xmath58 ( dalitz decay  @xcite ) and @xmath59 ( ads and glw  @xcite ) . the current combined direct @xmath6 measurements at the @xmath10 factories have large errors of about 30  @xcite . the expected lhcb @xmath6 sensitivities for 2  of data have been estimated for several of the above modes  @xcite ; we obtain @xmath60 . in addition , lhcb will use the decay   to measure @xmath6 with @xmath61 in 2  . by combining these measurements , lhcb will achieve a precision of @xmath62 in 2  . the theoretical error on the sm prediction is very small . note that some of the @xmath6 measurements will have a correlated systematic error , for example there is a dalitz model dependence . we take this into account when scaling to a 100  data sample and estimate that a precision of @xmath63 should be achievable . this would allow a very precise determination of the apex of the unitarity triangle . we presented a selection of key measurements that could be performed with an upgrade of the lhcb experiment , namely precision measurements of the weak mixing phase  in tree decays and @xmath8 transitions , and of the ckm angle  @xmath6 . many more measurements would also benefit from a ten - fold increase in statistics . with 100  of data , lhcb would be able measure @xmath64 in  with an estimated precision of @xmath65 . an estimated yield of 44k  events per year would allow probing right - handed currents by analysing forward - backward asymmetries ( @xmath66 ) and transversity angles . in addition , the lhcb upgrade would allow the probeing of np in charm decays and to study lepton flavour violation in @xmath67 decays .    in fig . [ fig : comparison ] we present a comparison of the expected sensitivities for the lhcb upgrade , scaled from 10 to 100  , with sensitivities at a super - b factory based on a @xmath68 data sample  @xcite . both approaches have comparable sensitivities for many measurements such as @xmath6 and @xmath69 . lhcb will be limited in reconstructing decays with ( more than one ) neutral final state particles ( photon , @xmath70 or @xmath36 ) . however only an upgrade of the lhcb experiment will allow an ultimate exploitation of  meson decays . while these estimates need to be confirmed with simulations at high luminosities , an excellent science case can be written for an lhcb upgrade . the trigger for the lhcb experiment comprises a first level trigger ( l0 ) implemented in hardware and a higher level trigger ( hlt ) running on a large cpu farm . the l0 trigger operates at 40 mhz and selects high transverse momentum objects in the electromagnetic and hadronic calorimeters ( @xmath71 , @xmath6 , charged hadrons ) and muons in the muon system ( @xmath72 , @xmath73 ) . pile - up events are vetoed . the output rate is limited to 1.1 mhz with @xmath74 latency . simulations show that the l0 muon trigger efficiency for reconstructable events is around 90% and that the output rate scales with luminosity up to @xmath75 . at the design luminosity , the l0 hadron trigger efficiencies are about 40% , and the output rate increases only slightly with luminosity . the vertex detector ( velo ) sensors  @xcite undergo radiation damage and it is expected that these will need to be replaced when 6 to 8  of luminosity has been collected . this illustrates that the existing trigger and at least some of the sub - detector systems will not allow operating the lhcb experiment at ten times the design luminosity . we have commenced studies which investigate how to upgrade the lhcb detector and triggers such that the experiment can operate at luminosities @xmath76 . two scenarios are under consideration . a step - by - step approach would foresee to replace the velo sensors with improved radiation - hard detectors which will be read out at 40 mhz . this would allow the addition of a first - level detached vertex trigger at l0 . the trigger tracker in front of the lhcb dipole magnet could also be included . this trigger could be implemented in fpgas , and it would need to be very fast and run within the l0 latency time . other lhcb detector systems will also need to be upgraded when increasing the luminosity , either due to unacceptably large occupancies or radiation effects . we are currently studying the luminosity limit of each lhcb sub - detector system . under discussion is the replacement of the central region of the rich1 photon detectors , an increase ( decrease ) of the inner ( outer ) tracker area , and a replacement of the inner region of the electromagnetic calorimeter .    another approach would be to change the read - out of all lhcb sub - detectors to 40  mhz . this has clear advantages as it would allow the implementation of a l0 displaced vertex trigger in a cpu farm . in fact all trigger decisions would be software - based which allows flexibility . however this approach requires a redesign of the front - end electronics which has implications for all sub - detector systems . besides the velo , all silicon sensors of the trigger and inner tracker and the rich photon detectors would need to be replaced . the r&d for the new front - end electronics will be able to profit from the developments of atlas and cms for super - lhc luminosities . r&d efforts have started on technologies for radiation - hard vertex detectors that will be able to operate in the radiation environments of the lhc and lhcb upgrade . the detector sensors will need to be able to operate at radiation doses of about @xmath77 . initial studies of czochralski and @xmath78-on-@xmath79 sensors irradiated up to @xmath80 24 gev protons/@xmath81 are promising and show that the charge collection efficiencies saturate at acceptable bias voltages  @xcite . three - dimensional sensors are another alternative that could be investigated .    two different vertex - detector geometries are envisaged . one is to shorten the strips , the other is to use pixels . removing the rf foil that separates the velo sensors from the primary beam - pipe vacuum would reduce the radiation length before the first measurement by 3% and improve the proper time resolution of @xmath10 meson decays . a preliminary study uses  decays simulated at a luminosity of @xmath82 . events with large numbers of interactions are employed to simulate larger effective luminosities up to @xmath83 . in a first step , the existing hlt velo trigger algorithm is added without modifications into the l0 trigger . we find that the minimum bias rate keeps rising with luminosity and saturates the bandwidth well below our target luminosity . a better approach is to combine the l0 trigger with a detached vertex trigger and to read out at a 40mhz rate . we require a transverse energy @xmath84 from the hadron trigger and combine this with a matched track that has a transverse momentum @xmath85 gev/@xmath86 and an impact parameter @xmath87 . in this combined trigger the minimum bias rate does not depend strongly on the luminosity and the triggered event yield scales linearly with the luminosity . in addition , the total trigger efficiency is 60% larger when compared with the existing baseline . we have presented a plan to upgrade the lhcb experiment after about five years of running . this programme is motivated by new physics beyond the standard model that can be studied with a ten times larger data sample . the lhcb upgrade would allow a precise measurement the weak mixing phase   at its standard model value and to probe new physics in @xmath10 meson decays at an unprecedented level of a few percent or better . initial estimates show that an excellent science case can be developed . the plan is to upgrade the lhcb detector such that it can operate at ten times the design luminosity , i.e. at @xmath76 . the upgraded experiment can run with or without an lhc luminosity upgrade . the lhcb upgrade will require a first - level detached vertex trigger . to achieve this , the vertex detector sensors will need to be replaced and many other sub - systems will need to be improved . feasibility studies for physics and trigger , and r&d for detectors and front - end electronics are now required . e. rodrigues , these proceedings . h. ruiz , these proceedings . l garrido , these proceedings . g. corti , these proceedings . a. carbone , these proceedings . y. xie , these proceedings . p. robbe , these proceedings . n. magini , these proceedings . m. smizanska , these proceedings . s. de capua , these proceedings . r. muresan , these proceedings . l. wolfenstein , phys . lett . * 51 * , ( 1983 ) 1945 . g. dambrosio _ et al . * b645 * ( 2002 ) 155 . _ , hep - ex/0701012 . l. fernandez , at physics at lhc , cracow , poland 2006 , cern / lhcb/2006 - 047 ( 2006 ) . z. ligeti _ * 97 * ( 2006 ) 101801 . heavy flavour averaging group , e. barberio _ et al . _ , hep - ex/0603003 , updated ( ichep06 ) : http://www.slac.stanford.edu/xorg/hfag/ m. raidal , phys . * 89 * ( 2002 ) 231803 . r antunes nobrega _ et al . _ , lhcb tdr 9 , cern / lhcc/2003 - 30 ( 2003 ) . j. libby , at ckm 2006 , nagoya ( 2006 ) . a. giri _ et al . , * d68 * ( 2003 ) 054018 . d. atwood , i. dunietz and a. soni , phys . * 78 * ( 1997 ) 3257 . m. gronau , phys . lett . , * b557 * ( 2003 ) 198 . m. hazumi , presented at flavour in the lhc - era workshop , cern ( 2005 ) . c. parkes , nucl . inst . and meth . * a569 * ( 2006 ) 115 .
Abstract:

Answer: the lhcb experiment will operate for about five years at a luminosity of @xmath0 and plans are to accumulate a data sample of @xmath1 here we present the physics programme and detector design for a future high luminosity phase of the lhcb experiment . an upgraded lhcb experiment would operate at ten times the design luminosity , i.e. at @xmath2 and aims to collect a data sample of @xmath3 over five years . this programme would allow the probe of new physics at an unprecedented level . key measurements include the @xmath4 mixing phase @xmath5 in  and  decays with a significant sensitivity to the small standard model prediction and a very precise measurement of the ckm angle  @xmath6 in tree diagram decays . initial studies of the modified lhcb trigger and detectors are presented . the upgraded lhcb experiment can run with or without an lhc luminosity upgrade .


Instruction: Write an abstract for the following acadamic paper.

Input:
Paper: we have recorded the flight of a fly during take off and landing using digital high speed photography . it is shown that the dynamics of flexible wings are different for these procedures . during this observation fly flew freely in a big box and it was not tethered .
Abstract: in this fluid dynamics video , we demonstrated take off and landing of a fly . the deformation of wings is in focus in this video .

Paper: in many scientific applications , we need to solve lower triangular problems and upper triangular problems , such as incomplete lu ( ilu ) preconditioners , domain decomposition preconditioners and gauss - seidel smoothers for algebraic multigrid solvers @xcite . the algorithms for these problems are serial in nature and difficult to parallelize on parallel platforms . gpu is one of these parallel devices , which is powerful in float point calculation and is over 10 times faster than latest cpu . recently , gpu has been popular in various numerical scientific applications . it is efficient for vector operations . researchers have developed linear solvers for gpu devices @xcite . however , the development of efficient parallel triangular solvers for gpu is still challenging @xcite . klie et al . ( 2011 ) investigated a triangular solver @xcite . they developed a level schedule method and a speedup of two was obtained . naumov ( 2011 ) from the nvidia company also developed parallel triangular solvers @xcite . he developed new parallel triangular solvers based on a graph analysis . the average speedup was also around two .    in this paper , we introduce our work on speeding triangular solvers . a new matrix format , hec ( hybrid ell and csr ) , is developed . a hec matrix includes two matrices , an ell matrix and a csr matrix . the ell part is in column - major order and is designed the way to increase the effective bandwidth of nvidia gpu . for the csr matrix , each row contains at least one non - zero element . this design of the csr part reduces the complexity of the solution of triangular systems . in addition , parallel algorithms for solving the triangular systems are developed . the algorithms are motivated by the level schedule method described in @xcite . our parallel triangular solvers can be sped up to seven times faster . based on these modified algorithms , ilu(k ) , ilut and domain decomposition ( restricted additive schwarz ) preconditioners are developed . numerical experiments are performed . these experiments show that we can speed linear solvers around ten times faster . the layout is as follows . in  2 , a new matrix format and algorithms for lower triangular problems and upper triangular problems are introduced . in  3 , parallel triangular solvers are employed to develop ilu preconditioners and domain decomposition preconditioners . in  4 , numerical tests are performed . at the end , conclusions are presented . for the most commonly used preconditioner ilu , the following triangular systems need to be solved : @xmath0 where @xmath1 and @xmath2 are a lower triangular matrix and an upper triangular matrix , respectively , @xmath3 is the right - hand side vector , @xmath4 is the unknown to be solved for , and @xmath5 is the intermediate unknown . the lower triangular problem , @xmath6 , is solved first , and then , by solving the upper triangular problem , @xmath7 , we can obtain the result @xmath4 . in this paper , we always assume that each row of @xmath1 and @xmath2 is sorted in ascending order according to their column indices . the matrix format we develop is denoted by hec ( hybrid ell and csr ) @xcite . its basic structure is demonstrated by figure  [ fig1 ] . an hec matrix contains two submatrices : an ell matrix , which was introduced in ellpack  @xcite , and a csr ( compressed sparse row ) matrix . the ell matrix has two submatrices , a column - indices matrix and a non - zeros matrix . the length of each row in these two matrices is the same . the ell matrix is in column - major order and is aligned when being stored on gpu . note that the data access pattern of global memory for nvidia tesla gpu is coalesced @xcite so the data access speed for the ell matrix is high . a disadvantage of the ell format is that it may waste memory if one row has too many non - zeros . in this paper , a csr submatrix is added to overcome this problem . a csr matrix contains three arrays , the first one for the offset of each row , the second one for column indices and the last one for non - zeros . for our hec format matrix , we store the regular part of a given triangular matrix @xmath1 in the ell part and the irregular part in the csr part . when we store the lower triangular matrix , each row of the csr matrix has at least one element , which is the diagonal element in the triangular matrix @xmath1 .          in this section , we introduce our parallel lower triangular solver to solve @xmath8 the solver we develop is based on the level schedule method @xcite . the idea is to group unknowns @xmath9 into different levels so that all unknowns within the same level can be computed simultaneously @xcite . for the lower triangular problem , the level of @xmath9 is defined as @xmath10 where @xmath11 is the @xmath12th entry of @xmath1 , @xmath13 is zero initially and @xmath14 is the number of rows . we define @xmath15 , which is the union of all unknowns that have the same level . here we assume that each set @xmath16 is sorted in ascending order according to the indices of the unknowns belonging to @xmath16 . define @xmath17 the number of unknowns in set @xmath16 and @xmath18 the number of levels . now , a map @xmath19 can be defined as follows : @xmath20 where @xmath21 is the position @xmath9 in the set @xmath22 when @xmath9 belongs to @xmath22 . with the help of the map , @xmath19 , we can reorder the triangular matrix @xmath1 to @xmath23 , where @xmath11 in @xmath1 is transformed to @xmath24 in @xmath23 . @xmath23 is still a lower triangular matrix . from this map , we find that if @xmath9 is next to @xmath25 in some set @xmath22 , the @xmath26th and @xmath27th rows of @xmath1 are next to each other in @xmath23 after reordering . it means that @xmath1 is reordered level by level , which implies that memory access in matrix @xmath23 is less irregular than that in matrix @xmath1 . therefore , @xmath23 has higher performance compared to @xmath1 when we solve a lower triangular problem . the whole algorithm is described in two steps , the preprocessing step and the solution step , respectively . the preprocessing step is described in algorithm [ alg1 ] . in this step , the level of each unknown is calculated first . according to these levels , a map between @xmath1 and @xmath23 can be set up according to equation ( [ equ - map ] ) . then the matrix @xmath1 is reordered . we should mention that @xmath1 can be stored in any kind of matrix format . a general format is csr . at the end , @xmath23 is converted to the hec format and as we discussed above each row of the csr part has at least one element .    . @xmath28 ;    start = level(i ) ; end = level(i + 1 ) - 1 ; solve the @xmath27th row ;    @xmath29 ;    the second step is to solve the lower triangular problem . this step is described in algorithm [ alg2 ] , where @xmath30 is the start row position of level @xmath26 . first , the right - hand side @xmath3 is permutated according to the map @xmath19 we computed . then the triangular problem is solved level by level and the solution in the same level is simultaneous . each thread is responsible for one row . at the end , the final solution is obtained by a permutation .    to solve the upper triangular problem @xmath31 , we introduce the following transferring map :    @xmath32 where @xmath14 is the number of rows in @xmath2 . with this map , the upper triangular problem is transferred to a lower triangular problem , and the lower triangular solver is called to solve the problem . the ilu factorization for a sparse matrix @xmath33 computes a sparse lower triangular matrix @xmath1 and a sparse upper triangular matrix @xmath2 . if no fill - in is allowed , we obtain the so - called ilu(0 ) preconditioner . if fill - in is allowed , we obtain the ilu(k ) preconditioner , where @xmath34 is the fill - in level . another method is ilut(p , tol ) , which drops entries based on the numerical values @xmath35 of the fill - in elements and the maximal number @xmath36 of non - zeros in each row @xcite . the performance of parallel triangular solvers for ilu(k ) and ilut is dominated by original problems . in this paper , we implement block ilu(k ) and block ilut preconditioners . when the number of blocks is large enough , we will have sufficient parallel performance . if the matrix is not well partitioned and the number of blocks is too large , the effect of ilu(k ) and ilut will be weakened . when we partition a matrix , the graph library metis @xcite is employed . as we discussed above , the effect of ilu preconditioners is weakened if the number of blocks is too large . the domain decomposition preconditioner is implemented , which was developed by cai et al . @xcite . the domain decomposition preconditioner we implement is the so - called restricted additive schwarz ( ras ) method . overlap is introduced , and , therefore , this preconditioner is not as sensitive as block ilu preconditioners . it can lead to good parallel performance and good preconditioning . in this paper , we treat the original matrix as an undirected graph and this graph is partitioned by metis @xcite . the subdomain can be extended according to the topology of the graph . then each smaller problem can be solved by ilu(k ) or ilut . in this paper , we use ilu(0 ) . assume that the domain is decomposed into @xmath37 subdomains , and we have @xmath37 smaller problems , @xmath38 , @xmath39 , @xmath40 , and @xmath41 . we do not solve these problems one by one but we treat them as one bigger problem : @xmath42 each @xmath43 is factorized by the ilu method , where we have @xmath44 then equation ( [ ddm ] ) is solved by our triangular solvers . in this section , numerical experiments are presented , which are performed on our workstation with intel xeon x5570 cpus and nvidia tesla c2050/c2070 gpus . the operating system is centos 6.2 x86_64 with cuda toolkit 4.1 and gcc 4.4 . all cpu codes are compiled with -o2 option . the data type of a float point number is double . the linear solver is gmres(20 ) . bilu , bilut and ras denote the block ilu(0 ) , block ilut and restricted additive schwarz preconditioners , respectively . the matrix used in this example is from a three - dimensional poisson equation . the dimension is 1,000,000 and the number of non - zeros is 6,940,000 . the ilu preconditioners we use in this example are block ilu(0 ) and block ilut(7 , 0.1 ) . the performance data is collected in table [ tab - ex1 ] . .performance of the matrix from the poisson equation [ cols= " < , < , < , < , < , < , < , < " , ]     for the block ilu(0 ) , we can speed it over 3 times faster . when the number of blocks increases , our algorithm has better speedup . the whole solving part is sped around 8 times faster . bilut is a special preconditioner , since it is obtained according to the values of @xmath1 and @xmath2 ; sometimes , the sparse patterns of @xmath1 and @xmath2 are more irregular than those in bilu . from table [ tab - ex1 ] , the speedups are a little lower compared to bilu . however , bilut reflects the real data dependence , and its performance is better in general . this is demonstrated by table [ tab - ex1 ] . the cpu version bilut always takes less time than the cpu bilu . but for the gpu versions , their performance is similar . the bilut is sped around 2 times faster . the whole solving part is sped around 6 times faster . the ras preconditioner is always stable . it also has better speedup than bilu and bilut . the triangular solver is sped around 6 times faster and the average speedup is around 9 . matrix atmosmodd is taken from the university of florida sparse matrix collection @xcite and is derived from a computational fluid dynamics ( cfd ) problem . the dimension of atmosmodd is 1,270,432 and it has 8,814,880 non - zeros . the ilu preconditioners we use in this example are block ilu(0 ) and block ilut(7 , 0.01 ) . the performance data is collected in table [ tab - ex2 ] . llllllll pre & blocks & cpu ( s ) & gpu ( s ) & speedup & pre cpu ( s)&pre gpu ( s ) & speedup +    bilu & 16 & 20.61 & 2.63 & 7.79 & 0.0248 & 0.0072 & 3.45 + bilu & 128 & 23.94 & 2.80 & 8.50 & 0.0244 & 0.0072 & 3.40 + bilu & 512 & 24.13 & 2.72 & 8.82 & 0.0241 & 0.0070 & 3.46 + bilut & 16 & 14.70 & 2.37 & 6.16 & 0.028669 & 0.0114 & 2.51 + bilut & 128 & 16.58 & 2.43 & 6.79 & 0.028380 & 0.0100 & 2.84 + bilut & 512 & 19.91 & 2.64 & 7.50 & 0.027945 & 0.0113 & 2.47 + ras & 256 & 23.45 & 2.66 & 8.75 & 0.0428 & 0.0072 & 5.96 + ras & 2048 & 25.75 & 3.28 & 7.81 & 0.0546 & 0.0083 & 6.59 +    the performance of bilu , bilut and ras is similar to that of the same preconditioners in example 1 . when bilu is applied , the triangular solvers are sped over 3 times faster and the speedup of the whole solving phase is about 7 . because of the irregular non - zero pattern of bilut , the speedup of bilut is around 2 . the average speedup of the solving phase is about 6 . the speedup increases when the number of blocks grows . ras is as stable as in example 1 . the triangular solvers are sped around 6 , and the average speedup of the solving phase is around 8 . a matrix from spe10 is used @xcite . spe10 is a standard benchmark for the black oil simulator . the problem is highly heterogenous and it is designed the way so that it is hard to solve . the grid size for spe10 is 60x220x85 . the number of unknowns is 2,188,851 and the number of non - zeros is 29,915,573 . the ilu preconditioners we use in this example are block ilu(0 ) and block ilut(14 , 0.01 ) . the performance data is collected in table [ tab - ex3 ] . llllllll pre & blocks & cpu ( s ) & gpu ( s ) & speedup & pre cpu ( s)&pre gpu ( s ) & speedup +    bilu & 16 & 92.80 & 12.78 & 7.25 & 0.0421 & 0.0118 & 3.54 + bilu & 128 & 86.22 & 12.05 & 7.14 & 0.0423 & 0.0119 & 3.56 + bilu & 512 & 92.82 & 12.87 & 7.20 & 0.0424 & 0.0119 & 3.56 + bilut & 16 & 32.00 & 7.17 & 4.46 & 0.0645 & 0.0747 & 0.86 + bilut & 128 & 42.51 & 7.82 & 5.42 & 0.0647 & 0.0747 & 0.86 + bilut & 512 & 47.44 & 8.80 & 5.37 & 0.0645 & 0.0747 & 0.86 + ras & 256 & 106.61 & 14.36 & 7.41 & 0.100 & 0.0198 & 5.07 + ras & 1024 & 110.36 & 16.36 & 6.73 & 0.124 & 0.0242 & 5.11 +    from table [ tab - ex3 ] , we can speed the whole solving phase 6.2 times faster when ilu(0 ) is applied . the speedup increases if we increase the number of blocks . the average speedup of bilu is about 3 and the average speedup for the whole solving stage is about 7 . in this example , bilut is the best , which always takes the least running time . however , due to its irregular non - zero pattern , we fail to speed the triangular solvers . the ras preconditioner is always stable , and the average speedup of ras is about 5 , while the average speedup of the solving phase is around 6.5 . we have developed a new matrix format and its corresponding triangular solvers . based on them , the block ilu(k ) , block ilut and restricted additive schwarz preconditioners have been implemented . the block ilu(0 ) is sped over three times faster , the block ilut is sped around 2 times , and the ras preconditioner is sped up to 7 times faster . the latter preconditioner is very stable and it can be served as a general preconditioner for parallel platform .
Abstract:

Answer: in this paper , we investigate gpu based parallel triangular solvers systematically . the parallel triangular solvers are fundamental to incomplete lu factorization family preconditioners and algebraic multigrid solvers . we develop a new matrix format suitable for gpu devices . parallel lower triangular solvers and upper triangular solvers are developed for this new data structure . with these solvers , ilu preconditioners and domain decomposition preconditioners are developed . numerical results show that we can speed triangular solvers around seven times faster . solver , gpu , parallel , linear system


Instruction: Write an abstract for the following acadamic paper.

Input:
Paper: we have recorded the flight of a fly during take off and landing using digital high speed photography . it is shown that the dynamics of flexible wings are different for these procedures . during this observation fly flew freely in a big box and it was not tethered .
Abstract: in this fluid dynamics video , we demonstrated take off and landing of a fly . the deformation of wings is in focus in this video .

Paper: the x - ray transient xte  j1752@xmath0223 was discovered by the _ rossi x - ray timing explorer _ ( _ rxte _ ) on 2009 october 23 . the 2010 outburst of xte  j1752@xmath0223 was monitored by the _ @xcite , the _ monitor of all - sky x - ray image _ ( _ maxi _ , * ? ? ? * ) , and _ swift _ @xcite in the x - ray domain . the outburst exhibited an evolving behavior in agreement with what would be expected from a stellar mass black hole binary .    in the preceding low - hard x - ray state , the radio counterpart was identified by the atca ( australia telescope compact array ) with a flux density of @xmath12  mjy and a flat spectrum consistent with that of a compact jet @xcite . xte  j1752@xmath0223 started to transit toward a high - soft x - ray state around 2010 january 20 ( e.g. * ? ? ? * ) . a 10-fold increase of radio flux density and a flat spectrum were observed by the atca on 2010 january 21 @xcite . the increased flux density suggested that a major ejection was imminent @xcite and , indeed , a transient jet was detected by the rapid - response evn ( european vlbi newtork ) observations on 2010 february 11 . the transitions from a low - hard x - ray state toward a high - soft x - ray state are in fact , commonly associated with discrete ejection events @xcite . together with follow - up vlba ( very long baseline array ) observations in 2010 february , the transient jet was found to show significant deceleration , possibly due to the jet interaction with the surrounding medium ( * ? ? ? * hereafter paper  i ) .     +    in this paper we present results from the additional three epochs that show a new transient ejection event and later a stationary feature which we identify as the core . the observations and data reduction are briefly described in section  [ sec2 ] . the very long baseline interferometry ( vlbi ) as well as the westerbork synthesis radio telescope ( wsrt ) imaging results are presented in section  [ sec3 ] . the proper motion of a transient ejecta detected in one of our earlier observations is analyzed in section  [ sec4 ] . we discuss our results on the compact jet region , the proper motion of the transient ejecta and their possible interaction with the large - scale medium in section  [ sec5 ] . all the vlbi experiments observed at 5  ghz during the 2010 outburst of xte  j1752@xmath0223 are listed in table  [ tab1 ] . the observations before march 2010 were already presented in paper  i. additional evn observations were carried out on 2010 march 22 , including effelsberg , the phased array wsrt , onsala , medicina , noto , toru , and the lovell telescope . there were two epochs of vlba observations on 2010 april 25 and 29 using the whole telescope array . all the experiments used the same phase - referencing source ( pmn  j1755@xmath02232 ) and the same duty - cycle . the vlbi calibration procedure was identical as well , and it is described in paper  i. we note that the new evn observations show a weak extended jet in the calibrator source , but since the radio core accounts for @xmath190% of the emission , this extension did not have a significant effect on the phase calibration in the first four epochs . in the wsrt data , we did self - calibration on pmn  j1755@xmath02232 and applied the solutions to the target . .the characteristics of vlbi images in figure  [ fig1 ] . [ cols="^,^,^,^,^,^,^ " , ]     223 component b. the straight line in the right panel shows the least - square fitting results : @xmath2  masd@xmath3 . note that the horizontal dotted line crossing each data point gives the bin length . the right panel shows the variation of the peak brightness after applying a series of proper motion corrections along the jet direction , as described in the text . , title="fig:",scaledwidth=23.0%,scaledwidth=30.0% ] 223 component b. the straight line in the right panel shows the least - square fitting results : @xmath2  masd@xmath3 . note that the horizontal dotted line crossing each data point gives the bin length . the right panel shows the variation of the peak brightness after applying a series of proper motion corrections along the jet direction , as described in the text . , title="fig:",scaledwidth=23.0%,scaledwidth=30.0% ] + component d was found to have no detectable proper motion with respect to the reference source at the last two vlbi epochs . the only explanation for this is that component d is the radio core of the system . we note that @xcite arrived to the same conclusion first : they independently analyzed data from the last two vlba epochs , and showed that component d was positionally consistent with the optical counterpart of xte  j1752@xmath0223 , within a @xmath4  mas error circle . they had only a 3.5@xmath5 detection on 2010 february 29 . in our analysis , we obtained a considerably higher significance of the detection reaching 5@xmath5 , firmly establishing that component d is related to the core region / compact jet in xte  j1752@xmath0223 and not related to a new transient ejecta .    compared with the most sensitive evn observations of 2010 march 22 during the x - ray soft state , the radio core re - brightened by a factor of @xmath610 . the non - detection in the earlier soft state is unlikely caused by the fast variability because the core was not detected on short time scales . the re - appearance provides a rare case of direct detection of the re - activation of radio emission in the compact jet during the soft to hard x - ray state transition , as described by the unified model for the coupling of accretion and ejection in black hole x - ray binaries @xcite . in fact , _ rxte _ observations on 2010 march 27 have revealed xte  j1752@xmath0223 in transition from the soft to the hard state @xcite . the compact jet appears highly variable . however , we note that the declination of the source is low and the strong apparent variation might be due to poor tropospheric conditions especially at the beginning and end of the observations , when the source elevation was very low . however the atca observations did show significant variations on timescales of days ( brocksopp et al . , in prep . ) . in literature , strong daily variations were also observed in xte  j1720@xmath0318 in the post low hard state @xcite . significant flux density variation on shorter ( @xmath7  day ) time scale were seen in grs  1915@xmath8105 @xcite and xte  j1650@xmath0500 @xcite .      according to the new proper motion measurements , component b is a transient jet component moving toward the same direction as component a ( cf . paper  i ) , in agreement with the core location discussed above . although we could not measure the proper motion of component c in the same way as component b , it lies between components a and b which indicates that there were at least three ejection events revealed by the vlbi observations during the 2010 outburst of xte  j1752@xmath0223 . the source shows a profound one - sided morphology , unprecedented for x - ray black hole binaries with transient ejecta . we searched for the receding jet in the vlbi images including the deep evn image . there were no receding ejecta found . if component a was ejected during the major x - ray outburst around mjd  55218 , a mean proper motion of @xmath128  masd@xmath3 is required to explain the large angular separation from the assumed core . considering the low proper motion of @xmath910  masd@xmath3 observed at the later stage , component a must have an extremely high proper motion at the early stage . the high proper motion observed in component b , an ejecta much closer to the core than component a , indirectly supports the association of component a with the major outburst . as a relativistic jet usually evolves with a decaying luminosity after the ejection , the non - detection of component b observed three days earlier is unlikely due to the insufficient sensitivity . assuming component b was ejected around mjd  55250.6 ( during another flare , brocksopp et al . in prep . ) , its average proper motion is @xmath158  masd@xmath3 , in agreement with our measurements . component b has a high apparent speed of @xmath10.7  c if it is at 3.5  kpc as estimated by @xcite via the spectral - timing correlation scaling technique ( model - dependent ) . however , if the source is located at a larger distance , e.g. near the galactic centre , then xte  j1752@xmath0223 might be another example of superluminal galactic jets , such as grs  1915@xmath8105 @xcite and gro  j1655@xmath040 @xcite . taking into account that component b might have already somewhat decelerated by the time it was detected , xte  j1752@xmath0223 is a promising galactic superluminal source candidate . the resulting doppler deboosting effect may be responsible for the non - detection of the receding transient jets , a rarely - seen property in stellar mass black hole binaries . there might be an extended emission in rough agreement with the position angle of the transient jet between xte  j1752@xmath0223 and g006.3904@xmath801.9844 in figure  [ fig3 ] . such the extended emission might be a large - scale diffuse relic emission associated with the transient . as the limited sensitivity @xmath10  mjy  beam@xmath3 , deep observations at lower frequencies are required to confirm its existence . a jet - blown bubble , @xmath15  pc in diameter , is found in  cygnus x-1 @xcite . it is inferred that there exists a cavity surrounding microquasars @xcite . as xte  j1752@xmath0223 showed strong jet decelerations , it is a promising candidate for the deep observations to hunt for the extended emission region . we investigated the proper motion of the closest ejecta to the assumed core among the detected transient jet components in the 2010 outburst of xte  j1752@xmath0223 . besides measuring proper motion through splitting the experiment into a few time bins , we looked for a proper motion via monitoring the variation of the image peak brightness of the proper motion - corrected visibility data . both ways give a very high proper motion and imply that the ejections are at least mildly relativistic . we also reported the detection of another new ejecta observed at the end of 2010 march . we provided another independent vlbi identification of the compact jet and have shown that it might be highly variable . moreover , we found a hint on the large - scale extended emission in the field of xte  j1752@xmath0223 with the wsrt observations . we thank the evn programme committee and the vlba proposal selection committee for prompt approval of our observations . e - vlbi developments in europe were supported by the ec dg - infso funded communication network developments project expres . the national radio astronomy observatory is a facility of the national science foundation operated under cooperative agreement by associated universities , inc . the evn is a joint facility of european , chinese , south african and other radio astronomy institutes funded by their national research councils . the wsrt is operated by astron with support from the netherlands foundation for scientific research .
Abstract:

Answer: the galactic x - ray transient xte  j1752@xmath0223 was shown to have properties of black hole binary candidates . as reported in our previous paper , we identified transient and decelerating ejecta in multi - epoch very long baseline interferometry ( vlbi ) observations with the european vlbi network ( evn ) and the nrao very long baseline array ( vlba ) . here we present new evn and vlba data in which a new transient ejection event and later a stationary component are identified . the latter is interpreted as a reappearance of the radio core / compact jet during the transition from soft to hard x - ray state . this component appears to be highly variable in brightness although effects of tropospheric instabilities might play a role too . we also re - analyze the earlier vlbi data and find that the transient ejecta closer to the core position has significantly higher proper motion , further strengthening the case for strongly decelerating ejecta on the scale of several hundred milli - arcsecond , never observed in x - ray binaries before . although the distance of the source is not well constrained , it is clear that these ejectas are at least mildly relativistic at the early stages . moreover , we show the large scale environment of the transient from the westerbork synthesis array data recorded in parallel during the evn run . [ firstpage ]    stars : individual : xte  j1752@xmath0223  stars : variable : others  ism : jets and outflows  radio continuum : stars  x - rays : binaries .


Instruction: Write an abstract for the following acadamic paper.

Input:
Paper: we have recorded the flight of a fly during take off and landing using digital high speed photography . it is shown that the dynamics of flexible wings are different for these procedures . during this observation fly flew freely in a big box and it was not tethered .
Abstract: in this fluid dynamics video , we demonstrated take off and landing of a fly . the deformation of wings is in focus in this video .

Paper: over the past decades magnetic thin films have been the focus of intense work . when their thickness is reduced surface and dimensionality effects become dominant , affecting e.g. magnetic anisotropy and thermal excitations . more recently the focus has been extended to systems of lower dimensionality , e.g. nanosized wires and dots . beyond the limits of chemical routes @xcite and lithography @xcite , self - organization is the most suited technique to achieve nanostructures down to the atomic size @xcite . self - organized dots @xcite and wires @xcite with a long range order could be grown on templates like surface reconstructions and vicinal surfaces . achieving self - organization on surfaces that spontaneously do not form a template is appealing , because it should be more versatile . concerning wires fabricated on _ non - vicinal _ surfaces , one route is to ion - etch continuous film under grazing incidence @xcite . here we propose an alternate route , consisting in growing magnetic wires on a non - magnetic nanotemplate of grooves , which is achived by a growth process out - of - equilibrium . the samples were grown in an ultra - high vacuum chamber using pulsed - laser deposition with a nd - yag laser at . the chamber is equipped with a quartz microbalance , sample heating and a translating mask for the fabrication of wedge - shaped samples . we use a 10 kev reflection high energy electron diffraction ( rheed ) setup with a ccd camera synchronized with laser shots , which permits operation during laser operation and thus deposition . an omicron room - temperature scanning tunneling microscope ( stm ) , used in the constant current mode , and an auger electron spectrometer ( aes ) are located in connected chambers . the metallic films are grown on commercial sapphire single crystals . a detailed description of the chamber and growth procedures can be found in @xcite . the magnetic measurements were performed ex - situ on capped samples , with a superconducting quantum interference device ( squid ) magnetometer . we start from mo(110 ) or w(110 ) buffer layers ( thickness @xmath3{nm}$ ] ) deposited on al@xmath4o@xmath5 ( 11@xmath60 ) substrates . the process consists of deposition at room temperature ( rt ) followed by annealing at @xmath7c , yielding a surface of quality similar to that of metal single crystals @xcite . in the case of w an underlayer of mo ( ) is deposited first at rt , to avoid twinning of the annealed film @xcite . it has been reported upon deposition at moderate temperature along the ( 110 ) plane , that fe @xcite and w @xcite develop a self - organized surface with uniaxial hills and grooves oriented along the [ 001 ] direction . the origin of this structure was explained by an anisotropic diffusion at the nucleation stage , followed by a kinetic effect related to the ehrlich - schwoebel barrier . the average angle of the grooves was reported and simulated to increase with nominal thickness and with lowering temperature . no limiting angle was discussed in @xcite , while in @xcite a limiting angle of @xmath8 associated with \{310 } facets was predicted . here we confirm the formation of grooves for w ( fig.1 - 2 ) and report it also for mo . we focus on deposition at moderate temperature for which the angle of the facets is significant , i.e. @xmath9c for mo and @xmath10c for w. monitoring rheed during growth we observe that the angle of the facets becomes stationary after a while , typically a few nanometers of nominal thickness . rheed and stm show in agreement that the stationary angle is close to @xmath11 , consistent with facets of type \{210 } ( fig.1 - 2 ) . thus the prediction of ref . @xcite for a stationary angle is confirmed , however the resulting facet is not the one that was foreseen .    the quality of the order is reflected on the rheed patterns by the occurrence of satellite streaks ( fig.1b ) , whose splitting agrees quantitatively with the period found by stm , e.g. for deposition at @xmath1c ( fig.1a ) . the average depth of the grooves is here . in @xcite it was argued that the period is selected during the nucleation stage in the sub - atomic - layer range of deposition . we thus deposited of w under decreasing temperature ( from @xmath12c to @xmath1c ) allowing us to reach a period of 10- and grooves of depth 2- . stm and rheed ( fig.2a - b ) reveal again the \{210 } faceting , which becomes stationary after 1- have been deposited . notice that the order is slightly decreased for this process , with respect to deposition at @xmath1c . we have used the surfaces with hills and grooves reported above as a template for the self - organization of magnetic wires . we focus on w as a material template because the faceting remains stable until a higher temperature than in the case of mo . we deposited fe at @xmath1c for which fe grows essentially layer - by - layer on _ flat _ w(110 ) @xcite . samples wedged from 0 ( w template ) to 3 atomic layers were investigated by stm . cross - sections across the wires ( fig.1c ) and calculation of the roughness on the stm images show that the depth of the grooves decreases with the amount of fe deposited . this suggests that fe grows inside the grooves and not uniformly or at the top of the hills , and thus forms wires . this conclusion is currently under investigation by quantitative aes peak analysis and tem . the process described here works equally well independently from the period of the template . we focus here on fe wires prepared on w(110 ) templates with a mean period of . for an fe nominal thickness of the wires are expected to have a mean width of and a mean height of from the topography observed by stm . the magnetic properties are examined _ ex - situ _ by squid using a -mo - capped sample . hysteresis loops were performed along two in - plane directions : parallel to fe[001 ] , i.e. along the wires and parallel to fe[1@xmath130 ] , i.e. across the wires . the easy magnetization axis is [ 001](fig.3a ) . the anisotropy energy is calculated along the [ 1@xmath130 ] loop like @xmath14 . this anisotropy originates as the sum of many contributions , among which only the dipolar energy can be estimated reliably : @xmath15{j / m^{3}}$ ] . other sources of anisotropy are surface ( nel - type ) expected to favor the [ 1@xmath130 ] direction for the ( bottom ) fe / w interface @xcite and [ 001 ] for the ( top ) mo / fe interface @xcite , step - edge for fe / w , expected to favor [ 1@xmath130 ] @xcite ; magneto - elastic  unlike thin films , investigated so far , stress is both in - plane and out - of plane owing to the occurrence of many steps , implying significant shear , so that no figure or even sign can be reliably predicted ; finally the bulk fe anisotropy is negligible . surprisingly @xmath16 despite this complex situation , similarly to ( fe , ag ) self - organized arrays of wires @xcite . + the loops measured with the external field h//[001 ] at different temperatures have a rather square shape , and at remanence full saturation is still observed ( fig.3b ) . the coercivity continuously decreases with temperature while m@xmath17 remains essentially unchanged suggesting a superparamagnetic behavior.the ultimate blocking temperature determined by the `` zero field cooling / field cooling '' process is @xmath18{k}$ ] ( fig.3c ) . mean _ value of t@xmath19 is around , estimated as the half - way of the zero - field cooling remagnetization curve . arrays of fe nanowires were obtained by deposition on self - organized kinetically faceted surfaces of w , with a period in the range 4- . arrays with a mean period of , prepared by deposition of an amount of of fe , display an in - plane uniaxial anisotropy along [ 001 ] , the direction of the wires . the wires are superparamagnetic at room temperature and the mean blocking temperature is . efforts to rise this blocking temperature are under way . we are grateful to ph . david and v. santonacci for technical support . fert a. , piraux l. , j. magn . magn . mater . * 200 * ( 1999 ) 338 - 358 . martn j.i . , nogus j. , kai liu , vicent j.l . , schuller ivan k. , j. magn . magn . mater . * 256 * ( 2003 ) 449 - 501 . fruchart o. , c. r. physique * 6 * ( 2005 ) 61 - 73 . rohart s. , repain v. , tejeda a. , ohresser p. , scheurer f. , bencok p. , ferre j. , rousset s. , phys . b * 73 * ( 2006 ) 165412 - 1 - 7 . voigtlnder b. , meyer g. , amer nabil m. , phys . b * 44 * ( 1991 ) 10354 - 10357 . gambardella p. , rusponi s. , veronese m. , dhesi s.s . , grazioli c. , dallmeyer a. , cabria i. , zeller p. , dederichs p.h . , kern k. , carbone c. , brune h. , science * 300 * ( 2003 ) 1130 - 1133 . valbusa u. , boragno c. , buatier de mongeot f. , j. phys . : condens * 14 * ( 2002 ) 8153 - 8175 . fruchart o. , jubert p.o . , eleoui m. , cheynis f. , borca b. , david p. , santonacci v. , linard a. , hasegawa m. , meyer c. , j. phys . condens . matter - submitted , cond - mat/0608261 fruchart o. , jaren s. , rothman j. , appl . * 135 * ( 1998 ) 218 - 232 . albrecht m. , fritzsche h. , gradmann u. , surf . * 294 * ( 1993 ) 1 - 9 . khler u. , jensen c. , wolf c. , schindler a.c . , brendel l. , wolf d.e . surf . sci . * 454 - 456 * ( 2000 ) 676 - 680 . gradmann u. , korecki j. , waller g. , appl . a * 39 * ( 1986 ) 101 - 108 . hauschild j. , elmers h.j . , gradmann u. , phys . rev . b * 57 * ( 1998 ) r677-r680 . fruchart o. , nozires j.p . , givord d. , j. magn . magn . mater . * 207 * ( 1999 ) 158 - 167 . borca b. , fruchart o. , meyer c. , j. appl . * 99 * ( 2006 ) 08q514 - 1 - 3
Abstract:

Answer: homoepitaxy of w(110 ) and mo(110 ) is performed in a kinetically - limited regime to yield a nanotemplate in the form of a uniaxial array of hills and grooves aligned along the [ 001 ] direction . the topography and organization of the grooves were studied with rheed and stm . the nanofacets , of type \{210 } , are tilted @xmath0 away from ( 110 ) . the lateral period could be varied from 4 to by tuning the deposition temperature . magnetic nanowires were formed in the grooves by deposition of fe at @xmath1c on such templates . fe / w wires display an easy axis along [ 001 ] and a mean blocking temperature @xmath2{k}$ ] . epitaxial growth , iron , tungsten , molybdenum , nanowires , self - organization , nanotemplate , rheed patterns , magnetic anisotropy


Instruction: Write an abstract for the following acadamic paper.

Input:
Paper: we have recorded the flight of a fly during take off and landing using digital high speed photography . it is shown that the dynamics of flexible wings are different for these procedures . during this observation fly flew freely in a big box and it was not tethered .
Abstract: in this fluid dynamics video , we demonstrated take off and landing of a fly . the deformation of wings is in focus in this video .

Paper: super - resolution ( sr ) is a promising technology that is expected to be applied to microscope time series images , satellite photographs , and so on . sr aims at reconstructing a spatially high - resolution ( hr ) image from multiple low - resolution ( lr ) images . when multiple lr images obtained for the same object have transformation between images , they will contain mutually complementary information , and it becomes possible to estimate the hr image . since the earliest work @xcite , sr has been realized through various methods . in this work , we handle sr in a bayesian framework @xcite . in a bayesian framework , what prior we use is quite important . for example , various markov random field ( mrf ) priors @xcite , the total variation ( tv ) prior @xcite , the huber prior @xcite , and patch - based priors @xcite have been used in image processing . these can represent image properties well and have good performance in sr , image restoration , and other applications .    as the sr estimator , we think posterior mean ( pm ) is suitable as an hr image estimator because we usually evaluate the accuracy of sr methods by l2-norm ( mean square error ) -based peak signal - to - noise ratio ( psnr ) , and pm is the optimal estimator when employing psnr as the objective function . to determine the exact pm of the hr image , all parameters other than the hr image should be marginalized out over the joint posterior distribution without using any point estimation . according to this meaning , the previous methods @xcite are not optimal . recently , the pm approach was proposed @xcite . it was the first method to employ the pm as the hr image estimator and does not use any point estimation of the registration parameters or the model parameters .    in this manuscript , we propose a new sr method that employs a compound gaussian mrf prior that can simultaneously represent smoothness and edge of image @xcite and utilizes variational bayes to calculate the optimal estimator , pm , with respect to the objective function of the psnr . this approach seems quite favorable , but possibly it was not proposed earlier because of an important limitation of variational bayes that a conjugate prior is needed , and because of the exponential - order calculation cost for a compound gaussian mrf prior . we solve these problems through simple taylor approximations . in section 2 , we show the formulation regarding the models and the estimator . in section 3 , we evaluate the proposed method . in section 4 , we discuss . first , we show the definitions of the gamma , bernoulli and gaussian distributions , the kullback - leibler ( kl ) divergence from distribution @xmath0 to @xmath1 , and psnr , used in this manuscript : @xmath2 here , @xmath3 is the gamma function , @xmath4 denotes the determinant of a given matrix , @xmath5 is the dimension of @xmath6 , and the angle brackets @xmath7 denote the expectation of @xmath8 with respect to a distribution @xmath9 . also , @xmath10 denotes a diagonal matrix . all the vectors in this manuscript are column vectors . here , these variables have absolutely nothing to do with the variables that appear later . we estimate an hr grayscale image @xmath11 from observed multiple lr grayscale images @xmath12 using an sr technique . the images @xmath13 and @xmath6 are regarded as lexicographically stacked vectors . the number of pixels for each lr image is assumed to be less than that of the hr image ; i.e. , @xmath14 . although we define the range of a pixel luminance value as infinite , we use @xmath15 for black , @xmath16 for white , and values between @xmath15 and @xmath16 for shades of gray . the hr image @xmath6 is geometrically warped , blurred , downsampled , and corrupted by noise @xmath17 to form the observed lr image @xmath13 . it is modeled as @xmath18_{k=1}^4      \equiv [ \theta_l , [ \vec{o}_l]_x , [ \vec{o}_l]_y , \gamma_l ] ^\top.\end{aligned}\ ] ] the @xmath19 is additive white gaussian noise ( awgn ) with precision ( inverse variance ) @xmath20 ( @xmath21 ) . @xmath22 is the transformation matrix that simultaneously applies warping , blurring , and downsampling . the details of this matrix are given in @xcite . @xmath23 is a four - dimensional vector consisting of the registration parameters : rotational motion parameter @xmath24 , translational motion parameter @xmath25 , and blurring parameter @xmath26 .        here , the `` causal '' gaussian mrf prior used in @xcite is defined as the joint distribution of @xmath6 and @xmath27 in the form of @xmath39 , and it differs from the compound one in that it is not simultaneously normalized about both @xmath6 and @xmath27 like eq . ( [ eqpriordistributionanotherxeta ] ) . a `` causal '' one is an approximation of the compound one , and it is easier to use than the compound one because simultaneous normalization of a compound one has an exponential - order calculation cost with respect to the dimensionality of the line process ; the calculation cost of a `` causal '' one is polynomial . additionally , though both @xmath6 and @xmath27 of a compound one has a markov property , in the `` causal '' one only @xmath6 has a markov property . however , in eq . ( [ eqpriordistributionanotherxeta ] ) , ignoring @xmath40 as in @xcite makes them take the same form and breaks either property . therefore , in section 3 , we propose a new approximation that does not ignore @xmath40 . the hyperparameter priors and the registration parameter priors are defined as @xmath41 where @xmath42 , @xmath43 , @xmath44 $ ] and @xmath45 $ ] . for a gamma distribution , the number of effective prior observations in the bayesian framework is equal to two times parameter @xmath46 . the above settings are considered sufficiently non - informative . here , we assume @xmath47 by setting its prior according to a gamma distribution similar to @xcite , resulting in an appropriate inference . we define the mean value for the gamma distribution as @xmath48 , @xmath49 , and @xmath50 . @xmath51 denotes an iterative step later introduced through variational bayes . the settings of the registration parameter priors are considered suitable for this sr task @xcite . note that the mean value @xmath52 of the @xmath26 prior is derived as the value equivalent to the anti - aliasing of the scale factor @xcite . first , we confirm that the joint distribution of all random variables can now be explicitly given as @xmath53 , \bphi ] .\end{aligned}\ ] ] we define the objective function using @xmath54 and optimal estimator as @xmath55      = \left\langle \bx \right\rangle_{p(\bx | \by)}.\end{aligned}\ ] ] since only lr images , @xmath56 , are available for the estimator , we sometimes explicitly express it as a function form , @xmath57 . we choose this objective function because we prefer good estimator performance on average over various hr images and the corresponding lr images . from eq . ( [ objectivefunction ] ) , the pm @xmath58 is the best estimator of the hr image . note that the posterior distribution of an hr image @xmath59 needs marginalization of all parameters other than @xmath6 over @xmath60 .      as stated above , we could derive the optimal estimator . however , we can not obtain the marginalized posterior distribution @xmath59 analytically . therefore , we solve this through approximation by using variational bayes @xcite . we impose the factorization assumption on the trial distribution @xmath61 . the optimal trial distribution is identified by minimizing the kl divergence between the trial and the true distributions as the best approximation of the true distribution : @xmath62 . in the common style of variational bayes @xcite , update equations are @xmath63 frequently , the application of variational bayes is difficult in practice because it requires a conjugate prior . we make the priors conjugate by using simple taylor approximations similar to @xcite . we also solve the exponential - order calculation cost problem of a compound gaussian mrf prior by using the same approach here . these approximations enable the analytical calculations in eq . ( [ eqrecurrenceinitiation ] ) . specifically , we apply the first - order taylor approximations for three non - linear terms . first , in eq . ( [ eqimageobservationprocess ] ) , @xmath22 is approximated around @xmath64 the same as in @xcite . second , the logarithm of the normalization term of eq . ( [ eqpriordistributionanotherxeta ] ) is approximated around @xmath65 . finally , in eq . ( [ eqpriordistributionanotherxeta ] ) , @xmath66 is approximated around @xmath67 = [ \bmu_\bmeta^{{(\!t\!+\!1\ ! ) } } , \ln \mu_\rho^{{(\!t\ ! ) } } , \ln \mu_\kappa^{{(\!t\!)}}]$ ] similar to @xcite . note that in @xcite , the third approximation is the key idea to solve the conjugate prior problem . in this work , we found it also enables us to solve the exponential - order calculation cost problem of a compound gaussian mrf prior . this makes it possible to calculate the normalization term .      from eq . ( [ eqrecurrenceinitiation ] ) and the taylor approximations , the trial distributions are obtained as the following distributions : @xmath68 the specific form of the trial distributions is omitted because of space limitations . for eq . ( [ eqrecurrenceinitiation ] ) , we do the following . first , we compute @xmath69 using @xmath70 . second , we compute @xmath71 using @xmath72 . last , we compute @xmath73 using @xmath74 and @xmath75 using @xmath76 . for the initial parameters of the trial distributions of @xmath27 and @xmath6 , we use non - informative values , @xmath77 . as the initial parameters for @xmath35 , @xmath37 , @xmath20 , @xmath38 , and @xmath78 we use the same values as their prior s values in eq . ( [ eqpriordistributionlambda ] ) , ( [ eqpriordistributionphi ] ) . we obtain the well approximated pm of @xmath6 as @xmath79 , for which the following convergence conditions hold for @xmath79 and each @xmath80 , @xmath81_{k } } & < 10^{-5 } \quad ( k=1,2,3,4),\end{aligned}\ ] ] where we defined @xmath82 $ ] as the scaling constant . 0    [ cols="^,^,^,^,^,^,^",options="header " , ]     we evaluated the proposed method using five gray - scale images with a size of @xmath83 pixels as shown in fig . [ figoriginal ] . from each image , @xmath84 images with a size of @xmath85 pixels were created by using eq . ( [ eqimageobservationprocess ] ) . the resolution enhancement factor was 4 . the transformation parameter @xmath78 was generated from the prior distribution in eq . ( [ eqpriordistributionphi ] ) , where it is similar to that in previous work @xcite . the noise level parameter @xmath20 was set for a signal - to - noise ratio ( snr ) of @xmath86 , @xmath87 , and @xmath88 db for each image . samples of the created images are shown in fig . [ figobserved ] . figure [ figinferred ] shows the estimated images under snr@xmath89db . we can see that the resolution of each image appeared to be better than that of the corresponding observed image . table [ tableresultspsnr ] shows the quantitative results compared to methods ( a ) , ( b ) , ( c ) , and ( d ) , where ( a ) is bilinear interpolation , ( b ) is the variational em approach with a causal gaussian mrf prior using the maximum a posteriori ( map ) estimator @xcite , ( c ) is the variational bayes approach with a tv prior using the map estimator @xcite , and ( d ) is the pm approach with a causal gaussian mrf prior using the pm estimator @xcite . note that we added a slight modification to methods ( b ) and ( c ) , the same as in @xcite , because they employ slightly different models . over @xmath90 experiments on each image and for each snr , we evaluated the results with regard to the expectation and the standard deviation of the psnr of the proposed method and the improvement in signal - to - noise ratio ( isnr ) defined as , @xmath91 where @xmath92 is the image estimated by the proposed method , and @xmath93 is the image estimated by the method to compare ( i.e. , ( a ) , ( b ) , ( c ) , and ( d ) in these experiments ) . we see that the isnrs of the images estimated by the proposed method were mostly better than those by the other methods . in the subjective visual comparison in fig . 4 , we also see that the edges are not overemphasized in the images estimated by the proposed method compared to those in the images estimated by other methods .    the calculation times with the proposed method and with methods ( b ) , ( c ) and ( d ) for each estimate , using an intel core i7 2600 processor , were almost the same , about @xmath86 minutes . [ figoriginal ]    [ figobserved ]    [ figinferred ]    [ figincompared ] in this manuscript , we have shown how we can adopt a compound gaussian mrf prior for pm sr . we got good results for almost all images and noise conditions . in the comparison with method ( c ) and ( d ) , the superiority of isnr was fair and the inferiority of isnr was small . regarding these some unfavorable results in table [ tableresultspsnr ] , we think that the reason is our numerical optimization method falls slightly short of optimization because the proposed method uses more approximations than method ( d ) and ( c ) . in addition , we consider the result compared to method ( c ) in the case of the pepper image in @xmath94 db noise is caused by unstable estimation of @xmath95 and @xmath37 , where method ( c ) fixed the value of @xmath95 to the true expected value in our implementation same to @xcite . therefore , we think compound gaussian mrf prior is considered preferable to a `` causal '' gaussian mrf prior as a natural image prior . regarding the estimator , we used the optimal estimator , the pm . from the experimental results , we see that the sr methods with the pm estimator ( i.e. , the proposed method and method ( d ) ) were more accurate than the sr methods with other estimators ( i.e. , method ( b ) and method ( c ) ) . this indicates that pm is an optimal estimator . regarding the calculation cost , our algorithm requires @xmath96 , and we must make our algorithm faster to apply this method to larger images . by using an idea similar to that of @xcite , we have developed a faster algorithm , but this algorithm causes obvious degradation of accuracy . we continue to search for a more effective way to reduce the calculation cost with less degradation of accuracy . we can say that the proposed method is an sr method with a favorable model and an optimal estimator . in addition , the proposed method does not need any parameter tuning , which is a favorable property in practice . we think our approach to the problem about the conjugate prior and the exponential - order calculation cost can be applied to many other problems , and we will try to do this in our future work .    99 tsai , r.y , and huang , t.s . , `` multiframe image restoration and registration , '' in _ advances in computer vision and image processing _ , ct : jai press , vol . 1 , pp . 317339 , greenwich , 1984 . m. e. tipping , and c. m. bishop , `` bayesian image super - resolution , '' in _ advances in nips 15 _ , pp . 12791286 , mit press , 2003 . hardie , r.c . , barnard , k.j . , and armstrong , e.e . , `` joint map registration and high resolution image estimation using a sequence of undersampled images , '' _ ieee trans . image process . _ , vol . 6 , no . 12 , pp . 16211633 , 1997 . a. kanemura , s. maeda , and s. ishii , `` hyperparameter estimation in bayesian image superresolution with a compound markov random field prior , '' _ ieee international workshop on machine learning for signal processing ( mlsp ) _ , pp . 181186 , 2007 . a. kanemura , s. maeda , and s. ishii , `` superresolution with compound markov random fields via the variational em algorithm , '' _ neural networks _ , vol . 22 , pp . 10251034 , 2009 . s. geman , and d. geman , `` stochastic relaxation , gibbs distributions , and the bayesian restoration of images , '' _ ieee trans . on pattern analysis and machine intelligence _ , vol . pami6 , no . 6 , pp . 721741 , 1984 . r. molina , j. mateos , a. k. katsaggelos , and m. vega , `` bayesian multichannel image restoration using compound gauss - markov random fields , '' _ ieee transactions on image processing _ , vol . 12 , pp . 16421654 , 2003 . r. chellappa , and a. jain , _ markov random fields : theory and application _ , san diego , ca : academic press , 1991 . f. -c . jeng , and j. w. woods , `` compound gauss - markov random fields for image estimation '' , _ ieee transactions on signal processing _ , vol . 3 , pp . 683697 , 1991 . t. katsuki , m. inoue , `` posterior mean super - resolution with a causal gaussian markov random field prior , '' _ arxiv:1109.5453v3 _ , 2012 s. d. babacan , r. molina , a. k. katsaggelos , `` variational bayesian super resolution , '' _ ieee transactions on image processing _ 4 , pp . 984999 , 2011 . s. villena , m. vega , s.d . babacan , r. molina , and a. k. katsaggelos , `` image prior combination in super - resolution image registration @xmath97 reconstruction , '' _ ieee international workshop on machine learning for signal processing(mlsp ) _ , pp . 355360 , 2010 . l. c. pickup , d. p. capel , s. j. roberts , and a. zisserman , `` bayesian image super - resolution , continued , '' in _ advances in nips 19 _ , mit press , 2007 . m. ranzato , v. mnih , g. hinton , `` generating more realistic images using gated mrf s '' . _ advances in nips 23 _ , mit press , 2010 . h. attias , `` inferring parameters and structure of latent variable models by variational bayes , '' in _ proc . , san francisco , ca , pp . 2130 , morgan kaufmann , 1999 . c. m. bishop , d. spiegelhalter , j. winn `` vibes : a variational inference engine for bayesian networks , '' in _ advances in nips 15 _ , pp . 777784 , mit press , 2003 .
Abstract:

Answer: this manuscript proposes a posterior mean ( pm ) super - resolution ( sr ) method with a compound gaussian markov random field ( mrf ) prior . sr is a technique to estimate a spatially high - resolution image from observed multiple low - resolution images . a compound gaussian mrf model provides a preferable prior for natural images that preserves edges . pm is the optimal estimator for the objective function of peak signal - to - noise ratio ( psnr ) . this estimator is numerically determined by using variational bayes ( vb ) . we then solve the conjugate prior problem on vb and the exponential - order calculation cost problem of a compound gaussian mrf prior with simple taylor approximations . in experiments , the proposed method roughly overcomes existing methods . super - resolution , fully bayesian approach , markov random field prior , variational bayes , taylor approximation .


Instruction: Write an abstract for the following acadamic paper.

Input:
Paper: we have recorded the flight of a fly during take off and landing using digital high speed photography . it is shown that the dynamics of flexible wings are different for these procedures . during this observation fly flew freely in a big box and it was not tethered .
Abstract: in this fluid dynamics video , we demonstrated take off and landing of a fly . the deformation of wings is in focus in this video .

Paper: binary bismuthide @xmath6-pdbi@xmath0 has attracted much interest recently as a promising candidate of topological superconductor ( ts)@xcite . topological superconductivity is a new state of matter possessing symmetry - protected surface states while the bulk states are fully gapped by superconducting pairing@xcite . the majorana fermions are believed to exist on the surface or vortex core in such tss , which may not only be of scientific importance , but also can lead to a wide - ranging applications in microelectronic devices and quantum computing . the centrosymmetric stoichiometric @xmath6-pdbi@xmath0 ( @xmath7@xmath8 5 k ) was claimed to be topologically nontrivial in view of the observation of the topologically - protected surface modes by spin- and angle - resolved arpes@xcite . however , no andreev bound states associated with majorana fermions are detectable through point - contact spectroscopy@xcite , in sharp contrast to the cases in cu - intercalated bi@xmath0se@xmath9@xcite and in - doped snte@xcite . on the other hand , it becomes the common wisdom that spin - orbit interaction ( soi ) in heavy elements is crucial for the topological states . it is therefore heuristic to ask what if we replace pd by heavier pt element with enhanced soi .    . ( a ) the primitive unit cell for hexagonal ptbi@xmath0 . the coordinates of pt(1 ) and pt(2 ) are ( @xmath10 and ( @xmath11 , respectively . ( b ) the structure as seen from a perspective along the @xmath1-axis . , width=340 ]    in this study , we substituted pt for pd in pdbi@xmath0 and found that this new material actually crystallizes in a distinct structure . unlike @xmath6-pdbi@xmath0 which has the tetragonal structure in an @xmath124/@xmath13 space group , ptbi@xmath0 crystallizes in space group p-3 with a hexagonal unit cell of @xmath2=@xmath3=6.553@xmath4 , @xmath1=6.165@xmath4@xcite it is also different from its homologue ptbi superconductor ( @xmath7=1 k ) with a monoclinic unit cell@xcite . the in - plane resistivity of ptbi@xmath0 shows metallic behaviors down to 2 k , the lowest temperature studied in this work . the intra - plane and inter - plane magnetization displays pronounced anisotropy , being diamagnetic with field aligned along the plane and paramagnetic when field is perpendicular to the plane . the magnetoresistance ( mr ) and hall resistivity measured on the same sample both show two types of carriers and the former one scales well to the semi - classical kohler s rule@xcite . ptbi@xmath0 single crystals were fabricated via a melt - growth method . the starting materials of high purity , bi powder(4n ) and pt powder ( 4n ) , were mixed thoroughly in the prescribed molar ratio of bi : pt = 2:1 ( 2 g in total weight ) . all these preparations were performed in a glove box filled with protective argon gas ( both h@xmath0o and o@xmath0 contents were limited below 0.1ppm ) . the mixtures were loaded and sealed in an evacuated quartz tube . this quartz tube was then heated to 700@xmath14c quickly in a sintering furnace and kept at this temperature for 48h , before being slowly cooled down to 450@xmath14c(3@xmath14c / h ) , and finally being quenched into cold water . large pieces of dark - gray plate - like ptbi@xmath0 single crystals of typical 7 - 8 mm in length were harvested . energy dispersive x - ray ( edx ) spectrometry confirms the stoichiometric ratio of the chemical composition ( 32.8 : 67.2 @xmath15 3.0% in molar percentage for pt : bi ) . the structure of crystals was characterized by powder x - ray diffraction ( xrd ) at room temperature using a rigaku diffractometer with cu @xmath16@xmath17 radiation and a graphite monochromator . lattice parameters were obtained by rietveld refinements . the magnetization was measured by vibrating sample magnetometry using a quantum design mpms-5 system . measurements of mr and hall effect were performed on the same sample by changing the field polarities . signal even in field was defined as mr and the odd component was calculated as hall resistivity . resistivity fits to @xmath18 very well below @xmath835k ( see the upper - left inset ) . ( b ) and ( c ) show the in - plane and inter - plane susceptibility under a field of 5koe , respectively.,width=347 ] the schematic view of the crystal structure of ptbi@xmath0 is shown in fig . it crystallizes in a hexagonal structure with the space group p-3 ( no.147 ) . its structure consists of alternate stacking of 2d pt layers and bismuth bilayers along the @xmath1-axis . in one primitive unit cell , there are three pt atoms , one being located at the corner of the polyhedron and the other two labelled as pt(1 ) and pt(2 ) in fig . 1 . the bi atoms are trigonally - coordinated . the xrd pattern of ptbi@xmath0 crystal is presented in fig.2 . a small trace of impurity phase , marked by the asterisks in panel ( a ) , was detectable in the powder x - ray pattern and only ( 00@xmath19 ) diffraction peaks were observed in the single - crystal x - ray , indicating good @xmath1-axis orientation of the as - grown samples . the calculated lattice parameters are @xmath2=@xmath3=6.553@xmath4 , @xmath1=6.165@xmath4 , in consistence with previous reported results@xcite .        zero - field in - plane resistivity is plotted in fig . the room temperature resistivity is about 0.12 m@xmath20 cm and it is metallic down to the lowest temperature we measured ( 2k ) . the residual resistivity ratio is approximately 50 for our samples , indicative of good sample quality . the sample is better characterized by the susceptibility measurements thereafter . remarkably , the magnetization of the sample shows large anisotropy with respect to the field orientations . as illustrated in fig . 3 , the in - plane magnetization @xmath21 is diamagnetic and varies little with @xmath5 down to 20k , below which it displays a significant upturn , whereas the inter - plane @xmath22 is paramagnetic instead and increases linearly with decreasing @xmath5 , followed by a downward trend below 20k . the origin of these intriguing magnetization behaviors is not clear . the magnetoresistive and hall response of a material can open a avenue for exploring the dispersion and dynamics of the charge carriers . first , in ptbi@xmath0 , it is noted that the absolute value of the mr , defined as @xmath23 , is rather large , reaching @xmath24400% at 2k in a magnetic field of 9 t . this large mr implies a rather large electron mean free path , hence a long relaxation time . however , this mr is damped very fast with increasing @xmath5 , as seen from the upper panels of fig . 4 . second , in single - band metals , the mr at small fields is usually quadratic and the hall resistivity varies linearly with field . however , in the two - band drude model , on the assumption of the field - independent carrier density and relaxation time , @xmath25 and @xmath26 can be written as@xcite    @xmath27    @xmath28    where @xmath29 and @xmath30 are electrical conductivity and hall coefficient for electron ( hole ) band , respectively . the mr and the hall signal for ptbi@xmath0 sample are exemplified in fig 4 at some selected temperatures . although the individual curves can be fitted with the above two - band equations reasonably well , plotted as the red solid line in each panel , we failed to model these two transport coefficients simultaneously with the same set of four parameters . these difficulties may arise from the simple assumption of the field independent charge carrier density and scattering , in analogy to the case in cuprates@xcite . nevertheless , given the quality of our fitting and the strong non - linearity of the hall resistivity , we strongly believe that the transport properties of this compound are governed by two - band charge carriers .    in standard metals , the mr @xmath31/@xmath32 at a certain temperature under a field @xmath33 has a general form known as the kohler s rule@xcite : @xmath31/@xmath32=@xmath34(@xmath35 ) . this rule can be derived from boltzmann transport theory , on the assumption of constant carrier number with @xmath5 and a single scattering rate on the fermi surface . from this rule , @xmath31/@xmath32 is literally independent of @xmath5 such that the plots of @xmath36 as a function of @xmath35 at distinct temperatures will collapse onto a single curve . interestingly , this rule , albeit its semiclassical origin , was found to be well obeyed in a large number of metals from conventional metals to some quantum matters . these involve the metals with two types of carriers@xcite , the pseudogap phase of the underdoped cuprates@xcite , quasi - one - dimensional metals@xcite as well as some topological semimetals@xcite . we examined this rule in ptbi@xmath0 ( fig . 5 ) and found that it is well obeyed in this material , over a wide field range ( up to 9 t ) and a broad @xmath5 window ( 2k-100k . above 100k , the mr tends to be negligible ) . moreover , the longitudinal mr in ptbi@xmath0 ( h @xmath37 i @xmath37 @xmath38 ) also shows two types of charge carriers and the validity of the kohler s rule ( data not shown ) . in recent work by sakano _ et al . _ , several topologically - protected surface states were observed by spin - resolved arpes in the ts candidate @xmath6-pdbi@xmath0@xcite . these non - trivial surface bands include one crossing the fermi level and the other one forming the dirac cone state 2 ev below the fermi level . it was noted that these topological surface states are _ all _ derived as a consequence of soi , although their respective microscopic details may be different . in ptbi@xmath0 , the soi ought to be stronger . owing to its good metallicity , however , the electrical transport is _ overall _ dominated by its bulk electrons and it looks more like a conventional good metal from transport perspective . in this material , the possible quantum linear mr arising from the degenerate dirac fermions in the quantum limit is not observed up to 9t@xcite . interestingly , this material was reported to superconduct below 150 mk@xcite , @xmath840 times lower than @xmath7 in pdbi@xmath0 . how the soi changes the electronic structure of ptbi@xmath0 , and induces the non - trivial surface states , if any , await more investigations , both theoretically and experimentally .    to summarize , we synthesized the single crystals of stoichiometric bismuthide ptbi@xmath0 by a solid - state reaction method . the samples were carefully characterized by combined procedures of xrd , ( magneto-)transport and susceptibility measurements . this compound shows prominent two - band transport behaviors with no clear signature from the possible surface states . however , the high - quality single crystals are now ready for prospective advanced experiments , especially for ones with more surface sensitivity . the authors would like to thank c. m. j. andrew , a. f. bangura for stimulating discussions . this work is sponsored by the national key basic research program of china ( grant no . 2014cb648400 ) , and by national natural science foundation of china ( grant no . 11474080 , u1432135 , 11611140101 ) . x.x . would also like to acknowledge the financial support from the distinguished young scientist funds of zhejiang province ( lr14a040001 ) and an open program from wuhan national high magnetic field center ( 2015kf15 ) . y. imai , f. nabeshima , t. yoshinaka , k. miyatani , r. kondo , s. komiya , i. tsukada , and a. maeda , j. phys . . jpn . * 81 * , 113708 ( 2012 ) . k. zhao , b. lv , y. xue , x. zhu , l. deng , z. wu , c. w. chu , phys . b * 92 * , 174404 ( 2015 ) . e. herrera , _ et al . b * 92 * , 054507 ( 2015 ) . m. sakano , k. okawa , m. kanou , h. sanjo , t. okuda , t. sasagawa , k. ishizaka , nat . comm . * 6 * , 8595 ( 2015 ) . j. kamark , _ et al . b * 93 * , 144502 ( 2016 ) . che , t. le , xiaofeng xu , xin lu , to appear in phys . hasan , c.l . kane , rev . phys . * 82 * , 3045 ( 2010 ) . hor , a.j . williams , j.g . checkelsky , p. roushan , j. seo , q. xu , h.w . zandbergen , a. yazdani , n.p . ong , and r.j . cava , phys . lett . * 104 * , 057001 ( 2010 ) . liang fu , c. l. kane , phys . b. * 76 * , 045302 ( 2007 ) . liang fu , c. l. kane , e. j. mele , phys . lett . * 98 * , 106803 ( 2007 ) . m. kriener , k. segawa , z. ren , s. sasaki , and y. ando , phys . * 106 * , 127004 ( 2011 ) . s. sasaki , m. kriener , k. segawa , k. yada , y. tanaka , m. sato , and y. ando , phys . 107 * , 217001 ( 2011 ) . m. novak , s. sasaki , m. kriener , k. segawa , and y. ando , phys . b * 88 * , 140502 ( 2013 ) . e. m. savitskii , v. v. baron , yu . v. efimov , m. i. bychkova , l. f. myzenkova , _ superconducting materials _ ( springer 1973 ) . matthias , phys . rev . * 92 * , 874 ( 1953 ) . m. kohler , ann . phys . * 32 * , 211 ( 1938 ) . n. luo , and g. h. miley physica c * 371 * , 259 ( 2002 ) . p. li , f. f. balakirev , and r. l. greene , phys . 99 * , 047003 ( 2007 ) . f. rullier - albenque , d. colson , a. forget , h. alloul , phys . * 103 * , 057001 ( 2009 ) . p. m. c. rourke , a. f. bangura , c. proust , j. levallois , n. doiron - leyraud , d. leboeuf , l. taillefer , s. adachi , m. l. sutherland , and n. e. hussey phys . b * 82 * , 020514(r ) ( 2010 ) . f. rullier - albenque , d. colson , a. forget , h. alloul , phys . rev 109 * , 187005 ( 2012 ) . chan , m.j . veit , c.j . dorow , y. ge , y. li , w. tabis , y. tang , x. zhao , n. barii and m. greven , phys . 113 * , 177005 ( 2014 ) . a. narduzzo , a. enayati - rad , s. horii , n.e . hussey phys . * 98 * , 146601 ( 2007 ) . x. xu , _ et al . _ , j. phys . : condens . matter * 27 * , 335701 ( 2015 ) . a. narayanan , m.d . watson , s.f . blake , n. bruyant , l. drigo , y.l . chen , d. prabhakaran , b. yan , c. felser , t. kong , p.c . canfield , a.i . coldea , phys . * 114 * , 117201 ( 2015 ) . abrikosov phys . b * 58 * , 2788 ( 1998 ) d.x . hor , j. xiong , r.j . cava , n.p . ong , science * 329 * , 821 ( 2010 ) . alkseevskii , yu . p. gaidkov , eksp . fiz * 25 * , 383 ( 1953 ) .
Abstract:

Answer: we report details of single crystal growth of stoichiometric bismuthide ptbi@xmath0 whose structure consists of alternate stacking of pt layer sandwiched by bi bilayer along the @xmath1-axis . the compound crystallizes in space group p-3 with a hexagonal unit cell of @xmath2=@xmath3=6.553@xmath4 , @xmath1=6.165@xmath4 . the magnetization data show opposite sign for fields parallel and perpendicular to the pt layers , respectively . the @xmath5-dependent resistivity is typical of a metal and the magnetic response shows clear two types of charge carriers and the validity of the semi - classical kohler s rule . its physical properties was discussed in comparison with recently proposed topological superconductor @xmath6-pdbi@xmath0 .