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: our understanding of the qcd matter at low temperatures ( @xmath0 ) and high baryon number densities ( @xmath1 ) has developed rapidly over the past few years ; for recent reviews , see @xcite . a likely scenario is that in this region , quarks form cooper pairs and new condensates develop . one of the most interesting questions today concerns the detailed mapping of the qcd phase diagram , especially in the above region where qcd may exhibit a colour superconducting phase . in this context , the light quark masses ( @xmath3 ) are expected to play an important role . figure 1 shows two likely scenarios corresponding to @xmath4 and @xmath5 . the transition between the quark - gluon plasma ( qgp ) phase and the colour superconducting phase is expected to be first order , in qcd @xcite .    in this paper , we study the first - order phase transition between the qgp phase and the colour superconducting phase ( fig . our approach is phenomenological and is similar in spirit to that adopted by ct and kharchenko @xcite . they , however , studied the bose - einstein condensation in a gas of trapped atomic hydrogen at temperatures @xmath6 and densities @xmath7 @xmath8 . temperatures and densities that we shall consider here are , of course , vastly different , and qgp has very little in common with the atomic hydrogen gas . _ interestingly , however , as we shall see , the same physics of activation energy barrier governs the formation of droplets in both cases . _ it is important to get a handle on the cold and dense quark matter because of its possible relevance to the neutron - star core . while the region of _ asymptotic _ densities is amenable to rigorous calculations , the same can not be said about the region of _ intermediate _ densities studied here . both numerical ( lattice ) and direct experimental studies of this region of the qcd phase diagram are out of reach at the moment . hence a phenomenological study such as the present one is appropriate and worthwhile . it has allowed us to make quantitative statements about the schematic phase diagram in fig . 1 .    alford et al . @xcite have discussed droplet formation in the transition between the hadronic matter phase and the 2sc phase ( fig . 1 ) , at @xmath9 , as a function of density . the chiral symmetry is restored inside their droplets ( which they call nucleons ) but is broken outside . neergaard and madsen @xcite also have studied the droplet formation but in the qgp to hadronic matter transition occurring at low baryon number densities . in this paper , we have considered the qgp to colour - superconducting phase transition occurring at high @xmath1 . _ to our knowledge , droplet formation in this transition has not been discussed previously in the literature . _ our droplets are the opposite of those in @xcite in the sense that the chiral symmetry is broken inside our droplet , but is restored outside . consider a gas of weakly interacting quarks , antiquarks and gluons in a volume @xmath10 . we take @xmath11 and @xmath12 to be the massless current quarks ( fig . 1 ) . we are interested in low temperatures ( @xmath0 ) and high number densities ( @xmath13 ) such that @xmath14 where @xmath15 is the thermal wave length ( see appendix ) . we also need the baryon number density @xmath16 to be large compared to that for the normal nuclear matter , namely 0.17 @xmath17 .    if the interactions among the particles are treated to the lowest order in @xmath18 , the strong coupling constant ( fig . 2 ) , the baryon number density is given by ( @xmath19 ) @xmath20,\ ] ] where @xmath21 2  ( spin ) @xmath22 3  ( colour ) @xmath22 3  ( flavour ) = 18 is the quark degeneracy factor , @xmath23 is the boltzmann constant and @xmath24 is the quark chemical potential @xcite . for a given baryon number density @xmath1 and temperature @xmath0 , eq . ( 1 ) uniquely fixes @xmath24 . this in turn determines the total energy and pressure of the gas as @xmath25,\\ p_{gas } & = & e_{gas}/3v.\end{aligned}\ ] ] here @xmath26 2  ( spin ) @xmath22 8  ( colour ) = 16 is the gluon degeneracy factor . it was possible to obtain the above analytic expressions because we have considered the gas of quarks as well as antiquarks , rather than a gas of quarks alone @xcite . at low @xmath0 and high @xmath1 , the @xmath24 terms in the pressure dominate and the gluonic pressure arising from the first term in eq . ( 2 ) , is quite negligible . we now consider the first - order phase transition from qgp to the superconducting phase as the temperature is lowered ( fig . 1 ) . in the latter phase , as a result of the attractive one - gluon - exchange interaction between quarks in the colour @xmath27 channel , quarks near the fermi surface tend to pair up as bosons . these are the cooper pairs which we assume to be massless . density fluctuations in the metastable qgp gas may provide suitable conditions for nucleation and growth of a droplet of the superconducting phase . the high - temperature phase ( qgp ) is present around the droplet and the low - temperature superconducting phase exists inside the droplet . we will presently explore the energetics for the formation of such a droplet by calculating the energy barrier between the two phases . we have given above the expressions for number density , energy and pressure of the qgp ; see eqs . ( 1)-(3 ) . we now derive corresponding expressions for the superconducting phase inside the droplet . before we begin , it is necessary to recall important differences between the bcs superconductor and the colour superconductor . ( a ) in the bcs theory of superconductivity the ratio of the energy gap ( @xmath28 ) to the fermi kinetic energy for a typical metal is of the order of @xmath29 . in contrast , in the problem under consideration , the gap @xmath28 is estimated to vary between several tens of mev to about 100 mev @xcite , and the quark chemical potential @xmath24 appearing in eq . ( 1 ) takes values between @xmath30 and @xmath31 mev for the @xmath1 and @xmath0 considered here . ( b ) in the bcs superconductor , the debye energy @xmath32 is much smaller than the chemical potential @xmath33 , while in the colour superconductor , the two may be of the same order of magnitude . ( c ) finally , if @xmath13 is the density of electrons and @xmath34 the pippard coherence length which measures the spatial extension of the cooper pair wave function , then in the bcs theory , @xmath35 . in the present case , on the other hand , @xmath36 fm for @xmath37 mev , and @xmath38 for the densities considered in this paper . this suggests that the @xmath39 pairs are rather compact objects unlike cooper pairs of electrons in a metal . for a good discussion of these and related issues , see a review article by kerbikov @xcite . we shall make use of the above points while deriving expressions for thermodynamics quantities of the droplet . for the sake of simplicity , we shall assume that all quarks in the energy interval @xmath40 form cooper pairs and all other quarks remain unpaired , and shall ignore other interactions . consider a spherical region of radius @xmath41 , volume @xmath42 , containing @xmath43 quarks out of which @xmath44 remain unpaired and @xmath45 undergo cooper pairing , when the pairing interaction is `` switched on '' . we have @xmath46 where @xmath47^{-1}$ ] with @xmath48 . the corresponding energies @xmath49 are @xmath50 in the last equation @xmath51 is the number of cooper pairs ( bosons ) in the droplet and @xmath52 is the energy of the ground state of the single - particle spectrum of a spherical cavity of radius @xmath41 . due to the small size of the droplet , even the lowest excited state has energy which is very large compared to the temperature . therefore all the @xmath39 pairs tend to occupy the ground state @xcite . next , entropies @xmath53 are given by @xmath54,\\ s_2 & = & s_1 + \frac{d_q v_d}{(2 \pi)^3 } \int^{\mu_q}_{\mu_q-\delta } 4 \pi p^2 dp \left [ f \ln f + ( 1-f ) \ln ( 1-f ) \right],\\ s_3 & = & 0.\end{aligned}\ ] ] @xmath55 can be taken to be zero for the reasons given above . finally , the pressures @xmath56 are calculated as follows : @xmath57 the thermodynamic potential @xmath58 for massless quarks is given by @xmath59.\ ] ] hence @xmath60 \nonumber \\ & = & \frac{e_1}{3v_d}-\frac{ktd_q}{2 \pi^2}\int_{\mu_q-\delta } ^{\mu_q } p^2 dp \ln\left[1 + \exp \beta ( \mu_q - p ) \right ] . \nonumber\end{aligned}\ ] ] performing the last integration by parts and rearranging the terms , one finally gets @xmath61 \}.\ ] ] the quantal pressure @xmath62 is given by @xmath63 it is easy to verify that if @xmath28 is set equal to zero , i.e. if no pairing is allowed to take place , @xmath45 , @xmath64 , @xmath62 vanish and @xmath44 , @xmath65 , @xmath66 , @xmath67 become equal to @xmath43 , @xmath68 , @xmath69 , @xmath70 , respectively . contribution of antiquarks to the above thermodynamic quantities is smaller by many orders of magnitude . although we have included it in our calculations its effect is negligible .    the droplet pressure @xmath71 is given by @xmath72 where @xmath67 and @xmath62 are as in eqs . ( 14)-(15 ) , @xmath73 and @xmath74 are the corresponding pressures due to antiquarks and @xmath75 is the pressure due to gluons in the droplet . consider once again the spherical region of radius @xmath41 . initially , i.e. before the pairing interaction is switched on , there are @xmath43 quarks , @xmath76 antiquarks and a certain number of gluons in it . in the final state , i.e. when the pairing interaction is on , there are @xmath44 quarks , @xmath51 quark pairs , @xmath77 antiquarks , @xmath78 antiquark pairs and the same number of gluons as before . the activation energy ( @xmath79 ) is defined as the difference between the helmholtz free energies ( @xmath80 ) in the final and initial states . since the gluonic contribution cancels out , one gets @xmath81 in the limit @xmath82 , this @xmath79 vanishes . the reader may wonder why surface tension of the droplet makes no appearance in our formalism . in the classical theory of nucleation , there are two competing contributions to the energy of a droplet : droplet formation reduces bulk free energy but only at the cost of surface energy . the former contribution wins if the droplet radius exceeds a certain critical value , while the surface tension ensures that a droplet with a smaller radius shrinks and disappears . in the formalism presented here , energy of the droplet is calculated quantum mechanically . it is clear from eq . ( 9 ) that small droplets have large quantal energy and tend to dissolve back into the qgp gas . this is similar to the effect of the surface tension in the classical theory . moreover , the thermal wavelength ( @xmath15 ) of the quarks is so large that a classification of particles as those localized on the surface and those in the interior breaks down . in short , we did not neglect the surface energy ; it is included in the total free energy of the quantal droplet . the precise boundaries of the various phases in fig . 1 are not known . since the baryon number density ( @xmath1 ) of the normal nuclear matter is @xmath83 @xmath17 , we have performed numerical calculations for @xmath84 0.5 to 3 @xmath17 and temperature @xmath85 0 to 100 mev . the resulting quark chemical potential @xmath24 varies over @xmath86 mev and the baryon chemical potential @xmath87 varies over @xmath88 mev . this fully covers the densities likely to be seen in the neutron - star cores . we have considered four different values of the superconducting gap parameter @xmath28 , namely , 1 , 10 , 50 and 100 mev .    since we are interested in a gas of weakly interacting quarks , antiquarks and gluons , we take @xmath18 to be zero . small values of @xmath18 do not affect the results in eqs . ( 1)-(3 ) significantly . we first choose a set of values for @xmath1 and @xmath0 , and solve eq . ( 1 ) numerically to determine @xmath24 . this @xmath24 is then used to evaluate @xmath89 as given by eq . also , once @xmath24 is determined , @xmath90 and @xmath91 can be evaluated separately by using @xmath92 ( above integrations can not be done analytically , and so it is not possible to get closed expressions like eq . ( 1 ) , for @xmath90 and @xmath93 . ) next , we choose a value for @xmath28 , and treating the droplet radius @xmath41 as a free parameter , calculate various other quantities ; see eqs . ( 4)-(17 ) . some of these results are shown in figs . 3 - 4 .    in our formalism , as in the classical theory , the activation energy when plotted as a function of the radius ( fig . 3 ) , initially increases from zero , attains a maximum , then starts decreasing monotonically eventually becoming negative for a sufficiently large value of the radius . in other words , the activation energy exhibits a barrier . _ the barrier originates as follows_. consider eq . as already stated , the antiquark terms are extremely small . the terms @xmath94 and @xmath95 are also negligible and in fact they tend to cancel each other . that leaves three terms , namely @xmath96 and @xmath64 . for small radii , @xmath65 and the quantal energy term @xmath64 dominate pushing the activation energy curve upwards . for large radii , the @xmath68 term dominates pulling the curve to the negative domain . as expected , the barrier becomes wider if @xmath0 is raised or if @xmath1 is lowered . the critical radii can be read off from fig . 3 . droplets with larger radii readily grow because that tends to reduce the activation energy . this in turn easily brings about the transition from qgp to the superconducting phase . figure 4 shows the ratio of the droplet pressure to the pressure of the surrounding gas , for @xmath97 mev . this ratio is significantly larger than unity for small droplets , which will tend to push the radii to larger values , thereby helping to bring about the qgp to superconducting phase transition . in the limit @xmath82 , this ratio tends to unity and is thus incapable of pushing the droplet radii to larger values . it is interesting to know the properties of the droplet when the pressure inside the droplet ( eq . 16 ) equals the pressure of the surrounding qgp gas ( eq . these results are shown in figs . the condition @xmath98 fixes @xmath41 . we solve the equation @xmath99 numerically to determine @xmath41 , for various values of @xmath1 , @xmath0 and @xmath28 . @xmath1 was varied from 0.5 to 3 @xmath17 in the steps of 0.5 @xmath17 . figure 5 shows the results for @xmath41 . as @xmath0 and/or @xmath1 are increased , @xmath41 decreases . however , @xmath41 is found to be rather insensitive to the variation in @xmath1 . although we have varied @xmath0 between 0 and 100 mev , temperatures relevant for neutron stars are close to zero . because of the small value of @xmath41 , the energy of even the lowest excited state inside the droplet turns out to be several hundreds of mev which is large compared to the temperature , effectively forcing all cooper pairs ( bosons ) to occupy the ground state . _ note also that at low temperatures , the droplet radius @xmath41 is larger than the pippard coherence length @xmath100 . _    figure 6 shows the activation energy @xmath79 . first see the upper half of fig . 6 . for a fixed @xmath1 , as @xmath0 increases , @xmath79 changes its sign from negative to positive . this is expected because , at higher @xmath0 , droplet formation is energetically unfavourable . secondly , as @xmath1 increases ( see curves labelled _ a _ to _ f _ ) , @xmath79 remains negative up to higher and higher @xmath0 . this is also expected because at higher @xmath1 , the transition to the superconducting phase is facilitated . now see the lower half of fig . 6 . it shows the @xmath101 contour , in the @xmath2 plane . this is also the phase transition curve separating the qgp phase at high @xmath0 and the superconducting phase at low @xmath0 : recall @xmath102 where @xmath80 are the free energies of the two phases . both @xmath103 and @xmath104 decrease as @xmath0 increases and cross each other at some temperature , say @xmath105 ( fig . since the system follows the path of minimum @xmath80 , there is a kink in @xmath106 . in other words , the slope @xmath107 has a discontinuity , indicating a first - order phase transition at @xmath105 . thus the @xmath101 contour serves to demarcate the region of the superconducting phase . note that as @xmath1 increases , the @xmath108 region in fig . 6 survives up to higher and higher @xmath0 . we now present an estimate of the latent heat ( @xmath109 ) associated with the first - order phase transition considered in this paper . the transition temperature @xmath105 for any given @xmath1 can be determined from fig . the latent heat is then estimated as @xmath110 . for @xmath97 mev , this turns out to be @xmath111 to @xmath30 mev and it increases with @xmath1 . in the limit @xmath82 , the two free - energy curves in fig . 7 coincide , @xmath79 vanishes , @xmath112 , there is no discontinuity in entropy , and the latent heat is zero . there is only one phase namely qgp , and there is no question of any phase transition . finally , fig . 8 shows 3-dimensional surface plots of the ( normalized ) activation energy @xmath113 defined as @xmath114 as a function of @xmath1 and @xmath0 , for @xmath115 and 100 mev . it is clear that in the limit @xmath116 , the activation energy surface is flat and there is no energy advantage for the qgp to superconducting phase transition . indeed there would be only one phase , namely qgp , which would extend up to the @xmath117 axis . on the other hand , for @xmath97 mev , the surface is such that the system would like to `` roll down '' the slope thereby bringing about the phase transition . it is interesting to compare the activation energy results for the qcd phase transition , obtained here with those for the bose - einstein condensation in a gas of trapped atomic hydrogen , presented in @xcite . in @xcite two gas densities were considered : @xmath118 @xmath8 and @xmath119 @xmath8 both at the same temperature @xmath120 . we have repeated their calculations for @xmath121 . their activation energy exhibits a barrier , and they denote the maximum value of the activation energy by @xmath122 . consider the ratio @xmath123 at the peak position and the corresponding number @xmath124 of particles in the droplet . in @xcite these two quantities are highly sensitive to the changes in the gas density and temperature . for the inset , it is @xmath125 a.u . ] for instance , for the above tiny change in the gas density , the ratio @xmath123 changes from 4.8 to 104.1 , and @xmath124 changes from 18 to 3862 . similar sensitivity was seen when @xmath0 was changed by a few @xmath126 . in contrast , the results presented here are relatively stable with respect to changes in @xmath1 and @xmath0 . in conclusion , we have studied hitherto unexplored first - order phase transition from qgp to a colour - superconducting phase , and have presented a simple model of the transition with very few assumptions . it provides numerical support to the prevalent ideas of phase structure of qcd in the low - temperature , high - baryon - number - density regime . we have studied properties of the droplets of the superconducting phase , formed at various points in the @xmath2 plane , for several values of the gap parameter . we have been able to estimate the extent of the superconducting phase in this plane and the latent heat of the transition . at low temperatures , droplet radii are found to be larger than the pippard coherence length which measures the spatial extent of the @xmath39 cooper pairs . rsb acknowledges the hospitality of the mcmaster university where this work was initiated . we thank mark alford and mustansir barma for their comments on this work . this research was partially supported by a grant from nserc ( canada ) .      if @xmath127 is the distribution function of a gas in equilibrium , the number density @xmath13 is given by @xmath128 . for a nonrelativistic classical gas at temperature @xmath0 , one has @xmath129,\ ] ] where @xmath130 and @xmath33 is the chemical potential . hence @xmath131 where @xmath15 is the thermal wave length .    for an extreme relativistic gas of massless classical particles , @xmath132,\ ] ] and hence @xmath133 we define the thermal wave length in this case as @xmath134 so that the fugacity @xmath135 once again equals @xmath136 . schematic phase diagrams of qcd : @xmath0 is the temperature and @xmath137 is the baryon chemical potential . @xmath138 is the critical point at which the line of first - order phase transitions ends . @xmath139 and @xmath140 ( colour - flavour - locked ) are 2-flavour - like and 3-flavour - like colour superconducting phases , respectively . the dashed line denotes the critical temperature at which quark - quark pairing vanishes . the thick arrow indicates the first - order phase transition considered in this paper . for the normal nuclear matter at @xmath117 , @xmath137 should be defined as @xmath141 where @xmath142 @xmath143 @xcite and @xmath144 is the nucleon mass , which gives @xmath145 mev . ( figure adapted from @xcite . ) activation energy @xmath79 vs droplet radius @xmath41 for various gap values @xmath28 . baryon number density @xmath1 ( @xmath17 ) and temperature @xmath0 ( mev ) are : ( a ) 0.5 and 5 ( b ) 0.5 and 100 ( c ) 3 and 5 ( d ) 3 and 100 , respectively . upper panel : activation energy @xmath79 vs temperature @xmath0 . baryon number densities @xmath1 ( @xmath17 ) are ( a ) 0.5 ( b ) 1 ( c ) 1.5 ( d ) 2 ( e ) 2.5 ( f ) 3 . lower panel : the @xmath101 contour or the phase transition line in the @xmath2 plane . sc : superconducting phase .
Abstract:

Answer: considering the low - temperature ( @xmath0 ) and high - baryon - number - density ( @xmath1 ) region of the qcd phase diagram , we present a model for the first - order phase transition between the quark - gluon plasma ( qgp ) and the recently proposed colour superconducting phase . we study nucleation of a droplet of the superconducting phase within the metastable qgp gas . numerical results for the activation energy , radius and other physical parameters of the droplets , at various temperatures , densities and gap parameters , are given . we have estimated the latent heat of the phase transition . in the @xmath2 plane , we are able to demarcate the region of the superconducting phase . pacs numbers : 12.38.mh , 12.38.aw , 26.60.+c , 64.60.qb    _ keywords _ : quark - gluon plasma , qcd phase transition , colour superconductivity , droplet formation , diquark    ' '' ''    0.3 cm e - mail : bhalerao@theory.tifr.res.in , bhaduri@smiley.physics.mcmaster.ca


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: exact solutions of the einstein equations are associated to highly idealized physical systems that have some exceptional geometrical properties . in some cases with a simple exact solution one can capture a significant part of the physical properties of non trivial systems . also , in nonlinear theories like general relativity and fluid dynamics the exact solutions play an important role in numerical analysis . these solutions can be used to test numerical codes and its outcomes . also they can be employed as initial conditions to describe more realistic situations , e.g. , a static solution can be used as part of the initial conditions for a full dynamical simulation . since the natural shape of an isolated self - gravitating fluid is axially symmetric , the solutions of einstein s field equations with this symmetry play a particularly important role in the astrophysical applications of general relativity . in particular disk like configurations of matter are of great interest , since they can be used as models of galaxies or accretion disks . also these disks can be used as starting point to represent more realistic models in which the bulge and halo of the galaxy are considered . solutions for static thin disks without radial pressure were first studied by bonnor and sackfield @xcite , and morgan and morgan @xcite , and with radial pressure by morgan and morgan @xcite . the first solution represent disks made of pressureless dust whereas the second disks with azimuthal pressure but without radial pressure . and the third a disk made of an anisotropic fluid with nonzero radial pressure . the bonnor - sackfield disk has a singular rim . these disks are finite . several classes of exact solutions of the einstein field equations corresponding to static @xcite and stationary @xcite thin disks have been obtained by different authors , with or without radial pressure . thin disks with radial tension @xcite , magnetic fields @xcite and magnetic and electric fields @xcite have been also studied . the non linear superposition of a disk and a black hole was first considered by lemos and letelier @xcite . this solution and its generalizations has been studied in some detail in @xcite . recently the stability of circular orbits of particles moving around black holes surrounded by axially symmetric structures have been considered in @xcite . for a recent survey on relativistic gravitating thin disks , see @xcite .    except for the pressureless disks all the other disks have as source matter with azimuthal pressure ( tension ) different from the radial pressure ( tension ) . however , in some cases these disks can be interpreted as the superposition of two counter - rotating perfect fluids . a detailed study of the counter - rotating model for the case of static thin disks is presented in @xcite . recently , more realistic models of thin disks and thin disks with halos made of perfect fluids were considered in @xcite .    in all the disks mentioned above an inverse style method was used to solve the einstein equations . the metric representing the disk is guessed and then it is used to compute the source ( energy - momentum tensor ) . this method was named by synge as the g - method @xcite in contrast to the t - method or direct method in which the source is given and the einstein equations are solved . the t - method , has been used to generate disks by the jena group @xcite . essentially , they are obtained by solving a riemann - hilbert problem . these solutions are highly non trivial , but they deserve special attention because of their clear physical meaning .    in the solutions obtained by the g - method the well known `` displace , cut and reflect '' method is used . the idea of the method is simple . given a solution of the vacuum einstein equations , a cut is make above all singularities or sources . the identification of this solution with its mirror image yields relativistic models of disks . in general , these disks are of infinite extension and finite mass . the aim of this paper is to consider disks beyond the thin disk limit to add a new degree of reality to these geometric models of galaxies . even though in first approximation the galactic disks can be considered to be very thin , e.g. , in our galaxy the radius of the disk is 10@xmath0 and its thickness is 1@xmath0 . in a more realistic model the thickness of the disk need to be considered . also it is well known in fluid mechanics that the addition of a new dimension can make dramatic changes in the dynamics of the fluid . in principle , this new dimension will also change the dynamical properties of the the disk source , e.g. , its stability . in this paper we generalize the `` displace , cut and reflect '' method in order to obtain thick disks models from vacuum solutions of einstein equations . we shall replace the surface of discontinuity of the metric derivatives with a `` thick '' shell in such a way that the matter content of the disk will be described by continuous functions with continuous first derivatives . this generalization can be named `` displace , cut , _ fill _ and reflect '' method . the disks obtained with this method , in general , will be of infinite extension and finite mass . also as in the case of thin disks , the matter that form the thick disks will not obey simple equations of state and in some regions of the disk the pressure can change sign given rise to tensions . although the energy condition will be fulfilled . the models of thick disks presented can be considered as generalizations of models of thin disks studied in references @xcite and @xcite . the article is divided as follows . in sec . ii we present , in some detail , the main idea of the `` displace , cut , _ fill _ and reflect '' method in newtonian gravity . the method is then applied , in sec . iii , to construct relativistic thick disks in weyl coordinates . we also study the general expression for the energy - momentum tensor of the disks . the method is illustrated by taking two simple weyl solutions that lead to thick disks in agreement with all the energy - conditions . in sec . iv we apply the method to the schwarzschild solution in isotropic cylindrical coordinates . the disk obtained also satisfy all the energy conditions , this disk have equal azimuthal and radial pressures and different vertical pressure . finally , in sec . v , we summarize our main results . @xmath1    @xmath2    @xmath3    the newtonian gravitational potential of a thin disk can be obtained by a simple procedure , the `` displace , cut and reflect method '' , due to kuzmin @xcite . the method can be divided in the following steps : first , choose a surface that divides the usual space in two parts : one with no singularities or sources an the other with the sources . second , disregard the part of the space with singularities . third , use the surface to make an inversion of the nonsingular part of the space . the result will be a space with a singularity that is a delta function with support on @xmath4 . this procedure is depicted in fig . [ fig : thind ] .    in order to obtain a thick disk we need to modify the above procedure . we essentially need to replace the surface of discontinuity with a `` thick '' shell in such a way that the matter content of the disk be described by continuous functions . now method has an additional step and can be named `` displace , cut , _ fill _ and reflect '' . after we disregard the part of the space with singularities , we put a thick shell below the surface . then we use the bottom surface of the shell to make the inversion . the procedure is illustrated in fig . [ fig : thickd ] . mathematically the method is equivalent to make the transformation @xmath5 , where @xmath6 is a constant and @xmath7 an even function of @xmath8 . in newtonian gravity , the potential @xmath9 is a solution of the laplace equation @xmath10 where @xmath11 are the usual cylindrical coordinates . after we make the transformation @xmath12 , the above equation leads to @xmath13 \phi_{,hh } , \label{eq : laph}\ ] ] where primes indicate differentiation with respect to @xmath8 . for the case of thin disks we take @xmath14 . note that @xmath15 and @xmath16 where @xmath17 is the heaveside function and @xmath18 the usual dirac distribution . by using ( [ eq : laph ] ) , the poisson equation leads to a mass density given by @xmath19 we have a surface distribution of matter located in the plane @xmath4 . in fig [ fig : figh ] , with dashed lines , we plot @xmath7 an its derivatives .    for the case of thick disks the function @xmath20 must be selected in such a way that @xmath21 and its first derivatives be continuous across the plane @xmath22 . let us take a function @xmath7 defined as @xmath23 hence , by taking the function @xmath20 above defined we can generate disks of thickness @xmath24 located in the region @xmath25 . in fig . [ fig : figh ] we plot @xmath7 an its derivatives ( solid lines ) . when @xmath26 the function @xmath7 is a linear function of @xmath8 such that @xmath27 . hence , its second derivative is zero . then the mass density vanish outside the disk . since the first derivative is continuous at @xmath28 ( see fig . [ fig : figh ] ) and the second derivative piecewise constant we have that the mass density , @xmath29 , will be well defined inside the disk , @xmath30 for @xmath31 and @xmath32 for @xmath33 .    as a simple example we can consider the usual potential for a mass point , written in cylindrical coordinates as @xmath34 by doing the transformation @xmath12 in the previous potential , we obtain the mass densities @xmath35 for a thin disk ( kuzmin - toomre disk @xcite ) , and @xmath36}{8 \pi a^2 { \rm r}^3 }   \nonumber \\ & &   \label{thicknden}\\ & + & \frac{3m(a^2 - z^2)(z^2 + 2ab)^2}{16 \pi a^4 { \rm r}^5 } , \nonumber \end{aligned}\ ] ] where @xmath37 , for a thick disk . when @xmath38 , the mass density will be positive everywhere . the function @xmath20 presented in ( [ eq : funh ] ) is the simplest function that have the desired properties . the part of the function in the domain @xmath39 can be changed by superposition of even functions of @xmath8 . in @xmath28 this new function need to be matched continuously with linear functions of @xmath8 such that @xmath27 , like in ( [ eq : funh ] ) . in this case the mass density can depend on the variable @xmath8 in a more general way . it is instructive to obtain the surface thin disk density ( say @xmath40 ) associated to ( [ thinnden ] ) , @xmath41 as a limit of the thick disk volume density ( [ thicknden ] ) . to obtain the surface density of the thin disk from the volume density ( [ thinnden ] ) we first do @xmath42 now to perform the thin disk limit we put @xmath43 and @xmath44 ( @xmath45 and @xmath46 artitrary constants ) in @xmath47 and we take the limit @xmath48 this limit is just ( [ sigma ] ) . when the matter is absent , the metric for a static axially symmetric spacetime can be cast without loosing generality as @xmath49 , \label{eq : metw}\ ] ] where @xmath21 and @xmath50 are functions of @xmath51 and @xmath8 only . the ranges of the coordinates @xmath52 are the usual for cylindrical coordinates ( weyl coordinates ) and @xmath53 . the einstein vacuum equations for this metric yield the weyl equations @xcite ,    @xmath54    [ eq : eqsw ]    from a solution of einstein vacuum equations corresponding to a weyl metric ( [ eq : metw ] ) , we can construct a thick disk model by means of the `` displace , cut , fill and reflect method '' using the transformation @xmath55 , with @xmath7 given by ( [ eq : funh ] ) . the energy - momentum tensor of the disk can be computed using the einstein equations in the matter , written as @xmath56 in units such that @xmath57 .    by using the einstein equations ( [ eq : eqsw ] ) , the nonzero components of @xmath58 are :    @xmath59 , \nonumber \\ \nonumber \\ t_\varphi^\varphi & = & \frac{e^{2(\phi - \lambda)}}{a^2 } [ a \lambda_{,h } + ( z^2 - a^2 ) ( \phi_{,h}^2 + \lambda_{,hh } ) ] , \\      \nonumber    \\ t_r^r & = & \frac{e^{2(\phi - \lambda)}}{a^2 } ( z^2 - a^2 ) \phi_{,h}^2 \ , \\ \nonumber    \\ t_z^z & = & \frac{e^{2(\phi - \lambda)}}{a^2 } ( a^2 - z^2 ) \phi_{,h}^2 \ , \end{aligned}\ ] ]    [ eq : emtw ]    valid for the region @xmath60 . outside this region we have @xmath61 .    for the above expressions we can see that the radial stress @xmath62 is negative ( we have radial tension ) . on the other hand , since @xmath63 , we have vertical pressure . defining the orthonormal tetrad @xmath64 , where    @xmath65    [ eq : tetrad ]    we can be cast the energy - momentum tensor in its canonical form @xmath66 here @xmath67 is the energy density , @xmath68 is the azimuthal stress , @xmath69 is the radial tension and @xmath70 is the vertical pressure .    from ( [ eq : emtw ] ) we get the `` effective newtonian '' density , @xmath71 @xmath72 . \label{eq : dnew}\ ] ] the strong energy condition requires that @xmath73 , whereas the weak energy condition imposes the condition @xmath74 . the dominant energy condition is equivalent to the requirement @xmath75 , @xmath76 and @xmath77 , see for instance @xcite . one can obtain the quantities associated to a general relativistic thin disk from the corresponding quantities associated to the thick disk using the same limit procedure described in the newtonian disk case .      as a first example we apply the `` displace , cut , fill and reflect '' method to obtain thick disks using the chazy - curzon solution @xcite , written in weyl coordinates as    @xmath78    where @xmath37 , @xmath79 and @xmath6 are positive constants and @xmath80 is given by ( [ eq : funh ] ) . now we rescale the variables and the parameters in terms of the disk thickness , @xmath81 . we @xmath82 , @xmath83 , @xmath84 , @xmath85 and @xmath86 . from ( [ eq : emtw ] ) and ( [ eq : dnew ] ) , we obtain    @xmath87 , \label{eq : dncc } \\ & &   \nonumber    \\ \tilde\epsilon & = & \frac{\tilde m \ e^{2(\phi - \lambda)}}{4 \tilde{\rm r}^8 } \left [ 4 ( 3 { \tilde z}^2 + 2 ( \tilde b - 1 ) ) ( \tilde{\rm r}^3 - \tilde m { \tilde r}^2 ) \tilde{\rm r}^2 + [ 6 ( \tilde{\rm r}^3 - 2 \tilde m { \tilde r}^2 ) + \tilde m \tilde{\rm r}^2 ) ( 1 - { \tilde z}^2 ) ( { \tilde z}^2 + 2 \tilde b)^2 \right ] , \label{eq : encc } \\ { \tilde p}_\varphi & = & \frac{{\tilde m}^2 e^{2(\phi - \lambda)}}{4 \tilde{\rm r}^8 } \left [ 4 { \tilde r}^2 [ ( 3 { \tilde z}^2 + 2 ( \tilde b - 1 ) ) \tilde{\rm r}^2 + 3 ( 1 - { \tilde z}^2 ) ( { \tilde z}^2 + 2 \tilde b)^2 ] + ( { \tilde z}^2 - 1 ) ( { \tilde z}^2 + 2 \tilde b)^2 \tilde{\rm r}^2 \right ] , \label{eq : p1cc } \\ & &   \nonumber    \\ { \tilde p}_r & = & \frac{{\tilde m}^2 e^{2(\phi - \lambda)}}{4 \tilde{\rm r}^6 } \left [ ( { \tilde z}^2 - 1 ) ( { \tilde z}^2 + 2 \tilde b)^2 \right ] , \label{p2cc } \\ & &   \nonumber    \\ { \tilde p}_z & = & \frac{{\tilde m}^2 e^{2(\phi - \lambda)}}{4 \tilde{\rm r}^6 } \left [ ( 1 - { \tilde z}^2 ) ( { \tilde z}^2 + 2 \tilde b)^2 \right ] , \label{p3cc}\end{aligned}\ ] ]    where @xmath88 , @xmath89 and @xmath90 .    equation ( [ eq : dncc ] ) and the condition @xmath91 imply @xmath92 . then , when @xmath91 the strong energy condition is satisfied . from ( [ eq : encc ] ) we conclude that @xmath93 whenever @xmath94 and @xmath91 . from the definition of @xmath95 we conclude that in order to have @xmath93 everywhere , we need @xmath96 and @xmath91 .    @xmath97    @xmath98    the behavior of the densities is better illustrated graphically . in order to have a disk in agreement with the weak and strong energy conditions , we take @xmath99 and @xmath100 . in fig . [ fig : dens1 ] we plot the effective newtonian density and the energy density in units of @xmath101 , @xmath102 and @xmath103 , as functions of @xmath104 and @xmath105 . we can see that the densities have a maximum at the center of the @xmath4 plane . then the densities decrease monotonously as @xmath51 increase and also the densities decrease for @xmath106 . we see that @xmath29 and @xmath107 have similar magnitudes .    the azimuthal stress , as we can see from ( [ eq : p1cc ] ) , is negative at the center of the disk ( @xmath108 ) , whereas it is positive for larges values of @xmath51 . the boundary between the region of negative stress and positive stress is the surface @xmath109 with @xmath110 given by ( [ eq : p1cc ] ) . the behavior of the stresses is also better illustrated graphically . in fig . [ fig : pres1 ] we plot the dimensionless azimuthal stress @xmath111 and vertical pressure @xmath112 , as functions of @xmath104 and @xmath105 . the radial pressure is given by @xmath113 . again we take @xmath99 and @xmath100 . the azimuthal stress is negative at the central region of the disk , then increases to have a positive maximum on the @xmath114 plane , for a value of @xmath115 . finally , it decreases monotonously for increasing @xmath51 and also for @xmath106 . the behavior of the vertical pressure is like of the densities , with a maximum at the center of the @xmath114 plane and then monotonously decreasing for increasing @xmath51 . also , @xmath116 for @xmath117 . from figs . [ fig : dens1 ] and [ fig : pres1 ] we can also see that the magnitude of the stresses is about a tenth of the magnitude of the densities . we have @xmath118 and so the disks are also in agreement with the dominant energy condition . thin disks based on the chazy - curzon metric were studied in ref . @xcite .      as a second example we take the weyl gamma metric also known as zipoy - voorhees solution @xcite , that in weyl coordinates can be cast as @xcite    @xmath119 , \\ & &   \nonumber    \\ \lambda & = & \frac{m^2}{2 k^2 } \ln \left [ \frac{({\rm r}_1 + { \rm r}_2)^2 - 4 k^2}{4 { \rm r}_1 { \rm r}_2 } \right],\end{aligned}\ ] ]    where @xmath120 , @xmath121 , @xmath79 , @xmath6 and @xmath122 are positive constants and @xmath7 is given by ( [ eq : funh ] ) . when @xmath123 this solution leads to the schwarzschild metric and when @xmath124 to the chazy - curzon solution of the previous section .    by using ( [ eq : emtw ] ) and ( [ eq : dnew ] ) we obtain , in terms of the dimensionless variables used in the previous section ,    @xmath125 \right\ } , \label{eq : dnzv } \\ & &   \nonumber    \\ \tilde\epsilon & = & \tilde\rho \left [ 1 - \frac{2 \tilde m { \tilde r}^2 ( \tilde{\rm r}_1 + \tilde{\rm r}_2)}{\tilde{\rm r}^2_1 \tilde{\rm r}^2_2 [ ( \tilde{\rm r}_1 + \tilde{\rm r}_2)^2 - 4 { \tilde k}^2 ] } \right ] + \frac{{\tilde m}^2 e^{2(\phi - \lambda)}(1 - { \tilde z}^2)(\tilde{\rm r}_1 - \tilde{\rm r}_2)}{4 { \tilde k}^2 \tilde{\rm r}^4_1 \tilde{\rm r}^4_2 }   \left[(\tilde{\rm r}_1 - \tilde{\rm r}_2 ) \tilde{\rm r}^2_1 \tilde{\rm r}^2_2 - 2 { \tilde r}^2 ( \tilde{\rm r}^3_1 - \tilde{\rm r}^3_2 ) \right ] , \label{eq : enzv } \\ & &   \nonumber    \\ { \tilde p}_\varphi & = & \frac{2 \tilde m \tilde\rho { \tilde r}^2 ( \tilde{\rm r}_1 + \tilde{\rm r}_2)}{\tilde{\rm r}^2_1 \tilde{\rm r}^2_2 [ ( \tilde{\rm r}_1 + \tilde{\rm r}_2)^2 - 4 { \tilde k}^2 ] } + \frac{{\tilde m}^2 e^{2(\phi - \lambda)}({\tilde z}^2 - 1)(\tilde{\rm r}_1 - \tilde{\rm r}_2)}{4 { \tilde k}^2 \tilde{\rm r}^4_1 \tilde{\rm r}^4_2 }   \left[(\tilde{\rm r}_1 - \tilde{\rm r}_2 ) \tilde{\rm r}^2_1 \tilde{\rm r}^2_2 - 2 { \tilde r}^2 ( \tilde{\rm r}^3_1 - \tilde{\rm r}^3_2 ) \right ] , \label{eq : p1zv } \\ & &   \nonumber    \\ { \tilde p}_r & = & \frac{{\tilde m}^2 e^{2(\phi - \lambda)}}{\tilde{\rm r}^2_1 \tilde{\rm r}^2_2 } \left [ \frac{({\tilde z}^2 - 1 ) ( { \tilde z}^2 + 2 \tilde b)^2}{(\tilde{\rm r}_1 + \tilde{\rm r}_2)^2 } \right ] , \\ \label{p2zv } & &   \nonumber    \\ { \tilde p}_z & = & \frac{{\tilde m}^2 e^{2(\phi - \lambda)}}{\tilde{\rm r}^2_1 \tilde{\rm r}^2_2 } \left [ \frac{(1 - { \tilde z}^2 ) ( { \tilde z}^2 + 2 \tilde b)^2}{(\tilde { \rm r}_1 + \tilde { \rm r}_2)^2 } \right ] , \label{p3zv}\end{aligned}\ ] ]    where @xmath126 , @xmath127 and @xmath128 .    from ( [ eq : dnzv ] ) and the condition @xmath129 we have @xmath130 , i.e. , the strong energy condition is satisfied . also from @xmath131 , with @xmath132 , we can show that @xmath133 also @xmath134 and @xmath135 give us @xmath136 then the condition @xmath137 can be cast as @xmath138 that leads also to @xmath139 . this last condition assure the nonsingular behavior of the energy density @xmath107 and the azimuthal stress @xmath140 .    from ( [ eq : enzv ] ) and @xmath108 we have that @xmath141 . when @xmath142 the condition @xmath143 } \leq 1 , \ ] ] gives us @xmath141 . since @xmath144 , @xmath145 , @xmath146 and @xmath147 we have , @xmath148 then the condition @xmath149 leads to @xmath150 this relation yields also @xmath151 as a condition to have @xmath152 . @xmath153    @xmath154    for any other value of @xmath51 and @xmath8 is not easy to obtain constraints over the parameters @xmath155 , @xmath156 and @xmath157 in order to have @xmath158 . the analysis is better done graphically . by considering different values of @xmath155 and @xmath156 that fulfill the condition ( [ eq : cond1 ] ) we find that @xmath141 everywhere in the disks only if we take for @xmath157 a value less than a tenth of the upper limit provided by the condition ( [ eq : cond2 ] ) . as an example , in fig . [ fig : dens2 ] we plot the dimensionless densities @xmath159 and @xmath160 for a disk with @xmath161 , @xmath162 and @xmath163 . we can see that , as in the chazy - curzon disk , the density have a maximum at the center of the @xmath4 plane and then it decreases monotonously as @xmath51 increase and also for @xmath106 . also @xmath29 and @xmath107 have similar magnitudes . the behavior of the stresses is also similar to the presented in the chazy - curzon disk . again , is better to do a graphical presentation . in fig . [ fig : pres2 ] we plot the dimensionless azimuthal stress @xmath111 and the vertical pressure @xmath112 for the disk with @xmath164 , @xmath162 and @xmath163 . again we have that @xmath165 . as with the chazy - curzon disk , the azimuthal stress is negative at the central region of the disk , then increase to have a positive maximum at the @xmath166 plane , for a value of @xmath167 , and finally decrease monotonously . the behavior of the vertical pressure is like of the densities , with a maximum at the center of the @xmath4 plane and then monotonously decreasing for increasing @xmath51 . also , @xmath168 for @xmath169 . from fig . [ fig : dens2 ] and [ fig : pres2 ] we have @xmath170 and so the disks are also in agreement with the dominant energy condition . thin disks based in the zipoy - voorhees metric were considered in ref . for a static spherically symmetric spacetime the metric in isotropic spherical coordinates @xmath171 can be cast as @xmath172 , \label{eq : mets}\ ] ] where @xmath21 and @xmath50 are functions of @xmath95 only . in isotropic cylindrical coordinates @xmath173 the metric ( [ eq : mets ] ) takes the form @xmath174 , \label{eq : meti}\ ] ] where now @xmath21 and @xmath50 depends on @xmath51 and @xmath8 . we will now apply the `` displace , cut , fill and reflect '' method to the schwarzschild solution in isotropic coordinates ,    @xmath175 , \\ & &      \nonumber    \\ \lambda & = & \ln \left [ 1 + \frac{m}{2{\rm r } } \right]^2 , \end{aligned}\ ] ]    [ eq : miso ]    with @xmath79 a positive constant and @xmath176 . now we put @xmath177 where @xmath6 is a positive constant and @xmath20 is given by ( [ eq : funh ] ) .    from the einstein equations ( [ eq : eqm ] ) and the orthonormal tetrad @xmath178 where    @xmath179    [ eq : tetrad2 ]    we find that the energy - momentum tensor of the disk can be written as @xmath180 here @xmath67 is the energy density , @xmath68 is the azimuthal stress , that is equal to the radial stress @xmath181 , and @xmath70 is the vertical stress . the effective newtonian density is given by @xmath182 .    from ( [ eq : miso ] ) we obtain , using the dimensionless variables previously defined ,    @xmath183}{(2\tilde{\rm r } - \tilde m)(2\tilde{\rm r } + \tilde m)^5 } , \label{eq : denis } \\ & &      \nonumber    \\ \tilde\epsilon & = & \frac{3 \tilde m \left[2 ( 3{\tilde z}^2 + 2(\tilde b - 1 ) ) \tilde{\rm r}^2 + 3 ( { \tilde z}^2 + 2\tilde b)^2 ( 1 - { \tilde z}^2)\right]}{2(2\tilde{\rm r } + \tilde m)^5 } , \label{eq : eneis } \\ & &      \nonumber    \\ { \tilde p}_\varphi & = & \frac{16 { \tilde m}^2 \left[(3{\tilde z}^2 + 2(\tilde b - 1 ) ) \tilde { \rm r}^2 + ( { \tilde z}^2 + 2\tilde b)^2 ( 1 - { \tilde z}^2)\right]}{(2\tilde{\rm r } - \tilde m)(2\tilde { \rm r } + \tilde m)^5 } , \label{eq : p1is } \\ & &      \nonumber    \\ { \tilde p}_r & = & \frac{16 { \tilde m}^2 \left[(3{\tilde z}^2 + 2(\tilde b - 1 ) ) \tilde{\rm r}^2 + ( { \tilde z}^2 + 2\tilde b)^2 ( 1 - { \tilde z}^2)\right]}{(2\tilde{\rm r } - \tilde m)(2\tilde{\rm r } + \tilde m)^5 } , \label{eq : p2is } \\ & &      \nonumber    \\ { \tilde p}_z & = & \frac{16 { \tilde m}^2 ( { \tilde z}^2 + 2\tilde b)^2 ( 1 - { \tilde z}^2)}{(2\tilde{\rm r } - \tilde m)(2\tilde{\rm r } + \tilde m)^5 } , \label{eq : p3is}\end{aligned}\ ] ]    from ( [ eq : eneis ] ) and @xmath184 we have @xmath141 . on the other hand , from ( [ eq : denis ] ) and @xmath185 we have @xmath92 . since @xmath186 , this last condition is equivalent to @xmath187 . therefore , when @xmath91 and @xmath188 we will have disks in agreement with the weak and strong energy conditions . also this values of @xmath157 assure the nonsingular behavior of @xmath29 , @xmath140 , @xmath189 and @xmath190 . we also have @xmath191 and @xmath192 . the vertical and horizontal stress are then pressures . @xmath193    @xmath194    as in the previous section , we perform a graphical analysis of the solution . in fig . [ fig : dens3 ] we plot @xmath159 and @xmath160 for a thick disk obtained from the schwarzchild isotropic solution with @xmath195 . the horizontal and vertical pressures are plotted in fig [ fig : pres3 ] . all the four quantities have a similar behavior , with a maximum at the center of the @xmath4 plane and then monotonously decreasing with increasing @xmath51 and @xmath8 . the relative magnitudes of the densities and pressures are such that @xmath196 . we have @xmath197 and so the disks are in agreement with all the energy conditions . thin disks based on the schwarzschild solution in isotropic coordinates were studied in @xcite , whereas thin disks based in the schwarzschild metric in weyl coordinates were studied in @xcite . for the disks in isotropic coordinates we have matter with radial pressure equal to the azimuthal pressure ( isotropic matter ) an for the disks in weyl coordinates we have zero radial pressure . we presented a method to obtain exact general relativistic thick disks as a generalization of the `` displace , cut and reflect '' method commonly used to obtain newtonian and relativistic thin disks . the generalization was be done by means of the transformation @xmath12 , where @xmath7 is an even function of @xmath8 and @xmath6 is a positive constant . the function @xmath20 must be selected in such a way that the metric tensor and and its first derivatives will be continuous across the plane @xmath4 . all the cases considered leads to thick disks with similar behavior of the energy and newtonian effective densities : a maximum at the center of the central plane of the disks , the @xmath4 plane , and then monotonously decreasing for increasing @xmath51 and @xmath8 . we found that , when the method is applied to vacuum weyl spacetimes , the thick disks presents radial tension and vertical pressure . the azimuthal stress is negative at the central region of the disks , then have a positive maximum and finally decrease monotonously for large values of @xmath51 and @xmath106 , where @xmath24 is the thickness of the disk . the disks obtained are in full agreement with all the energy conditions .    on the other hand , when the method is applied to the schwarzschild isotropic metric all the stresses are pressures and have a behavior like the densities . the disks obtained are also in full agreement with all the energy conditions . we plan to extend the models of thick disk presented along this lines by considering more elaborate functions @xmath20 and by the incorporation of new properties like rotation , either electric or magnetic fields or both . also we believe that the study of stability in this disks can produce some non trivial results . o. semerk , _ towards gravitating disks around stationary black holes_. in _ ` gravitation : following the prague inspiration ' to celebrate the 60th birthday of jiri bick _ , eds . o. semerk , j. podolsky and m. zofka , ( world scientific , singapore , 2002 ) , p. 111 .
Abstract:

Answer: a method to construct exact general relativistic thick disks that is a simple generalization of the `` displace , cut and reflect '' method commonly used in newtonian , as well as , in einstein theory of gravitation is presented . this generalization consists in the addition of a new step in the above mentioned method . the new method can be pictured as a `` displace , cut , _ fill _ and reflect '' method . in the newtonian case , the method is illustrated in some detail with the kuzmin - toomre disk . we obtain a thick disk with acceptable physical properties . in the relativistic case two solutions of the weyl equations , the weyl gamma metric ( also known as zipoy - voorhees metric ) and the chazy - curzon metric are used to construct thick disks . also the schwarzschild metric in isotropic coordinates is employed to construct another family of thick disks . in all the considered cases we have non trivial ranges of the involved parameter that yield thick disks in which all the energy conditions are satisfied .


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 study of the equation of state ( eos ) , transport properties , and mixing rules of hydrogen ( h ) and helium ( he ) under extreme condition of high pressure and temperature is not only of fundamental interest but also of essential practical applications for astrophysics @xcite . for instance , giant planets such as jupiter and saturn require accurate eos as the basic input into the respective interior models in order to solve hydrostatic equation and investigate the solubility of the rocky core @xcite . on the other side , the evolution of stars and the design of thermal protection system is assisted by high precision transport coefficients of h - he mixtures at high pressure @xcite . in addition , the viscosity and mutual diffusion coefficients are also important input properties for hydrodynamic simulations in modelling the stability of the hot spot - fuel interfaces and the degree of fuel contamination in inertial confinement fusion ( icf ) @xcite .    since direct experimental access such as shock wave experiments is limited in the mbar regime @xcite , the states deep in the interior of jupiter ( @xmath2 mbar ) and saturn ( @xmath3 mbar ) @xcite can not be duplicated in the laboratory . as a consequence , theoretical modelling provides most of the insight into the internal structure of giant planets . the eos of h - he mixtures have been treated by a linear mixing ( lm ) of the individual eos via fluid perturbation theory @xcite and monte carlo simulations @xcite . recently , several attempts have been made to calculate eos of h - he mixtures by means of quantum molecular dynamic ( qmd ) simulations . @xcite applied local density approximation of density functional theory ( lda - dft ) calculations for solid h - he mixtures , implying demixing for jupiter and saturn at 15000 k for a he fraction of @xmath4 . @xcite , lorenzen _ et al . _ @xcite , and militzer @xcite performed qmd simulations by using generalized gradient expansion ( gga ) instead of the lda in order to evaluate the accuracy of the lm approximation and study the demixing of h - he at mbar pressures . @xcite introduced car - parrinello molecular dynamics ( cpmd ) simulations to calculate the excess gibbs free energy of mixing at a lower temperature compared to that set in the work of lorenzen _ et al . _ @xcite and militzer @xcite .    considering the transport properties , qmd @xcite and orbital - free molecular dynamics ( ofmd ) @xcite simulations have been introduced to study hydrogen and its isotropic deuterium ( d ) and tritium ( t ) . self - diffusion coefficients in the pure h system and mutual diffusion for d - t mixtures were determined for temperatures @xmath51 to 10 ev and equivalent h mass densities 0.1 to 8.0 g / cm@xmath1 @xcite . qmd and ofmd simulations of self - diffusion , mutual diffusion , and viscosity have recently been performed on heavier elements ( fe , au , be ) @xcite and on mixtures of li and h @xcite . the present work selects h - he mixture as a representative system and examines some of the standard mixing rules with respect to the eos and transport properties ( viscosity , self and mutual diffusion coefficients ) in the warm dense regime that covers standard extreme condition as reached in the interiors of jupiter and saturn . the thermophysical properties of the full mixture and the individual species have been derived from qmd simulations , where the electrons are quantum mechanically treated through finite - temperature ( ft ) dft and ions move classically . in the next section , we present the formalism for qmd and for determining the static and transport properties . then , eos , viscosity , and diffusion coefficients for h - he mixtures are presented , and the qmd results are compared with the results from reduced models and qmd based linear mixing models . finally , concluding remarks are given . in this section , a brief description of the fundamental formalism employed to investigate h - he mixtures is introduced . the basic quantum mechanical density functional theory forms the basis of our simulations . the implementation of schemes in determining diffusion and viscosity is discussed . mixing rules that combine pure species quantities to form composite properties is also presented . qmd simulations have been performed for h - he mixtures by using vienna _ ab initio _ simulation package ( vasp ) @xcite . in these simulations , the electrons are treated fully quantum mechanically by employing a plane - wave ft - dft description , where the electronic states follow the fermi - dirac distribution . the ions move classically according to the forces from the electron density and the ion - ion repulsion . simulations have been performed in the nvt ( canonical ) ensemble where the number of particles @xmath7 and the volume are fixed . the system was assumed to be in local thermodynamic equilibrium with the electron and ion temperatures being equal ( @xmath8 ) . in these calculations , the electronic temperature has been kept constant according to fermi - dirac distribution , and ion temperature is controlled by noe thermostat @xcite .    at each step during md simulations , a set of electronic state functions [ @xmath9 for each * k*-point are determined within kohn - sham construction by @xmath10 with @xmath11 in which the four terms respectively represent the kinetic contribution , the electron - ion interaction , the hartree contribution and the exchange - correlation term . the electronic density is obtained by @xmath12 then by applying the velocity verlet algorithm , based on the force from interactions between ions and electrons , a new set of positions and velocities are obtained for ions . all simulations are performed with 256 atoms and 128 atoms for pure species of h and he , and as for the case of the h - he mixture , a total number of 245 atoms ( 234 h atoms and 11 he atoms ) for a mixing ratio @xmath13 ( corresponding to the h - he immiscibility region determined in the work of morales _ et al . _ @xcite ) has been adopted , where a cubic cell of length @xmath14 ( volume @xmath15 ) is periodically repeated . the simulated densities range from 1.0 to 4.0 g / cm@xmath1 for pure h system . as for pure he and h - he mixture , the size of the supercell is chosen to be the same as that for pure h to secure a constant electron number density ( in the range @xmath0/m@xmath1 ) . the temperature from 4000 k to 20000 k has been selected to highlight the conditions in the interiors of jupiter and saturn . the convergence of the thermodynamic quantities plays an important role in the accuracy of qmd simulations . in the present work , a plane - wave cutoff energy of 1200 ev is employed in all simulations so that the pressure is converged within 2% . we have also checked out the convergence with respect to a systematic enlargement of the * k*-point set in the representation of the brillouin zone . in the molecular dynamic simulations , only the @xmath16 point of brillouin zone is included . the dynamic simulation is lasted 20000 steps with time steps of 0.2 @xmath17 0.7 fs according to different densities and temperatures . for each pressure and temperature , the system is equilibrated within 0.5 @xmath17 1 ps . the eos data are obtained by averaging over the final 1 @xmath17 3 ps molecular dynamic simulations . the self - diffusion coefficient @xmath18 can either be calculated from the trajectory by the mean - square displacement @xmath19 or by the velocity autocorrelation function @xmath20 where @xmath21 is the position and @xmath22 is the velocity of the @xmath23th nucleus . only in the long - time limit , these two formulas of @xmath18 are formally equivalent . sufficient lengths of the trajectories have been generated to secure contributions from the velocity autocorrelation function to the integral is zero , and the mean mean - square displacement away from the origin consistently fits to a straight line . the diffusion coefficient obtained from these two approaches lie within 1 % accuracy of each other . here , we report the results from velocity autocorrelation function . we have also computed the mutual - diffusion coefficient @xmath24 from the autocorrelation function @xmath25 with @xmath26 where the concentration and particle number of species @xmath27 are denoted by @xmath28 and @xmath29 , respectively , and the total number of particles in the simulation box @xmath30 . the quantity @xmath31 is the thermodynamic factor related to the second derivation of the gibbs free energy with respect to concentrations @xcite . in the present simulations , @xmath31 value has been adopted equal to unity since studies with leonard - jones and other model potentials have shown that for dissimilar constituents the @xmath31-factor departs from unity by about 10% @xcite . the viscosity @xmath32 has been computed from the autocorrelation function of the off - diagonal component of the stress tensor @xcite @xmath33 the results are averaged from the five independent off - diagonal components of the stress tensor @xmath34 , @xmath35 , @xmath36 , @xmath37 , and @xmath38 . different from the self - diffusion coefficient , which involves single - particle correlations and attains significant statistical improvement from averaging over the particles , the viscosity depends on the entire system and thus needs very long trajectories so as to gain statistical accuracy . to shorten the length of the trajectory , we use empirical fits @xcite to the integrals of the autocorrelation functions . thus , extrapolation of the fits to @xmath39 can more effectively determine the basic dynamical properties . both of the @xmath18 and @xmath40 have been fit to the functional in the form of @xmath41 $ ] , where @xmath42 and @xmath43 are free parameters . reasonable approximation to the viscosity can be produced from the finite time fitting procedure , which also serves to damp the long - time fluctuations . the fractional statistical error in calculating a correlation function @xmath44 for molecular - dynamics trajectories @xcite can be given by @xmath45 where @xmath43 is the correlation time of the function , and @xmath46 is the length of the trajectory . in the present work , we generally fitted over a time interval of [ 0 , @xmath47 .      here , we examine two representative mixing rules . the first , termed density - matching rule ( mrd ) with the inspiration of a two - species ideal gas . the second , termed pressure - matching rule ( mrp ) , which follows from two interacting immiscible fluids .    in the mrd , the volume of the individual species is set equal to that of the mixture ( @xmath48 ) , and qmd simulations are performed for h at a density of @xmath49 and he at @xmath50 at a temperature @xmath51 . then , pressure predicted by mrd is determined by simply adding the individual pressures from the pure species h and he simulations . other transport coefficients , such as mutual diffusion and viscosity , follow the same prescription and are summarized as @xmath52 the superscript is used to denote values predicted from mrd . the derived pressure based on density mixing rule generally follows from the ideal noninteracting h and he gas in a volume @xmath53 . the mrp has a more complicated construction compared to mrd . mrp can be characterized as the following prescription : @xmath54 in this case , we have performed a series of qmd simulations on the individual species h and he , where the volumes change under a constraint @xmath55 until the individual pressures equal to each other @xmath56 . the total pressure becomes the predicted value . here , we use the excess or electronic pressure @xmath57 to evaluate this mrp mixing rule . composite properties such as mutual diffusion and viscosity are evaluated by combining the individual species results via volume fractions @xmath58 . finally , we also derive properties of the mixture from a slightly more complex mixing rule @xcite , as so - called binary ionic mixture ( bim ) : @xmath59 with @xmath60 the predicted mutual - diffusion coefficient or the viscosity . the subscript @xmath61 denotes the mixture and @xmath23 the pure species . in this section , the wealth of information derived from qmd calculations are mainly presented through figures , and the general trends of the eos as well as transport coefficients are concentrated in the text . it is , therefore , interesting to explore not only to get insight into the interior physical properties of giant gas planets but also to examine a series of mixing rules for hydrogen and helium . additionally , one can consider the influence of helium on the eos and transport coefficients of mixing . high precision eos data of hydrogen and helium are essential for understanding the evolution of jupiter and target implosion in icf . experimentally , the eos of hydrogen and helium in the fluid regime have been studied through gas gun @xcite , chemical explosive @xcite , magnetic driven plate flyer @xcite , and high power laser @xcite . since these experiments were limited by the conservation of mass , momentum , and energy , the explored density of warm dense matter were limited within @xmath62 times of the initial density . recently , a new technique combined diamond anvil cell ( dac ) and high intensity laser pulse has successfully been proved to provide visible ways to generate shock huguniot data of hydrogen over a significantly broader density - temperature regime than previous experiments @xcite . however , the density therein was still restricted within 1 g / cm@xmath1 .        in our simulations , wide range eos for h , he , and h - he mixtures have been determined according to qmd method . the eos can be divided into two parts that is , contributions from the noninteracting motion of ions ( @xmath63 ) and the electronic term ( @xmath64 ) , @xmath65 where @xmath64 is calculated directly through dft . in fig . [ fig_eos1 ] ( a ) , we have compared our results of @xmath64 with that of holst _ et al . _ @xcite , where the electronic pressure is expressed as a smooth function in terms of density and temperature , and the results agree with each other with a very slight difference ( accuracy within 5% ) . in the simulated density and temperature regime , we do not find any signs indicating a liquid - liquid phase transition ( @xmath66 ) or plasma phase transition ( @xmath67 ) , which are characterized by molecular dissociation and ionization of electrons , respectively . with considering the mixing of he into h , the electronic pressure is effectively reduced , as has been shown in fig . [ fig_eos1 ] ( b ) . mrd accounts for contribution from noninteracting h and he subsystems in the volume of the mixtures . in mrd , the pressure contributed from noninteracting ions is the same as that of the mixture , but the electronic pressure is much lower due to the low electronic density in the pure species simulations . it is indicated that the electronic pressure is underestimated by mrd model at about 8% @xmath17 9% ( see fig . [ fig_eos2 ] ) . for mrp model , we have firstly performed a series of pure species simulations at a wider density ( temperature ) regime compared to h - he mixtures . then , the simulated eos data are fitted into smooth functions in terms of density and temperature . under the constraint of @xmath68 , we have predicted the electronic pressures @xmath69 according to mrp model by solving pure species eos function at certain densities and temperatures , as shown in fig . [ fig_eos2 ] . it is indicated that the mrp model agrees better with direct qmd simulations ( accuracy within 3% ) , the difference mainly come from the ionic interactions between h and he species after mixing .       and 8 g / cm@xmath1 , respectively . ( a ) self - diffusion coefficient ; ( b ) viscosity . the direct qmd simulated results are presented by black crosses , while the fitted results are denoted by blue dashed lines . ] qmd simulations have been performed within the framework of ft - dft to benchmark the dynamic properties of h , he , and h - he mixture in the wdm regime . illustrations for the self - diffusion coefficients and viscosity ( for h and he at densities of 2.0 g / cm@xmath1 and 8.0 g / cm@xmath1 , respectively ) at a temperature 12000 k , as well as their fits are shown in fig . [ fig_dandeta ] . the trajectory of the present simulations lasts 4.0 @xmath17 14.0 ps , and correlation times between 1.0 and 15.0 fs . as a consequence , the computational error for the viscosity lies within 10% . after accounting for the fitting error and extrapolation to infinite time , a total uncertainty of @xmath17 20% can be estimated . the uncertainty in the self - diffusion coefficients is smaller than 1% , due to the additional @xmath70 advantage given by particle average .    . for helium , the density is 8 g / cm@xmath1 . the electron number density for pure species ( h and he ) is @xmath71/m@xmath1 . ]    dynamic properties of wdm are generally governed by two dimensionless quantities , namely , ionic coupling ( @xmath16 ) and electronic degenerate ( @xmath72 ) parameter . the former one is defined by the ratio of the potential to kinetic energy @xmath73 , with @xmath74 the ionic charge , and @xmath75 the ion - sphere radius ( @xmath76 is the number density ) . the latter one @xmath77 , where @xmath78 is fermi temperature . it has been reported that dynamic properties such as diffusion coefficients and viscosity can be represented purely in terms of ionic coupling parameter @xmath16 according to molecular dynamics or monte carlo simulations based on one component plasma ( ocp ) model @xcite , where ions move classically in a neutralizing background of electrons . for instance , hansen _ @xcite introduced a memory function to analyze the velocity autocorrelation function , and obtain the diffusion coefficient in terms of @xmath79 with the plasma frequency @xmath80 . based on classical molecular dynamic simulations , bastea @xcite has fitted the viscosity into the following form @xmath81 with @xmath82 , @xmath83 , @xmath84 , and @xmath85 . since ocp model is restricted to a fully ionized plasma , we use @xmath74=1.0 ( or 2.0 ) for hydrogen ( or helium ) to compute the self - diffusion coefficient and viscosity . in fig . [ fig_qmd_ocp ] ( a ) , we show comparison between qmd and ocp model @xcite for hydrogen and helium at densities of 2 g / cm@xmath1 and 8 g / cm@xmath1 . the general tendency for the self - diffusion coefficient with respect to temperature is similar for qmd and ocp model , however , the difference up to @xmath1760 % is observed between the two results . for the viscosity [ fig . [ fig_qmd_ocp ] ( b ) ] , ocp @xcite predicts smaller ( larger ) values for hydrogen ( helium ) compared to qmd simulations . the viscosity is governed by interactions between particles and ionic motions , contribution from the former one decrease with the increase of temperature , while , it increases for the latter one . as a consequence , the viscosity may have local minimum along temperature . for hydrogen , the local minimum locates around 10000 k and 14000 k indicated by qmd and ocp model , respectively . while in the case of helium , we do not observe any signs for the local minimum in the simulated regime . /m@xmath1 and temperature from 4000 @xmath17 20000 k. ]    in fig . [ fig_mix ] , we have shown the mutual diffusion coefficient @xmath86 and viscosity @xmath87 for h - he mixture with an electron number density of @xmath88/m@xmath1 , and results from mixing models are also provided . the transport coefficients predicted by mrd can be directly evaluated through eq . ( [ eq_mrd ] ) . for mrp model , we have firstly fitted the self - diffusion coefficients and viscosity in terms of density and temperature , after determining the volume for each species under the constraint of @xmath89 , the transport coefficients are then obtained . in bim model , we have used @xmath90 and @xmath91 , then , the transport coefficients are determined by eq . ( [ eq_bim ] ) . here we would like to stress that in some mixture studies based on average atom models , the properties of pure species are derived from perturbed - atom models , where boundary conditions are introduced from the surrounding medium by treating a single atom within a cell . in the present work , the dynamic properties of different mixing rules originate from qmd calculations of the individual species . despite divorced of the h - he interactions , the pure species calculations still contain complex intra - atomic interactions based on large samples of atoms . the mutual diffusion coefficient of h - he mixture shows a linear increase with respect of temperature , as indicated in fig . [ fig_mix](a ) . the data from mrp and bim models have a better agreement with qmd simulations compared with that of the mrd model , where ion densities are reduced and results in a larger diffusion coefficient . the viscosity of h - he mixture has a more complex behavior than pure species under extreme condition . as shown in fig . [ fig_mix ] ( b ) , mrd is valid at low temperature , while mrp works at higher temperature . bim rule moves the results into better agreement with the h - he mixture , leaving within 30% or better for the simulated conditions . in summary , we have performed systematic qmd simulations of h , he , and h - he mixture in the warm dense regime for electron number density ranging from @xmath0/m@xmath1 and for temperatures from 4000 to 20000 k. the present study concentrated on thermophysical properties such as the eos , diffusion coefficient , and viscosity , which are of crucial interest in astrophysics and icf . various mixing rules have been introduced to predict dynamical properties from qmd simulations of the pure species and compare with direct calculations on the fully interacting mixture . we have shown that mrd and mrp rules produce pressures within about 10 % of the h - he mixture , however , the mutual diffusion coefficients are as different as 75 % and it is 50 % for the viscosity . bim rule generally gives better agreement with the mixture results . we have also compared our qmd results with ocp model for the pure species . this work was supported by nsfc under grants no . 11275032 , no . 11005012 and no . 51071032 , by the national basic security research program of china , and by the national high - tech icf committee of china . 99 w. b. hubbard , science * 214 * , 145 ( 1981 ) . h. f. wilson and b. militzer , astrophys . j. * 745 * , 54 ( 2012 ) . h. f. wilson and b. militzer , phys . lett . * 108 * , 111101 ( 2012 ) . d. bruno , c. catalfamo , m. capitelli , g. colonna , o. de pascale , p. diomede , c. gorse , a. laricchiuta , s. longo , d. giordano , and f. pirani , phys . plasmas * 17 * , 112315 ( 2010 ) . s. atzeni and j. meyer - ter - vehn , _ the physics of inertial fusion : beam plasma interaction , hydrodynamics , hot dense matter _ , international series of monographs on physics ( clarendon press , oxford , 2004 ) .    j. d. lindl , _ inertial confinement fusion : the quest for ignition and energy gain using indirect drive _ , ( springer - verlag , new york , 1998 ) . m. d. knudson and m. p. desjarlais , phys . lett . * 103 * , 225501 ( 2009 ) . w. lorenzen , b. holst , and r. redmer , phys . lett . * 102 * , 115701 ( 2009 ) . d. j. stevenson , phys . b * 12 * , 3999 ( 1975 ) . w. b. hubbard and h. e. dewitt , astrophys . j. * 290 * , 388 ( 1985 ) . j. e. klepeis , k. j. schafer , t. w. barbee iii , and m. ross , science * 254 * , 986 ( 1991 ) . j. vorberger , i. tamblyn , b. militzer , and s. a. bonev , phys . b * 75 * , 024206 ( 2007 ) . b. militzer , phys . b * 87 * , 014202 ( 2013 ) . o. pfaffenzeller , d. hohl , and p. ballone , phys . lett . * 74 * , 2599 ( 1995 ) . l. a. collins , j. d. kress , d. l. lynch , and n. troullier , j. quant . transf . * 51 * , 65 ( 1994 ) . i. kwon , j. d. kress , and l. a. collins , phys . b * 50 * , 9118 ( 1994 ) . l. collins , i. kwon , j. kress , n. troullier , and d. lynch , phys . e * 52 * , 6202 ( 1995 ) . j. clrouin and s. bernard , phys . e * 56 * , 3534 ( 1997 ) . i. kwon , l. collins , j. kress , and n. troullier , epl * 29 * , 537 ( 1995 ) . j. clrouin and j .- f . dufrche , phys . e * 64 * , 066406 ( 2001 ) . j. d. kress , j. s. cohen , d. a. horner , f. lambert , and l. a. collins , phys . e * 82 * , 036404 ( 2010 ) . f. lambert , j. clouin , and s. mazevet , epl * 75 * , 681 ( 2006 ) . f. lambert , j. clouin , and g. zrah , phys . e * 73 * , 016403 ( 2006 ) . c. wang , y. long , x .- he , and p. zhang , phys . e * 87 * , 043105 ( 2013 ) . d. a. horner , f. lambert , j. d. kress , and l. a. collins , phys . b * 80 * , 024305 ( 2009 ) . m. a. morales , c. pierleoni , e. schwegler , and d. m. ceperley , proc . * 106 * , 1324 ( 2009 ) . j. p. hansen and i. r. mcdonald , _ the theory of simple liquids _ , 2nd ed . ( academic press , new york , 1986 ) , chap . 8 . m. schoen and c. hoheisel , mol.phys . * 52 * , 33 ( 1984 ) ; * 52 * , 1029 ( 1984 ) . m. p. allen and d. j. tildesley , _ computer simulation of liquids _ ( oxford university press , new york , 1987 ) . j. d. kress , james s. cohen , d. p. kilcrease , d. a. horner , and l. a. collins , phys . e * 83 * , 026404 ( 2011 ) . r. zwanzig and n. k. ailawadi , phys . rev . * 182 * , 280 ( 1969 ) . s. bastea , phys . e * 71 * , 056405 ( 2005 ) . w. j. nellis , rep . * 69 * , 1479 ( 2006 ) . g. v. boriskov , a. i. bykov , r. ilkaev , v. d. selemir , g. v. simakov , r. f. trunin , v. d. urlin , a. n. shuikin , and w. j. nellis , phys . b * 71 * , 092104 ( 2005 ) . m. d. knudson , d. l. hanson , j. e. bailey , c. a. hall , j. r. asay , and c. deeney , phys . b * 69 * , 144209 ( 2004 ) . g. w. collins , l. b. da silva , p. celliers , d. m. gold , m. e. foord , r. j. wallace , a. ng , s. v. weber , k. s. budil , and r. cauble , science * 281 * , 1178 ( 1998 ) . t. r. boehly , d. g. hicks , p. m. celliers , t. j. b. collins , r. earley , j. h. eggert , d. jacobs - perkins , s. j. moon , e. vianello , d. d. meyerhofer , and g. w. collins , phys . plasmas * 11 * , l49 ( 2004 ) . d. g. hicks , t. r. boehly , p. m. celliers , j. h. eggert , s. j. moon , d. d. meyerhofer , and g. w. collins , phys . b * 79 * , 014112 ( 2009 ) . p. loubeyre , s. btygoo , j. eggert , p. m. celliers , d. k. spaulding , j. r. rygg , t. r. boehly , g. w. collins , and r. jeanloz , phys . b * 86 * , 144115 ( 2012 ) . b. holst , r. redmer , and m. p. desjarlais , phys . b * 77 * , 184201 ( 2008 ) . f. lambert , ph.d . thesis , universit paris xi . commissariat l nergie atomique , 2007 . b. bernu and p. vieillefosse , phys . a * 18 * , 2345 ( 1978 ) . z. donk , b. nyri , l. szalai , and s. holl , phys . lett . * 81 * , 1622 ( 1998 ) ; z. donk and b. nyri , phys . plasmas * 7 * , 45 ( 2000 ) . j. daligault , phys . lett . * 96 * , 065003 ( 2006 ) ; * 103 * , 029901(e ) ( 2009 ) . j. p. hansen , i. r. mcdonald , and e. l. pollock , phys a * 11 * , 1025 ( 1975 ) . r. m. more , adv . at . * 21 * , 305 ( 1985 ) .
Abstract:

Answer: thermophysical properties of hydrogen , helium , and hydrogen - helium mixtures have been investigated in the warm dense matter regime at electron number densities ranging from @xmath0/m@xmath1 and temperatures from 4000 to 20000 k via quantum molecular dynamics simulations . we focus on the dynamical properties such as the equation of states , diffusion coefficients , and viscosity . mixing rules ( density matching , pressure matching , and binary ionic mixing rules ) have been validated by checking composite properties of pure species against that of the fully interacting mixture derived from qmd simulations . these mixing rules reproduce pressures within 10% accuracy , while it is 75 % and 50 % for the diffusion and viscosity , respectively . binary ionic mixing rule moves the results into better agreement . predictions from one component plasma model are also provided and discussed .