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: accurate prediction of finite amplitude acoustic waves traveling through biological media are essential in developing new therapy and imaging techniques for medical ultrasound applications . numerous experimental studies show that attenuation @xmath0 of biological media exhibits a power law dependence on frequency @xmath1 over the frequency range used in medical applications @xcite @xcite . in this sense , considering an initial monochromatic plane wave traveling through nonlinear media , two opposite effects govern the final wave amplitude : on one hand , higher harmonics appear and their amplitude grows as a consequence of nonlinear progressive wave steepening ; and on the other hand , the damping for each harmonic is different , following the above power law . therefore , differences in the frequency power law model can lead to huge differences in attenuation of higher harmonics , and hence the inclusion of frequency dependence attenuation is critical for correctly predict nonlinear propagation . classical thermo - viscous @xcite sound attenuation exhibit a squared frequency dependence , @xmath2 , however the empirically fitted power law for soft tissues typically ranges between @xmath3@xcite @xcite @xcite . for many tissues the attenuation shows a dependence that can be modeled by @xmath4 close to the unity , where the value varies from different tissues . moreover , for some individual examples the local value of @xmath4 exhibit lower values at low ultrasound frequencies and tends to 2 in the high frequency limit @xcite . despite the study of the physical mechanism besides this complex frequency dependence is out of the scope of this work , there exist numerous phenomenological approaches for including the observed losses in the acoustic equations@xcite @xcite . on the other hand , is common in literature to describe the losses in soft tissues as multiple - relaxation processes @xcite @xcite @xcite , where the relaxation frequencies can be associated to either a tissue specific physical mechanism or empirically optimized to fit the observed tissue attenuation @xcite , @xcite . moreover , fractional partial differential operators has been demonstrated the ability to describe frequency power law attenuation @xcite @xcite . these operators can be included in the modeling by means of time @xcite , space @xcite or combined time - space fractional derivatives @xcite @xcite . the latter operators can be biological motivated and derived from tissue micro - structure using fractal ladders based on networks of springs and dashpots @xcite . however , these fractional loss operators can be also derived from a continuum of relaxation process @xcite , that suggest that the two approaches can be equivalent under certain conditions @xcite .    in order to solve these models , many time - domain numerical methods have been developed . attenuation modeled by relaxation processes can be solved by means of finite - differences in time - domain ( fdtd ) solvers in linear regime @xcite and in nonlinear regime applied to augmented burger s equation @xcite , khokhlov - zabolotskaya - kuznetsov ( kzk ) @xcite and westervelt @xcite nonlinear wave equations . on the other hand , time - dependent fractional derivatives can be solved in nonlinear regime @xcite by convolutional operators . this approach requires the memory storage of certain time history , and although the memory can be strongly reduced compared to direct convolutions , this algorithm employs up to ten auxiliary fields and a memory buffer of three time steps . in order to overcome this limitation , time - space fractional derivatives or fractional laplacian @xcite @xcite can be used to model frequency power laws without time - domain convolutional operators . recently k - space and pseudo - spectral methods have been applied in order to solve fractional laplacian operators efficiently in nonlinear regime @xcite . the aim of this work is to present a numerical method that solves the complete set of equations ( continuity , momentum and state equations ) including nonlinear propagation and frequency power law attenuation based in multiple relaxation processes . the numerical method presented here is based on the fdtd method . the inclusion of relaxation processes in the presented formulation avoids convolutional operators so only one extra field per relaxation process is needed and no memory buffer is needed . the paper is organized as follows : in sec . [ s : model ] the model equations that describes the problem are exposed , sec . [ s : numerical ] describes the computational method presented in this work and in sec . [ s : validation ] the method is validated comparing the numerical results with known analytic solutions for linear and nonlinear regimes . the principles of mass and momentum conservation lead to the main constitutive relations for nonlinear acoustic waves , which for a fluid can be expressed as @xcite @xmath5 and @xmath6 where @xmath7 is the total density field,@xmath8 is the particle velocity vector,@xmath9 is the pressure,@xmath10 and @xmath11 are the coefficients of shear and the bulk viscosity respectively . the acoustic waves described by this model exhibit viscous losses with squared power law dependence on frequency . in order to include a power law frequency dependence on the attenuation , a multiple relaxation model will be added into the time domain equations . the basic mechanism for energy loss in a relaxing media is the appearance of a phase shift between the pressure and density fields . this behavior is commonly modeled as a time dependent connection at the fluid state equation , that for a fluid retaining the nonlinear effects up to second order an be expressed as @xcite@xcite : @xmath12 where @xmath13 is the density perturbation over the stationary density @xmath14 , @xmath15 is the nonlinear parameter , @xmath16 is the small amplitude sound speed , and @xmath17 is the kernel associated with the relaxation mechanism . the first two terms describe the instantaneous response of the medium and the convolutional third term accounts for the  memory time " of the relaxing media . thus , by choosing an adequate time function for the kernel @xmath17 the model can present an attenuation and dispersion response that fits the experimental data of the heterogeneous media . however , the direct resolution of the constitutive relations ( [ eq : continuity]-[eq : st_cont ] ) in this integral form is a complex numerical task due to the convolutional operator . thus , instead of describe @xmath17 with a specific time domain waveform , the response of the heterogeneous medium can be alternatively described by a sum of @xmath18 relaxation processes with exponential time dependence as : @xmath19 with the @xmath20-th order kernel expressed as @xmath21 where @xmath22 is the heaviside piecewise function @xmath23 , @xmath24 , @xmath25 is the characteristic relaxation time and @xmath26 the relaxation parameter for the @xmath20-th order process . this last dimensionless parameter controls the amount of attenuation and dispersion for each process as @xmath27 , where @xmath28 is the sound speed in the high frequency limit associated to @xmath20-th order relaxation process , also known as the speed of sound in the  frozen " state @xcite . in order to describe relaxation without the need of including a convolutional operator , we shall define a state variable @xmath29 for each process as @xmath30 thus , using the convolutional property @xmath31 , the time derivative of the relaxation state variable obeys the following relation for the n - th order process : @xmath32 where @xmath33 is the dirac delta function . using the eq . ( [ eq : sn ] ) this relation becomes a simple ordinary differential equation for each process as : @xmath34 using again convolutional properties , we can substitute eq . ( [ eq : dsdt ] ) into ( [ eq : relax_define ] ) , and the relaxing nonlinear state eq . ( [ eq : st_cont ] ) becomes : @xmath35 moreover , if  frozen " sound speed for @xmath18 mechanisms is defined as @xmath36 , eq . ( [ eq : st_final_c0 ] ) leads to : @xmath37 due to the smallness of the relaxation parameter , @xmath26 , i.e. when weak dispersion is modeled , the sound speed in the high frequency limit reduces to @xcite : @xmath38 note eq . ( [ eq : st_final_cinf ] ) for a mono - relaxing media is equivalent to that can be found in literature @xcite @xmath39 thus , the constitutive relations to solve by means of the numerical method in the nonlinear regime are the continuity eq . ( [ eq : continuity ] ) , the motion eq . ( [ eq : momentum ] ) and the second order fluid state relaxing eq . ( [ eq : st_final_cinf ] ) , where the state variable @xmath29 obeys the relation ( [ eq : dsdt ] ) for the @xmath20-th order relaxation process . although the aim of this work is to model biological media , the generalized formulation presented here can be used to describe the attenuation and hence the dispersion observed in other relaxing media , as the relaxation processes of oxygen and nitrogen molecules in air or the relaxation associated with boric acid and magnesium sulfate in seawater @xcite .      on the other hand , if small amplitude perturbations are considered , an equivalent derivation of this model can be expressed for multiple relaxation media @xcite . thus , for an homogeneous inviscid relaxing fluid the linearized continuity and motion eq . ( [ eq : continuity]-[eq : momentum ] ) reduces to @xmath40 and @xmath41 and linearizing the fluid state eq . ( [ eq : st_final_cinf ] ) we obtain : @xmath42 these equations can be solved directly in this form , however , if expressed in pressure - velocity formulation the density field is no longer necessary and computational effort can be reduced . thereby , assuming a linear  instantaneous " compressibility @xmath43 , and substituting eq . ( [ eq : lin_st ] ) into eq . ( [ eq : lin_cont ] ) yields @xmath44 then , taking the time derivative of the state variable eq . ( [ eq : dsdt ] ) we get @xmath45 finally , substituting again the linearized state eq . ( [ eq : lin_st ] ) and arranging terms the linearized continuity equation leads to @xmath46 on the other hand , the state evolution equation can be expressed as a function of the acoustic pressure as @xmath47 thus , the linearized governing eq . ( [ eq : lin_momen],[eq : lin_pconti ] ) for a relaxing media are expressed in a pressure - velocity formulation and can be solved together with the coupled state evolution equation ( [ eq : lin_dsndt ] ) by means of standard finite differences numerical techniques @xcite . in this way , lossless linear acoustics equations can be obtained by setting @xmath48 or in the limit when the relaxation times @xmath49 . the relaxation behavior described by this linearized model is achieved too by the formulation described in ref . @xcite , where the relaxation coefficients @xmath26 and the relaxation variable @xmath29 are defined in a different , but analogous way . equations ( [ eq : continuity ] , [ eq : momentum ] , [ eq : st_final_cinf ] , [ eq : dsdt ] ) are solved by means of finite - differences in time - domain ( fdtd ) method . cylindrical axisymmetric @xmath50 coordinate system is considered in this work , however , the method can be derived in other coordinate system . as in the standard acoustic fdtd method @xcite , the particle velocity fields are discretized staggered in time and space respect to the density and pressure fields , fig . uniform grid is considered , where @xmath51 , @xmath52 , @xmath53 , with @xmath54 and @xmath55 as the radial and axial spatial steps , and @xmath56 is the temporal step .     and the @xmath20-th order relaxation process state fields @xmath57 are evaluated at same discrete location as the density @xmath58 . particle velocity fields are discretized staggered in both space and time respect to the density , pressure and the @xmath20-th order state fields.,width=302 ]    thus , centered finite differences operators are applied over the partial derivatives of the governing equations . in order to fulfill the conservation principles over each discrete cell of the domain interpolation is needed over the off - center grid variables @xcite . the @xmath59 component of eq . ( [ eq : momentum ] ) is expressed in a cylindrical axisymmetric system as @xmath60 in the present method , each term of the above expression is approximated by centered finite differences evaluated at @xmath61 , @xmath62 , @xmath63 . this equation can be solved obtaining an update equation for @xmath64 . in the same way , the @xmath65 component of the motion equation ( [ eq : momentum ] ) is expressed as @xmath66 each term of this expression is approximated by centered finite differences and evaluated at @xmath67 , @xmath68 , @xmath69 . an update equation is obtained solving this equation for @xmath70 . the equation ( [ eq : continuity ] ) in cylindrical axisymmetric coordinate system is expressed as @xmath71 following the same procedure , each term of the above expression is approximated by centered finite differences and evaluated at @xmath72 , @xmath73 , @xmath74 , and the update equation is obtained solving this expression for @xmath75 . finally , the equation ( [ eq : dsdt ] ) is solved for @xmath76 as @xmath77      in order to fulfill the conservation principles for all the discrete cells an averaging of the off centered values is proposed @xcite , and the temporal averaging of the fields leads to an implicit system of equations . however , to obtain an explicit solution , the equations can be solved using an iterative solver . in this way , the future field values depend not only in the past values but also on unknown future values : @xmath78 @xmath79 @xmath80 the fields in the above expressions with temporal index @xmath81 are the temporal averaged values and are calculated by the iterative solver as : @xmath82 @xmath83 @xmath84 the initial value for the iterative solver is obtained by means of the explicit fdtd linear solution . this value is injected into the implicit system ( [ eq : solver1]-[eq : solver3 ] ) . due to the fact that the linear value is very similar to the nonlinear solution for a small time step , the iterative algorithm converges fast . in all the simulations presented in the validation section one iteration is used . finally , a leap - frog time marching is applied to solve each time step until the desired simulation time is reached . perfectly matched layers ( pml ) @xcite were placed in the limits of the domain ( @xmath85 and @xmath86 ) to avoid spurious reflections from the limits of the integration domain . inside the pml domains linearized acoustics equations were solved using the complex coordinate screeching method @xcite . these absorbent boundary conditions have reported an attenuation coefficient of 55.2 db for a layer of 30 elements and broadband wave with 1 mhz central frequency and non - normal incidence angle . the staggered grid is terminated on velocity nodes , so the coupling to the pml layers is implemented by only connecting these external nodes . this allows us to prevent the singularity of the cylindrical coordinate system : due to the staggered grid , the only variable discretized at @xmath87 is @xmath88 , and axisymmetric condition @xmath89 is applied there . to solve spatial differential operators on boundaries some  phantom " nodes must be created with the conditions : @xmath90      the stability for the lossless algorithm follows the courant - friedrich - levy ( cfl ) condition , so for uniform grid @xmath91 the maximum duration of the time step is limited by @xmath92 where @xmath93 is the number of dimensions ( in our case , cylindrical axisymmetric coordinate system , @xmath94 ) . however , numerical instabilities have been observed when @xmath95 . due to this empirical relation , the maximum values for relaxation frequencies are limited too by the chosen spatial discretization by the simple relation @xmath96 where @xmath97 is the maximum relaxation frequency for all processes , @xmath98 is the number of elements per wavelength and @xmath99 the frequency of the propagating wave . at higher amplitudes nonlinear effects induce the growing of harmonics of the fundamental frequency of the initial wave . the diffusive viscosity term in the motion equation , eq . ( [ eq : momentum ] ) , attenuates the small wavenumbers , damping the  node to node " numerical oscillations and ensuring numerical stability in weakly nonlinear regime . however , when discontinuities are present in the solution , extra numerical techniques must be employed to guarantee convergence . in this work artificial viscosity @xcite was added when shock waves are present in the solution . thus , by choosing the artificial viscosity as @xmath100 the algorithm shows consistency when @xmath101 , so if stability is achieved by the cfl condition , the convergence is guaranteed . thereby , the viscosity term ensures the entropy condition is fulfilled @xcite , and the method is able to solve shock waves in viscous fluids . in order to validate the frequency dependent attenuation and dispersion of the numerical method a simulation was done in linear regime including a single relaxation process . a homogeneous medium was considered , where the parameters were set to the typical values for water at 20 c : @xmath102 m / s , @xmath103 kg / m@xmath104 , @xmath105 , @xmath106 pa@xmath107s . a single relaxation process was included , with a characteristic relaxation time of @xmath108 and relaxation modulus of @xmath109 that leads to a frozen sound speed of @xmath110 m / s . in this case , the numerical parameters were set to @xmath111 m and @xmath112 s. a plane wave front traveling in @xmath113 direction was considered . thus , at @xmath114 the media was excited with a negative second derivative of a gaussian function : @xmath115 here , the central frequency of the broadband signal was set to @xmath116 mhz , and the amplitude was set to @xmath117 @xmath118 pa , small enough to neglect the nonlinear propagation terms . as the wave propagates , a single relaxation process change the complex amplitude as @xmath119 , where @xmath120 can be expressed as @xcite @xmath121 thus , the theoretical attenuation for the relaxation processes was estimated from the real part of @xmath120 , and including viscous losses leads to @xmath122 in order to compute the attenuation of the numerical method , simulated pressure was recorded at two locations @xmath123 and @xmath124 , and attenuation and phase velocity were estimated from the spectral components over the bandwidth of the input signal . the numerical attenuation was calculated as @xmath125 where @xmath126 is the fourier transform of the measured pressure waveforms at points @xmath123 and @xmath124 . the results , plotted in fig . 2 , show excellent agreement between the attenuation obtained by numerical solution and the predicted by the theory for a single relaxation process .     and high @xmath127 frequency limit attenuation are plotted as dash - dotted lines.,width=302 ]    for small amplitude acoustic perturbations in a weakly dispersive media , the phase velocity due to the relaxation processes is related to the attenuation through the kramers - kronig relations @xcite . thus , the theoretical phase velocity can be predicted for relaxing media as @xcite @xmath128 on the other hand , in order to study the propagation speed for each spectral component of the numerical method , the phase velocity is computed as @xmath129 using the same simulation parameters , fig . 3 shows the retrieved phase velocity and the theoretical dispersive response of the single relaxation process . the agreement between the model solution and the analytic prediction is again excellent , demonstrating that the inclusion of relaxation processes in nonlinear equations by means of the proposed model exhibit attenuation and dispersion in a correct way . however , typical fdtd numerical dispersion can be observed in the high frequency limit ( @xmath130 mhz ) , where the cumulative phase error lead to phase velocity mismatching . this error can be mitigated in the present numerical algorithm by increasing the number of elements per wavelength . an optimization algorithm has been used to fit the numerical attenuation response due to multiple relaxation , eq . ( [ eq : alpha_t ] ) , to a frequency power law @xmath131 . in order to find the proper relaxation coefficients , this algorithm uses the @xmath132 function in the optimization toolbox in matlab v7.13 . thus , an optimization of the relaxation times @xmath25 and relaxation modulus @xmath26 values has been done in order to minimize the relative error between the target power law and the computed attenuation . similar results and accuracy , but longer computational times , where obtained by using a genetic algorithm for the relaxation coefficients optimization . plane wave propagation was considered and simulation parameters were @xmath102 m / s , @xmath103 kg / m@xmath104 , @xmath105 , @xmath133 pa@xmath107s , @xmath116 mhz , @xmath134 m , @xmath135 s ; that leads to 26 elements per wavelength and a cfl number of 0.9 . only two independent relaxation processes were employed in this section to obtain the target frequency power laws .    following the above procedure , the relaxation times @xmath25 and relaxation modulus @xmath26 were optimized for different frequency power laws covering the range of that observed in tissues @xmath136 $ ] . the attenuation coefficient @xmath137 was chosen for each fit to present an attenuation @xmath138 db / cm at 10 mhz . the fitting was developed over the typical frequency range for medical ultrasound applications , i. e. 1 to 20 mhz . the results for the attenuation curves are plotted in fig . 4 , where the theoretical and the numerical predictions agree over the frequency range used for the fitting .    $ ] and the attenuation retrieved.,width=302 ]    the fitted exponents for each curve and the goodness of the fit are listed in table [ table:1:fit ] , showing that two relaxation process are enough to correctly exhibit a power law in a relative wide frequency range . .results for the fitting procedure for different power frequency attenuation laws [ cols= " < , < , < , < , < " , ]     [ table:3:deltac ]      in order to study the convergence of the numerical calculations to an analytical solution of the model in the nonlinear regime , a medium with frequency squared dependence attenuation is implemented using the adequate relaxation times and relaxation modulus as explained above . the solution for the frequency squared absorption is compared with the analytical solution for a plane wave traveling through a thermo - viscous fluid proposed by mendousse @xcite : @xmath139 where @xmath140 is the goldberg number as @xmath141 . here the diffusion coefficient @xmath142 , @xmath143 is the shock formation distance , the acoustic mach number @xmath144 , @xmath145 is the particle velocity excitation at @xmath146 and the wavenumber @xmath147 . ( white circles ) , @xmath148 ( black circles ) and @xmath149 ( asterisk).,width=302 ]    in this case , a @xmath116 mhz sinusoidal plane wave propagation was considered with an amplitude of 6 mpa , leading to a shock propagation distance of @xmath150 . thus , defining the distance normalized by the shock formation distance @xmath151 , the pressure waveforms @xmath152 for the numerical and the analytic solutions are shown in figure 6 for @xmath153 . the numerical solution agrees with the predicted by the analytical solution ( eq . [ eq : medousse ] ) at different distances .    in order to study the accuracy of the algorithm , the amplitude of the first three harmonics has been extracted for the numerical and analytic solutions and plotted versus @xmath154 in figure 7 . this result shows that the wave steepening due to the nonlinear processes are well resolved by the numerical method presented here . however , as a consequence of the numerical dispersion of the fdtd algorithm , the numerical sound speed does not exactly match the theoretical sound speed , so a phase error is always present on the solution . the observed relative error for the first three harmonic amplitude decreases due to grid coarsening by a square law ( i. e. the numerical method is second order accuracy ) . thus , the magnitude of the error mainly depends on the traveled propagated distance and the number of elements per wavelength . in the case of a traveled distance of 85 @xmath155 , a grid of 32 elements per wavelength is needed to obtain a relative error below 1% for the third harmonic . in this case the relative error of the lower harmonics is always lower , where the fundamental harmonic error was 0.072 % . the grid can be reduced if shorter path length are considered or if tolerance requirements for the higher spectral components are less restrictive , saving computational resources in higher dimensional reference systems ( i.e. 3d . cartesian or cylindrical coordinate systems ) . for the three first harmonic components of the analytical ( continuous lines ) and numerical solutions , where the evolution of the fundamental component @xmath99 ( white circles ) , @xmath156 ( black circles ) , @xmath157 ( asterisks ) and @xmath158 ( crosses ) agrees with the analytical solution.,width=302 ] in the present work a model for nonlinear acoustic waves in relaxing media is presented in time - domain formulation which does not require convolutional operators . a numerical solution by means of finite differences in time domain have been obtained , showing that the theoretical attenuation and dispersion due to relaxation processes can be achieved by the numerical method with accuracy . these result can be used to model typical relaxation process ( e. g. the processes observed in air , associated to the molecules of oxygen and nitrogen , or in seawater , associated to the relaxation of boric acid and magnesium sulfate ) .    moreover , in this work a method for model a frequency power law on attenuation by means of multiple relaxation was described . after an optimization of the relaxation coefficients , numerical results show that a pair of relaxation processes are enough to obtain frequency power law with accuracy ( relative error below one percent ) for medical ultrasound predictions . therefore , only two auxiliary fields are necessary to model a frequency power law attenuation over a typical medical ultrasound frequency band .    besides , the proposed method can achieve local variations of the power in the frequency power law , so an arbitrary attenuation curve can be modeled by means of the proper optimization of the relaxation coefficients . the only observed limitation in that the local exponent of the frequency power law must be lower than 2 . this feature of the presented method is an advantage when compared with most fractional derivatives methods , where the attenuation follows an exact but unique frequency power law over the entire frequency range . on the other hand , due to the numerical behavior of the finite differences in time domain method , an increase on the dispersion has been observed , mismatching the phase velocity in the higher frequency limit . this effect leads to cumulative phase errors on large propagation distances , so a convenient grid refinement must be employed , increasing the computational cost of the algorithm . moreover , if path length is on the order of thousand wavelengths and phase error is crucial for a specific prediction , correct dispersion can be accurately achieved by _ k_-space or pseudospectral numerical methods saving computational time and memory . [ acknowledgments ] the work was supported by spanish ministry of science and innovation through project fis2011 - 29731-c02 - 01 . n. jimnez acknowledges financial support from the universitat politcnica de valncia ( spain ) through the fpi-2011 phd grant .
Abstract:

Answer: this paper describes a model for nonlinear acoustic wave propagation through absorbing and weakly dispersive media , and its numerical solution by means of finite differences in time domain method ( fdtd ) . the attenuation is based on multiple relaxation processes , and provides frequency dependent absorption and dispersion without using computational expensive convolutional operators . in this way , by using an optimization algorithm the coefficients for the relaxation processes can be obtained in order to fit a frequency power law that agrees the experimentally measured attenuation data for heterogeneous media over the typical frequency range for ultrasound medical applications . our results show that two relaxation processes are enough to fit attenuation data for most soft tissues in this frequency range including the fundamental and the first ten harmonics . furthermore , this model can fit experimental attenuation data that do not follow exactly a frequency power law over the frequency range of interest . the main advantage of the proposed method is that only one auxiliary field per relaxation process is needed , which implies less computational resources compared with time - domain fractional derivatives solvers based on convolutional operators .


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: one of the simplest and most common ways to categorize an entity as living or non - living is to stipulate that if we were to decompose that entity into its constituent parts and could then reconstruct the whole by assembling the parts , it is non - living ; otherwise , it is living . that is , a living thing is more than just the collection of its parts , but a non - living thing is not . then how is the gap between living and non - living things expressed ? r. rosen  @xcite claims that a central feature of living things is _ complexity _ , where a system is said to be _ complex _ if its behavior can not be captured by models of that system ; otherwise , that system is said to be _ simple_. ( that is , the complexity of a system depends primarily on the models of the system , and only secondarily on the system per se . ) in rosen s description , use of the term _ complex _ in reference to a system is equivalent to use of the term _ incomputable _ or _ notwell - formed _ in reference to its models .    in his model of living things , rosen considers a logical paradox ( such as russell s paradox  @xcite ) as a _ metabolism - repairsystem _ ( _ mrsystem _ )  @xcite . an mr system consists of two sets ( a set @xmath6 of _ rawmaterials _ and a set @xmath7 of _ behaviors _ ) and three functions ( a _ metabolicfunction_@xmath8@xmath9 , a _ repairfunction_@xmath10@xmath11 , and a _ replicationfunction_@xmath12@xmath13 ) , where @xmath14 and @xmath15 are onto functions and @xmath16 . ( for sets @xmath17 denotes the set of all morphisms from @xmath18 to @xmath19  @xcite . ) one remarkable feature of rosen s model is the following :    1 . there do not exist such onto functions @xmath20 and @xmath21 . if they are assumed to exist , we obtain a contradiction .    one way to treat such a contradiction is to invoke hyperset theory  @xcite . a _ hyperset _ is defined as a graphable set  in particular , one that can be represented by a hyperset diagram , which is a digraph ( directed graph ) with the property that every node is either a set contained in the hyperset or an element of such a set , and every edge is directed from some set  @xmath22 contained in the hyperset to one of the elements of  @xmath22  @xcite . thus edge @xmath23 expresses the relationship @xmath24 , i.e. , @xmath25 is an element of  @xmath26 . for example , the hyperset diagram of @xmath27 is shown in fig . [ fighyperset1_p5 ] . there are two directed edges : @xmath28 and @xmath29 . a relationship of the form @xmath28 leads to the well - known contradiction known as russell s paradox  @xcite ; however , such a relationship is interpreted as just a loop structure in hyperset theory .    , width=132,height=43 ]    property ( i ) can be expressed in hyperset theory by its way . let us represent @xmath30 by the set @xmath31 ; thus @xmath32 and @xmath33 are represented by the sets @xmath34 and @xmath35 respectively . then we obtain the hyperset diagram of mr as shown in fig . [ fighyperset2_p5 ]  @xcite . accordingly , a system can be said to be complex if it can not be well - formed in standard set theory but it can be well - formed in hyperset theory . matsuno and gunji , in their theory of internal measurement  ( i m )  @xcite , propose a different way , called a _ perpetualequilibratingmechanism_. they address property ( i ) as follows :    1 . any solution to a problem ( paradox ) is inevitably a pseudo - solution , and the pseudo - solution at any given time step `` always '' triggers a problem that that is to be solved at the next time step . this leads to an evolutionary process that goes on forever and thereby allows for emergent properties . therefore , the paradox is treated as the force that drives the time evolution of the system .    in i m , the gap between the whole and the aggregation of the parts drives the time evolution of the dynamical system . therefore , i m can be understood as a proposal of a dynamical system which is based on the idea of a perpetual loop between two different logical layers , such as the parts and the whole , with a decomposition @xmath36 of the whole into the parts , and an integration @xmath37 ( as a pseudo - solution ) of the parts into the whole . thus we may express an i m model in terms of a triad : the two different logical layers and a mediator / interface to adjust the two layers in an _ inconsistent _ manner  @xcite . i m claims that the relationship between the two layers ( or the relationship between @xmath36 and @xmath37 ) is not consistently determined , and for this reason they are perpetually changing relative to each another , where the term inconsistent means that we obtain a logical paradox if we assume consistency between them . any finite boolean network ( bn ) has a fixed transition matrix , hence a finite bn ultimately reaches an attractor and stays there ( i.e. , the computation ends ) . in section  [ secvbn_p5 ] , we define a virtual boolean network ( vbn ) as a bn in which the set of incoming nodes for every node is the entire set of nodes of the bn , and we show that every finite bn can be treated as a vbn . in section  [ secexpansions_p5 ] , we introduce a labeling function to label the states of a vbn , and we discuss the probabilistic boolean network ( pbn ) and the hidden markov model ( hmm ) as extensions of a vbn . the purpose of an hmm is to find an optimal solution ; in the case of a living thing , however , it is not always possible to find such a solution . thus , in section  [ secdbn_p5 ] , we abandon the attempt to find an optimal solution , and we introduce a dynamical boolean network ( dbn ) , which is a vbn whose set of nodes is fixed but whose transition matrix can change from one time step to another . thus even a finite dbn does not necessarily have any attractors . in our dbn , a labeling function plays a role in maintaining a perpetual loop between the of the vbn at different time steps .    in section  [ secdbn_p5 ] , we present a construction of a dbn  a construction that incorporates stochasticity in a number of respects . in section [ secsimulation_p5 ] , we present the results of a number of simulations of a dbn with two nodes . the purpose of each simulation was to investigate the set of rule vectors actually visited by the dbn , which is equivalent to the set of transition matrices actually visited by the dbn . in subsection [ subsectype1_p5 ] , we present the results of a simulation of a dbn based on the construction in section  [ secdbn_p5 ] , which we call a type  1 simulation . the simulations presented in the remaining three subsections of section  [ secsimulation_p5 ] use dbn constructions which are similar to the one in section  [ secdbn_p5 ] , the only difference being that , for @xmath38 , we chose the transition matrix @xmath0 of the dbn for time step  @xmath1 to be @xmath39 , where @xmath3 is the transition matrix in our original construction and @xmath40 is the matrix representation of some permutation @xmath41 of the states of the dbn . those three simulations correspond to three different classes of permutations . in our type  2 simulation , @xmath41 was randomly chosen from the set of permutations that consist of either the identity or a single 2-cycle or the product of two disjoint 2-cycles . in our type  3 simulation , @xmath41 was randomly chosen from the set of all 24 permutations of the four states of a dbn with two nodes . in our type  4 simulation , @xmath41 was constructed from a labeling function . there are 81 rule vectors ( out of a total of 256 possible rule vectors for a dbn with two nodes ) that were not visited at all after the fifth time step of any trial in the type  1 and type  2 simulations . however , all 256 rule vectors were visited in every trial of the type  3 simulation , as well as in slightly more than half of the trials of the type  4 simulation . indeed , six of the 81 aforementioned rule vectors were the ones visited most frequently of all in the type  4 simulation . a total of six single - node transition rules are associated with those six rule vectors ( each rule vector corresponds to an ordered pair of single - node rule vectors ) . three of those six single - node transition rules are transition rules of a linear vbn ( lvbn ) ; the other three are the negations of the first three . thus our conclusion is that the type  4 simulation presents new insight into linearity of a vbn . a boolean network ( bn )  @xcite is a digraph in which each vertex has an internal state and a transition rule . the vertices of a bn are hereinafter referred to as nodes . each internal state has a value of either 0 or 1 at any given time . the transition rule for a given node specifies the next state of that node as a function of the current internal states of its incoming ( input ) nodes .    a bn can be finite or infinite . the discussion of bn s in this paper is limited to finite bn s ( i.e. , those that have finitely many nodes ) . some of the assertions made herein about bn s do not extend to the infinite case . an example of a bn is shown in fig . [ kauffman1_p5 ] . there are three nodes ( @xmath42 ) , so the bn has a total of @xmath43 ( @xmath44 ) possible states:@xmath45 the set of incoming nodes for each node of a bn is specified by the digraph . here , the sets of incoming nodes for @xmath46 , and @xmath47 are @xmath48 , and @xmath49 , respectively . since nodes @xmath26 and @xmath25 have just one incoming node apiece , each of them has just two possible input states : @xmath50 and @xmath51 . node @xmath47 has two incoming nodes , hence it has four possible input states : ( 0,0),(0,1),(1,0),(1,1 ) . given the initial state of a bn , we obtain the time evolution of the system ( the sequence of states visited ) by repeated application of the transition rules . since the number of nodes is finite , the number of possible states is finite , so for every initial state the system evolves to some finite cyclic sequence of states , i.e. , a finite sequence of states of the bn which is visited _ ad infinitum from some time on , with no other states intervening_. such a sequence is called an attractor . clearly , each state evolves to one and only one attractor . thus the set of states of a bn can be partitioned in such a way that states in the same block of the partition evolve to the same attractor , and states in different blocks evolve to different attractors . all possible states and attractors of the bn in fig . [ kauffman1_p5 ] are depicted in the transition diagram in fig . [ kauffman2_p5 ] . in this particular bn , only four of the states are visited infinitely often , and there are two attractors : one of length 3 ( @xmath52 ) and one of length 1 ( @xmath53 ) . ( the length of an attractor is the number of states in the cycle . ) the set of states can be partitioned as @xmath54 , where @xmath55 and@xmath56 and every state in  @xmath57 evolves to the attractor of length 3 .    , width=415,height=113 ]    a digraph with node set @xmath58 is specified by assigning to every node a set of incoming nodes . assigning a set of incoming nodes to a single node can be regarded as selecting an element of @xmath59 , where @xmath60 denotes the set of all morphisms from @xmath18 to @xmath19  @xcite . ( here , 2 denotes the set @xmath61 . ) the term _ morphism _ as used here is synonymous with _ function _ ; hence @xmath59 is the set of all binary functions ( functions whose values are elements of the set @xmath62 ) with @xmath58 as domain . each such binary function can be represented as a binary sequence ( a sequence of 0 s and 1 s ) in which the @xmath63th term is the value of the function at node @xmath63 . the set of incoming nodes for node @xmath64 is represented by the binary sequence in which the @xmath63th term is 1 if node @xmath63 is an incoming node for @xmath64 , and 0 otherwise . thus assigning to every node in @xmath58 a set of incoming nodes ( i.e. , selecting an element of @xmath59 for every node ) can be regarded as selecting an element of @xmath65 . alternatively , specification of a digraph can be regarded as selecting an element of @xmath66 , since @xmath67 via the mapping of @xmath68 to the element @xmath69 of @xmath65 which is defined by @xmath70 . for every @xmath71 , denote @xmath72 by @xmath73 . then @xmath74 and , for every @xmath75 , @xmath76 note that @xmath77    the internal states of the nodes of a bn with node set  @xmath58 can be regarded as the components of a state  @xmath78 , where @xmath79 is the number of nodes in  @xmath58 and @xmath80 is the internal state of node @xmath64 ( i.e. , @xmath81 ) . there are two ways to specify the transition rules for a bn . the first method is completely general . [ deffirstruledescription_p5 ] the transition rule @xmath82 for node @xmath64 assigns the next internal state ( either 0 or 1 ) to node  @xmath64 . the assignment is a function of the current internal states of the nodes in the set @xmath83 of incoming nodes for @xmath64 . let @xmath84 , where @xmath85 is the number of nodes in  @xmath83 . then@xmath86 where @xmath87    the transition rules for the nodes of a bn can be regarded as the components of a rule vector  @xmath88:@xmath89    the second method of specifying the transition rules applies only to a bn for which @xmath90 ( hence @xmath91 ) for every node @xmath64 , where @xmath92 denotes @xmath93 . we will refer to such a bn as a virtual boolean network ( vbn ) . [ defsecondruledescription_p5 ] if @xmath90 for every @xmath64 , then the transition rules @xmath82 for the individual nodes can be incorporated into a single transition rule @xmath94 : @xmath95    for node @xmath64 , denote @xmath96 by @xmath97 . for @xmath98 , we will say that node @xmath63 is _ virtually disconnected _ from node @xmath64 if @xmath99 is invariant with respect to substitution of input @xmath100 with @xmath101 , where@xmath102 we will say that the nodes which are not virtually disconnected from node @xmath64 are _ virtual incoming nodes _ for node  @xmath64 , and we will denote the set of virtual incoming nodes by @xmath103 ; hence the set of nodes which are virtually disconnected from node @xmath64 is @xmath104 . any bn can be regarded as a vbn , by defining @xmath97 in such a way that it extends @xmath82 and all the nodes in @xmath92 are virtually disconnected from @xmath64 . [ ex1_p5 ] let @xmath105 , and fix a node @xmath64 . suppose that @xmath106 and @xmath107 , and that the transition rule @xmath82 for node @xmath64 is @xmath108\ ] ]    the domain of rule @xmath82 is just the set of states of node @xmath26 , but @xmath82 can be embedded into rule @xmath97 defined by @xmath109\begin{array}{lc|cccc } \emph{in}&(x_{a},x_{b},x_{c})&(1,0,0)&(1,0,1)&(1,1,0)&(1,1,1)\\\hline \emph{out}&(x_{i})&0&0&0&0\end{array } \end{array}\right]\ ] ]    the domain of @xmath97 is the set of states of all the nodes . note that , for all @xmath110 , the embedding satisfies @xmath111 hence nodes @xmath25 and @xmath47 are virtually disconnected from node  @xmath64 . we can translate a specification of a vbn into a specification of a bn , though the resulting bn is not necessarily unique ( since @xmath83 can be any superset of @xmath103 ) . for node  @xmath64 in example  [ ex1_p5 ] , we could define @xmath83 as @xmath49 , and the transition rule @xmath112 for node @xmath64 as    @xmath113\ ] ]    note that , for all @xmath114,@xmath115 furthermore , both @xmath82 and @xmath112 embed into the transition rule @xmath97 given in example  [ ex1_p5 ] . the transition diagram of a bn is a digraph in which the vertices are the states of the bn and every vertex has just one out - going edge . if the set of nodes of the bn is  @xmath58 , then the bn has a total of @xmath116 states . thus the transition diagram can be represented as a @xmath117 _ adjacency matrix _ . the adjacency matrix @xmath3 of any digraph is a binary matrix ( i.e. , a matrix of 0 s and 1 s ) such that @xmath118 if there is an out - going edge from the vertex labeled  @xmath64 to the vertex labeled  @xmath63 ( and @xmath119 otherwise ) . since each state of a bn has just one out - going edge , exactly one of the entries in each row of the adjacency matrix of a bn is a 1 . in the case of a bn , the adjacency matrix is also called the _ transition matrix_.    [ proprules2trans1_p5 ] the transition matrix @xmath3 of a bn with vertex set  @xmath58 can be obtained from the rule vector @xmath88 ; however , @xmath88 can not be recovered from @xmath3 . * proof * ( @xmath120 ) : label the rows of  @xmath3 with the states of the bn , and use the same labels ( in the same order ) for the columns of  @xmath3 . let @xmath78 be a state of the bn . then , in the row of  @xmath3 which is labeled with state  @xmath121 , the unique 1 is located in the column labeled with state  @xmath122 , where , for every  @xmath123  is computed by applying transition rule @xmath82 to  @xmath121 . + ( @xmath124 ) : to determine transition rule  @xmath82 , its domain ( namely , @xmath83 ) must be known ; however , @xmath83 can not be derived from  @xmath3 . @xmath125    the following proposition follows from the fact that , for every node in a vbn with node set  @xmath58 , the set of incoming nodes is @xmath58 . [ proprules2trans2_p5 ] in a vbn , the transition matrix @xmath3 can be derived from the transition rule  @xmath94 , and @xmath94  can be derived from  @xmath3 . * proof * by proposition  [ proprules2trans1_p5 ] , it suffices to prove that @xmath126 . assume that the rows and columns of  @xmath3 are labeled with the states of the vbn , as indicated in the proof of proposition  [ proprules2trans1_p5 ] , and let @xmath127 . the output of transition rule  @xmath97 , which is the next internal state of node  @xmath64 , is a function of the current state of the bn , so let @xmath121 be a state of the vbn . furthermore , let @xmath63  be such that , in the row of  @xmath3 which is labeled with state  @xmath121 , the unique 1 is located in column  @xmath63 , and let @xmath122 be the state that labels column  @xmath63 . then @xmath128 is the output of @xmath97 that corresponds to input  @xmath129    from this point on , we assume that every bn is a vbn , so we need not specify the set of incoming nodes for any node . every state @xmath121 of a vbn with node set  @xmath58 is an element of @xmath59 , and the transition rule @xmath94 for such a vbn is an element of @xmath130 .    [ deflbn_p5 ] a vbn is said to be a _ linear vbn ( lvbn ) _ if , for every node  @xmath64 , transition rule @xmath97 is of the form @xmath131 where @xmath121 is the current state of the vbn and @xmath132 is a vector whose components are 0 s and 1 s . ( clearly , the number of components of  @xmath132 must be equal to the number of nodes in the vbn . ) [ propmbn1_p5 ] in an lvbn , node  @xmath64 is virtually disconnected from node  @xmath63 if and only if the @xmath63th component of the vector @xmath132 in definition  [ deflbn_p5 ] is 0 .    * proof*the output of transition rule  @xmath97 ( namely , @xmath133 ) is invariant under the transformation that consists of replacing the current internal state  @xmath134 of node  @xmath63 with  @xmath135 ( and leaving the current internal states of the other nodes unchanged ) if and only if the @xmath63th component of @xmath132 is 0 . @xmath125    for a vbn with set of states  @xmath22 and node set  @xmath58 , the number of possible single - node transition rules is @xmath136 . in the case of a vbn with two nodes , this is @xmath137 . the output of each of the 16 possible single - node transition rules for a vbn with two nodes is given in table  [ tablerule_p5 ] as a function of the input state of the vbn ; tabulated along with the output of each transition rule is the number @xmath138 of virtual incoming nodes for each of the node(s ) to which that rule would be applied . there are only four possible transition rules for a node of an lvbn with two nodes ( namely , the rules which are presented in boldface in the table ) . those four transition rules ( 1 , 7 , 11 , and 13 ) correspond to the vectors @xmath139 , and @xmath140 , respectively , where @xmath141 is as in definition  [ deflbn_p5 ] . if nodes 1 and 2 of an lvbn are @xmath26 and @xmath25 ( in that order ) , then the sets of virtual incoming nodes for each of the nodes to which rules 1 , 7 , 11 , and 13 would be applied are @xmath142 , and @xmath143 , respectively . .outputs of the transition rules for a vbn with two nodes . [ cols="^,^,^,^,^,^,^,^,^,^,^,^,^,^,^,^,^,^",options="header " , ]     the 1,000 trials of the type  4 simulation can be grouped into two classes :    * trials in which a handful of rule vectors were visited very frequently ( visited at roughly 1,000 or more different time steps ) and the remaining rule vectors were visited very infrequently ( visited at roughly no more than 30 different time steps ) * trials in which the variation in the number of visits per rule vector was very small    about 70@xmath144 of the trials of the type  4 simulation fall into class ( i ) ; the remaining 30@xmath144 are in class ( ii ) . all 94 trials with less than 90@xmath144 coverages are in class ( i ) . the numbers of time steps at which rule vector @xmath145 was visited in two specific trials of the type  4 simulation with at least 90@xmath144 coverage apiece are shown in the maps in fig . [ figresulttype4_3_p5 ] . the map in ( a ) is for a typical class ( i ) trial , and the map in ( b ) is for a typical class ( ii ) trial . only six rule vectors ( namely , @xmath146 , and @xmath147 ) were visited very frequently ( visited at roughly 1,000 or more different time steps ) in the class  ( i ) trial associated with fig . [ figresulttype4_3_p5](a ) . white is used for the remaining 250 rule vectors , each of which was visited no more than 22 times . in about 20@xmath144 of the class ( i ) trials , at least one rule vector other than those six was among the rule vectors visited very frequently ( visited at roughly 1,000 or more different time steps ) .    the map shown in fig . [ figresulttype4_3_p5](b ) , for a class ( ii ) trial , bears a much greater resemblance to the maps for the trials of the type  3 simulation than to those for the class ( i ) trials of the type  4 simulation . was visited in two trials of the type  4 simulation with at least 90% coverage apiece : ( a ) a typical class ( i ) trial and ( b ) a typical class ( ii ) trial , width=491,height=229 ]    all six of the rule vectors enumerated earlier are from fig . [ figresultiti3_p5](a ) . each of those six rule vectors was visited a total of more than 300,000 times in the class  ( i ) trials of the type  4 simulation , while each of the 175 rule vectors not from fig . [ figresultiti3_p5](a ) had a total of fewer than 50,000 visits in the class  ( i ) trials . the total number of visits ( in units of @xmath148 ) of rule vector @xmath145 in the class ( i ) trials of the type  4 simulation ( summed over all 1,000 trials ) is shown in the map in fig . [ figbc_p5](a ) . only 54 of the 256 rule vectors were visited very frequently in at least one of the class  ( i ) trials : 41 of the 81 rule vectors from fig . [ figresultiti3_p5](a ) , including the six rule vectors enumerated earlier ; and 13 of the other 175 rule vectors . those 54 rule vectors are depicted by red squares in fig . [ figbc_p5](b ) .    ) of rule vector @xmath149 in the class ( i ) trials of the type  4 simulation . ( b ) the 54 rule vectors ( indicated by red squares ) that were visited very frequently in at least one of the class ( i ) trials.,width=491,height=245 ]    the maximum time step at which rule vector @xmath145 was visited in the type  4 simulation ( where the maximum is taken over all 1,000 trials ) is shown in the map in fig . [ figmaxtimestep4_p5 ] . of the 256 rule vectors , 203 were visited at the very last ( 10,000th ) time step of at least one trial ; those rule vectors are depicted in white in the figure . red is used for the other 53 rule vectors , each of which was visited at some time step in the range 9,9929,999 but was never visited at the last time step . was visited in the type  4 simulation ( where the maximum is taken over all 1,000 trials),width=302,height=238 ]    fig . [ figresulttype4_4_p5 ] shows a graph of the cumulative percentage of rule vectors visited in our type 4 simulation , averaged over the 550 trials in which all 256 rule vectors were visited , vs.  the number of time steps . again , the blue curve is the result of our simulation , and the green curve is the parametric curve given by @xmath150 . the number of time steps actually taken to visit all 256 rule vectors ( in the 550 trials of our simulation with 100@xmath144 coverage ) was about 4,500 , which is about 2.9 times as large as the theoretical value and 2.7 times the number of steps taken in the type  3 simulation .    , where @xmath151 varies from 1 to 256 and @xmath152 is given by ( [ eqvisitratingfunction2_p5]).,width=283,height=223 ] at each time step in our simulations of a dbn , we incorporated randomness into several processes : in generating pseudo - transition diagrams from the output digraph , in choosing a pseudo - transition diagram and then transforming it to the transition diagram of a vbn , in choosing a labeling function for the states of the dbn , and ( in types  2 , 3 , and  4 ) in choosing a permutation of the states of the dbn .    in every trial of the type  1 simulation , fewer than 70@xmath144 of the rule vectors were visited . furthermore , the 256 rule vectors split into two groups : the set of 81 rule vectors that were not visited after the fifth time step in any trial of the first type  1 simulation ( and not visited at all in the two additional type  1 simulations ) , and the set consisting of the other 175 rule vectors , each of which was visited considerably more frequently than those in the set of 81 . those features of the type  1 simulation were largely preserved under the kinds of permutations of the states of the dbn which were chosen at random in the type  2 simulation ( namely , the identity permutation , the 2-cycles , and the products of disjoint pairs of 2-cycles ) . the salient features of the type  1 and type  2 simulations were not preserved under the kinds of permutations chosen in the type  3 and type  4 simulations . ( any permutation could be selected in type  3 ; however , the permutation @xmath41 selected at time step  @xmath1 in the type  4 simulation depended in part on the labeling function @xmath153 . ) a remarkable contrast was found between the results of the type  3 and type  4 simulations on one hand , and those for type  1 and type  2 on the other . all 256 rule vectors were visited in every trial of the type  3 simulation , and six of the 81 rule vectors which were not visited after the sixth time step of any trial in the type  1 and type  2 simulations were the rule vectors visited the most frequently of all  by far  in the type  4 simulation . the results of the one simulation we ran in which @xmath41 was randomly chosen from the complement of the set of permutations used in type  2 were similar to those for our type  3 simulation , in that the total number of visits per rule vector varied little from one rule vector to another . fully 41 of the aforementioned 81 rule vectors were visited very frequently ( visited at roughly 1,000 or more different time steps ) in at least one class  ( i ) trial of our type  4 simulation . six of those 41 rule vectors had a total of over 300,000 visits in the type  4 simulation . each component of every one of those six rule vectors is a single - node transition rule from the set @xmath154 . three of those ( rule numbers 7 , 11 , and 13 ) are transition rules of an lvbn ; each of the other three ( 4 , 6 , and 10 ) is the negation of a single - node transition rule of an lvbn . the only single - node transition rules of an lvbn which are not a component of any of those six rule vectors are rule numbers 1 and 16and both of those transition rules correspond to a node of an lvbn that has no virtual incoming nodes . thus it appears that we can regard the class ( i ) trials of the type  4 simulation as a restriction of the single - node transition rules of a vbn to the single - node transition rules of an lvbn and their negations , hence that linearity may have played a significant role in our type  4 simulation . therefore , linearity may be a worthwhile topic for further study .
Abstract:

Answer: we propose a dynamical boolean network ( dbn ) , which is a virtual boolean network ( vbn ) whose set of states is fixed but whose transition matrix can change from one discrete time step to another . the transition matrix @xmath0 of our dbn for time step  @xmath1 is of the form @xmath2 , where @xmath3 is a transition matrix ( of a vbn ) defined at time step  @xmath1 in the course of the construction of our dbn and @xmath4 is the matrix representation of some randomly chosen permutation @xmath5 of the states of our dbn . for each of several classes of such permutations , we carried out a number of simulations of a dbn with two nodes ; each of our simulations consisted of 1,000 trials of 10,000 time steps each . in one of our simulations , only six of the 16 possible single - node transition rules for a vbn with two nodes were visited a total of 300,000 times ( over all 1,000 trials ) . in that simulation , linearity appears to play a significant role in that three of those six single - node transition rules are transition rules of a linear virtual boolean network ( lvbn ) ; the other three are the negations of the first three . we also discuss the notions of a probabilistic boolean network and a hidden markov model  in both cases , in the context of using an arbitrary ( though not necessarily one - to - one ) function to label the states of a vbn . [ thm]note    internal measurement ; boolean network ; probabilistic boolean network ; hidden markov model ; linearity


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: epitaxy is an important concept for the fabrication of films with good crystalline quality like overlayers , multilayers , compound materials , and ordered alloys . such films are technologically important materials with adjustable electronic , magnetic , and optical properties  @xcite . epitaxial growth also allows to study fundamental aspects of low - dimensional structures and interfacial effects , and it allows to produce structures that are unstable in the bulk  @xcite . for example , body - centered cubic ( bcc ) co does not exist as a bulk phase but was stabilized by epitaxial growth on a gaas substrate  @xcite .    in coherent or pseudomorphic film growth , the in - plane film lattice parameters are determined by the substrate - in - plane lattice parameters . this allows to expose a material to static , non - isotropic but homogeneous strain which is a valuable mean to influence its intrinsic properties  @xcite . experiments on ferromagnetic bulk - like films showed , that the magnetic moments , the curie temperature , and the magnetic anisotropy can be tuned in a quasi - continuous manner by varying the substrate  @xcite . using piezoelectric substrates , even a continuous variation of the lattice parameters and of the related film properties was achieved  @xcite .    here , we predict by density - functional calculations scanning 30 non - magnetic metallic elements , that bulk - like epitaxial films of at least six metals show first - order martensitic transitions upon variation of the epitaxial strain . for example , we find jumps in the lattice parameter of la by @xmath3  and in the superconducting critical temperature of v by 20% . the examples considered in the present work include only the simplest case of tetragonal non - magnetic monoatomic films . we however emphasize that our main result , the possibility of isostructural phase transitions  @xcite under epitaxial strain and a related discontinuous dependence of intrinsic film properties on the interface lattice parameters , is not restricted to certain structures , compositions , or ground states . we model the structures of the considered films by the epitaxial bain path ( ebp )  @xcite . the ebp model applies to the bulk part ( i.e. , the interior ) of a coherently grown thick film  @xcite . it assumes that the film can be described by a bulk structure with periodicity in all three dimensions . due to epitaxial coherency , the in - plane film lattice parameters are rigidly tied to the substrate lattice parameters in the substrate - film interface via the underlying epitaxial relationship . on the other hand , the film is free to relax its geometry perpendicular to the interface due to the absence of forces . summarizing , the model neglects all surface and interface effects except the determination of the in - plane film lattice parameters by the substrate .    in the past , the ebp model was used to predict metastable states in transition metals  @xcite , to investigate the magnetic order in strained overlayers  @xcite , and to identify the mother phases of strained bulk - like films  @xcite . in the same context , we recently suggested ru , os , and u films to be ferromagnetic in certain ranges of epitaxial strain  @xcite . focusing on the important case of substrates with four - fold surface symmetry and body - centered tetragonal ( bct ) films with ( 001 ) orientation , we denote the in - plane film lattice parameter by @xmath4 and the out - of - plane lattice parameter by @xmath5 , see . the mentioned absence of forces perpendicular to the substrate - film interface allows for a relaxation of @xmath5 minimizing the total film energy @xmath6 at any fixed value of @xmath4  @xcite . hence , the ebp of a certain bct material is the set of points @xmath7 , for which @xmath8 the side condition determines the relaxation of @xmath5 for a given @xmath4 and defines the relation @xmath9 . quantities traced _ along _ the ebp only depend on @xmath4 , e.g. , the total energy @xmath10 . the difference between @xmath11 and the ground state energy of the film material , @xmath12 , corresponds to the general definition of an epitaxial strain energy  @xcite . a strict implementation of the notion ` relaxation of @xmath5 minimizing @xmath6 ' leads to relevant , until now disregarded consequences for all materials to which strained coherent epitaxy applies . as we show , the function @xmath13 can have two equal - valued global minima at @xmath14 for certain values of @xmath15 . thus , @xmath9 can be discontinuous at @xmath16 in a first - order transition . this includes , by epitaxial coherency , also the dependence of @xmath2 on the substrate lattice parameters in the interface plane . in this way , the possibility to tune overlayer - material specific properties can be restricted to a discontinuous manner . however , one can switch these properties on purpose , if one succeeds to control the lattice parameter @xmath4 in a narrow interval around @xmath16 , e.g. , by a piezoelectric substrate .    in the following , we verify the existence of degenerate global minima in @xmath13 by high - precision density - functional theory ( dft ) calculations , carried out with the full - potential local - orbital scheme fplo-7.00 - 28  @xcite . using the local - density approximation according to ref . @xcite ( pw92 ) and a scalar - relativistic mode for elements with atomic number @xmath17 , a full - relativistic mode otherwise , we scanned the ebp of 30 elements with the atomic numbers 20 - 23 , 29 , 30 , 38 - 48 , 56 , 57 , 71 - 80 , and 92 in a wide range of parameters @xmath4  @xcite . the convergence of numerical parameters and the basis set were carefully checked  @xcite . in order to converge the total energy per atom at a level smaller than @xmath18 mev , linear - tetrahedron integrations with blchl corrections were performed on a @xmath19 mesh in the full brillouin zone , apart from the elements cu , ru , ir , and pt , for which a denser @xmath20 mesh was required . all detected degenerate minima in @xmath13 were cross - checked by additional calculations employing the perdew - burke - ernzerhof ( pbe96 ) parameterization of the exchange - correlation functional  @xcite ( results listed in table  [ table : discontinuities ] ) . delineation of the bct lattice ( closed symbols ) in the fcc lattice ( @xmath21 , open symbols , two unit cells shown ) and in the bcc lattice ( @xmath22 ) . bct lattice and bcc lattice are on top of each other.,title="fig:",height=113]delineation of the bct lattice ( closed symbols ) in the fcc lattice ( @xmath21 , open symbols , two unit cells shown ) and in the bcc lattice ( @xmath22 ) . bct lattice and bcc lattice are on top of each other.,title="fig:",height=113 ] vanadium is one of the elements which ebp was studied in several investigations  @xcite . this provides us with the following , as we will argue incomplete , picture . the ebp of vanadium includes two energy minima : one corresponds to the stable bcc phase ; the other one at @xmath23  , @xmath24 , is a local minimum , see . this local minimum corresponds to an unstable bct bulk state  @xcite , since its energy is lowered by a shear deformation . et al . _ @xcite however succeeded to stabilize coherently grown vanadium - films on cu\{001 } and ni\{001 } which structures were determined to be strained derivatives of this bct state . the face - centered cubic ( fcc ) structure was found at a saddle point of @xmath6 and assigned to a maximum of @xmath25  @xcite . total energy of vanadium with respect to the energy of the stable bcc phase , @xmath26 , ( left hand ordinate ) and the electronic density of states ( dos ) at the fermi level , @xmath27 , ( right hand ordinate ) along the ebp . the two discontinuities of the ebp are marked by solid bars on both the @xmath4 and the @xmath28 scale . the first discontinuity at @xmath29 results in a pronounced drop in @xmath30 . the second discontinuity at @xmath31 is related to a peak in the dos . the energies are spline interpolated to guide the eye . the pw92 equilibrium lattice parameters of the cubic structures are @xmath32 ; @xmath33 . the brown broken line denotes the bcc structure . ]    by symmetry , the bcc structure ( @xmath22 ) and the fcc structure ( @xmath21 ) are stationary points of @xmath6  @xcite , but the side condition in definition determines , _ whether they are located at the ebp or not_. at variance with previous reports  @xcite , we find that the fcc structure does _ not _ belong to the ebp of vanadium , because it is unstable against a tetragonal distortion at fixed @xmath34 . we note , that in terms of linear elasticity theory this instability of fcc vanadium with respect to a relaxation of @xmath5 corresponds to a negative elastic constant @xmath35 , since @xmath36 . this is obvious from the plot of @xmath37 in . there , @xmath37 has a local maximum at @xmath38 and the global minimum is situated at @xmath39 , @xmath40 . the characteristic feature of this instability is the occurrence of a double well in @xmath37 vs.@xmath41 @xmath5 with the global minimum at @xmath2 and a local minimum . according to , double wells also occur for other values of @xmath4 in the vicinity of @xmath42 . the energies of both minima vary with @xmath4 such that the two minima are degenerate at a critical value @xmath43 . all in all , the behavior shown in is typical for a first - order phase transition . one can readily read off the energy barrier of somewhat less than 10 mev per atom , an upper limit for the hysteresis of @xmath44  , and the discontinuity of @xmath2 , @xmath45  .    in general , a saddle point in @xmath6 ( here , at @xmath46 ) and a discontinuity in the ebp do not have the same @xmath4-coordinates . in the case of vanadium , the already described ( first ) discontinuity occurs at @xmath47 . the associated change of @xmath28 amounts to @xmath48 . surprisingly , the ebp of vanadium shows yet another discontinuity at @xmath49 with @xmath50 .     cuts of the total energy @xmath6 of vanadium for several fixed in - plane lattice parameters @xmath4 ( legend ) near the fcc saddle point . for each @xmath4 , the global minimum at @xmath2 is indicated by a filled symbol . solid lines are spline interpolations to the data ( symbols ) and guide the eye . ] a depicts a contour plot of @xmath51 of vanadium together with the ebp . the two minima of @xmath51 correspond to the two minima of @xmath25 from and the fcc structure is the saddle point as demanded by symmetry . the contour diagram ( see inset ) reveals the following : there are pairs of lines of constant energy , @xmath52 , where each line of the pair encloses one of the two minima and the two lines overlap in their @xmath4-coordinates in the vicinity of the saddle point . the double wells in arise from this overlap . a similar peculiarity of the energy landscape at @xmath53  causes the second discontinuity and the related kink of @xmath25 visible in . contour plots in the tetragonal parameter space : ( a ) total energy landscape and ( b ) @xmath30 of vanadium . the two minima of @xmath51 are indicated by @xmath54 and the saddle point at the fcc structure is indicated by @xmath55 . the inset shows the saddle point region in more detail . the global minimum corresponds to the stable bcc phase , the local minimum is a bct structure . the ebp ( black dot - dashed line ) has two discontinuities , indicated by black dashed lines . brown hatched lines signal axial ratios of the cubic structures . the energy and the density of states were calculated on a uniform @xmath56 grid with a distance of @xmath57 between two adjacent points . ] contour plots in the tetragonal parameter space : ( a ) total energy landscape and ( b ) @xmath30 of vanadium . the two minima of @xmath51 are indicated by @xmath54 and the saddle point at the fcc structure is indicated by @xmath55 . the inset shows the saddle point region in more detail . the global minimum corresponds to the stable bcc phase , the local minimum is a bct structure . the ebp ( black dot - dashed line ) has two discontinuities , indicated by black dashed lines . brown hatched lines signal axial ratios of the cubic structures . the energy and the density of states were calculated on a uniform @xmath56 grid with a distance of @xmath57 between two adjacent points . ] while the specific energy landscape in the vicinity of the saddle point leads to a discontinuous ebp , the mere presence of a saddle point  be it symmetry dictated or not  does not automatically lead to a discontinuous ebp . for example , if abscissa and ordinate axes in a were interchanged , determining the ebp would mean to search for minima along horizontals , which would result in a continuous ebp including the fcc structure .    in order to analyze the origin of the discontinuities , we consider the electronic density of states ( dos ) at the fermi level , @xmath30 , see and b. the transition at @xmath58 is accompanied by a pronounced change in @xmath30 due to a significant change of the unit cell volume . this change depends monotonously on @xmath5 , see b. further analysis of the total energy shows that the discontinuity at @xmath58 is mainly driven by the kinetic energy contribution . since there is no obvious peculiarity in @xmath30 , the double well has to be assigned to an integral effect of the whole dos . at the position of the second transition @xmath59 , however , there is a distinct peak of @xmath60 in @xmath30 along the ebp ( ) . the same characteristic peak appears as tooth - like feature in the large-@xmath4 region at approximately @xmath61 in b. the ebp of vanadium , which crosses this area , becomes discontinuous in order to avoid the large @xmath30 ( jahn - teller behavior ) . at the local maximum of the related double well of @xmath62 ( not shown ) , @xmath63 . cr*6c element & xc & @xmath16 & @xmath64 & @xmath65 & @xmath66 & @xmath67 & @xmath68 + & & [ @xmath69 & [ @xmath69 & [ @xmath69 & [ @xmath69 & % & [ mev / atom ] + v & pw92 & 2.637(3 ) & 3.97 & 3.39 & 0.58 & 16 & 313 + & pbe96 & 2.698(2 ) & 3.96 & 3.51 & 0.44 & 12 & 280 + v & pw92 & 3.345(5 ) & 2.54 & 2.37 & 0.17 & 7 & 323 + & pbe96 & 3.447(2 ) & 2.59 & 2.46 & 0.13 & 5 & 330 + nb & pw92 & 2.933(2 ) & 4.38 & 3.76 & 0.62 & 15 & 390 + & pbe96 & 2.988(2 ) & 4.41 & 3.87 & 0.54 & 13 & 354 + ru & pw92 & 3.010(1 ) & 3.26 & 3.03 & 0.23 & 7 & 685 + & pbe96 & 3.062(2 ) & 3.31 & 3.05 & 0.26 & 8 & 627 + la & pw92 & 4.043(1 ) & 4.50 & 3.67 & 0.83 & 20 & 178 + & pbe96 & 4.232(2 ) & 4.65 & 3.91 & 0.74 & 17 & 103 + os & pw92 & 3.063(1 ) & 3.28 & 3.02 & 0.26 & 8 & 907 + & pbe96 & 3.104(1 ) & 3.39 & 3.03 & 0.36 & 11 & 900 + ir & pw92 & 3.077(2 ) & 3.23 & 2.89 & 0.24 & 11 & 630 + & pbe96 & 3.118(1 ) & 3.31 & 2.93 & 0.38 & 12 & 818 +    we searched the ebps of the 30 elements listed above for the occurrence of first - order transitions in the same way as described for the case of vanadium  @xcite . the search space was restricted to @xmath70 and had a resolution @xmath71 ( @xmath72 was subsequently reduced to pinpoint @xmath16 ) . seven transitions were found in this way , see table  [ table : discontinuities ] . apart from vanadium , also the group-5 element niobium ( but not tantalum ) and the elements ruthenium , lanthanum , osmium , and iridium exhibit discontinuous ebps . for la and ir , we identified a negative @xmath73 of the bcc structure as the driving force of the instability , while in nb the fcc structure has @xmath74 . bcc ru and os possess small positive @xmath73 but the discontinuities of the ebp occur in the vicinity of the bcc structure at a lattice parameter being slightly smaller and larger than @xmath75 , respectively . no second instability was found for any element except vanadium . all seven first - order transitions were also found in additional calculations relying on pbe96 , see table  [ table : discontinuities ] . compared to the pw92 results , @xmath16 is shifted to larger values in pbe96 . the present work focuses on non - magnetic elements , where we found first - order martensitic transitions under epitaxial strain due to either a large peak in @xmath30 or a more subtle interplay of all contributions to the total energy in the vicinity of a saddle point . in order to find out whether the epitaxial relation of a particular substance is continuous or not , one usually has to calculate this relation in detail . an fcc or bcc structure with negative @xmath73 always results in a discontinuous ebp but a ( small ) positive @xmath73 does not exclude a discontinuity .    in the case of magnetic systems , high - spin to low - spin transitions  @xcite or a change of the ground state magnetic order  @xcite are usually accompanied by discontinuities of the lattice geometry . an example might be the @xmath4-dependence of the magnetic ground state of bct fe found by qiu _ et al . _ these authors did however not discuss the consequence of a transition into a different magnetic state for the lattice parameter @xmath5 of coherent films . we now return to the case of vanadium to demonstrate , how strongly the structural discontinuity may alter another particular intrinsic property , namely the superconducting critical temperature @xmath76 . vanadium is a well - known superconductor and possesses one of the highest @xmath76-values among all elements , @xmath77k at normal pressure  @xcite . a significant change of @xmath76 may be expected at @xmath58 , since @xmath30 jumps by @xmath78 , see . the calculations for @xmath76 of bct vanadium were carried out with the abinit - dft package in the version 6.12  @xcite employing a pseudopotential ( pp ) generated with the fhi98pp code  @xcite for pw92 . with this pp , we found slightly shifted values of @xmath79 and of @xmath80 compared with the fplo results . the phonon properties were obtained through density - functional perturbation theory  @xcite as implemented in abinit  @xcite . to compute the superconducting properties of v , we used a plane - wave kinetic - energy cutoff of 650ev and a @xmath81 mesh in the full brillouin zone for linear tetrahedron integrations . the phonon quantities are sampled by eight wave vectors in the full brillouin zone . we verified the accuracy of our approach for bcc vanadium for which we obtained an electron - phonon coupling constant @xmath82 and a logarithmic - averaged phonon frequency @xmath83 . for comparison , tunneling experiments found @xmath84  @xcite . we derived a reasonable value of @xmath85k for bcc vanadium according to the allen - dynes modified mcmillan formula  @xcite . this equation contains one parameter , the effective coulomb repulsion @xmath86 , which we set to a value @xmath87 commonly used for transition metals  @xcite . @xmath76 across the discontinuity @xmath58 was computed for the points ( @xmath88 coordinates ) ( @xmath89 ) and ( @xmath90 ) . @xmath91 changes by @xmath92 and @xmath76 by @xmath93 or 20% across @xmath58 ( from 3.5k at @xmath94 to 4.1k at @xmath95 ) . the actual value of @xmath76 may depend on the film thickness as this was found in experiments for bcc vanadium films  @xcite , but we expect that the relative change will not substantially be affected by this dependence . we have shown that coherently grown epitaxial films can exhibit a discontinuous relation between in - plane and out - of - plane lattice parameters in a strain - driven first - order phase transition . in particular , we predict that bct bulk - like films of the elements v , nb , ru , la , os , and ir undergo martensitic lattice transformations induced by epitaxial strain with related volume changes up to 20% . as a consequence , the intrinsic film properties can be switched , if one succeeds to control the substrate lattice parameter appropriately . we demonstrated this possibility by the example of vanadium , where the superconducting critical temperature changes by about 20% at a critical strain value . since the mechanisms of the discussed instabilities may be present in arbitrary metallic systems , our results are not limited to monoatomic overlayers , bct films or specific ground states . verification of the predicted transitions should be possible by advanced diffraction methods applied to coherently grown overlayers on a series of substrates with quasi - continuous lattice spacings . an indirect verification could be achieved by measuring any other structure - dependent physical property of the film material . we thank k drr , s gemming , s haindl , s fhler , and a mbius for discussions . 10 url # 1#1urlprefix[2][]#2 ohring m 1992 _ the material science of thin films _ ( san diego : academic press )                                                              gonze x , rignanese g  m , verstraete m , beuken j  m , pouillon y , caracas r , jollet f , torrent m , zerah g , mikami m , ghosez p , veithen m , raty j  y , olevano v , bruneval f , reining l , godby r , onida g , hamann d  r and allan d 2005 _ zeit . kristallogr . _ * 220 * 558562    gonze x , amadon b , anglade p  m , beuken j  m , bottin f , boulanger p , bruneval f , caliste d , caracas r , cte m , deutsch t , genovese l , ghosez p , giantomassi m , goedecker s , hamann d  r , hermet p , jollet f , jomard g , leroux s , mancini m , mazevet s , oliveira m  j  t , onida g , pouillon y , rangel t , rignanese g  m , sangalli d , shaltaf r , torrent m , verstraete m  j , zerah g and zwanziger j  w 2009 _ computer phys . commun . _ * 180 * 25822615
Abstract:

Answer: coherent epitaxial growth allows to produce strained crystalline films with structures which are unstable in the bulk . thereby , the overlayer lattice parameters in the interface plane , @xmath0 , determine the minimum - energy out - of - plane lattice parameter , @xmath1 . we show by means of density - functional total energy calculations that this dependence can be discontinuous and predict related first - order phase transitions in strained tetragonal films of the elements v , nb , ru , la , os , and ir . the abrupt change of @xmath2 can be exploited to switch properties specific to the overlayer material . this is demonstrated for the example of the superconducting critical temperature of a vanadium film which we predict to jump by 20% at a discontinuity of @xmath2 .