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: shape - memory alloys ( smas ) are active crystalline materials with desirable mechanical properties , exploited in a number of applications from engineering to medicine , as actuators , sensors and dampeners @xcite . in these substances , effects such as super- ( or pseudo-)elasticity and shape memory originate from a reversible martensitic transformation triggered by either stress or temperature . these phenomena are strongly affected by the formation of austenite / martensite microstructures , typically including twins and habit planes in various arrangements which may act as proper avenues for the reversible phase change  @xcite . a number of theoretical approaches have been developed for the analysis of such morphologies , going back to the classical ` crystallographic theory of martensite ' @xcite , based on considerations of kinematic compatibility between phases and variants . this viewpoint has been much developed in recent decades  @xcite .    in parallel , the experimental investigation of sma microstructure has adopted a variety of means , primarily optical and electron microscopy , originating an extensive literature , see for instance the references in @xcite . developments in the techniques of full - field measurements have recently also enabled the experimental analysis of the microstructural strain field over extended regions of a sample s surface . moir interferometry and digital image correlation are typical examples of such techniques , which collect information on the response of the tested sma samples by distinguishing the different phases through their strain amplitude  @xcite . infrared thermography has also been employed to measure the inhomogeneous temperature distributions associated with the phase change in smas , giving relevant information on the features of the martensitic transformation  @xcite . the spatial distribution of the phases and variants over a sample can also be deduced through the grid method @xcite , wherein this information is derived by analyzing the deformation of a periodic grid deposited on the sample , obtaining the strain levels which are generally different among the phases and their variants . a main advantage of the grid method is that it offers a good compromise between strain resolution and spatial resolution @xcite . however , this full - field measurement technique have not yet been employed to investigate phase - change intermittency in memory materials . this inhomogeneity is a fundamental feature of the martensitic transformation in smas , wherein under a smooth thermal or mechanical driving , the strain progresses through a set of discrete avalanche - like events , corresponding to transitions between neighboring metastable states . such effects have been previously studied through the measurement of calorimetric and/or acoustic signals  @xcite . these techniques have revealed many aspects of the bursty character of the martensitic transformation , with events exhibiting , in many cases , power - law size distributions , possibly after training  @xcite . this links reversible sma martensites to other material systems exhibiting intermittent dynamics and scale - free behavior , in magnetism , superconductivity , brittle fracture , crystal plasticity @xcite , in turn framing them within the background of a wide variety of complex systems in nature exhibiting avalanche - mediated behavior , as in turbulence , earthquakes , computer or social networks , and financial markets @xcite . while recent investigations have examined quantitatively some spatial aspects in the evolution of martensitic transition phenomena @xcite , there is a lack of systematic quantitative information on the strain events derived from the analysis of evolving strain maps during the phase change in smas . in the present study we used the grid method to perform the first full - field investigation of the intermittent progress of the microstructural strain field in a cu - al - be single crystal across the mechanically - induced martensitic transformation . this was done through uniaxial force - controlled loading tests , during which we have recorded and analyzed the super - elastic and hysteretic stress - strain behavior of the sample , as well as surface strain data in the forward and reverse transition . inspired by earlier work on the stress - controlled loading of smas in  @xcite ( see also @xcite ) , an experimental apparatus based on gravity was used , which applied to the sample a low - rate perfectly monotonic force , controlled by very slow flow water pumps , at almost isothermal conditions . the loading device was designed to be capable of very small force increments , enabling us to highlight the intermittency in the alloy , through the quantitative analysis of the strain field features measured over the sample . in contrast to conventional testing machines , as used for instance in  @xcite , in the present apparatus the crystal may freely adjust its orientation in relation to the vertical loading direction . a main reason for this was obtaining relatively less complex microstructures , developed in the presence of minimal extraneous effects such as external friction , lattice plastification , or thermal inhomogeneity , to enable us to investigate the strain transformation intermittency occurring in the crystal in its most elementary and basic form . this paper is organized as follows . [ section_2 ] presents the sample under study , the experimental setup , the loading conditions annd the procedure employed for image processing to retrieve strain maps . [ section_3 ] presents an analysis of the hysteresis and microstructures deduced from the strain evolution recorded during a full loading test . [ section_4 ] discusses the observed intermittent character of the phase transformation . [ section_5 ] gives some conclusive remarks . further details on the experimental apparatus and on data acquisition and treatment are given in the appendix . the supplemental material @xcite comprises four video files showing different aspects of the experimental results . we tested a cual@xmath0be@xmath1 ( wt . % ) single crystal sample , the same used in the loading experiments described in  @xcite , to which we refer for additional information on the crystallographic properties and orientation of the lattice in the sample ( see the appendix ) . this sma , under uniaxial tensile loading , undergoes a reversible martensitic transformation , from a cubic austenite to a long - period stacking monoclinic martensite with twelve variants , which are all compatible with the austenitic phase @xcite . the two - fold axes in the monoclinic variants originate from the cubic axes in the austenite @xcite . [ figure_1](a ) shows the sample geometry . aluminum tabs were clamped in two pairs of jaws by means of screws onto the ends of the sample , to avoid damage or sliding once positioned in the testing machine . a bidirectional grid ( pitch @xmath2 mm ) was transferred using the technique described in  @xcite on the monitored gauge region ( dimensions @xmath3 mm @xmath4 38 mm respectively along the @xmath5- and @xmath6-directions ) covering nearly all the exposed surface of the sma sample , with thickness @xmath7 mm . see the appendix for more details .     at various stages of the loading - unloading cycle ( austenite is in blue , martensite in red ) . the stress - strain hysteresis regards the average strain @xmath8 vs.  the average stress @xmath9 . the green dots mark points on the transformation plateaus for reference in the figures below . also shown are maps of all the in - plane strain components and of the angle @xmath10 , corresponding to the average strain value @xmath11 on loading ( top ) , and to @xmath12 on unloading ( bottom ) . notice that in the @xmath13 map the austenite is in red and the martensite in blue . see the supplemental video  # 1  @xcite . ] the loading apparatus is schematically shown in fig . [ figure_1](b ) .  the sample , which before testing was austenitic with some possible residual martensite , was suspended at its top end through a ball joint which allows for any rotations . a tank was hung through a cable to the jaws at the bottom end of the sample . the tank was filled with water during the experiments , enabling us to apply a gravity - load in the vertical direction . loading and unloading of the tank were controlled by electronic pumps producing a constant and slow flux of water . this allowed for a careful stress control of both the forward and reverse austenite - to - martensite phase change upon loading and unloading , respectively . as transformation stresses depend on the temperature , we ensured the stability of the ambient temperature @xmath14c during the entire duration of the test ( more than 45h ) . due to the very slow driving , the testing conditions were almost isothermal . before using the pumps , a pre - load was applied to the sample to reduce the test duration shortening the first part of the elastic climb where martensitic transformation events are negligible . this gave the reference configuration for the ensuing strain maps . the forward and reverse transformation plateaus were spanned through a slowly - varying , monotonic loading controlled by adjusting the water flow through the pumps .    in what follows we specifically report on a test in which we have observed the whole forward and reverse transformations under these conditions : ( a ) pre - load of 34.37mpa ; ( b ) constant loading rate 1.055mpa / h up to 57.29mpa ( duration : about 22h ) ; ( c ) constant unloading rate @xmath15mpa / h down to 35.95mpa ( duration : about 23h 15min ) . due to the inherent feedback loop in the control , very slow monotonic rates such as the ones above , obtained with the present device , are not possible with conventional testing machines . as a comparison , in earlier stress - controlled tests , such as in  @xcite , the rate was about 37mpa / h in the gauge region of the sample . we have gathered about 20,000 images of the sample over the 45-hour duration of the test , from which we have determined the spatial distribution of the strain levels on the sample surface by applying suitable grid - image processing . this enabled us to obtain four maps from each image during the loading - unloading cycle , giving the in - plane linear strain fields @xmath13 , @xmath16 , @xmath17 , and the rotation - angle field @xmath10 about the @xmath18-direction ( see fig . [ figure_1 ] for the axes orientation ) . more details on the data acquisition and treatment are given in the appendix . fig .  [ hysteresis ] presents the hysteretic stress - strain relation recorded during the loading - unloading cycle . further information is given in the supplemental video  # 1  @xcite . [ hysteresis ] also shows the evolving strain maps through the forward and reverse transformations , where we see that , expectedly , @xmath19 and @xmath20 . [ evolution_nu ] shows the associated hysteresis in the evolution of the martensitic volume fraction @xmath21 , defined as the percentage of the sample surface where  @xmath16 exceeds the threshold  @xmath22 , which is about half of the maximum local strain value recorded during the test .     vs.  the average stress @xmath9 during the stress - driven forward and reverse transformation in fig . [ hysteresis ] . the fact that @xmath21 does not reach 1 originates from the possible presence of both residual martensite at the lowest loading as well as some residual austenite at the highest loading . ] we observe in fig . [ hysteresis ] the development of the strain microstructure as the phase transformation progresses from austenite ( blue region in the @xmath16 map ) towards a unique variant of martensite ( red region in the @xmath16 map ) , which is selected among the twelve possible variants due to the suitable alignment of its lattice with the imposed load . the presence of a single variant in this test is indicated by the very homogeneous values for the field @xmath16 in the red martensitic domain of the sample ( with average approaching @xmath23 near the end of the transformation ) , and by the fact that the shear component @xmath17 has constant sign . the occurrence of a single variant is also suggested by fig . [ strain_clouds ] , showing how the material evolves between strain values mostly clustering at two levels , corresponding to ( strained ) austenite and to a single ( strained ) variant of martensite . at the intermediate loading configurations in fig . [ strain_clouds ] a trail is noticed between the pure - phases strain values . this is due to the elastic deformation of the phases lattices under while the load is present , as well as to the fact that any small scale phase domains in the sample are below the spatial resolution of the apparatus ( see the appendix ) , and thus give rise to averaged strain values in the corresponding pixels . the evolution of the material from the austenitic energy well to the martensitic one along the upper transformation plateau can also be viewed in the supplemental video  # 2  @xcite . are shown . the color bar indicates the fraction of pixels with the given strain values , with grid pitch on the @xmath13 and @xmath16 axes equal respectively to 0.005 and 0.007 . see also the supplemental video  # 2  @xcite.,width=283,height=264 ]    the @xmath16 strain maps in fig . [ hysteresis ] show that the forward transformation develops through needle - like bands which progressively enlarge their width and whose length eventually crosses the entire sample . their inclination agrees with the theoretical estimates @xcite based on the compatibility equations giving the habit - planes for this transformation @xcite . the snapshots in fig . [ hysteresis ] also confirm the existence of relative lattice rotations between regions occupied by different phases , with jumps of the angle @xmath10 from + 0.03 to @xmath240.03 rad across the austenite - martensite interfaces , again in good agreement with the theoretically computed values @xcite . the @xmath16-field in fig . [ hysteresis ] appears more complex along the reverse transformation , possibly involving also a second habit plane @xcite ( see more details below , and the supplemental videos  @xcite ) . the evolution on a loading cycle is thus hysteretic not only in terms of the stress - strain relation , but also in the distinct microstructural morphologies developed in loading vs.  unloading . this might be a factor contributing to a slight mechanical irreversibility in the phase change .     along the mid - section  ab of the sample . the detail in the upper left corner of the figure shows the evolution of @xmath16 in the interval [ 0h , 14h ] at a magnified strain scale.,height=340 ]    fig . [ time_evolution_e ] shows , as an example , the time evolution of the strain component  @xmath16 along the central section ab of the gauge region . most of the phase transformation occurs within the time window [ 14h , 19h ] on loading ( the upper plateau in the stress - strain curve of fig . [ hysteresis ] ) , and [ 38h , 43h ] on unloading , with a total stress increase / decrease of about 5mpa on each plateau . this gives the overall slant of the hysteresis loop in the stress - strain plane in the present test . the upper left corner of fig . [ time_evolution_e ] shows the evolution of  @xmath16 in the interval [ 0 h , 14 h ] at a magnified strain scale . the strain profiles in fig . [ time_evolution_e ] show the evolution of the phase microstructures in fig . [ hysteresis ] , where we also notice that the end pixels a and b and their adjacent regions do not undergo significant strain . it is interesting to contrast these results with those recorded in earlier tests performed on the same cu - al - be sample by means of a conventional uniaxial loading machine , which also exploited the grid method @xcite . [ comparison ] reports the @xmath25-strain maps in the two cases , near the end of the transformation when comparable small fractions of austenite remain in the sample ( blue regions ) . the martensitic ( red ) regions show the different microstructures developed in the two tests . as noted above , the present loading device produced a homogeneous @xmath16-field ( fig . [ comparison ] , left ) , while the conventional machine produced a martensitic zone with two alternating strain levels ( fig . [ comparison ] , right ) , which we interpret as martensite twinning . the development of twin bands is indeed expected in the conventional machine , as the latter does not allow for rotations of the sample , forcing the material to activate more than one martensite variant to accommodate the imposed load . conversely , a single variant can accommodate the loads in the present testing device , as the latter allows for an overall rotation to accompany the transformation strain in the sample . this helped in producing less complex microstructures , for this first full - field study of strain intermittency in smas .     and @xmath16 denote respectively the hencky and linear strain components . the faint horizontal lines in the left snapshot are artifacts due to slight pitch variations in the grid , not completely eliminated by the data treatment . the small dots originate from local grid defects . , height=313 ] recorded at each consecutive image for the average strain @xmath8 during loading ( a ) , and unloading . the reference green dots mark the transformation plateaus . reported in red is the corresponding time evolution of @xmath8 , see also fig . [ hysteresis ] . the corresponding heavy - tailed probability distributions for @xmath26 are also shown ( logarithmic binning).,scaledwidth=47.0% ]    we study here various aspects of the non - continuous progress of the observed forward and reverse martensitic transformation . we see in fig . [ hysteresis ] that the global variable @xmath8 displays a rather continuous evolution at the time scale of the entire test ( the same holds for @xmath21 in fig . [ evolution_nu ] ) . however , fig . [ spikes_delta_epsilon_mean ] shows that the evolution in time of the increments @xmath27 obtained by comparing each subsequent image taken during the test , produces a very spiky signal , evidencing a pronounced strain intermittency in the phase - transforming sample under the slowly - varying monotonic load . [ spikes_delta_epsilon_mean ] shows the corresponding heavy - tailed probability distributions @xmath26 for the @xmath27 amplitudes over almost two decades ( we threshold the values of @xmath8 and @xmath27 at @xmath28 , see the appendix ) . yet another form of hysteresis in the forward vs.  the reverse transformation is highlighted in this way , as we see from fig . [ spikes_delta_epsilon_mean ] that most of the jumps @xmath27 are systematically smaller on unloading than on loading . however , the largest events are more numerous in the reverse than in the forward transformation . analogous behavior is observed for the spiky signal given by the phase - fraction jumps @xmath29 ( not shown ) . [ spikes_delta_epsilon_mean ] also evidences that the recorded strain intermittency is non - stationary along the transformation cycle , whereby there is more activity in the evolution of @xmath27 near the center of both the plateaus of the hysteresis as opposed to their ends ( more details on this below ) .    , evaluated at each consecutive image , along the mid - section of the sample . see also fig . [ time_evolution_e ] . ( b ) time evolution of @xmath30 at two pixels p1 and p2 . the time interval is [ 13h ; 18h ] , where most of the phase transformation occurred on loading . the reference green dots delimit the transformation strain plateaus . ]    the intermittent evolution of the global variables such as ( @xmath21 and ) @xmath8 originates from the intermittency in the underlying strain field in the sample . the latter is highlighted in fig . [ local_bursts ] , which reports the bursty behavior of the strain increments @xmath30 at some typical locations on the surface ( the increments again refer to each consecutive image ) . the main panel shows the space - time evolution of @xmath30 along the sample mid - section on loading . the lower panels plot the temporal evolution of @xmath30 at two pixels p1 and p2 of the same section . strain events occur throughout the forward and reverse phase change in the crystal , with the larger bursts proportionally likely to include phase transformation in the lattice as opposed to purely elastic local deformation . we notice high - intensity ( yellow - to - red ) localized strain bursts occur throughout the sample , indicating the local evolution of the phase microstructures in the crystal . the largest events in fig .  [ local_bursts ] are of the order of @xmath31 , but the scale is limited to less than @xmath32 to better show the evolution of @xmath30 in the intensity interval where most of the strain bursts occur for the two considered pixels p1 and p2 .    to quantify statistically the degree of pixel - level burstiness in the strain - field evolution across the sample , we consider the amplitudes of @xmath30 recorded at all pixels along the entire transformation cycle . we threshold the values of @xmath30 at @xmath33 , to stay above the noise level of the local strain measurements ( see the appendix ) . it results that the pixel values of @xmath30 throughout the sample span about an order of magnitude , with heavy - tailed size distributions which are rather abruptly truncated , as shown in fig . [ single_pixel_statistics ] separately for the forward and the reverse phase change . each single pixel thus gives to the overall strain change a contribution that is bounded by the transformation strain between the austenitic and martensitic configurations . the large transformation events must thus necessarily exhibit some spatial structure , obtained as the sum of many smaller local events .     of the strain jumps @xmath30 ( logarithmic binning ) computed for all pixels in the sample during the forward ( squares ) and reverse ( crosses ) transformation.,width=302 ] we now characterize the long - range aspects of the avalanche - mediated phase transformation . to take into account all the components of the evolving strain field , we consider the norm @xmath34 of the strain jumps rather than the sole @xmath30 component . we then check for any pixel p on which @xmath35 is , at any time , higher than @xmath33 . given any such p , we track a strain avalanche as the collection of pixels neighboring p whereon @xmath35 is also above threshold . in this way , about 14,000 events are obtained from the strain evolution data collected along the hysteresis cycle . we recall that the time interval between images in the present tests is almost ten seconds ( see the appendix ) , which is many orders of magnitude larger than the temporal scales of avalanche durations as recorded for instance through acoustic emission @xcite . thus the avalanches evidenced here largely pertain only to a single image , and can give no indication on the actual time dynamics of the transformation , as each avalanche merges a large and variable number of bursty microscale events . we first consider the number @xmath36 of avalanches detected per image in the sample . on loading , we record a total of almost 8,000 events , and about 6,000 events in the reverse transformation . in both case avalanching produces a non - stationary spiky signal , as shown for instance in fig . [ avalanche_statistics](a ) for the direct transformation . the corresponding heavy - tailed statistics in fig . [ avalanche_statistics](b ) evidence the very fluctuating character of these co - operative strain events , with @xmath36 spanning two orders of magnitude . the supplemental video  # 3  @xcite highlights the strain avalanching in the sample on loading , which helps visualize the bursty development of the martensite bands mentioned earlier . the supplemental video  # 4  @xcite shows the corresponding evolution of the strain avalanches , with their more complex microstructures , along the reverse transformation .     of avalanches detected per consecutive image ( black ) and its running sum ( red ) . the green dots mark the reference points on the transformation strain plateaus . ( b ) log - log plot of the corresponding probability distribution @xmath37 ( logarithmic binning).,width=302 ]     and @xmath38 of the strain avalanche magnitude @xmath39 and the avalanche size @xmath40 during the transformation cycle ( logarithmic binning ) . ( c ) log - linear plot of the correlation between @xmath39 and the ratio @xmath41 for the recorded strain avalanches . ( d)-(e ) values of the exponents @xmath42 respectively for @xmath43 and @xmath38 , determined by the maximum likelihood method as a function of the lower cutoff imposed to the data.,width=340,height=510 ]    for each avalanche we then consider its _ size _ @xmath40 given by the total number of involved pixels , and its _ magnitude _ @xmath39 given by the integral of @xmath35 over the avalanche . the quantity @xmath39 measured for each strain avalanche gives the contribution of such event to the intermittent deformation detected in the sample at each consecutive image as the loading was applied , and , in particular , to the sample elongation detailed in fig . [ spikes_delta_epsilon_mean ] . the probability distributions for both @xmath40 and @xmath39 , computed over the full cycle are shown in fig . [ pow_law ] . while none of the heavy - tailed statistics for strain intermittency considered so far ( figs .  [ spikes_delta_epsilon_mean]-[avalanche_statistics ] ) seemingly exhibit power - law character , the log - log plots for @xmath38 and @xmath43 in fig . [ pow_law](a)-(b ) appear in remarkable agreement with the emergence of scale - free behavior during the transition process , which is well documented by the power - law size statistics for acoustic - emission events reported in many cases for the martensitic transformation of smas  @xcite . we notice that the power - law character of the distributions of figs . [ pow_law](b)-(c ) is affected by a slightly larger frequency of events with sizes of the order of @xmath44 , which are mainly due to spurious avalanches due to grid defects . also , fig . [ pow_law](c ) shows there is a positive correlation between @xmath40 and @xmath39 , with a super - linear tendency for the avalanches with large @xmath40 to have larger @xmath39 . this leads to an increase ( to about six ) of the number of decades recorded for the avalanche magnitudes @xmath39 in fig . [ pow_law](a ) , as compared to those exhibited by the avalanche size @xmath40 in fig . [ pow_law](b ) . furthermore , the defect - related events in the small plume of fig .  [ pow_law](c ) for @xmath44 , widen the corresponding range of @xmath39 , and , at the corresponding scales in figs .  [ pow_law](a)-(b ) , locally slightly perturb both the distributions @xmath38 and @xmath43 . this affects the stability of these distributions exponents , which , as we see in figs . [ pow_law](d)-(e ) , can both evaluated only as being in the range 1.5 - 1.8 . this is compatible with the exponent 1.6 recently proposed for the size statistics of the optically - observed microstructural transformation events recorded in  @xcite at the surface of a ni - mn - ga polycrystal undergoing a temperature - driven martensitic transformation .     of the avalanches magnitude distribution @xmath43 respectively for the forward ( a ) , and for the reverse ( b ) transformation . the exponent values are determined by the maximum likelihood method as a function of the lower cutoff imposed to the data.,height=105 ]    , @xmath45 ( logarithmic binning ) , of the avalanche magnitudes @xmath39 during three consecutive equal - duration time intervals on the _ forward_-transformation plateau . the inset indicates each interval , with , in parenthesis , the corresponding number of events :  1265 , 6024 , 374 . ( b ) values of the exponent @xmath42 for each distribution @xmath46 , determined by the maximum likelihood method as a function of the lower cutoff imposed to the data . ( c - d ) data analogs of ( a - b ) for three consecutive equal - duration time intervals on the plateau for the _ reverse _ transformation , with corresponding number of events :  213 , 4798 , 614 . ]    a further hysteretic aspect of the strain evolution in our material is manifested by the different exponents of the distributions @xmath43 recorded for the strain avalanches on loading vs.  unloading , see fig . [ separate_exponents ] . we notice that the exponent in the reverse transformation , near 2 , is higher than the one on the forward transformation , near 1.5 . an exponent asymmetry of this type was also recorded in @xcite for the temperature - driven phase change in a nimnga sma . also acoustic emission data indicate a forward - reverse asymmetry in some sma polycristals @xcite . we make a final brief analysis of the non - stationarity of the transformation process , already remarked from figs . [ spikes_delta_epsilon_mean ] and [ avalanche_statistics ] . we consider six consecutive time sub - intervals of equal duration over the transformation cycle ( three for the forward and three for the reverse transformation ) . the number of events in each interval are indicated in the insets of figs . [ non - stationarity](a)-(c ) . we see that the non - stationarity of the process is marked by a much higher strain activity in the central portions of both plateaus , as opposed to their end intervals . this is seemingly analogous to the non - stationary signal recorded in other studies of sma phase change @xcite , or in the behavior of some hysteretic magnetic systems  @xcite . the present data also show there is a higher avalanching activity when the alloy is mostly austenitic than when martensite is predominant . the non - stationary character of the strain intermittency across the transformation is also evidenced through the probability distributions @xmath46 of the avalanche magnitudes @xmath39 pertaining to each one of the six time intervals , see figs . [ non - stationarity](a ) and ( c ) . [ non - stationarity](b ) and ( d ) show that the exponents of the @xmath46 do not all have the same value , with the two distributions relative to the end intervals of the reverse transformation actually not possessing power - law character . as a further aspect of the asymmetry in the two phase change directions , we have from figs . [ spikes_delta_epsilon_mean ] and  [ non - stationarity ] that the recorded above - threshold avalanches not only are more numerous , but are also typically larger on loading vs.  unloading . we have performed a full - field study , via the grid method , of the progress of the reversible cubic - monoclinic martensitic transformation in a cu - al - be shape - memory single crystal . the phase change was driven by a stress - controlled loading device based on gravity capable of very small force increments , which applied to the sample a low - rate , monotonic load at largely isothermal conditions . we have thus recorded the stress - strain hysteresis in the forward and reverse transformation , and of the corresponding strain evolution in the microstructure , i.e.  in the spatial assemblies of the phases / variants in the sample . this was done by following the evolution in time of the surface strain maps in the sample obtained from the deformed grid data .    in general , due to the breaking of symmetry in the phase transition , the martensite exhibits multiple variants , with different orientations with respect to the parent lattice , and the detailed understanding of the microstructure is of wide interest in mechanics and materials science , as it is a key element influencing the macroscopic properties and behavior of a multiphase crystal @xcite . however , the design of the loading device sought to obtain for this study the simplest microstructures in a sma single crystal , involving only the parent austenite and one variant of martensite . the very slow driving allowed us to highlight for the first time in a quantitative way the strain intermittency that marks the martensitic transformation in smas . we have characterized the non - smooth progress of the transformation by first studying the intermittency for both the global ( body - averaged ) and local ( pixel - level ) strain measures . we have thereafter identified the organization of possibly long - range events , i.e.  strain avalanches , and studied their statistical properties ( see also the supplementary videos @xcite ) . many earlier experiments based on the acoustic emission accompanying the martensitic transformation in smas have shown some decades of avalanche - size scaling in the austenite - martensite dynamics @xcite . the statistics that we have now uncovered for our strain events indicate that power - law behavior pertains also to a broad distribution of strain - avalanche intensities , covering almost six decades well into the macroscopic sample scales . these results extend and make more quantitatively precise the results recently reported in @xcite about the possible scale - free behavior of the surface optical activity in sma martensitic transformation . we interpret this strain phenomenology as deriving from the evolution of the crystalline substance within a complex energy landscape whose local minimizers represent distinct microstructural morphologies @xcite . the strain avalanches evidenced here mark the progress of the material through such minimizers under the loading @xcite . the associated display of bursty scale - free behavior , as we presently observe , may confirm the indications about the critical nature of the reversible martensitic transformation in smas  @xcite , in analogy to what is observed also in the plastic flow of many crystals @xcite . we have also described a number of other aspects of the observed phase change , evidencing especially the non - stationarity of the material response along the transformation plateaus , and the asymmetric behavior characterized in the forward vs.  the reverse transformation process . further experimental work along the lines developed in the present study can be envisaged in various directions . this should advance our understanding , analysis , and possible control of phase coexistence and its evolution under changing external conditions , and may aid in the design of materials with targeted properties .    . we thank e. vives for several conversations . we acknowledge funding of the italian prin contract 200959l72b004 . nb acknowledges support from the italian group of mathematical physics gnfm - indam . the online supplemental material @xcite accompanying this paper comprises four video files illustrating a number of different aspects of the phase transformation . hereafter we give some details on the experimental methods . [ figure_1 ] shows schematically the cual@xmath0be@xmath1 ( wt .  % ) single crystal sample and the testing apparatus . the sample dimensions are @xmath47 mm@xmath48 respectively along the @xmath5- , @xmath6- , and @xmath18-directions in fig . [ figure_1 ] , the vertical distance between the clamped aluminum tabs being @xmath49 mm . the monitored gauge region has dimensions @xmath3 mm @xmath4 38 mm respectively along the @xmath5 and @xmath6 axes . the rotation matrix giving the orientation of the austenitic cubic single crystal of cu - al - be with respect to the reference axes was measured by x - ray diffraction , and is given by , in the same reference frame of fig . [ figure_1 ] ( see also @xcite ) :      we measured the in - plane displacement , the linear strain components and the local rotation on the gauge region of the sample through the grid method . this technique derives these quantities from the images of a grid captured by a camera as the sample deforms under the loading . grids are transferred onto the sample surface using the procedure described in  @xcite . the grid pattern ( square , with pitch @xmath51  mm ) was printed on a thin polymeric substrate using a 12,000-dpi high - definition printer and then laid on a thin adhesive layer ( e@xmath52 glue from epotecny , whose white color optimizes the visual contrast with the black grid lines . after curing for about @xmath53h at @xmath54c , the polymeric substrate is carefully removed , but the black ink of the grid lines remain bonded on the thin adhesive layer and the sample becomes ready for testing . the deformation of the grid can reasonably be assumed to be the same as that of the surface of the sample as the glue layer is very thin ( some tenths of mm ) , and has negligible stiffness compared to that of the sample , making this a non - intrusive measurement technique . the maximum strain that can be sustained by the grid without debonding or cracking is about @xmath55 , which can well accommodate any strain reached during the martensitic transformation in our alloy .    in our experiment , the images of the deformed grid during the tests were captured by a sensicam qe camera featuring a @xmath56-bit/@xmath57 pixel sensor and a @xmath58-mm sigma lens . the camera was secured on a mounting plate whose position was adjusted so that the lines of the grid were parallel to the pixels of the sensor along both the horizontal and vertical axes . an appropriate magnification allowed to clearly distinguish the black lines of the grid . this magnification was adjusted so that one grid pitch was encoded with 5  pixels . three movable light guides fed by a kl @xmath59 lcd cold light source were used during image acquisition to obtain a nearly - uniform lighting of the grid . about 20,000 images of the deforming sample were obtained over the whole duration of the test , one about every 8 sec , with 15-min pauses about every 100min for data recording and refilling the water reservoirs which feed the pumps . as the grids are the superimposition of two perpendicular straight line patterns , each of them can be considered as a periodic marking , whose change due to the deformation during the test is proportional to the sought in - plane displacement  @xcite : @xmath60 where @xmath61 and @xmath62 are the local phases along @xmath5 and @xmath6 , respectively . in practice , @xmath61 and @xmath62 are deduced from the grid images by calculating at any pixel the windowed fourier transform ( wft ) , see  @xcite , for which a gaussian window @xcite was employed in the present study . the possible appearance of local fictitious strains induced by slight grid printing defects was strongly reduced by using a motion compensation technique described in @xcite . the above measurement and data - treatment methods provide strain maps with a good compromise between strain resolution and spatial resolution , a key factor for distinguishing details in the strain maps . the spatial resolution of the technique , i.e.  the shortest distance between two independent measurements , is considered here to the equal to the width of the gaussian envelope used in the wft , which we consider equal to @xmath63 according to the classical 3-@xmath64 rule . as we set @xmath65 , which is the smallest value that can be chosen for @xmath64  @xcite , and as 5  pixels encode one grid pitch  @xmath66 , we obtain a spatial resolution equal to @xmath67 pixels in this study , that is , 1.2  mm on the sample . random noise propagating from the camera sensor is one of the main sources of disturbance in the strain maps . considering only the latter , and taking into account that the camera was set to average 128 frames to provide one grid image , a threshold equal to @xmath68 for @xmath35 has been estimated @xcite , and has been used to separate the avalanches within the noise floor in sect . [ section_4 ] .  the same value has also been considered for the strain components @xmath13 , @xmath16 and @xmath17 . correspondingly , for the average strain components and their variations , such as @xmath8 , @xmath30 , etc . , computed over a whole strain map , we have a lower threshold , as we rely on  @xmath69 measurements to calculate a single global quantity for each field . the obtained value was then rounded up to be equal @xmath70 , which was used as threshold for the analysis in sect . [ section_2 ] . the strain - increment maps show in some cases the presence of parasitic fringes , probably originating from a slight moir effect between grid and camera sensor . such quasi - periodic noise , for which there is at present no automatic elimination procedure available in the case of strain distributions @xcite , may shade some actual strain - map features , thus spuriously increasing the number of apparently independent strain avalanches . this undesirable effect constitutes one of the limits of the present analysis . we have however heuristically performed avalanche identification also by using thresholds for @xmath35 lower than the value @xmath68 pertaining to the sole sensor noise . we checked directly that this only minimally affects the statistical results reported in figs . 10 - 13 , confirming their robustness . see supplemental material online : # 1 evolution of the strain maps in the forward transformation https://dl.dropboxusercontent.com/u/24282280/video_1.asf ; # 2 evolution of the strain clustering in the forward transformation https://dl.dropboxusercontent.com/u/24282280/video_2.asf ; # 3 strain avalanching in the forward transformation https://dl.dropboxusercontent.com/u/24282280/video_3.asf ; # 4 strain avalanching in the reverse transformation https://dl.dropboxusercontent.com/u/24282280/video_4.asf . rosinberg and e.  vives , in `` disorder and strain induced complexity in functional materials '' ( t.  kakeshita , t.  fukuda , a.  saxena , a.  planes , eds . ) , springer series in materials science * 148 * , springer verlag ( 2012 ) .
Abstract:

Answer: we study experimentally the intermittent progress of the mechanically - induced martensitic transformation in a cu - al - be single crystal through a full - field measurement technique : the grid method . we utilize an in - house especially designed gravity - based device , wherein a system controlled by water pumps applies a perfectly monotonic uniaxial load through very small force increments . the sample exhibits hysteretic superelastic behavior during the forward and reverse cubic - monoclinic transformation , produced by the evolution of the strain field of the phase microstructures . the in - plane linear strain components are measured on the sample surface during the loading cycle , and we characterize for the first time the strain intermittency in a number of ways , showing the emergence of power - law behavior for the strain avalanching over almost six decades of magnitude . we also describe the non - stationarity and the asymmetry observed in the forward vs.  the reverse transformation . the present experimental approach , which allows for the monitoring of the reversible martensitic transformation both locally and globally in the crystal , proves useful and enhances our capabilities in the analysis and possible control of transition - related phenomena in shape - memory alloys .


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: with the progress in short  pulse laser technology @xcite high intensity electromagnetic fields are becoming broadly available . in particular , they make it possible to create strongly correlated quantum plasmas under extreme nonequilibrium conditions which opens a broad range of applications , e.g. @xcite . at the same time , optical techniques for time  resolved diagnostics are improving remarkably @xcite . these developments create the need for a quantum kinetic theory of dense nonideal plasmas in intense laser fields . nonequilibrium properties of dense plasmas in which collisions are important are usually studied on the basis of kinetic equations of the boltzmann type . however , in spite of their fundamental character , boltzmann  like kinetic equations have a number of shortcomings , in particular in view of their application to dense plasmas in intense laser fields :    i ) : :    they are valid only for times larger than the correlation ( or    collision ) time @xmath0 , ii ) : :    they conserve only the mean kinetic energy instead of the sum of    kinetic and potential energy , iii ) : :    they are valid only in the weak field limit since the corresponding    collision integrals are independent of the electromagnetic field , iv ) : :    they are not applicable to high  frequency processes ( fields ) , where    @xmath1 , cf . i ) .    obviously , in the case of strong correlations , high  frequency electromagnetic fields and/or short  time phenomena generalizations are necessary . _ generalized kinetic equations for correlated plasmas _ have been derived already in the 60ies by prigogine @xcite , zwanzig @xcite , kadanoff and baym @xcite , balescu @xcite , silin @xcite , klimontovich @xcite and others . in recent years , the increasing interest in ultrafast processes has revived the theoretical activities , e.g. @xcite , acompanied by progress in numerical solutions , e.g. @xcite ; for textbook overviews , see @xcite . furthermore , kinetic equations for classical plasmas in _ high  frequency fields _ have been derived for the first time in papers of silin @xcite . among other problems , he computed the high  frequency conductivity of a plasma . we mention that , for the weak  field limit , this problem has been studied by many authors , including oberman et al . @xcite and dubois et al . an essential further development of the theory has been given by klimontovich and co  workers @xcite . klimontovich used his powerful technique of second quantization in phase space @xcite to investigate the density  density and microfield fluctuations in low and high  frequency fields . this allowed him to derive collision integrals for classical plasmas in strong fields which take into account dynamical screening and to derive a complete theory of transport processes @xcite . a result of central importance is an expression for the collisional heating rate and the electron  ion collision frequency in strong fields in terms of the imaginary part of the inverse dielectric function @xmath2 @xcite . recently , expressions of the same form were derived again @xcite . for a recent overview on the collision frequency in laser plasmas , we refer to mulser et al . @xcite . the above kinetic theories for plasmas in electromagnetic fields were limited to classical plasmas . a first extension to quantum plasmas was given by silin and uryupin @xcite . more recently , a kinetic equation for dense quantum plasmas in strong static fields has been derived @xcite , whereas a systematic quantum kinetic theory for plasmas in strong fields of arbitrary time dependence was presented in ref . @xcite . there , electron  electron and electron  ion quantum collision integrals were derived within the static born approximation ( landau collision integrals ) .    in this paper , we extend this theory to the case of full dynamical screening for dense quantum plasmas under arbitrary nonequilibrium conditions by using the random phase approximation ( polarization approximation ) for the particle  particle scattering processes . our approach is based on the nonequilibrium green s functions formalism which allows for the most straightforward derivation and for an explicit solution of the gauge problem . we derive a kinetic equation which is a generalization of klimontovich s classical result @xcite and , on the other hand , generalizes previous quantum results for the case of zero field @xcite and static field @xcite . furthermore , we derive results for the polarization and screening properties and longitudinal field fluctuations in a strong field of arbitrary time  dependence . we consider the time evolution of a dense charged particle system under the influence of a strong time  dependent electromagnetic field and inter  particle correlations . it is instructive to recall first the motion of classical free charges in an external field . from integrating newton s equation , @xmath3 , we obtain the velocity change of a particle with charge @xmath4 and mass @xmath5 in the field @xmath6 during a time interval @xmath7 $ ] , @xmath8 and the field induced displacement @xmath9 in the equations above , we have dropped contributions from the acceleration and velocity at the initial moment @xmath10 since they are not related to the field . a further important quantity is the average kinetic energy @xmath11 where , in case of a periodic field , @xmath12 is the oscillation period . below we will also need the change of the relative velocity of a particle pair @xmath13 gained in the field @xmath14 @xmath15 and the change of the inter  particle distance @xmath16 @xmath17    the most important special case is that of a _ harmonic time dependence _ @xmath18 which leads to the following explicit results for the quantities introduced above @xmath19 , \label{va - h } \\ \delta{\bf r}_a(t , t ' ) & = & { \bf r}_a^0 \left[\omega ( t - t')\sin \omega t + \cos \omega t - \cos \omega t'\right ] . \label{va - ra - h}\end{aligned}\ ] ] here , we introduced two important quantities , the `` quiver '' velocity @xmath20 and the so  called excursion amplitude @xmath21 furthermore , the cycle averaged ( @xmath22 ) kinetic energy gain of a charged particle ( [ e - kin - av ] ) is the so  called ponderomotive energy @xmath23    obviously , the above results may be extended to electromagnetic fields with arbitrary time dependence , for example by expanding the field in terms of harmonic components . nevertheless , it is useful to explicitly consider a second situation frequently encountered in modern applications : pulsed fields , e.g. those produced by femtosecond lasers . we will consider pulses of the following form @xmath24 and @xmath25 otherwise . typically , @xmath26 , although modern femtosecond laser pulses may be as short as a few periods of the main frequency @xmath27 . for the field ( [ e - pulse ] ) , we obtain the velocity change and displacement @xmath28 , \quad { \bf v}_a^{0\pm } \equiv \frac{e_a { \bf e}_0}{m_a \omega^{\pm } } , \quad \omega^{\pm } \equiv \omega_p \pm \omega , \\ \delta{\bf r}_a(t , t ' ) & = &   \sum_{s=\pm}{\bf r}_a^{0s } \left[-\omega^s ( t - t')\cos \omega^s t + \sin \omega^s t - \sin \omega^s t'\right ] , \ ] ] with @xmath29 , whereas the change of relative velocity and two  particle distance , @xmath30 and @xmath31 follow from @xmath32 and @xmath33 by replacing @xmath34 by @xmath35 and @xmath36 by @xmath37 , respectively . the above results trivially include the case of a time - independent electric field which is recovered by letting in eq . ( [ efeld ] ) @xmath38 . the corresponding results are @xmath39    it is instructive to consider a number of parameters which characterize the state of the plasma , field strength , quantum properties etc :    1 . the field strength can be characterized by the ratio @xmath40 of the amplitude of the oscillation velocity ( quiver velocity ) @xmath41 , eq . ( [ va - h ] ) , to the thermal velocity @xmath42 . the relative importance of the field and of particle  particle interaction is characterized by @xmath43 , where @xmath44 is the amplitude of the field induced displacement ( [ va - ra - h ] ) , and @xmath45 is the debye radius ( interaction range ) . 3 .   the frequency of the field has to be compared to the eigenfrequencies of the plasma , most importrantly , the electron langmuir ( plasma ) frequency , @xmath46 , which reflects competition between field frequency and plasma density effects . the relevance of collisional processes depends on the ratio @xmath47 , where @xmath48 is the total collision frequency of electrons in the plasma . the photon energy is characterized by its ratio to the thermal energy , @xmath49 . modern lasers easily produce strong fields which satisfy the inequalities @xmath50 and @xmath51 . in high  frequency fields and/or plasmas of moderate density , @xmath52 , which allows to treat collisions perturbatively , see below . coulomb interaction between the charged carriers as well as quantum effects , obviously , may drastically modify the free particle behavior . scattering of two particles with charges @xmath53 in quantum states @xmath54 and @xmath55 on the coulomb potential @xmath56 leads to a transfer of momentum @xmath57 between them , so after the collision time @xmath58 the particles are in momentum states @xmath59 and @xmath60 . while conventional kinetic approaches treat collision as instantaneous , @xmath61 , this is not appropriate for correlated plasmas as well as in the presence of rapidly varying fields with @xmath62 not being small . in this case , during the collision time , the scattering partners will be accelerated by the external field which is called _ intra  collisional field effect _ , which essentially modifies the scattering process . using a quantum language , during the collision time , the particles may absorb photons of the electromagnetic field which is the familiar _ inverse bremsstrahlung _ , or re  emit them ( bremsstrahlung ) .    the kinetic treatment of two  particle scattering on the coulomb potential ( [ coulomb ] ) leads to the well  known divergencies at short and long wavelengths . while the first is naturally cured by a quantum theoretical approach , the origin of the latter is the long range of the coulomb interaction . the familiar solution lies in the replacement of the bare coulomb potential ( [ coulomb ] ) by a screened one @xmath63 where @xmath64 is the retarded dielectric function . this is not only of fundamental interest but has also practical relevance . the dielectric function includes collective plasma oscillations and instabilities which , especially in nonequilibrium situations , may strongly enhance scattering , transport and energy exchange with the electromagnetic field ( anomalous transport ) . the simplest approximation for the dielectric function is the random phase approximation ( rpa ) being the quantum generalization of the vlasov dielectric function @xmath65 where @xmath66 denotes the single  particle energy . interestingly , this result was derived by klimontovich and silin @xcite two years before lindhard @xcite . while this approximation is applicable to nonequilibrium situations in which the wigner distributions @xmath67 are weakly time  dependent ( when the time scale @xmath68 is much longer than @xmath69 ) , on short times or for fast processes , such as in high frequency fields , generalizations are necessary . such generalizations avoid the assumption of separation of the two time scales and lead to an explicit dependence of the functions @xmath70 and @xmath71 on two times . moreover , the presence of an electromagnetic field may be expected to modify the dielectric and screening properties , leading to a very complex problem of coupled particle , screening and field dynamics . the appropriate theoretical concept to tackle this problem is provided by quantum field theory . numerous concepts have been developed to describe the mentioned above dynamics of particles and fields . among them , the most systematic and powerful is the theory of nonequilibrium green s functions . it is based on the method of relativistic quantum field theory , where charged particles and the longitudinal and transverse electromagnetic field are described on equal footing by field operators @xcite . from the equations of motion for the field operators - the dirac equation and maxwell s equations , one can derive equations of motion for all quantities of interest . among them , the most important are two  time correlation functions ( green s functions ) which allow for systematic and far  reaching generalizations of traditional kinetic theory . in this paper , we focus on nonrelativistic particle dynamics and start our derivations from the familiar equations of motion for the particle correlation functions @xmath72 and @xmath73 while the electromagnetic field is treated classically .      the field theoretical description of plasmas is based on the creation and annihilation operators @xmath74 and @xmath75 @xcite which are defined to guarantee the spin statistics theorem , @xmath76 where @xmath77 has been assumed . the upper ( lower ) sign refers to bosons ( fermions ) , @xmath78 , and @xmath79 labels the particle species . below , we will drop the spin index and assume fermions . the nonequilibrium state of a correlated plasma is described by the two  time correlation functions which are statistical averages ( with the initial density operator of the system ) of field operator products @xmath80 where @xmath81 and @xmath82 are , in nonequilibrium , independent from one another . they contain the complete dynamical and statistical information . the latter follows from their elements _ along the time diagonal _ : the one - particle density matrix is immediately obtained from the function @xmath73 according to @xmath83 whereas the dynamical information ( e.g. the single  particle spectrum and the correlations ) follows from the function values _ across the diagonal _ in the @xmath84plane , in particular , from the _ spectral function _ @xmath85 , @xmath86 where @xmath87 are the retarded and advanced green s functions , defined below in eq . ( [ gra - def ] ) . in the following , it will often be convenient to use microscopic and macroscopic time and space variables being defined as @xmath88 in particular , in cases where the microscopic variables vary on much smaller scales than the macroscopic ones , it is advantageous to perform a fourier transformation with respect to @xmath89 and/or @xmath90 which leads to the frequency and momentum variables @xmath91 and @xmath92 , respectively . in particular , eq . ( [ f - def ] ) then yields the familiar wigner distribution function @xmath93    the time evolution of the correlation functions in an electromagnetic field is determined by the kadanoff  baym equations @xcite @xmath94g_a^{\gtrless}(1 , 1 ' ) - \int d{\bar { \bf r}}_1 \ , \sigma_a^{\rm hf}(1,{\bar { \bf r}}_1t_1 ) g_a^{\gtrless}({\bar{\bf r}}_1t_1,1 ' ) \nonumber\\ & & = \int_{t_0}^{t_1 } d{\bar 1 } \ , \left[\sigma_a^>(1,{\bar 1})-\sigma_a^<(1,{\bar 1})\right ] \ , g_a^{\gtrless}({\bar 1 } , 1 ' ) - \int_{t_0}^{t'_1 } d{\bar 1 } \ , \sigma_a^{\gtrless } ( 1,{\bar 1})\,\left[g_a^>({\bar 1 } , 1')- g_a^<({\bar 1 } , 1')\right ] , \label{kb_eq}\end{aligned}\ ] ] which have to be fulfilled together with the adjoint equations . here , @xmath95 denotes the initial time where the system is assumed to be uncorrelated ( otherwise , the equations have to be supplemented with an initial correlation contribution to @xmath96 , cf . @xmath97 is the hartree  fock selfenergy ( mean  field energy with exchange ) , @xmath98 and @xmath99 are the correlation selfenergies which will be discussed below .    for the following derivations , it is useful to introduce , in addition , the retarded and advanced green s functions @xmath100\left\{g_a^>(1,1')-g_a^<(1,1')\right\ } , \label{gra - def}\ ] ] which obey the simpler equations @xmath101 g_a^{r / a}(1,1 ' ) \nonumber\\ & & -\int d2 \,\sigma_a^{r / a}(1,2)g_a^{r / a}(2,1')\,\ , = \,\ , \delta(1 - 1 ' ) . \label{gra_eq}\end{aligned}\ ] ] in equations ( [ kb_eq ] ) and ( [ gra_eq ] ) , the electromagnetic field is given by the vector and scalar potentials @xmath102 and @xmath103 and will be treated classically . @xmath102 is the full vector potential ( external plus induced ) which obeys maxwell s equations , whereas @xmath103 is understood as to be due to external sources only , the induced longitudinal field is fully accounted for in the screened coulomb potential @xmath71 which enters the selfenergies @xmath104 and @xmath105 , see below .    although one can directly analyze and solve the two  time kadanoff  baym equations ( [ kb_eq ] ) , e.g. @xcite , for an overview , see @xcite , it is easier to consider the kinetic equation for the wigner distribution function ( [ f - def2 ] ) , which we will be concerned with here . this equation is immediately obtained from the equal time limit , @xmath106 , of eq . ( [ kb_eq ] ) plus its adjoint . introducing further the variables @xmath107 and @xmath90 , eqs . ( [ r - r - def ] ) and ( [ t - tau - def ] ) , we obtain , after fourier transformation with respect to @xmath90 , for the spatially homogeneous case , @xmath108 this is an exact equation and , therefore , well suited for deriving generalized kinetic equations . however , this equation is not closed yet since it contains under the collision integral functions depending on two times . therefore , to obtain explicit expressions for the collision integral , one has :    1 . to find appropriate approximations for the self energy . for this , the green s functions approach provides powerful approximation schemes based on feynman diagrams which allow for a very systematic development of the theory . here , we are interested in the plasma dynamics with screening effects properly included . therefore , the appropriate choice for the selfenergy will be the random phase approximation ( rpa ) ; 2 .   to express the correlation functions @xmath109 as functionals of the wigner functions @xmath67 ( reconstruction problem ) . this problem can be solved approximately on the basis of the generalized kadanoff  baym ansatz ( gkba ) of lipavsk et al . @xcite , see sec . [ gkba_ss ] . it is well known that the electromagnetic field can be introduced in various ways ( gauges ) which may lead to essentially different explicit forms of the resulting kinetic equations . although alternative derivations are successfully applied too , gauge invariance becomes a particular problem if the resulting kinetic equations are treated by means of approximations , such as retardation or gradient expansions . a critical issue is that the result of these approximations maybe essentially different in different gauges , see e.g. @xcite for examples . to avoid these difficulties , we will formulate the theory in terms of correlation functions which are made explicitly gauge  invariant . while the main results have been presented in ref . @xcite , here we provide some additional details .    in this section , we use a co - variant 4-vector notation as it makes the following transformations more compact and symmetric . the corresponding definitions are @xmath110 and the convention @xmath111 is being used . one readily proofs that the kadanoff  baym equations ( [ kb_eq ] ) remain covariant under gauge transformations , i.e. , under the following transformations of the potentials and field operators @xmath112 the corresponding gauge transform of the green s functions leads to @xmath113 } g_a(x , x)\ , . \ ] ] following an idea of fujita @xcite , we now introduce a gauge  invariant green s function @xmath114 which is given by the modified fourier transform @xmath115\right\ } g_a(x , x ) , \label{gauge - ft}\end{aligned}\ ] ] where use has been made of the identity @xmath116 indeed , one readily confirms that under any gauge transform ( [ gauge - a ] ) , the phase factors cancel , and @xmath117 , @xcite .    in the following , we focus on spatially homogeneous electric fields and use the vector potential gauge @xmath118 in this case , relation ( [ gauge - ft ] ) simplifies to    @xmath119 g_a({\bf r},\tau;{\bf r},t ) , \label{gauge - ftv}\end{aligned}\ ] ]    what means that the gauge  invariant green s function @xmath120 follows from the wigner transformed function @xmath121 by replacing the canonical momentum @xmath92 by the gauge  invariant kinematic momentum @xmath122 according to @xmath123 let us illustrate this for the examples studied in sec . [ sec_bas ] . for a _ harmonic electric field _ given by eq . ( [ efeld ] ) , the vector potential and the momentum relation become , according to eq .  ( [ v - gauge ] ) , @xmath124 similarly , for a _ pulsed field _ of the form eq . ( [ e - pulse ] ) , the result for the vector potential and the momenta is @xmath125 finally , for a _ static field _ we obtain , @xmath126    for the derivations below , we will need the gauge invariant fourier transform of the convolution of two functions which , in the homogeneous case , is given by @xmath127 after straightforward manipulations which involve the back transform of ( [ gauge - ftv ] ) , we arrive at @xmath128 \nonumber \\ & & \qquad \times\,\ , c \left[{\bf k}+ \frac{e_a}{c}\int_{t'_1}^{t_1 } dt''\ , \frac{{\bf a } ( t'')}{t_1-t'_1 } - \frac{e_a}{c } \int^{\bar{t}}_{t'_1 } dt''\ , \frac{{\bf a}(t '' ) } { \bar{t}-t'_1};\,\bar{t},t'_1\right ] . \quad \label{conv - ft}\end{aligned}\ ] ] in particular , the derivation of the collision integral in the kinetic equation for the wigner function , requires the equal  time limit of this expression , @xmath106 , @xmath129 \nonumber \\ & & \qquad \times\,\ , c \left[{\bf k}+\frac{e_a}{c}{\bf a}(t ) - \frac{e_a}{c}\int_{\bar t}^{t } dt '' \ , \frac{{\bf a}(t'')}{t-{\bar t } } ; \,\bar{t},t\right ] , \quad \label{conv - ftt}\end{aligned}\ ] ] where @xmath130 and @xmath131 will be replaced by @xmath132 and @xmath133 . notice that in this case the momentum arguments of @xmath130 and @xmath131 are equal . to simplify the notation below , we introduce the field induced momentum shift @xmath134 which has the important property @xmath135 where @xmath136 is nothing but minus the momentum gain of a free particle in the field , @xmath137 , cf . ( [ va ] ) . another important relation follows from multiplication by the time interval : @xmath138 where @xmath139 is just the field induced displacement of a free particle , eq . ( [ ra ] ) , i.e. @xmath140 .    definition ( [ ka ] ) allows us to rewrite eq . ( [ conv - ftt ] ) as @xmath141 \ , c \left[{\bf k}+{\bf k}^a_a(t,{\bar t});\ , t,{\bar t}\right ] . \label{conv - ftt2}\end{aligned}\ ] ]      as noted in sec .  [ kbe_ss ] , for the derivation of the collision integral in the kinetic equation , we need to express the functions @xmath72 and @xmath73 in terms of the wigner function . in addition , such a reconstruction ansatz involves the retarded and advanced green s functions @xmath142 for which suitable expressions have to be found . we first determine these quantities for free particles in an electromagnetic field which allows to simplify eq . ( [ gra_eq ] ) to @xmath143 g_{a}^{r / a}({\bf p};t_1,t'_1)=\delta ( t_1-t'_1),\ ] ] which is solved immediately by @xmath144 ^ 2/2m_a\right ] , \label{gra_free}\end{aligned}\ ] ] and @xmath145 is obtained from the symmetry relation @xmath146^*$ ] . from this result , we can calculate the spectral function @xmath147 , eq . ( [ a - def ] ) , @xmath148 ^ 2/2m_a\right ] . \label{a_free}\end{aligned}\ ] ] obviously , the results ( [ gra_free ] ) and ( [ a_free ] ) are gauge - dependent since @xmath149 and @xmath150 are functions of the canonic momentum @xmath92 . but one can easily obtain the corresponding gauge  invariant results by applying the transform ( [ gauge - ft ] ) , with the result @xmath151 } , \label{gra_fg } \\ \nonumber\\ \mbox{where}\quad s_a({\bf a};\tau , t ) & = & \frac{e_a^2}{2 m_a c^2 } \left [ \int_{t-\frac{\tau}{2}}^{t+\frac{\tau}{2 } } dt ' a^2(t ' ) + \frac{1}{\tau } \left(\int_{t-\frac{\tau}{2}}^{t+\frac{\tau}{2 } } dt ' { \bf a}(t')\right)^2\right]\,.\end{aligned}\ ] ] this result has a simple physical interpretation . for a free particle without field , the spectral function shows free undamped oscillations along @xmath89 ( i.e. perpendicular to the time diagonal ) with the one  particle energy @xmath152 , and its fourier transform is @xmath153 . \label{afree}\ ] ] this clearly underlines the meaning of the functions @xmath87 and , in particular , the spectral function @xmath79  they contain the full information on the single  particle energy spectrum . furthermore , in a correlated system , the single  particle spectrum is affected by interactions with other particles . this leads to a shift of the oscillation frequency and to damping of the oscillations , i.e. to finite life time effects , and it is reasonable to call the corresponding single  particle excitations _ quasi  particles_. on the other hand , the result ( [ gra_fg ] ) reflects the influence of an electromagnetic field on the particle spectrum , while correlation effects have been neglected . ( [ gra_fg ] ) shows that the field causes a time  dependent shift of the single  particle energy which , obviously , reflects the well  known fact that the proper eigenstates of the system contain the electromagnetic field and are given by volkov states @xcite . the spectrum may even contain additional peaks which becomes particularly transparent in the limiting case of a harmonic time dependence : for the field ( [ efeld ] ) , the time integrations in @xmath154 can be performed , and simple trigonometric relations lead to @xcite @xmath155 , \label{gra_fgh}\end{aligned}\ ] ] where @xmath156 is the ponderomotive potential which was introduced in eq . ( [ e - pond - h ] ) . the first term in the brackets leads to a shift of the single - particle energy , the average kinetic energy of the particles increased by @xmath156 . the remaining terms modify the spectrum qualitatively giving rise to additional peaks which are related to _ photon sidebands _ @xcite . now we turn to the solution of the reconstruction problem . the simplest solution is the common _ kadanoff  baym ansatz _ , @xmath157 where @xmath158 and @xmath159 , and the upper ( lower ) sign refers to @xmath72 ( @xmath73 ) . indeed , the two  time functions @xmath160 and @xmath161 are now expressed in terms of one  time wigner distribution functions and a known spectral function . however , due to the expected retardation effects , this ansatz is not applicable here . as mentioned above , a more general solution which properly takes into account retardation ( memory ) effects is the _ generalized kadanoff  baym ansatz _ proposed by lipavsk et al . @xcite which reads @xmath162 where for the functions on the time diagonal @xmath163 , cf . ( [ f - def2 ] ) . within the quasiparticle approximation and with static selfenergies , eq .  ( [ gkba ] ) is exact . in more complex situations , it is an approximation to the exact reconstruction solution , which has prooved extremely successful in many applications . in particular , it has been used for more general selfenergies and also with more general propagators @xmath164 , e.g. @xcite . we , therefore will use this ansatz below . ( [ gkba ] ) is written in terms of the momentum @xmath92 and is , therefore , gauge  dependent . to transform this relation into a gauge  invariant form , we use its coordinate representation , @xmath165 and apply the gauge  invariant fourier transform ( [ gauge - ftv ] ) together with the back transforms of @xmath164 and @xmath166 which leads to the _ gauge  invariant generalization of the gkba _ @xmath167 - f_a^{\gtrless}\left[{\bf k}-{\bf k}^a_a(t , t ' ) ; t_1\right ] \ , g_a^a({\bf k};t_1,t'_1 ) , \label{gkba_gi}\end{aligned}\ ] ] where the definition ( [ ka ] ) for @xmath168 has been used . as in the field  free case , the first term is nonzero only for @xmath169 and the second in the opposite case . notice the difference of the time arguments in the two distribution functions . we now come back to the time  diagonal limit of the kadanoff  baym equations , cf . ( [ kb_eq ] ) and ( [ f_eq ] ) , and derive the quantum kinetic equation for a plasma in a laser field thereby fully taking into account dynamical screening . again , it is advantageous to derive this equation for the gauge  invariant wigner distribution . to this end , we take the fourier transform ( [ gauge - ftv ] ) of the time  diagonal kadanoff  baym equation @xmath170 where the full arguments are , according to the convolution relation ( [ conv - ftt2 ] ) , given by @xmath171 \ , g_a^{\lessgtr}\left[{\bf k}_a+{\bf k}_a^a(t,{\bar t});{\bar t},t\right ] . \nonumber\end{aligned}\ ] ] this expression is valid for arbitrary approximations for the selfenergies @xmath172 and @xmath173 . in our previous paper @xcite , we used the simple static born approximation . here , we are interested in a fully selfconsistent _ inclusion of dynamical screening _ , so the appropriate choice is the random phase approximation ( rpa ) .      starting from the familiar expression in coordinate representation , application of the gauge  invariant fourier transform ( [ gauge - ftv ] ) straightforwardly leads to the following gauge  invariant result @xmath174 which transforms the collision integral of eq.([f_eqgi ] ) into @xmath175 \nonumber\\ & \times & i\hbar v_{aa}^{s\,>}({\bf q};t,{\bar t})\ , g_a^<\left[{\bf k}_a+{\bf k}_a^a(t,{\bar t});{\bar t},t\right]\ , \quad - \quad [ \;>\ ; \longleftrightarrow \;<\ ; ] \bigg\}\end{aligned}\ ] ]    in eq . ( [ sigma - rpa ] ) we introduced the correlation functions of the screened potential ( plasmon correlation functions ) @xmath176 which contain the whole screening problem and are directly related to the correlation function of the _ longitudinal field fluctuations _ ( microfield fluctuations ) @xcite @xmath177 . \label{dede}\end{aligned}\ ] ] while in the classical case , the contributions from @xmath178 and @xmath179 are equal , in the quantum case a symmetrization is useful which is indicated by the bar over the fluctuation term . @xmath180 and @xmath181 can be related to the retarded and advanced screened potentials via the optical theorem @xmath182 where @xmath183 and @xmath184 obey the following equation of motion ( dyson equation ) @xmath185 and @xmath186 are related to the nonequilibrium inverse dielectric function according to @xmath187^{-1}. \label{idf}\end{aligned}\ ] ] in the above equations , @xmath188 is the bare coulomb potential ( [ coulomb ] ) and @xmath189 the retarded and advanced longitudinal polarization functions ( plasmon selfenergies ) . to close this system of equations , the polarization functions have to be expressed in terms of the particle correlation functions for which the simplest approximation is provided by the rpa , @xmath190 and @xmath191 follows from the relation @xmath192^*$ ] . this set of equations completely defines the non  markovian polarization approximation ( rpa ) for a quantum plasma in a strong transverse field . what is left to obtain a closed expression for the collision integral in equation ( [ f_eqgi ] ) is to apply the gauge  invariant gkba ( [ gkba_gi ] ) together with the free  particle approximation ( [ gra_fg ] ) to all two  time functions . this leads to the following results for the optical theorem : @xmath193 } \nonumber\\ & & \times f_b^{\gtrless}\left[{\bf k}+{\bf q}+{\bf q}_b(t_3,t_4);t_4\right]\ , f_b^{\lessgtr}\left[{\bf k}+{\bf q}_b(t_3,t_4);t_4\right ] \nonumber\\ & & + \theta(t_4-t_3)\,e^{-\frac{i}{\hbar}\left [ \left(\epsilon^b_{{\bf k}+{\bf q}}-\epsilon^b_{{\bf k}}\right)(t_3-t_4 ) + { \bf q}{\bf r}_b(t_4,t_3 ) \right ] } \nonumber\\ & & \times f_b^{\gtrless}\left[{\bf k}+{\bf q}+{\bf q}_b(t_4,t_3);t_3\right]\ , f_b^{\lessgtr}\left[{\bf k}+{\bf q}_b(t_4,t_3);t_3\right ] \bigg\}. \label{vs - gkba}\end{aligned}\ ] ] similarly , the result for @xmath194 can be transformed to @xmath195- f_a\left[{\bf k+q}+{\bf q}_a(t , t');t'\right ] \right\ } , \label{pira - gkba}\end{aligned}\ ] ] where the momentum shift @xmath136 and field  induced displacement @xmath139 were defined above in eqs . ( [ ka - dif ] ) and eqs . ( [ ka - dt ] ) , respectively . in the absence of the electromagnetic field , ( @xmath196 ) , eq . ( [ pira - gkba ] ) reduces to the well  known nonequilibrium rpa  polarization function . the effect of the field is two  fold : first , it introduces an additional retardation @xmath136 in the distributions ( intra  collisional field effect ) and second , it leads to a modification of the one  particle energies in the exponent given by @xmath139 which we discussed in detail in ref . @xcite . we now can transform the collision integral eq . ( [ coli1 ] ) by applying the gkba ( [ gkba_gi ] ) to @xmath197 , and using the result for @xmath198 and @xmath199 , eq . ( [ vs - gkba ] ) . after straightforward calculations , we obtain the final result @xmath200}\ , f_a^>\left[{\bar { \bf k}}_a+{\bf q}_a(t,{\bar t});{\bar t}\right ] f_a^<\left[{\bf k}_a+{\bf q}_a(t,{\bar t});{\bar t}\right]\ , \nonumber\\ & \times&\int_{t_0}^t d t_3 \;\int_{t_0}^{\bar t } d t_4 v^{s\,r}({\bf k}_a-{\bar { \bf k}}_a;t_1,t_3)\ , v^{s\,a}({\bf k}_a-{\bar { \bf k}}_a;t_4,t_2)\,\times \nonumber\\ & & \bigg [ \theta(t_3-t_4)\,e^{-\frac{i}{\hbar}\left [ \left({\bar \epsilon}_b-\epsilon_b\right)(t_3-t_4 ) + ( { \bf k}_b-{\bar { \bf k}}_b){\bf r}_b(t_3,t_4 ) \right ] } \nonumber\\ & & \times f_b^>\left[{\bar { \bf k}}_b+{\bf q}_b(t_3,t_4);t_4\right]\ , f_b^<\left[{\bf k}_b+{\bf q}_b(t_3,t_4);t_4\right ] \nonumber\\ & & + \theta(t_4-t_3)\,e^{-\frac{i}{\hbar}\left [ \left({\bar \epsilon}_b-\epsilon_b\right)(t_3-t_4 ) -({\bf k}_b-{\bar { \bf k}}_b){\bf r}_b(t_4,t_3 ) \right ] } \nonumber\\ & & \times f_b^>\left[{\bar { \bf k}}_b+{\bf q}_b(t_4,t_3);t_3\right]\ , f_b^<\left[{\bf k}_b+{\bf q}_b(t_4,t_3);t_3\right ] \bigg ] \quad - \quad [ \;>\ ; \longleftrightarrow \;<\ ; ] \bigg\}.\end{aligned}\ ] ] this is the general non  markovian collision integral for a homogeneous weakly coupled _ dynamically screened plasma in an electromagnetic field_. it is a generalization of numerous results which were previously obtained by various authors , including our static screening result for strong time  dependent fields @xcite and the result of silin and uryupin @xcite , the rpa result for a static field of morawetz @xcite and the field - free rpa results of kuznetsov @xcite and haug and ell @xcite . furthermore , it generalizes previous results obtained for classical plasmas by silin , oberman et al . , klimontovich and others @xcite . in particular , as we will see below , the classical dynamical screening result of klimontovich and puchkov @xcite is straightforwardly recovered from the collision integral ( [ coli3 ] ) . this collision integral is the basis for computing electron  ion collision frequencies , plasma heating inverse bremsstrahlung effects etc . , thereby fully taking into account dynamical screening , plasma instabilities and anomalous transport    despite the complicated structure of the integral ( [ coli3 ] ) , a direct numerical integration of the kinetic equation ( [ f_eqgi ] ) appears to be within reach , as recently solutions of non  markovian rpa - type equations ( without longitudinal field ) for semiconductors have been reported @xcite . on the other hand , to gain deeper insight in the physical processes contained in the collision term ( [ coli3 ] ) , it is useful to consider analytical simplifications . if the collision frequency is low compared to the oscillation frequency of the field , i.e. if the parameter @xmath52 , one may follow an idea of silin @xcite and solve the kinetic equation ( [ f_eqgi ] ) with a perturbation ansatz @xmath201 where @xmath202 obeys the collisionless equation @xmath203 , \label{silin_an}\end{aligned}\ ] ] where @xmath204 is an arbitrary function depending on the initial conditions . the equation for @xmath205 reads @xmath206 \right\}. \label{silin_an1}\end{aligned}\ ] ] with this scheme , there follow essential simplifications of the above results because it effectively eliminates the time retardation of the distribution functions in the collision integrals @xcite . indeed , one easily verifies that the arguments of the distributions which appear in the formulas above , now become @xmath207 \longrightarrow f_{a0}\left[{\bf k}+\frac{e_a}{c}{\bf a}(t)\right ] , \label{farg - silin}\end{aligned}\ ] ] and do not depend on the time @xmath10 anymore . this simplification allows to compute the transport , screening and fluctuation properties quite efficiently . we first consider the modification of the longitudinal polarization @xmath208 . straightforward transformations of eq . ( [ pira - gkba ] ) including a change of the momentum integration variable lead to @xmath209 using this result and the adiabatic approximation ( neglecting the ion contribution to the polarization ) , we obtain from the dyson equation , eq .  ( [ vdyson ] ) , for the retarded screened potential @xmath210 where again , @xmath211 denotes the screened potential in the zero field limit and without retardation in the distribution functions . simplifications are possible also for the collision integrals . in particular , we obtain for the electron  ion scattering term , @xmath212\tau } \nonumber\\ & & \ ; \times \int_{0}^{\infty}d\tau_1\int_{-\infty}^{0}d\tau_2\frac{1 } { \epsilon_{0}^{r}({\bf q},\tau_1)}\frac{1}{\epsilon_{0}^{a}({\bf q},\tau_2 ) } \ , e^{-\frac{i}{\hbar}\frac{e_e}{m_e c}{\bf q } \int\limits_{t-\tau-\tau_2}^{t-\tau_1 } d{\bar t}\,{\bf a}({\bar t } ) } \left\{f_{e0}({\bf k - q})-f_{e0}({\bf k } ) \right\ } , \quad \label{iei - silin}\end{aligned}\ ] ] where @xmath213 are the field  free inverse dielectric functions . ( [ iei - silin ] ) is the generalization of klimontovich s result @xcite to quantum plasmas and . finally , we consider the simplifications introduced by the ansatz ( [ silin_an1 ] ) to the field fluctuations . the quantity of central importance is the correlation function of the temporal and spatial microfield fluctuations @xmath214 , eq . ( [ dede ] ) . a lengthy but straightforward calculation leads to the following result @xmath215(t_1-t_2 ) } } { |\epsilon_{0}({\bf q},\epsilon^{e}_{{\bf k+q}}-\epsilon^{e}_{{\bf k}})|^2 } \left\{f_{e0}({\bf k+q})[1-f_{e0}({\bf k})]+ f_{e0}({\bf k})[1-f_{e0}({\bf k+q})]\right\ } \nonumber\\ & & \qquad + \;e_i^2 \ , e^{-\frac{i}{\hbar}[\epsilon^{i}_{{\bf k+q}}-\epsilon^{i}_{{\bf k}}](t_1-t_2 ) } \int_{0}^{\infty } d\tau_1 \ , \frac{e^{-\frac{i}{\hbar } [ \epsilon^{i}_{{\bf k+q}}-\epsilon^{i}_{{\bf k}}]\tau_1 } } { \epsilon_{0}^{r}({\bf q},\tau_1 ) } \nonumber\\ & & \qquad\quad\times \int^{0}_{-\infty } d\tau_2 \ , \frac{e^{\frac{i}{\hbar } [ \epsilon^{i}_{{\bf k+q}}-\epsilon^{i}_{{\bf k}}]\tau_2 } } { \epsilon_{0}^{a}({\bf q},\tau_2 ) } \ ; e^{-\frac{i}{\hbar}\left(\frac{e_e}{m_e}-\frac{e_i}{m_i}\right){\bf q } \int\limits_{t_2-\tau_2}^{t_1-\tau_1 } d{\bar t}\,\frac{{\bf a}}{c}({\bar t } ) } \left[f_{i0}({\bf k+q})+f_{i0}({\bf k})\right ] \bigg\ } , \nonumber \label{de - silin}\end{aligned}\ ] ] where the first term in parantheses ( second line ) is the electron contribution , and the second ( third and fourth line ) results from the ions . again , this is a generalization of klimontovich s remarkable result @xcite who considered the classical limit and the equal time fluctuations , @xmath77 . our result fully includes the two  time fluctuations which are directly measurable quantities . from the above fluctuation spectrum , all major observables of dense quantum plasmas in a strong laser field can be computed . the corresponding analysis will be presented in a forthcoming paper . in this paper , we have presented a gauge  invariant derivation of the quantum kinetic equation for dense plasmas in a laser field . our main result , eq . ( [ coli3 ] ) , generalizes previous work to quantum systems . this equation can be used to calculate the transport properties of a dense plasma in a laser field on arbitrary time scales , i.e. , over the whole frequency range . the use of the random phase approximation allows for a highly consistent treatment of the combined effect of internal longitudinal fields ( dynamical screening ) and transverse electromagnetic fields , including intense laser pulses . in particular , it allows to investigate the influence of the electromagnetic field on the two  particle scattering process and the screening properties of the plasma and on the screening buildup in the presence of a strong field .    besides , the presented gauge  invariant approach is completely general and can be extended straightforwardly to more complex situations , including strong coupling effects , bound states , impact and field ionization . moreover , it can be directly generalized to relativistic systems und ultra  intense fields . it is a great pleasure to dedicate this paper to yuri lvovich klimontovich on his seventyfith birthday . all of us have , at different times , come in contact with the remarkable papers and monographs of klimontovich . we have very much enjoyed countless discussions with him as well as rewarding scientific collaborations on problems of kinetic theory of nonideal gases and plasmas in general , in particular on dielectric properties @xcite , kinetic equations with bound states @xcite , foundation of reaction  diffusion equations for nonideal plasmas @xcite , selforganization and structure formation @xcite and many others . the authors acknowledge stimulating discussions with g. kalman , p. mulser and r. sauerbrey . this work is supported by the deutsche forschungsgemeinschaft ( schwerpunkt `` wechselwirkung intensiver laserfelder mit materie '' ) , and by the european commission through the tmr network silasi . it is interesting to note that klimontovich s phase space densities @xmath216 @xcite are essentially the classical limit of the field operators . therefore , his kinetic theory is of similar generality as modern field theory , cf . his text books @xcite .
Abstract:

Answer: a quantum kinetic theory for correlated charged  particle systems in strong time dependent electromagnetic fields is developed . our approach is based on a systematic gauge  invariant nonequilibrium green s functions formulation . extending our previous analysis @xcite we concentrate on the selfconsistent treatment of dynamical screening and electromagnetic fields which is applicable to arbitrary nonequilibrium situations . the resulting kinetic equation generalizes previous results to quantum plasmas with full dynamical screening and includes many body effects . it is , in particular , applicable to the interaction of dense plasmas with strong electromagnetic fields , including laser fields and x - rays . furthermore , results for the modification of the plasma screening and the longitudinal field fluctuations due to the electromagnetic field are presented .