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: understanding the energy transfer mechanisms during thermal non - equilibrium between electrons and the lattice is critical for a wide array of applications . non - equilibrium electron properties on time scales of 10 - 100s of femtoseconds are most efficiently observed with pulsed laser measurement techniques.@xcite laser irradiation of a metal film or nanostructure with an ultrashort laser pulse pushes the electron gas out of equilibrium ; describing the evolution of this non - equilibrium distribution has been the subject of intense research for two decades . a majority of investigations so far employ various approximate models , typically based on free - electron models and empirical electron - phonon interactions , to calculate the energy absorption , electron - electron thermalization and electron - phonon relaxation . @xcite however , a complete _ ab initio _ description of the time evolution and optical response of this non - equilibrium electron gas from femtosecond to picosecond time scales has remained elusive , especially because of the empirical treatment of electron - phonon interactions.@xcite    the initial electron thermalization via electron - electron scattering is qualitatively described within the landau theory of fermi liquids.@xcite the subsequent relaxation of the high temperature electron gas with the lattice is widely described by the two - temperature model ( ttm),@xcite given by coupled differential equations for the electron and lattice temperatures , @xmath0 and @xmath1 ,    c_e(t_e ) & = ( _ e t_e ) - g(t_e ) ( t_e - t_l ) + s(t ) + c_l(t_l ) & = ( _ p t_l ) + g(t_e ) ( t_e - t_l ) . [ eqn : ttm ]    here , @xmath2 and @xmath3 are the thermal conductivities of the electrons and phonons , @xmath4 is the electron - phonon coupling factor , @xmath5 and @xmath6 are the electronic and lattice heat capacities , and @xmath7 is the source term which describes energy deposition by a laser pulse . in nanostructures , the temperatures become homogeneous in space rapidly and the contributions of the thermal conductivities drop out . a vast majority of studies , both theoretical and experimental , treat the remaining material parameters , @xmath4 , @xmath5 and @xmath6 , as phenomenological temperature - independent constants.@xcite    figure  [ fig : schematic ] schematically shows the time evolution of the electron and lattice temperatures in a plasmonic metal like gold , and the role of the _ temperature - dependent _ material properties . the electronic density - of - states and the resultant electronic heat capacity @xmath5 determine the peak electron temperature @xmath0 reached after electron - electron thermalization . the electron - phonon matrix elements and the resulting coupling strength @xmath4 determine the rate of energy transfer from the electrons to the lattice , which along with @xmath5 determines the rate of relaxation of @xmath0 . finally , the phonon density of states and the resulting lattice heat capacity @xmath6 determine the rise in lattice temperature @xmath1 . a key challenge in the quantitative application of ttm models is the determination of these _ temperature - dependent _ material parameters . with pulsed lasers , it is possible to absorb sufficient energy in plasmonic nanostructures to melt the metal once the electrons and lattice have equilibrated.@xcite the highest electron temperature , @xmath8 accessible in repeatable measurements is therefore limited only by the equilibrated lattice temperature being less than the melting temperature @xmath9 of the metal,@xcite which yields the condition @xmath10 . starting at room temperature @xmath11  k and using our calculations of the electron and lattice heat capacities , @xmath5 and @xmath6 , we find @xmath12 5700 , 8300 , 7500 and 6700  k respectively for aluminum , silver , gold and copper . for gold and copper in particular , these temperatures are sufficient to change the occupations of the @xmath13-bands @xmath14  ev below the fermi level . consequently , it is important to derive the temperature dependence of these material parameters from electronic structure calculations rather than free - electron like models.@xcite    therefore to accurately predict the transient optical response of metal nanostructures , we account for the electron - temperature dependence of the electronic heat capacity , electron - phonon coupling factor and dielectric functions . these properties , in turn , require accurate electron and phonon band structures as well as electron - phonon and optical matrix elements . we recently showed that _ ab initio _ calculations can quantitatively predict optical response , carrier generation and electron transport in plasmonic metals in comparison with experiment , with no empirical parameters.@xcite in this article , we calculate @xmath5 , @xmath4 and the temperature and frequency - dependent dielectric function , @xmath15 from first principles . these calculations implicitly include electronic - structure effects in the density of states and electron - phonon interaction matrix elements , and implicitly account for processes such as umklapp scattering . we show substantial differences between our predictions and those from simplified models due to the energy dependence of the electron - phonon matrix elements , especially at high electron temperatures . the paper is organized as follows . we start with the theoretical background and computational methods used in the calculations of the electron heat capacity , phonon coupling and temperature dependent dielectric function of plasmonic materials ( section  [ sec : compdetails ] ) . in section  [ sec : dos ] , we show calculations of the electron heat capacity and its dependence on the electron temperature due to the electronic density of states . analogously , section  [ sec : phonondos ] presents the lattice - temperature dependence of the lattice heat capacity due to the phonon density of states . next , in section  [ sec : geph ] we show a key result of the paper : temperature dependence of the electron - phonon coupling strength accounting for energy dependence of the electron - phonon matrix elements . finally , section  [ sec : dielectric ] presents the temperature and frequency dependence of the dielectric function , including direct ( interband ) , phonon - assisted and drude intraband contributions . section  [ sec : conclusions ] summarizes our results and discusses their application to plasmonic nanostructures in various experimental regimes . we perform density - functional theory ( dft ) calculations of the electronic states , phonons , electron - phonon and optical matrix elements , and several derived quantities based on these properties , for four plasmonic metals , aluminum , copper , silver and gold . we use the open - source plane - wave density - functional software named ` jdftx'@xcite to perform fully relativistic ( spinorial ) band structure calculations using norm - conserving pseudopotentials at a kinetic energy cutoff of 30  hartrees , and the ` pbesol ' exchange - correlation functional ( perdew - burke - ernzerhof functional reparametrized for solids)@xcite with a localized ` + @xmath16 ' correction@xcite for the @xmath13-bands in the noble metals . ref .   shows that this method produces accurate electronic band structures in agreement with angle - resolved photoemission ( arpes ) measurements within 0.1  ev . we calculate phonon energies and electron - phonon matrix elements using perturbations on a @xmath17 supercell . in our calculations , these matrix elements implicitly include umklapp - like processes . we then convert the electron and phonon hamiltonians to a maximally - localized wannier function basis,@xcite with @xmath18 @xmath19-points in the brillouin zone for electrons . specifically , we employ 24 wannier centers for aluminum and 46 spinorial centers for the noble metals which reproduces the density functional theory ( dft ) band structure exactly to at least 50  ev above the fermi level .    using this wannier representation , we interpolate the electron , phonon and electron - phonon interaction hamiltonians to arbitrary wave - vectors and perform dense monte carlo sampling for accurately evaluating the brillouin zone integrals for each derived property below . this dense brillouin zone sampling is necessary because of the large disparity in the energy scales of electrons and phonons , and directly calculating dft phonon properties on dense @xmath19-point grids is computationally expensive and impractical . see ref .   for further details on the calculation protocol and benchmarks of the accuracy of the electron - phonon coupling ( eg . resistivity within 5% for all four metals ) . the electronic density of states ( dos ) per unit volume @xmath20 where @xmath21 are energies of quasiparticles with band index @xmath22 and wave - vector @xmath23 in the brillouin zone bz , directly determines the electronic heat capacity and is an important factor in the electron - phonon coupling and dielectric response of hot electrons . above , the band index @xmath22 implicitly counts spinorial orbitals in our relativistic calculations , and hence we omit the explicit spin degeneracy factor . figure [ fig : dos ] compares the dos predicted by our relativistic pbesol+@xmath16 method with a previous non - relativistic semi - local estimate@xcite using the pbe ( perdew - burke - ernzerhof ) functional,@xcite as well as a free electron model @xmath24 for which @xmath25 . the free electron model is a reasonable approximation for aluminum and the pbe and pbesol+@xmath16 density - functional calculations also agree reasonably well in this case . ( @xmath26 for aluminum . ) the regular @xmath27 @xmath19-point grid used for brillouin zone sampling introduces the sharp artifacts in the dos from ref .  , compared to the much denser monte carlo sampling in our calculations with 640,000 @xmath19-points for au , ag , and cu , and 1,280,000 @xmath19-points for al .    for the noble metals , the free electron model and the density functional methods agree reasonably near the fermi level , but differ significantly @xmath14  ev below the fermi level where @xmath13-bands contribute . the free electron models ignore the @xmath13-bands entirely , whereas the semi - local pbe calculations predict @xmath13-bands that are narrower and closer to the fermi level than the pbesol+@xmath16 predictions . the @xmath16 correction@xcite accounts for self - interaction errors in semi - local dft and positions the @xmath13-bands in agreement with arpes measurements ( to within @xmath28  ev).@xcite additionally , the dos in the non - relativistic pbe calculations strongly peaks at the top of the @xmath13-bands ( closest to the fermi level ) , whereas the dos in our relativistic calculations is comparatively balanced between the top and middle of the @xmath13-bands due to strong spin - orbit splitting , particularly for gold . below , we find that these inaccuracies in the dos due to electronic structure methods previously employed for studying hot electrons propagates to the predicted electronic heat capacity and electron - phonon coupling . the electronic heat capacity , defined as the derivative of the electronic energy per unit volume with respect to the electronic temperature ( @xmath0 ) , can be related to the dos as @xmath29 where @xmath30 is the fermi distribution function . the term @xmath31 is sharply peaked at the fermi energy @xmath32 with a width @xmath33 , and therefore the heat capacity depends only on electronic states within a few @xmath34 of the fermi level . for the free electron model , taylor expanding @xmath35 around @xmath32 and analytically integrating ( [ eqn : ce ] ) yields the sommerfeld model @xmath36 , which is valid for @xmath37 ( @xmath38  k ) . above , @xmath39 , @xmath40 and @xmath41 are respectively the number density , fermi energy and fermi wave - vector of the free electron model .    at temperatures @xmath37 , the electronic heat capacities are much smaller than the lattice heat capacities,@xcite which makes it possible for laser pulses to increase @xmath0 by @xmath42  kelvin , while @xmath1 remains relatively constant.@xcite figure [ fig : ce ] compares @xmath5 from the free - electron sommerfeld model with predictions of ( [ eqn : ce ] ) using dos from pbe and pbesol+@xmath16 calculations . the free - electron sommerfeld model is accurate at low temperatures ( up to @xmath43  k ) for all four metals .    with increasing @xmath0 , @xmath31 in ( [ eqn : ce ] ) is non - zero increasingly further away from the fermi energy , so that deviations from the free electron dos eventually become important . for aluminum , the dos remains free - electron - like over a wide energy range and the sommerfeld model remains valid throughout . for the noble metals , the increase in dos due to @xmath13-bands causes a dramatic increase in @xmath5 once @xmath0 is high enough that @xmath31 becomes non - zero in that energy range . copper and gold have shallower @xmath13-bands and deviate at lower temperatures compared to silver . additionally , the @xmath13-bands are too close to the fermi level in the semilocal pbe calculations of ref .   which results in an overestimation of @xmath5 compared to our predictions based on the more accurate relativistic pbesol+@xmath16 method . similarly , the phonon dos per unit volume @xmath44 where @xmath45 are energies of phonons with polarization index @xmath46 and wave - vector @xmath47 , directly determines the lattice heat capacity , @xmath48 where @xmath49 is the bose occupation factor .    within the debye model , the phonon energies are approximated by an isotropic linear dispersion relation @xmath50 up to a maximum debye wave vector @xmath51 chosen to conserve the number of phonon modes per unit volume . this model yields the analytical phonon dos , @xmath52 , where @xmath53 are the speeds of sound for the one longitudinal and two degenerate transverse phonon modes of the face - centered cubic metals considered here.@xcite    figure  [ fig : phonondos ] compares the dft - calculated phonon dos with the debye model predictions , and shows that the debye model is a good approximation for the dos only up to 0.01  ev . however , figure  [ fig : cl ] shows that the corresponding predictions for the lattice heat capacities are very similar , rapidly approaching the equipartition theorem prediction of @xmath54 at high temperatures , which is insensitive to details in the phonon dos . in fact , the largest deviations of the debye model are below 100  k and less than 10  % from the direct calculations for all four metals . we therefore find that a simple model of the phonons is adequate for predicting the lattice heat capacity , in contrast to the remaining quantities we consider below which are highly sensitive to details of the phonons and their coupling to the electrons . in section  [ sec : dos ] we have shown that the electronic heat capacity , which determines the initial temperature that the hot electrons equilibrate to , is sensitive to electronic structure especially in noble metals at high @xmath0 where @xmath13-bands contribute . now we analyze the electron - phonon coupling which determines the subsequent thermalization of the hot electrons with the lattice . we show that details in the electron - phonon matrix elements calculated using dft also play a significant role , in addition to the electronic band structure , and compare previous semi - empirical estimates of the @xmath0-dependent phonon coupling to our direct calculations . the rate of energy transfer from electrons at temperature @xmath0 to the lattice ( phonons ) at temperature @xmath1 per unit volume is given by fermi s golden rule as    & g(t_e ) ( t_e - t_l ) + & = @xmath55 _ nn ( b - _ - , ) + & _ - , |g^-,_n,n |^2 s_t_e , t_l(_n , _ n , _ - , )    with @xmath56 here , @xmath57 is the unit cell volume , @xmath45 is the energy of a phonon with wave - vector @xmath58 and polarization index @xmath46 , and @xmath59 is the electron - phonon matrix element coupling this phonon to electronic states indexed by @xmath60 and @xmath61 . above , @xmath62 is the difference between the product of occupation factors for the forward and reverse directions of the electron - phonon scattering process @xmath63 , with @xmath64 and @xmath65 being the fermi and bose distribution function for the electrons and phonons respectively . using the fact that @xmath66 for an energy - conserving process @xmath67 by detailed balance , we can write the electron - phonon coupling coefficient as @xmath68 this general form for _ dft - based _ electronic and phononic states is analogous to previous single - band / free electron theories of the electron - phonon coupling coefficient , see for example the derivation by allen et al.@xcite note that unlike previous empirical models , here the coupling coefficient depends on the lattice temperature @xmath1 as well , but we omit the @xmath1 label in @xmath4 to keep the notation consistent with previous approaches,@xcite and present results below for @xmath69  k ( ambient temperature ) . the direct evaluation of @xmath4 using ( [ eqn : g ] ) requires a six - dimensional integral over electron - phonon matrix elements from dft with very fine @xmath19-point grids that can resolve both electronic and phononic energy scales . this is impractical without the recently - developed wannier interpolation and monte carlo sampling methods for these matrix elements,@xcite and therefore our results are the first parameter - free predictions of @xmath4 , derived entirely from dft . previous theoretical estimates of @xmath4 are semi - empirical , combining dft electronic structure with empirical models for the phonon coupling . for example , wang et al.@xcite assume that the electron - phonon matrix elements averaged over scattering angles is independent of energy and that the phonon energies are smaller than @xmath34 , and then approximate the electron - phonon coupling coefficient as @xmath70 where @xmath71 is the electron - phonon mass enhancement parameter and @xmath72 is the second moment of the phonon spectrum.@xcite lin et al.@xcite treat @xmath73 as an empirical parameter calibrated to experimental @xmath4 at low @xmath0 obtained from thermoreflectance measurements , and extrapolate it to higher @xmath0 using ( [ eqn : glin ] ) . see refs .   and for more details .    for clarity , we motivate here a simpler derivation of an expression of the form of ( [ eqn : glin ] ) from the general form ( [ eqn : g ] ) . first , making the approximation @xmath74 ( which is reasonably valid for @xmath0 above room temperature ) allows us to approximate the difference between the electron occupation factors in the second line of ( [ eqn : g ] ) by @xmath75 ( using energy conservation ) . additionally , for @xmath76 , the third line of ( [ eqn : g ] ) simplifies to @xmath77 . with no other approximations , we can then rearrange ( [ eqn : g ] ) to collect contributions by initial electron energy , @xmath78 with @xmath79 therefore , the primary approximation in previous semi - empirical estimates@xcite is the replacement of @xmath80 by an energy - independent constant @xmath73 , used as an empirical parameter .    fig . [ fig : h ] compares our calculations of this energy - resolved electron - phonon coupling strength , @xmath80 , with previous empirical estimates of @xmath73 , and fig . [ fig : g ] compares the resulting temperature dependence of the electron - phonon coupling , @xmath4 , from ( [ eqn : g ] ) and semi - empirical methods([eqn : glin ] ) . for noble metals , @xmath4 increases sharply beyond @xmath81  k because of the large density of states in the @xmath13-bands . however , @xmath80 is smaller by a factor of @xmath82 in the @xmath13-bands compared to near the fermi level . therefore , assuming @xmath80 to be an empirical constant@xcite results in a significant overestimate of @xmath4 at high @xmath0 , compared to the direct calculations . additionally , the shallowness of the @xmath13-bands in the semi - local pbe band structure used in ref .   lowers the onset temperature of the increase in @xmath4 , and results in further overestimation compared to our predictions . our predictions agree very well with the experimental measurements of @xmath4 available at lower temperatures for noble metals.@xcite in fact , the semi - empirical calculation based on @xmath83 underestimates the room temperature electron - phonon coupling for these metals ; the significant overestimation of @xmath4 seen in fig . [ fig : g ] is in despite this partial cancellation of error . this shows the importance of detailed dft electron - phonon matrix elements in calculating the coupling between hot electrons and the lattice . experimental measurements of the electron - phonon coupling in noble metals are reliable because of the reasonably clear separation between a fast electron - electron thermalization rise followed by a slower electron - phonon decay in the thermoreflectance signal . in aluminum , these time scales significantly overlap resulting in strong non - equilibrium effects and making experimental determination of the equilibrium electron - phonon coupling @xmath4 difficult . consequently , the value of @xmath4 for al is not well agreed upon.@xcite using a simplified single - band free - electron - like model of the electrons , ref . estimates @xmath84  w / m@xmath85k for thermalized electrons at 2000  k , which is 1.5x larger than @xmath86  w / m@xmath85k for non - thermalized electrons with the same amount of energy . in figure [ fig : g](a ) , our predictions using ( [ eqn : g ] ) which assumes equilibrium are 2x larger than the experimental estimates@xcite which implicitly include the non - equilibrium effects . on the other hand , the semi - empirical model of ref . assumes thermalized electrons , but fits to experimental data that includes non - thermal effects ( and matches experiment by construction ) . the single - band - model non - equilibrium predictions do not match experiment because it assumes a simple model for electron - phonon matrix elements that ignores umklapp processes.@xcite ultimately , quantitative agreement with experiments for aluminum ( for the right reasons ) therefore requires an extension of our non - empirical dft approach ( [ eqn : g ] ) to include non - equilibrium effects , a subject of current work in our group . the final ingredient for a complete theoretical description of ultrafast transient absorption measurements is the temperature - dependent dielectric function of the material . we previously showed@xcite that we could predict the imaginary part of the dielectric function @xmath87 of plasmonic metals in quantitative agreement with ellipsometric measurements for a wide range of frequencies by accounting for the three dominant contributions , @xmath88 we briefly summarize the calculation of these contributions and focus on their electron temperature dependence below ; see ref .   for a detailed description .    the first term of ( [ eqn : imepstot ] ) accounts for the drude response of the metal due to free carriers near the fermi level , with the zero - frequency conductivity @xmath89 and momentum relaxation time @xmath90 calculated using the linearized boltzmann equation with collision integrals based on dft.@xcite the second and third terms of ( [ eqn : imepstot ] ) ,    @xmath91    @xmath92    capture the contributions due to direct interband excitations and phonon - assisted intraband excitations respectively . here @xmath93 are matrix elements of the momentum operator , @xmath94 is the electric field direction ( results are isotropic for crystals with cubic symmetry ) , and all remaining electron and phonon properties are exactly as described previously . the energy - conserving @xmath95-functions are replaced by a lorentzian of width equal to the sum of initial and final electron linewidths , because of the finite lifetime of the quasiparticles . the dielectric function calculated using ( [ eqn : imepstot]-[eqn : imepsphonon ] ) depends on the electron temperature @xmath0 in two ways . first , the electron occupations @xmath96 directly depend on @xmath0 . second , the phase - space for electron - electron scattering increases with electron temperature , which increases the lorentzian broadening in the energy conserving @xmath95-functions in ( [ eqn : imepsdirect ] ) and ( [ eqn : imepsphonon ] ) . we calculate electron linewidths from dft using fermi golden rule calculations for electron - electron and electron - phonon scattering at room temperature , as detailed in ref .  . these calculations are computationally expensive and difficult to repeat for several electron temperatures ; we instead use the linewidths at room temperature with an analytical correction for the @xmath0 dependence . the electron - phonon scattering rate depends on the lattice temperature , but is approximately independent of @xmath0 because the phase space for scattering is determined primarily by the electronic density - of - states and electron - phonon matrix elements , which depend strongly on the electron energies but not on the occupation factors or @xmath0 . the phase space for electron - electron scattering , on the other hand , depends on the occupation factors and @xmath0 because an electron at an energy far from the fermi level can scatter with electrons close to the fermi level . the variation of this phase - space with temperature is primarily due to the change in occupation of states near the fermi level , and we can therefore estimate this effect in plasmonic metals using a free electron model .    within a free electron model , the phase - space for electron - electron scattering grows quadratically with energy relative to the fermi level , resulting in scattering rates @xmath97 at zero electron temperatures , as is well - known.@xcite we can extend these derivations to finite electron temperature to show that the energy and temperature - dependent electron - electron scattering rate @xmath98 \label{eqn : tautedep}\ ] ] for @xmath99 and @xmath100 . within the free electron model , the constant of proportionality @xmath101 , where the background dielectric constant @xmath102 and the thomas - fermi screening energy scale @xmath103 are typically treated as empirical parameters.@xcite here , we extract @xmath104 by fitting ( [ eqn : tautedep ] ) to the electron - electron scattering rates at room temperature @xmath105 calculated using dft.@xcite the resulting fit parameters are listed in table  [ tab : de ] . we then estimate the total scattering rates at other temperatures by adding @xmath106 to the total dft - calculated results ( including electron - phonon scattering ) at @xmath105 . note that we could have equivalently fit the dft - calculated scattering rates at zero temperature , but the fermi golden rule results at room temperature are less noisy at finite @xmath19-point sampling , and moreover these rates do not differ appreciably for electron energies more than @xmath107  ev away from the fermi level anyway .    . coefficient of the temperature dependence of the electron - electron scattering rate as given by ( [ eqn : tautedep ] ) , extracted from fits to the energy dependence of dft - calculated electron - electron scattering rates at room temperature.@xcite [ tab : de ] [ cols="<,^,^,^,^",options="header " , ]     finally , we use the kramers - kronig relations to calculate @xmath108 from @xmath109 . figure  [ fig : eps300k ] compares the dft - predicted dielectric functions with ellipsometry measurements@xcite for a range of frequencies spanning from near - infrared to ultraviolet . note that we scale the @xmath110-axis by @xmath111 , where @xmath112 is the free - electron plasma frequency , in order to display features at all frequencies on the same scale . we find excellent agreement for aluminum within 10  % of experiment over the entire frequency range , including the peak around 1.6  ev due to an interband transition . the agreement is reasonable for noble metals with a typical error within 20  % , but with a larger error @xmath113  % for certain features in the interband @xmath114 transitions due to inaccuracies in the @xmath13-band positions predicted by dft ( especially for silver ) . in the present work , the pbesol+@xmath16 band structure is typically accurate to @xmath28  ev,@xcite compared to errors @xmath115  ev in @xmath13-band positions predicted by semi - local dft functionals@xcite and qualitative inadequacies of free - electron - like models that ignore @xmath13 bands entirely . consequently , our chosen method has the potential to provide the most reliable predictions of metal dielectric functions , especially for the electron temperature dependence that we discuss next . ( empirical fits such as drude - lorentz models can be more accurate by construction at one temperature,@xcite but do not predict temperature dependence . )    figures  [ fig : epschange400k ] , [ fig : epschange1000k ] and [ fig : epschange5000k ] show the change of the dft - calculated complex dielectric function ( solid lines ) upon increasing the electron temperature @xmath0 from room temperature to 400  k , 1000  k and 5000  k respectively , while the lattice remains at room temperature.@xcite for all four metals , the response from infrared to ultraviolet frequencies is dominated by ` sharp ' features due to interband transitions that broaden with increasing temperature . in the remainder of this section , we analyze these sharp interband features in greater detail using a simpler analytic model of the ( @xmath116 ) transitions ( shown in dashed lines in the aforementioned figures ) .    the strongest temperature dependence in noble metals results from transitions between the highest occupied @xmath13-band to the fermi level near the @xmath117 point , as shown in figure  [ fig : transitionmodel](a ) . assuming a parabolic dispersion and a constant transition matrix element , this temperature dependence can be modeled as@xcite @xmath118 the denominator captures the joint density of states for transitions between the bands , and the numerator counts unoccupied states near the fermi level , which introduces the temperature dependence . above , @xmath119 fills in the real part of the dielectric function , given the imaginary part using the kramers - kronig relation . lccc & ag & au & cu +   + @xmath120  [ ev/@xmath121 & 8.98 & 9.01 & 10.8 + @xmath122  [ ev/@xmath121 & 0.0175 & 0.0240 & 0.0268 +   + @xmath123 [ ev@xmath124 & 70 & 22 & 90 + @xmath125  [ ev ] & 0.31 & 0.96 & 0.98 + @xmath126  [ ev ] & 3.36 & 1.25 & 1.05 + @xmath127 & 5.4 & 3.4 & 16.1 +    table  [ tab : epsparameters ] lists the parameters for the parabolic band approximation obtained from the pbesol+@xmath16 band structures . figure  [ fig : epschange400k ] shows that this approximation ( dashed lines ) captures the correct shape of @xmath128 for small changes in @xmath0 . however , this model underestimates the @xmath0 dependence for higher @xmath0 because it ignores the quadratic increase in broadening of the electronic states due to increased electron - electron scattering , as figures  [ fig : epschange1000k ] and [ fig : epschange5000k ] show . aluminum exhibits a sharp change in the dielectric function around @xmath129  ev , which results from several transitions to / from the fermi level near the w point as figure  [ fig : transitionmodel](b ) shows . additionally two of the involved bands are not parabolic , making it difficult to construct a simple model like ( [ eqn : epschangemodel ] ) . therefore , simplified models are adequate for qualitative analysis of lower temperature excitation experiments in noble metals,@xcite but dielectric functions from first - principles dft calculations are necessary for a quantitative analysis of higher temperature experiments and a wider range of materials and probe frequencies . our parameter - free dft calculations of electron - phonon coupling , electron and lattice heat capacities , and dielectric functions show qualitative differences from free - electron and previous semi - empirical estimates because of the substantial energy dependence of electron - phonon matrix elements and electronic density of states . these changes are particularly important for gold and copper at transient electron temperatures greater than 2000  k because of the change in occupations of the @xmath13-bands situated @xmath14  ev below the fermi level in these metals . the temperature dependence of the optical response is , in particular , important for a wide range of applications beyond understanding ultrafast measurements . we show that while simple models can account for some of the qualitative features of the change in dielectric function for small changes in temperature , an electronic structure treatment is essential to quantitatively account for the complete frequency and temperature dependence , including effects such as carrier linewidth broadening and transitions between multiple non - parabolic bands . given the dearth of published temperature - dependent dielectric functions in the literature , we include detailed tables of our predictions for electron temperatures up to 8000  k , and spanning frequencies from the infrared to the ultraviolet , in the supplementary information.@xcite    this work has direct implications for analysis of experimental pump - probe studies of metal nanostructures and is the subject on ongoing work in our group . with the predicted material properties we anticipate a parameter - free description of the spectra obtained in transient absorption studies since we implicitly account for all the microscopic processes in the non - equilibrium dynamics of electrons in plasmonic metals . this material is based upon work performed by the joint center for artificial photosynthesis , a doe energy innovation hub , supported through the office of science of the u.s . department of energy under award number de - sc0004993 . this research used resources of the national energy research scientific computing center , a doe office of science user facility supported by the office of science of the u.s . department of energy under contract no . de - ac02 - 05ch11231 . the authors acknowledge support from ng next at northrop grumman corporation . p. n. is supported by a national science foundation graduate research fellowship and by the resnick sustainability institute . b. is supported by a national science foundation graduate research fellowship , a link foundation energy fellowship , and the doe ` light - material interactions in energy conversion ' energy frontier research center ( de - sc0001293 ) . 55ifxundefined [ 1 ] ifx#1 ifnum [ 1 ] # 1firstoftwo secondoftwo ifx [ 1 ] # 1firstoftwo secondoftwo `` `` # 1''''@noop [ 0]secondoftwosanitize@url [ 0 ]  + 12$12  & 12#1212_12%12@startlink[1]@endlink[0]@bib@innerbibempty @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) http://link.aps.org/doi/10.1103/physrevb.43.4488 [ * * , ( ) ] link:\doibase    10.1063/1.4906553 [ ( ) ] link:\doibase 10.1021/cr1002547 [ * * , ( ) ] @noop * * ,   ( ) @noop _ _ ,  vol . ( ,  ) link:\doibase    10.1364/ol.40.000340 [ * * ,   ( ) ] link:\doibase 10.1021/jp9917648 [ * * ,   ( ) ] link:\doibase 10.1039/c2cp43771f [ * * , ( ) ] link:\doibase    10.1063/1.1517187 [ * * ,   ( ) ] \doibase    http://dx.doi.org/10.1016/j.apsusc.2013.12.074 [ * * ,   ( ) ] @noop * * ,   ( ) @noop * * ,   ( ) link:\doibase 10.1103/physrevb.65.214303 [ * * ,   ( ) ] link:\doibase    10.1103/physrevb.50.15337 [ * * ,   ( ) ] \doibase http://dx.doi.org/10.1063/1.4904355 [ * * ,   ( ) ] http://dx.doi.org/10.1063/1.4870040 [ * * ,   ( ) ] link:\doibase 10.1103/physrevb.77.075133 [ ( ) ] @noop _ _ ,  vol .  ( ,  ) @noop _ _ ( ,  ) @noop _ _ ( ,  ) @noop _ _ ,  vol .  ( ,  ) link:\doibase 10.1103/physrevb.57.11334 [ * * ,   ( ) ] @noop * * ,   ( ) http://dx.doi.org/10.1103/physrevb.79.195405 [ * * , ( ) ] @noop * * ,   ( ) @noop * * ,   ( ) http://dx.doi.org/10.1021/cr050161r [ * * , ( ) ] http://dx.doi.org/10.1038/nnano.2015.165 [ * * ,   ( ) ] @noop * * ,   ( ) link:\doibase 10.1021/jp983141k [ * * ,   ( ) ] @noop _ _ ,  crc handbook of chemistry and physics ( ,  ) @noop * * ,   ( ) @noop `` , ''   ( ) link:\doibase    10.1103/physrevlett.100.136406 [ * * ,   ( ) ] link:\doibase    10.1103/physrevb.57.1505 [ * * ,   ( ) ] @noop * * ,   ( ) link:\doibase 10.1103/physrevb.65.035109 [ * * ,   ( ) ] @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) link:\doibase 10.1103/physrevb.76.165108 [ * * ,   ( ) ] @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) link:\doibase 10.1103/physrevb.87.035139 [ * * ,   ( ) ] link:\doibase 10.1103/physrevb.70.235125 [ * * ,   ( ) ] @noop _ _  ( ,  ,  ) link:\doibase 10.1364/ao.37.005271 [ * * , ( ) ] \doibase http://dx.doi.org/10.1016/0039-6028(73)90359-2 [ * * ,   ( ) ]
Abstract:

Answer: ultrafast laser measurements probe the non - equilibrium dynamics of excited electrons in metals with increasing temporal resolution . electronic structure calculations can provide a detailed microscopic understanding of hot electron dynamics , but a parameter - free description of pump - probe measurements has not yet been possible , despite intensive research , because of the phenomenological treatment of electron - phonon interactions . we present _ ab initio _ predictions of the electron - temperature dependent heat capacities and electron - phonon coupling coefficients of plasmonic metals . we find substantial differences from free - electron and semi - empirical estimates , especially in noble metals above transient electron temperatures of 2000  k , because of the previously - neglected strong dependence of electron - phonon matrix elements on electron energy . we also present first - principles calculations of the electron - temperature dependent dielectric response of hot electrons in plasmonic metals , including direct interband and phonon - assisted intraband transitions , facilitating complete theoretical predictions of the time - resolved optical probe signatures in ultrafast laser experiments .


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: reliable prediction of durability of structural components and assemblies is a fundamental requirement in various branches of engineering : transport , power generation , manufacturing and many others . this requirement led to the development of various analytical approaches to structural integrity , including elasto - plastic fracture mechanics , low cycle and high cycle fatigue analysis , creep and damage analysis , and to the introduction and standardization of appropriate methods of material property characterization . once the material properties are determined from series of controlled laboratory experiments , numerical models of complex assemblies are used to predict the location and time of failure . underlying this methodology is the assumption that material properties and damage accumulation models validated in the laboratory setting can be successfully transferred to the prototype ( full scale object ) , provided certain requirements of scale independence are fulfilled . one of the difficulties arising in the way of applying this methodology is the fact that in the process of assembly materials undergo additional processing operations that modify their internal structure ( e.g. grain size and texture , composition ) and internal load distribution ( residual stress ) . residual stress appears to be a particularly difficult parameter to account for . unlike e.g. microstructure and composition , residual stress is associated with the entire assembly , and often undergoes significant changes if a testing piece is extracted from it for investigation . yet residual stress is often the most crucial parameter that determines the durability of an assembled structure . welding is a joining and fabrication technique that relies on local melting or softening of the parent material with the purpose of allowing it to fuse together with the filler material and the opposing piece to which a joint is being made . in the process of welding the material undergoes a complex thermal and mechanical history , and the resulting assembly inherits a memory of this process in the form of microstructural and compositional changes , and residual stress distribution . in the recent decades significant efforts have been devoted by many researchers to the development of detailed numerical models of the welding process ; yet at the present time reliable prediction of material structure and residual stress state at the outcome of a welding process remains elusive , its use being primarily limited to qualitative identification of , for example , the most favourable conditions for producing a weld with lower residual stress and distortion . however , it remains necessary to predict durability of assemblies in the presence of welding - induced residual stresses , since this joining method remains in widespread use e.g. in the aerospace industry . in this situation experimental methods of residual stress determination come to the fore , since they provide information about the central link in the chain _ processing _ - _ residual stress _ - _ structural integrity_. methods of residual stress evaluation can be notionally split into _ mechanical _ and _ physical_. mechanical methods of residual stress determination rely on detecting the deformation ( strain ) or distortion in the test piece that arises following some cutting or material removal . perhaps the most well - known of such techniques is hole drilling , that involves plunging a sharp fast drill into the surface of material , and detecting strain change from the original state at the surface of the material using a specially designed strain rosette . the results of a hole drilling experiment are interpreted using an elastic numerical model . another method developed recently is known as the contour method @xcite in which the test piece is carefully sliced using spark erosion and a two - dimensional map of displacement in the direction normal to the cutting plane is collected . this map is then used as the input for an elastic finite element numerical model of the piece , and allows the residual stress to be determined in the plane of the section . physical methods of residual stress evaluation rely on the determination of some parameter of the system that is known to correlate with the residual stress . perhaps the most well - known of the physical methods of residual stress determination is x - ray diffraction . in a diffraction experiment a beam of particles ( e.g. photons or neutrons ) is directed at a particular location within a polycrystalline sample , and an intensity profile of the scattered particles is collected , either as a function of scattering angle for fixed energy beam ( monochromatic mode ) , or as a function of photon energy for fixed scattering angle ( white beam mode ) . in both cases the pattern displays distinct peaks that correspond to distances between crystal lattice planes that are prevalent within the sample . if strain - free distances are known for the same sets of planes , then the measurements allow the calculation of residual direct elastic strains referring to specific orientations within the crystal and the sample . the most widespread laboratory implementation of the x - ray diffraction method for the determination of residual stress is known as the @xmath0 technique . in this technique a series of measurements is carried out to collect the data for elastic lattice strains for a set of directions that deviate from the sample normal by a varying angle @xmath1 . an assumption is then made that these measured strains correspond to a consistent strain state within homogeneous isotropic linear elastic solid , that allows the stress state within the sample to be deduced . an important observation that needs to be made at this point concerns the fact that residual stress is not _ measured _ in such an experiment , but merely deduced on the basis of certain assumptions , including that ( i ) that the material is uniform , isotropic and continuous , ( ii ) that strain values measured at different angles of tilt , @xmath1 , are obtained from the same group of grains within the same gauge volume ; ( iii ) that the component of stress normal to the sample surface vanishes within the gauge volume ; etc . the above assumptions are in fact approximations whose validity may or may not be readily proven , or which are , in the worst case , simply wrong . the diffraction of neutrons and high energy synchrotron x - rays provides a unique non - destructive probe for obtaining information on strains deep in the bulk of engineering components and structures , e.g. @xcite . this method has become a mature tool for the determination of residual strain states in small coupons , and developments are under way to establish the facilities for performing high resolution measurements directly on larger engineering components @xcite . a particular feature of high energy x - ray diffraction is that the radiation is primarily scattered forward , i.e. in directions close to that of the incident beam @xcite . therefore small diffraction angles have to be used , usually @xmath2 . two difficulties follow . firstly , it is difficult to measure strains in directions close to that of the incident beam . this is due to the fact that the scattering vector is always the bisector of the incident and diffracted beams . hence for high energy x - rays the strain measurement directions form a shallow cone . for a scattering angle of @xmath3 this cone has the angle of @xmath4 . in practice this means that strain directions accessible for the high energy x - ray diffraction techniques are close to being normal to the incident beam . this situation is in stark contrast with that encountered in laboratory x - ray diffraction where near backscattering geometry is used , and measured strains are in directions close to being parallel with the incident beam . secondly , it is difficult to achieve good spatial resolution in the direction of the incident beam , due to the elongated shape of the scattering volume . although rotating the sample may help to overcome these difficulties , in practice this is often limited by the sample shape and absorption , and means that often only two components of strain ( in the sample system ) are known with sufficient accuracy and resolution . neutron diffraction strain analysis has the advantage that it is more often possible to measure the lattice parameter in three mutually perpendicular directions . it is therefore sometimes claimed that this is the only method capable of true residual stress measurement . however , it must be noted even then that the residual stress evaluation involves making certain assumptions : that indeed three principal directions were chosen ; that the strain - free lattice parameters have been correctly determined for all three directions ; that the correct values of young s modulus and poisson s ratio were used in the calculation . in other words , stress evaluation relies on calculations based on certain assumptions . furthermore , the conventional interpretation procedures remain point - wise , i.e. make no attempt to establish correlation between measurements at different locations , and to check whether the results conform to the basic requirements of stress and moment balance within each arbitrarily chosen sub - volume , and that strain compatibility and traction - free surface conditions are satisfied . the purpose of the foregoing discussion was to establish the basic fact that residual stress state is never measured directly , be it by mechanical or physical methods , but always deduced by manipulating the results of the measurements in conjunction with certain models of deformation . it is therefore correct to say that residual stress measurement is one area of experimental activity where the development and validation of numerical models needed for the interpretation of data occupies a particularly important place : without adopting some particular model of deformation is it impossible to present measurement results in terms of residual stress .    to give a very general example , when a ruler is pressed against the sample to determine its length , the implication is that the sample and ruler surfaces are collocated all along the measured length ; and that the length of the ruler between every pair of markers is exactly the same . only if that is so then the reading from the ruler is correct . the approach adopted in this study rests on the explicit postulate that it is necessary to make informed assumptions about the nature of the object ( or distribution ) that is being determined in order to regularise the problem . interpretation of any and every measurement result relies on a model of the object being studied .    in the present paper we are concerned with a model of residual stress generation to which we refer as the eigenstrain technique . the term eigenstrain and notation @xmath5 for it were introduced by toshio mura @xcite to designate any kind of permanent strain in the material arises due to some inelastic process such as plastic deformation , crystallographic transformation , thermal expansion mismatch between different parts of an assembly , etc . in some simple cases , e.g. in the analysis of residually stressed bent beams , it is possible to derive explicit analytical solutions for residual stresses as a function of an arbitrary eigenstrain distribution @xcite . in the more general context it is apparent that eigenstrains are responsible for the observed residual stresses , and therefore can be thought of as the source of residual stress ; yet eigenstrain is not a priori associated with any stress , nor can it be measured directly , as it does not manifest itself in the change of crystal lattice spacing or any other physical parameter . in fact , if a uniform eigenstrain distribution is introduced into any arbitrarily shaped object , then no residual stress is produced . in other words , eigenstrains are associated with misfit and mismatch between different parts of an object . this conclusion is particularly interesting on the context of the foregoing discussion about engineering assemblies and the nature of residual stresses in them . the following discussion is based on the analysis of the fundamental equations describing the generation of residual elastic stresses and strains from a distribution of eigenstrains . most often the problem has to be addressed in the inverse sense : residual stresses or residual elastic strains are somehow known in a number of locations , while the unknown underlying distribution of eigenstrain sources of residual stresses needs to be found . while the direct problem of residual stress determination from eigenstrain distribution can be classed as easy ( linear , elastic ) , the inverse problem is more difficult . however , it is important to emphasize that once the eigenstrain distribution is found , it can be used to solve the easy direct problem , so that the residual stress distribution becomes known in full . the procedure of finding the underlying eigenstrain distribution and reconstructing the complete residual stress state is entirely similar in principle to any other method of residual stress determination discussed above : the experimental data are treated using some suitable numerical model , and the residual stress state is deduced . there are some distinct advantages offered by the eigenstrain approach . firstly , the solution obtained in this way is forced to satisfy the requirements of total strain compatibility and stress equilibrium , that often become violated if less sophisticated methods of data interpretation are used . secondly , once the eigenstrain distribution is deduced it allows the prediction of the object s distortion and residual stress re - distribution during any subsequent machining operation , such as sectioning or surface layer removal .    in the following section we present a formulation of the direct problem of residual stress determination from the known eigenstrain distribution . we then formulate an efficient non - iterative variational approach for solving the inverse problem , and describe briefly its possible numerical implementations . we then apply the method to the analysis of experimental data for residual elastic strains in a single pass electron beam weld in a plate of aerospace nickel superalloy in718 , obtained using high energy synchrotron x - ray diffraction . we demonstrate how the eigenstrain distribution can be found that minimizes the disagreement between the measurements and the model prediction , and also how the method allows the refinement of the strain - free lattice parameter across the weld . we show reconstructions of the complete residual stress field within the plate . finally , we carry out sensitivity analysis to determine whether the solution obtained in terms of the eigenstrain distribution ( and hence in terms of the reconstructed residual stress state ) is stable with respect to the amount of experimental residual elastic strain data available . , and transverse direction , @xmath6 , were carried out along each line , allowing macroscopic residual stresses to be calculated . ] figure [ fig : one ] illustrates the dimensions of the experimental specimen used in the present study , and the location of the measurement points . electron beam welding was used to manufacture a flat rectangular plate by joining two elongated strips ( 3 mm thick , 200 mm long and approximately 25 and 35 mm wide ) .    the sample used for the experiment was made from in718 creep resistant nickel superalloy used in the manufacture of aeroengine components , such as combustion casings and liners , as well as disks and blades . the composition of the alloy in weight percent is approximately given by 53% ni , 19% fe , 18% cr , 5% nb , and small amounts of additional alloying elements ti , mo , co , and al . apart from the matrix phase , referred to as @xmath7 , the microstructure of the alloy may show additional precipitate phases , referred to as @xmath8 , @xmath9 , and @xmath10 . the primary strengthening phase , @xmath9 , has the composition @xmath11 and a body - centred tetragonal structure , and forms semi - coherently as disc - shaped platelets within the @xmath7 matrix . it is highly stable at @xmath12c , but above this temperature it decomposes to form the @xmath8 @xmath13 phase ( between @xmath14c and @xmath15c ) , and @xmath10 , having the same composition as @xmath9 ( between @xmath16c and @xmath17c ) . at large volume fractions and when formed continuously along grain boundaries , the @xmath10 is thought to be detrimental to both strength and toughness @xcite . the @xmath10 phase that forms is more stable than the @xmath9 phase , and has an orthorhombic structure @xcite . welding is known to give rise to residual tensile stresses along the weld line and in the immediately adjacent heat affected zones ( haz ) , while compressive residual stress is found further away from the seam . the most significant residual stress component is the longitudinal stress that we denote by symbols @xmath18 or @xmath19 that have the same meaning throughout . high energy synchrotron x - ray diffraction measurements were carried out on the id11 and id31 beamlines at the esrf using monochromatic mode and a scanning diffractometer . the energy of about 70kev was selected by the monochromator , and the 111 reflection of the @xmath7 matrix phase was used . the beam spot size on the sample was approximately 1 mm ( horizontal ) by 0.25 mm ( vertical ) . line scans were performed with the scan step size of 1 mm along the four lines indicated in figure [ fig : one ] , lying 0 mm , 2 mm , 10 mm and 50 mm above the lower edge of the weld plate . both the horizontal ( transverse ) strain component , @xmath6 , and the vertical ( longitudinal ) strain component @xmath20 were evaluated . assuming that the state of plane stress persists in the thin welded plate , this allowed the stress components @xmath21 and @xmath19 to be calculated . , along each of the four lines , calculated from two mutually orthogonal components of strain under the assumption of plane stress . note that the curve corresponding to the edge of the plate ( 0 mm ) does not show uniform zero stress . ] figure [ fig : two ] illustrates the results of the measurement interpreted in this straightforward way , plotted as the longitudinal residual stress @xmath19 as a function of horizontal position measured from the nominal centre of the weld line . it is apparent from the plots that the stress profiles evolve both in terms of the magnitude and shape away from the edge of the plate . originally the results were interpreted by assuming a constant value of @xmath22 , the unstrained lattice spacing , everywhere within the plate . however , this led to the physically unacceptable result of non - zero longitudinal stress existing at the bottom edge of the plate . the calculation was then repeated imposing the stress - free condition at the bottom edge , and allowing @xmath22 to vary only as the function of the horizontal coordinate @xmath23 along the bottom edge of the weld plate so as to produce stress free condition at that edge . distributions of inelastic strains contained in the sample act as sources of residual stresses . indeed , in the absence of inelastic ( permanent ) strain of some origin ( or indeed , when such inelastic strain is uniform throughout the sample ) , then any specimen is stress free in the absence of external loading . for a known non - uniform eigenstrain distribution @xmath24 the elastic residual strains ( and hence residual stresses ) in the body can be found by the formulae @xcite @xmath25 the above formula is in principle applicable to bodies of arbitrarily complex shape , provided the elastic constants @xmath26 are known , together with the corresponding green s function @xmath27 . in practice green s functions can be found only for bodies of simple geometry , e.g. infinitely extended two - dimensional or three - dimensional solid . the fundamental value of the above formula , however , lies in the statement that for known eigenstrain distribution the elastic response of the body containing it can be readily found .    for convoluted geometries the finite element method provides a suitable method of solving the above direct problem of finding residual elastic strains from given eigenstrains . we are interested here in the problem that often arises in residual stress measurement and interpretation . let there be given a set of measurements ( with certain accuracy ) of strains and stresses collected from a finite number of points ( sampling volumes ) within a bounded specimen . we would like to solve the inverse problem about the determination of unknown eigenstrains from this incomplete knowledge of elastic strains or residual stresses . the limited accuracy and lack of completeness of measurements suggest that direct inversion of ( 1 ) may not be the preferred solution . in fact the method chosen must be sufficiently robust to furnish approximate solutions even in this case . the incorporation of eigenstrain into the finite element method framework can be accomplished via the use of anisotropic pseudo - thermal strains . in the present case we concentrated our attention on the determination of a single eigenstrain component , @xmath28 , longitudinal with respect to the extent of the weld joint . it is clear that in practice this is not the only eigenstrain component likely to be present . however , it is also apparent that this component is the most significant in terms of its effect on the longitudinal stress , @xmath18 . it is worth noting that the procedure outlined below is not in any way restricted to the determination of single component eigenstrain distributions , but is in fact entirely general .     along the line 50 mm above the lower edge of the welded plate . stresses computed from measurements are shown as markers , while the model prediction is shown by the continuous curve . the stress values are plotted as a function of horizontal position @xmath29 with respect to the nominal centre of the weld line . ] the eigenstrain distribution was introduced in the form of a truncated series of basis functions @xmath30 each of the basis functions @xmath31 is chosen in the variable separable form as @xmath32 and index @xmath33 provides numeration for the entire pair set @xmath34 . functions @xmath35 and @xmath36 are chosen to reflect the nature of the eigenstrain dsitribution in the present problem . it was found that for the functions of the horizontal coordinate , @xmath35 , a suitable choice is provided by the chebyshev polynomials @xmath37 , with the argument @xmath38 scaled from the canonical interval @xmath39 $ ] to the plate width . for the functions of the vertical coordinate , @xmath36 , powers @xmath40 , were used .    as stated earlier , the solution of the direct eigenstrain problem can be readily obtained for any eigenstrain distribution by an essentially elastic calculation within the fe model . this task is easily accomplished for the basis functions @xmath41 . furthermore , due to the problem s linearity , the solution of the direct problem described by a linear combination of individual eigenstrain basis functions @xmath41 with coefficients @xmath42 is given by the linear superposition of solutions with the same coefficients . this observation provides a basis for formulating an efficient variational procedure for solving the inverse problem about the determination of underlying eigenstrain distribution .    the problem that we wish to address here stands in an inverse relationship to the direct eigenstrain problem . this is the situation most commonly encountered in practice : the residual elastic strain distribution is known , at least partially , e.g. from diffraction measurement . the details of the preceding deformation process need to be found , such as distribution of eigenstrains within the plastic zone . alternatively , in the absence of non - destructive measurements of residual elastic strain , changes in the elastic strain may have been monitored , e.g. using strain gauges , in the course of material removal . questions arise immediately regarding the invertibility of the problem ; its uniqueness ; the regularity of solution , i.e. whether the solution depends smoothly on the unknown parameters . although do not attempt to answer these questions here , we present a constructive inversion procedure and also provide a practical illustration of its stability . denote by @xmath43 the distribution of the longitudinal stress component @xmath19 arising from the eigenstrain distribution given by the @xmath44th basis function @xmath31 . evaluating @xmath43 at each of the @xmath45th measurement points with coordinates @xmath46 gives rise to predicted values @xmath47 . due to the problem s linearity the linear combination of basis functions expressed by equation ( [ eq : one ] ) gives rise to the stress values at each measurement point given by the linear combination of @xmath48 with the same coefficients @xmath42 , i.e. @xmath49 . denote by @xmath50 the values of the same stress component @xmath51 at point @xmath46 deduced from the experiment . in order to measure the goodness of the prediction we form a functional @xmath52 given by the sum of squares of differences between actual measurements and the predicted values , with weights : @xmath53 where the sum in @xmath54 is taken over all the measurement points . the choice of weights @xmath55 remains at our disposal and can be made e.g. on the basis of the accuracy of measurement at different points . in the sequel we assume for simplicity that @xmath56 , although this assumption is not restrictive . the search for the best choice of model can now be accomplished by minimising @xmath52 with respect to the unknown coefficients , @xmath42 , i.e. by solving @xmath57    due to the positive definiteness of the quadratic form ( [ eq : three ] ) , the system of linear equations ( [ eq : four ] ) always has a unique solution that corresponds to a minimum of @xmath52 . the partial derivative of @xmath52 with respect to the coefficient @xmath42 can be written explicitly as @xmath58    we introduce the following matrix and vector notation @xmath59 the entities appearing in equation ( [ eq : six ] ) can be written in matrix form as : @xmath60 hence equation ( [ eq : five ] ) assumes the form @xmath61 the solution of the inverse problem has thus been reduced to the solution of the linear system @xmath62 for the unknown vector of coefficients @xmath63 .     along the line 10 mm above the lower edge of the welded plate . stresses computed from measurements are shown as markers , while the model prediction is shown by the continuous curve . the stress values are plotted as a function of horizontal position @xmath29 with respect to the nominal centre of the weld line . ] once the coefficients @xmath42 have been determined the eigenstrain distribution in equation ( [ eq : one ] ) is introduced into the finite element model , and the complete stress - strain field is reconstructed by solving the direct problem . by construction the corresponding stress field satisfies the conditions of equilibrium within arbitrary sub - volume of the model , and traction - free boundary conditions are enforced . the total strain field composed of the residual elastic strains and eigenstrains satisfies the conditions of compatibility . the optimal agreement with the experimental measurements is achieved in the least squares sense over all eigenstrain distributions spanned by the functional space of equation ( [ eq : one ] ) . figure [ fig : three ] shows the comparison between the experimental values shown by the markers and the reconstructed stress profile ( continuous curve ) along the line 50 mm above the lower edge of the weld plate . figure [ fig : four ] shows a similar comparison for the line at 10 mm from the edge , and figure [ fig : five ] for the line 2 mm above the lower edge of the plate . note the difference in the scales used for the three figures that explains the greater apparent scatter in the last plot . it is also worth recalling that as a result of adjustment of the @xmath22 values longitudinal stress @xmath19 is equal to zero along the line 0 mm lying at the edge , and hence the plot is not shown . along the line 2 mm above the lower edge of the welded plate . stresses computed from measurements are shown as markers , while the model prediction is shown by the continuous curve . the stress values are plotted as a function of horizontal position @xmath29 with respect to the nominal centre of the weld line . ]     along the four lines lying 0 mm , 2 mm , 10 mm and 50 mm above the lower edge of the plate , respectively . ] figure [ fig : six ] illustrates the nature of the eigenstrain distribution used in the variation model for residual stress reconstruction . the four curves indicate the variation of ( compressive ) eigenstrain as a function of the horizontal coordinate @xmath29 along the four lines lying 0 mm , 2 mm , 10 mm and 50 mm above the lower edge of the plate , respectively . figure [ fig : seven ] shows a two - dimensional contour representation of the underlying eigenstrain field determined using the variational approach , shown for the lower half of the welded plate . recall that symmetry is implied with respect to the upper edge of the plot . figure [ fig : eight ] shows a contour plot of the reconstructed von mises stress field in the lower half of the welded plate . figure [ fig : ten ] shows a contour plot of the reconstructed longitudinal @xmath19 stress field in the lower half of the welded plate .     in the lower half of the welded plate . ]     along the line at @xmath64 from the lower edge of the plate from three models using the data for 50,10,2 mm ; 50,10 mm ; 50,2 mm , respectively . ]     along the line at @xmath65 from the lower edge of the plate from three models using the data for 50,10,2 mm ; 50,10 mm ; 50,2 mm , respectively . ]     along the line at @xmath66 from the lower edge of the plate from three models using the data for 50,10,2 mm ; 50,10 mm ; 50,2 mm , respectively . ]     along the weld line @xmath67 obtained from three models using the data for 50,10,2 mm ; 50,10 mm ; 50,2 mm , respectively . ] obtained from three models using the data for 50,10,2 mm ; 50,10 mm ; 50,2 mm , respectively . ] obtained from three models using the data for 50,10,2 mm ; 50,10 mm ; 50,2 mm , respectively . ] obtained from three models using the data for 50,10,2 mm ; 50,10 mm ; 50,2 mm , respectively . ] ( the lower edge of the weld plate ) obtained from three models using the data for 50,10,2 mm ; 50,10 mm ; 50,2 mm , respectively . ] it may appear at first glance that the proposed reconstruction procedure is akin to some kind of a trick , since the amount of information presented in figures [ fig : eight ] and [ fig : ten ] seems significantly greater than that originally avalable from the measurements in figure [ fig : two ] . in fact , the reconstructions shown in figures [ fig : eight ] and [ fig : ten ] are not just interpolations , and do contribute additional information to the analysis . by the very nature of the reconstruction process the possible fields included in the consideration are only those that satisfy all the requirements of continuum mechanics . this amounts to a very significant additional constraint being placed on data interpretation . provided the analysis of the experimental data is carried out in terms of eigenstrain , all of the predicted fields necessarily conform to these constraints , furnishing additional insight into the residual stress field being studied . the question that we would like to tackle further concerns the sensitivity of the solution , i.e. the determined eigenstrain distribution and the reconstructedelastic fields , to the selection of experimental data included in the analysis . instead of attempting to provide a general analytical answer to this question at this point we perform some tests on the data set available within this study , as follows .    in the results shown in the previous section all of the experimental data available was used in the reconstruction . in the discussion that follows below this model will be referred to as the y=50,y=10,y=2 model , since the data along these lines was used in the reconstruction . we now perform variational eigenstrain analysis on two further models , the y=50,y=10 and y=50,y=2 models , i.e. omitting line y=2 and y=10 respectively . figure [ fig : twelve ] shows that the reconstructed residual stress @xmath18 plot along the line at 50 mm from the lower edge of the weld plate is quite insensitive to the omission of some data . figure [ fig : thirteen ] represents the plot of @xmath19 along the line 10 mm from the lower edge of the weld plate . clearly the greatest deviation from the complete model results is found in the y=50,y=2 model , in which the data along the line y=10 itself was omitted . however , it s encouraging to note that the qualitative nature of the residual stress distribution is reconstructed correctly , although quantitative difference in the magnitude is observed . note , however , that this is in fact the result of prediction of the residual stress from measurements conducted remotely from the point of interest made without recourse to any process model whatsoever . figure [ fig : fourteen ] represents the plot of @xmath19 along the line 2 mm from the lower edge of the weld plate . comparison between predictions made by different models once again demonstrates considerable stability of the prediction with respect ot data omission . figure [ fig : fifteen ] presents the plot of @xmath19 along the line @xmath67 , i.e. the weld line . it is found that the agreement between the models y=50,y=10,y=2 and y=50,y=10 is remarkably good . however , comparison between the y=50,y=10,y=2 and y=50,y=2 models shows that omitting y=10 exerts a significant influence on the predicted residual stress distribution . note , however , that the data along the @xmath68direction is very sparse in the first place , so perhaps this result is not entirely unexpected .    figure [ fig : sixteen ] shows the comparison plot for the eigenstrain variation perpendicular to the weld line at @xmath64 obtained from three models using the data for 50,10,2 mm ; 50,10 mm ; 50,2 mm , respectively . remarkable stability of eigenstrain determination with respect to data omission is observed here . figure [ fig : seventeen ] shows the comparison plot for the eigenstrain variation perpendicular to the weld line at @xmath65 obtained from three models using the data for 50,10,2 mm ; 50,10 mm ; 50,2 mm , respectively . the two models using the data from the 10 mm line show a very close agreement , while the model y=50,y=2 shows some deviation , although even in that case the agreeement remains good . figure [ fig : eighteen ] shows the comparison plot for the eigenstrain variation perpendicular to the weld line at @xmath66 obtained from three models using the data for 50,10,2 mm ; 50,10 mm ; 50,2 mm , respectively . once again , agreement between the models remains good evenwhen the data from the line in question is omitted . finally , figure [ fig : nineteen ] shows the comparison plot for the eigenstrain variation perpendicular to the weld line at @xmath69 ( edge of the plate ) obtained from three models using the data for 50,10,2 mm ; 50,10 mm ; 50,2 mm , respectively , confirming the stability of the eigenstrain determination procedure .    the above analysis does not aim to provide a rigorous proof of the regularity or stability of the proposed inversion procedure . however , it does serve to illustrate that the removal of some data ( or its absence in the first place ) does not appear to lead to any significant artefacts that might raise doubts in the utility of the proposed approach . it is the authors belief that the variation approach to eigenstrain determination and residual stress reconstruction introduced in the present paper has the potential to provide very significant improvement in the quality of experimental data interpretation for the purpose of residual stress assessment . the scope of the newly proposed approach is very wide : it can be used with some success to study the data form hole drilling , slitting , sachs boring and many either destructive and non - destructive technqiues . furthermore , the eigenstrains introduced into the finite element model in the way described here provide an excellent framework for considering subsequent inelastic deformation mechanism accompanying various processing operations and in situ loading , including creep and/or relaxation during heat treatment , residual stress evolution in fatigue , etc . these research directions are being pursuded by the authors . the authors would like to acknowledge the support of uk department of trade and industry and rolls - royce plc under project adam - darp . a.m.korsunsky , s.p.collins , r.a.owen , m.r.daymond , s.achtioui , and k.e.james . fast residual stress mapping using energy dispersive synchrotron x - ray diffraction on station 16.3 at the srs . , 9:7781 , 2002 . j.liu , k.kim , m.golshan , d.laundy , and a.m.korsunsky . energy calibration and full - pattern refinement for strain analysis using energy - dispersive and monochromatic x - ray diffraction . , 38:661667 , 2005 .
Abstract:

Answer: we present the formulation for finding the distribution of eigenstrains , i.e. the sources of residual stress , from a set of measurements of residual elastic strain ( e.g. by diffraction ) , or residual stress , or stress redistribution , or distortion . the variational formulation employed seeks to achieve the best agreement between the model prediction and some measured parameters in the sense of a minimum of a functional given by a sum over the entire set of measurements . the advantage of this approach lies in its flexibility : different sets of measurements and information about different components of the stress - strain state can be incorporated . we demonstrate the power of the technique by analysing experimental data for welds in thin sheet of a nickel superalloy aerospace material . very good agreement can be achieved between the prediction and the measurement results without the necessity of using iterative solution . in practice complete characterisation of residual stress states is often very difficult , due to limitations of facility access , measurement time or specimen dimensions . implications of the new technique for experimental analysis are all the more significant , since it allows the reconstruction of the entire stress state from incomplete sets of data .