Instruction: Write an abstract for the following acadamic paper.

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

Paper: the last decade has witnessed a significant advancement of single molecule manipulation and visualization techniques  @xcite providing access to the distribution of physical properties across many individual molecules and not just average properties as was the case of traditional biochemical techniques . optical tweezers ( ot ) and atomic force microscopy ( afm ) in particular , enable the study of mechanical properties of proteins such as those in the extracellular matrix , in the cytoskeleton and in the muscle , that _ in vivo _ are exposed to stretching forces . the mechanical properties of macromolecules obtained in these experiments may be directly connected to corresponding thermodynamical quantities  @xcite with a bit of caution , since mechanical experiments are usually out - of - equilibrium , as discussed in this paper . on the other hand , the mechanical properties of biomolecules may be of direct importance for what concerns bendability ( dna ) , translocation of nucleic acids and proteins across cellular membranes , rigidity and elasticity ( structural proteins ) .    in a force - measuring afm experiment , the tip of a microscopic cantilever is pressed against a flat , gold - covered substrate , coated with a thin layer of protein . upon retraction of the substrate , a protein molecule that may have been adsorbed to the cantilever tip is then stretched . finally , the force the protein opposes to the stretching is computed from the cantilever deflection , and a force - extension plot is drawn . in a typical stretching experiment the force - extension curve shows a saw - tooth pattern , each peak corresponding to the unfolding of a single domain . however , if the protein is composed by several different modules , it is difficult to associate each peak to the corresponding domain . moreover , in order to neglect the tip - substrate interaction , it is necessary to have a sufficiently long protein . a common strategy proposed in the literature to address these problems is to use an engineered protein composed by several tandem repeats of the same kind of domain . a typical choice is to use domains belonging to the immunoglobulin  @xcite , fibronectin iii  @xcite or cadherin  @xcite superfamilies , characterized by @xmath0-sandwich structures . alpha - helical domains , such as those of the cytoskeletal protein spectrin  @xcite , have also been used . many data about the mechanical properties of modular proteins can be extracted from the force - extension profiles . the amplitude of each peak , in fact , represents the mechanical stability of the corresponding domain , while the spacing between peaks reflects the length of the unfolded domain and is thus proportional to the number of amino - acids it comprises . the experimental data also show a logarithmic dependence of the height of the peaks on the pulling velocity , so that higher forces are required for domain unraveling when the pulling occurs very quickly . the features of the force - extension curves ( figure  [ figvassalli ] ) are accounted for by the entropic elasticity of polymer chains in solution . the entropy of a polymer , in fact , is maximal in the random coil state , whereas it tends to zero in the fully extended conformation . the force required to stretch a polymer thus reflects the entropy loss and it grows in a non - linear way as the molecule is lengthened . the entropic elasticity of biopolymers is currently modeled through the worm - like chain ( wlc )  @xcite and freely - jointed chain ( fjc )  @xcite models . the wlc model describes a molecular chain as a deformable rod whose stiffness is determined by the persistence length @xmath1 ( the length scale over which the polymer loses orientation order ) . the functional relation between the external force and the fractional extension @xmath2 ( @xmath3 is the end - to - end distance and @xmath4 is the contour length ) is approximatly given by the interpolating formula @xmath5 where @xmath6 is the boltzmann constant and @xmath7 is the temperature .    in the fjc model , the polymer is portrayed as a chain of rigid segments linked by frictionless joints so that each segment can point in any direction irrespective of the orientations of the others . a measure of the stiffness of the chain is represented by the kuhn segment length @xmath8 ( the average length of the segments ) . the analytic relation between the average end - to - end distance @xmath9 and the stretching force @xmath10 is @xmath11    in our model , the polymer is described by an array of binary variables representing native contacts that can be in either of two states : formed or broken . the cantilever , on the other hand , is modeled as a harmonic spring in series with the molecule . the energy of the system is the sum of a harmonic term and a term of long - range interaction modeling in a bulk , coarse - grained way , the chemical interactions stabilizing the native conformation of the protein . in fact , rather than providing a detailed description , we account for chemical interactions through a folding prize attributed to a domain when the fraction of intact contacts is above a threshold . we assume any contact breakdown to produce an equal increment in the molecule length . let us consider first the force - extension characteristic of our model in the absence of folding prize . let @xmath10 be a constant pulling force , @xmath12 the number of broken contacts and @xmath13 the length increment per cleaved contact . if no folding prize is attributed to the molecule in a native - like state , the hamiltonian of the system is @xmath14 and the partition function is @xmath15 notice that this is also the partition function of an ising - like model in one dimension without coupling among spins . the average end - to - end distance can thus be computed as @xmath16    , @xmath17 , @xmath18 ; fjc with @xmath19 , ising with @xmath20 @xmath21 , hook ( harmonic potential ) with @xmath22 . ]    the variation of length @xmath23 versus @xmath10 in the wlc , fjc and our ( ising ) models is shown in figure  [ models ] , where we have adjusted the constant @xmath13 ( corresponding to persistence length @xmath1 so to make the curve coincide for small elongations . despite the extreme simplicity of our assumptions , the three curves are qualitatively similar . this similarities could be increased by adding contact - contact interactions or considering different elongations for the various contact breaking events .    although the wlc and fjc models more accurately represent the physics of a polymer , their statistical mechanics treatment is so complex that one generally employes the average formulas and for each domain , complemented with an independence hypothesis of domains and an ad - hoc treatment of the unfolding event  @xcite , based on an extension of bell s expression  @xcite or full kramer s theory  @xcite for the rupture rate coefficient in the presence of a time - dependent external force .    , @xmath24 , @xmath25 , @xmath26 , @xmath27 and @xmath28 ( see section  [ model ] for the illustration of the model ) . ]    the independence assumption is questionable , since all domains are coupled by the presence of the cantilever . this difference may be explicited by comparying computations  @xcite in which the position of the cantilever is observed ( isometric ) while the force may fluctuate , with computations in which the force is maintained constant ( isotensional ) . we can obtain exact comparisons of the two different set - up for our model ( see section  [ model ] ) , as shown in figure  [ iso ] . it can be noticed that the isotensional model does not show any peaks ( the peaks are here smoothed due to the small length of the single module ) .    moreover , the peaks in the force - extension profile ( fig .  [ figvassalli ] ) are a signature of a first - order transition . as we shall show in the following , this transition naturally arise in our model due to the long - range coupling ( the folding prize ) .    despite its extreme simplicity , our approach captures some important aspects of the physics of ig domain stretching . steered molecular dynamics simulations performed in schulten s group  @xcite , in fact , showed that the mechanical unfolding of the i27 module occurs only after the breakdown of a patch of six hydrogen bonds bridging the @xmath29 and @xmath30 @xmath0-strands . the rupture of these critical bonds was shown to be the key event allowing the full unraveling of the molecule under an external force . even if this pattern was originally observed in a specific protein , it could be hypothesized a more widespread distribution . makarov and coworkers  @xcite performed monte carlo simulations of titin forced unfolding . during these simulations the number of hydrogen bonds at time @xmath31 , @xmath32 , undergoes a random walk . it was concluded that a critical value @xmath33 of the number of hydrogen bond does exist , such that when @xmath34 , the domain unfolds very rapidly . makarov also showed that the number of bonds is roughly one , when the force is low enough , whereas for very large pulling rates ( and thus large pulling forces ) , it is likely to be equal to six , recovering the findings by lu and schulten  @xcite . the accuracy of the model developed by makarov and coworkers , allows quantitative comparisons with experimental data at the expense of very long simulation times and the need to assume the knowledge of the free energy profile . their model is also based on the hypothesis of independence of domains which , however , might be incompatible with the coupling introduced by the cantilever . since the transitions shown in afm experiments are out of equilibrium  @xcite , thermal fluctuations may play a fundamental role . our model , conversely , is so simplified that we can compute exactly the free energy of a multi - domain protein for the equilibrium case , and perform long simulations in the out of equilibrium case . the model we propose shows interesting similarities to a g model with force rescaling . recent theoretical studies  @xcite show that the ability of g models to simulate the cooperativity of the folding process can be enhanced by imparting an extra energetic stabilization to the native state so as to simulate specific interactions appearing only after the assembly of native - like structures . a rigorous approach to simulate the stabilizing interactions peculiar of the native state , would be the use of two different analytic expressions of the force - field inside and outside the native basin . the same purpose can be pursued through a much simpler strategy  @xcite , by rescaling the conformational force when the fraction of native contacts crosses a pre - chosen threshold . this approach appears to be similar to the one employed in our model . the energy function of our model comprises an elastic term and a contact term . the harmonic term accounts for the elasticity of the cantilever , while the protein can be thought of as a soft spring so that its contribution to the elasticity of the system is negligible . the contact term of our energy function , on the other hand , is a stabilizing contribution that the protein receives only when the fraction of intact contacts crosses a threshold . this approach is thus equivalent to a force rescaling occurring whenever the polymer enters the native basin , and the contact term of the energy function appears to be a square well . finally , from a purely formal point of view , our description of the polymer as an array of binary variables , bears some similarity with the model proposed by galzitskaya and finkelstein ( gf )  @xcite . in the gf model , in fact , each residue of the polypeptide chain can be either in an ordered ( native ) or disordered ( non - native ) state , encoded by the two possible values of a binary variable . this approach , similarly to ours , significantly narrows the conformational space , that consists of @xmath35 conformations only , for a polymer with @xmath36 residues . this approach , while drastically reducing the computation time , is in agreement with the zimm - bragg model  @xcite , widely employed to describe the helix - coil transition in heteropolymers . the outline of the paper is as follows . in section  [ model ] we describe the model and the simulation procefdure ; in section  [ single ] we study the equilbrium behavio r of a single domain ; in section  [ fluctuations ] we investigate the role of fluctuations in the reciprocal influence between successive unfolding events ; in section  [ pulling - rate ] we investigate the dependence of the unfolding force on the pulling rate and we sketch a simple analytic treatment of the unfolding probability in the limit of an extremely high pulling rate ; finally , in section  [ conclusions ] we draw the conclusions of our work . the afm in its simplest arrangement is just a spring ( the cantilever ) whose deflection , and thus the applied force , is measured as a function of its position . the system , as shown in figure  [ setup ] , can therefore be modeled as a harmonic spring in series with a protein composed by @xmath37 tandem repeats of the same domain . each domain @xmath38 is simply portrayed as a sequence of contacts ( @xmath39 , @xmath40 ) that can be either intact ( @xmath41 ) or broken ( @xmath42 ) , where @xmath36 is the total number of native contacts inside each domain . the length of the polymer chain can be simply computed as @xmath43 , where @xmath13 represents the incremental elongation associated to each contact breakdown and @xmath44 is the number of broken contacts in domain @xmath38 . for the sake of simplicity , in all our computations we set @xmath45 . the length of the spring , on the other hand , is just @xmath46 where @xmath47 is the extension of the spring - polymer system , using the rest position as the reference point . the energy of a domain configuration @xmath48 is modeled by the sum of two contributions : a harmonic term @xmath49 , accounting for the presence of the spring , and the sum over all domains @xmath38 of a term @xmath50 , related to the long - range interactions among the monomers , @xmath51    the harmonic term is expressed as @xmath52 where @xmath53 is the harmonic constant .    in our simplified representation , if the fraction of intact contacts @xmath54 in a domain @xmath38 is below a given threshold @xmath55 , the domain receives a folding prize @xmath56 proportional to the number of possible contacts . the interaction term is thus computed as @xmath57 where @xmath58 is the heaviside step function @xmath59    in summary , the energy of a configuration is just a function of the number of broken contacts in each domain @xmath60^{2 }     -an\theta(\theta n - n_{j}),\ ] ] where @xmath61 and @xmath62 this assumption speeds up the computations , that may be performed in terms of the @xmath63 .    a stretching simulation starts from a completely folded initial structure , where no contact is broken . the protein - spring length @xmath47 , chosen as the control parameter , is linearly increased from @xmath64 to @xmath65 in @xmath66 steps during the simulation ; @xmath67 where @xmath68 . for each value of @xmath47 , we compute the average length of the spring @xmath69 in an equilibrium simulation as @xmath70 where @xmath71 is the partition function @xmath72 and the multeplicity factor @xmath73 is given in terms of the number of possible microscopic configurations containing @xmath63 cleaved contacts .    in real stretching experiments , however , the polymer is subjected to a finite pulling velocity , so that the molecule can not be considered in an equilibrium condition . in order to consider this effect , we use monte carlo simulations . for each value of @xmath47 , @xmath7 monte carlo steps are performed , each involving @xmath74 elementary steps . the elementary step consists in the random choice of a domain and in the attempt to increase or decrease by one the number of contacts with probability equal to the fraction of broken or intact contacts respectively . the trial move is then accepted or rejected with a probability derived from the heat - bath criterion , @xmath75 where @xmath0 is the inverse temperature . the average length of the spring @xmath69 , corresponding to the current position @xmath47 is computed averaging over @xmath7 monte carlo steps . we first analyze the entropic effects related to temperature through the analysis of computations in the absence of a folding prize , and then investigate the role of long - range interaction by setting a non - zero prize on a single - domain polymer . after that , we discuss a set of computations on a three - domain protein , showing the importance of the coupling due to the cantilever in the mechanism of mechanical unfolding and , in particular , they explain how the first unfolding event affects the following ones . in the last part , we discuss the relation between pulling rate and unfolding force , finding a logarithmic law . the section is completed with a analytic treatment of the unfolding probability valid in the limit of high pulling rate . .legend of the symbols appearing in the figures . [ cols="<,<",options="header " , ]     the features of the saw - tooth @xmath76 profile are affected by several simulation parameters . the folding prize @xmath77 is related to the enthalpic component of the free energy and stabilizes the folded conformation of the protein . as a consequence , when @xmath77 is large , the polymer tends to remain in the compact conformation so as to retain the significant folding prize and the increase in @xmath47 leads to a stretching of the cantilever spring . only when @xmath69 is very large it becomes enthalpically favourable for the polymer to unfold because the decrease in harmonic energy due to the cantilever relaxation more than compensates the loss of the folding prize . thus , as @xmath77 is increased , higher and higher values of @xmath69 are required for the elastic energy to compensate the folding prize and the peak of the @xmath76 plot becomes accordingly higher and higher . the harmonic constant @xmath53 of the cantilever spring plays a role basically opposite to that of the folding prize . in fact , when @xmath53 is large , smaller average extensions @xmath69 are required for the harmonic energy to balance the folding prize so that larger @xmath53s result in a less pronounced peak in the @xmath76 plot . the role of @xmath77 and @xmath53 is exemplified by the simulations portrayed in figure  [ ak ] .    . other computation parameters : @xmath78 ; @xmath79 ; @xmath80 ; @xmath81 ; @xmath82 ; @xmath83 ; @xmath84 . ]    as just discussed , the role of the folding prize and of the harmonic constant is related to the enthalpic term of the free energy . by contrast , the folding threshold @xmath55 influences the entropic contribution to the free energy . a small value of @xmath55 in fact , implies that only a small fraction of the contacts can be broken without loss of the folding prize . if @xmath85 represents the number of broken contacts in a moment preceding the unfolding event , then the number of microscopic conformations allowed will be @xmath86 , and the entropy will be @xmath87 $ ] . if @xmath88 , the unraveling of the domain will increase the degeneracy @xmath89 and thus the entropy , so that a moderate extension of the spring will be sufficient for the enthalpic term to be more than compensated by the entropic one . by contrast , when @xmath90 , the breakdown of the domain and the resulting increase in the number of disrupted contacts , will bring the degeneracy and the entropy further away from the maximum thus disfavouring the unfolding event and causing the peak of the @xmath76 profile to become higher . the pattern of increase in height of the unfolding peak as @xmath55 takes on higher values is shown in figure  [ theta ] .     and the inverse temperature @xmath0 . for small @xmath55s ( a ) higher temperatures favour the unfolding , whereas the opposite is true for large @xmath55s ( b ) . common computation parameters : @xmath78 ; @xmath79 ; @xmath80 ; @xmath81 ; @xmath82 ; k = 0.1 . ] the relevance of the entropic contribution on free energy computation strongly depends on temperature that may amplify the role of the threshold @xmath55 . as already noticed , in fact , when @xmath88 , the breakdown of the domain is entropically favoured as it brings the degeneracy closer to its maximum . since this gain in entropy becomes larger and larger as the temperature is increased , for small @xmath55 the unfolding event becomes more and more favourable as @xmath0 is decreased , and the height of the peak of the @xmath76 plot will also decrease accordingly . for @xmath90 the opposite pattern can be observed . in fact , the entropy loss due to the decrease of the degeneracy resulting from the unfolding event , is magnified as @xmath0 is decreased leading to higher and higher peaks in the @xmath76 plot . the interplay between the parameters @xmath55 and @xmath0 is shown in figure  [ thetabeta ] .      ) causes the first unfolding event to favour the following ones . panels ( c ) and ( d ) : the first unfolding event does not affect the second and third ones as a result of the high threshold @xmath91 . computation parameters : @xmath92 ; @xmath93 ; @xmath94 ; @xmath95 ; @xmath96 ; @xmath97 ; @xmath98 . ]    . simulation parameters : @xmath92 ; @xmath93 ; @xmath94 ; @xmath95 ; @xmath99 ; @xmath100 ; @xmath97 ; @xmath84 . ]    before studying the mutual influence of the unfolding events , let us illustrate the features of the free energy landscape of a single module near the unfolding transition . in figure  [ freeenergy ] we show the evolution of the free energy landscape in correspondence of the unfolding transition in a simulation with @xmath101 and @xmath100 ( see also top left panel of fig . [ fig : fluctuations ] ) . the profile is characterized by a cusp - like , narrow well corresponding to the folded state , and a wide , smooth well related to the unfolded state . the width of the two wells , in fact , depends on the number of conformations that the system can explore : in the folded state , the system is stretched due to the action of the spring and no fluctuations are allowed so that only one conformation will be populated ; after the transition , the residues of the collapsed domain become free to fluctuate and many conformations will be explored thus determining a very wide well in the free energy profile . the figure shows that for low values of @xmath47 ( and thus low values of the elastic force ) , the free energy of the reference folded state @xmath102 is lower than that of the unfolded state @xmath30 , thus forbidding the breakdown of the domain ; as @xmath47 is increased , the free energy of the folded state increases , until it finally becomes higher than that of the folded conformation and the stretching transition occurs . this is a typical example of a first - order phase transition . we now consider a polymer composed by @xmath93 tandem repeats with @xmath78 contacts each . in figure [ fig : fluctuations ] we compare two simulations performed with the same parameter set except for a different threshold @xmath55 . in the simulation with @xmath100 the peaks corresponding to the second and third unfolding events are less pronounced than the first one thus suggesting that the unfolding of a domain actually favours further unfolding events . conversely , in the simulation with @xmath91 , the three peaks feature almost the same height showing that the first unfolding event has little or no influence at all on the following ones . each unfolding event in figure  [ fig : fluctuations ] corresponds to a peak in the variance plot because the unraveling of a domain increases the fluctuations of the polymer and spring length . the variance peaks are characterized by a high and narrow spike followed by a smoother region that decreases more slowly . the shape of the variance peak is related to the regions of the free energy landscape explored by the system during the unfolding transition . the spikes in the variance plots ( that are truncated for the sake of graphical clarity ) correspond to the situation with @xmath103 when both wells are explored by the system and the variance @xmath104 is related to the distance between the two wells . on the other hand , the smooth regions of the variance plots refer to the case with @xmath105 when the system only explores the unfolded region of the free energy landscape whose width correlates with the variance . figure  [ fig : fluctuations ] shows that the height of the smooth region of the variance peaks increases with the order of the unfolding event . this trend is due to the fact that , with each unfolding event , @xmath106 new monomers become free to fluctuate and the number of accessible conformations increases accordingly . figure  [ fig : fluctuations ] also shows that the value of the variance @xmath104 after each unfolding event , gradually decrease as @xmath47 is increased , because the disruption of the contacts of the domain just collapsed allow an extension of the molecule in the direction of the pulling force so as to avoid as far as possible a further stretching of the spring that would cause an increase in energy . as a result , narrower and narrower regions of the conformation space become accessible to the polymer and the variance is lowered . it is worthwhile noticing , however , that the extension of the unfolded domain is hindered by the subsequent decrease in entropy so that the number of contacts actually broken before each unfolding event is smaller than the maximum number allowed by the loss of the folding prize in the previous domain breakdown .    when @xmath100 , 26 contacts ( a number of the order of @xmath107 ) are broken before the first unfolding event . this value is consistent with the number of contacts that can be disrupted without loss of the folding prize . after the first collapse event , the number of contacts broken in the simulation rises to 115 , thus determining an increase of the fluctuations and favouring the breakdown of another domain . this explains why the second peak of the @xmath76 plot is less pronounced than the first one . the occurrence of the second unfolding determines a further increase of the number of cleaved contacts to 205 . this results in a easy breakdown of the last domain of the protein and the last peak of the @xmath76 plot is therefore less high than the second one . this scenario is significantly different for @xmath91 . for large values of @xmath55 , in fact , only a small number of residues can be recruited for fluctuations after each unfolding event . as a consequence , the variance @xmath104 rapidly goes to zero after each unfolding event and the fluctuations are extinguished before the force threshold for unfolding can be crossed . the following unfolding events , similarly to the first one , will occur when the increase in harmonic energy balances the loss of the folding prize and therefore the height of the unfolding @xmath69-peaks will be roughly the same . , if the temperature is sufficiently high ( @xmath84 , panels ( a ) and ( b ) ) , the first unfolding event favours the following ones due to the role of fluctuations . at low temperatures however ( @xmath108 , panels ( c ) and ( d ) ) , the fluctuations rapidly become negligible and the heights of the peaks become roughly the same due to the identical energetic features of the domains . computation parameters : @xmath92 ; @xmath93 ; @xmath94 ; @xmath95 ; @xmath96 ; @xmath100 ; @xmath109 . ]    the role of fluctuations in determining the relative heights of the peaks of the @xmath76 plot , and thus the coupling among domains , is confirmed through simulations performed at different temperatures . at very low temperatures , in fact , the entropic term of the free energy becomes negligible and the height of the peaks of the @xmath76 plot only depends on the energetic term . as we are considering a protein composed by identical domains , we expect the @xmath110 peaks to be identical if the temperature is sufficiently low . this pattern can be observed in figure  [ temperature ] where we compare two simulations performed at different temperatures , namely @xmath84 and @xmath108 .      ; @xmath111 ; @xmath112 ; @xmath81 ; @xmath113 ; @xmath96 ; @xmath114 ; @xmath115 ; @xmath84 . ]    up to now , we discussed equilibrium stretching computations , _ i.e. _ we assumed an infinitely slow pulling . however , at the typical time scales of an afm stretching experiment , the polymer is pulled at a finite velocity and the unfolding is a non - equilibrium process , as testified by the differences between the unfolding and the refolding force - extension profiles ( see e.g.figure  [ montecarlo ] ) . in fact , while the unfolding profile features the typical saw - tooth pattern , upon relaxation of the unfolded polymer , the trace exhibits no discontinuities that would indicate refolding . a consequence of the irreversibility of the stretching process is that the unfolding force depends on the pulling speed . actually , if the loading rate @xmath116 ( with @xmath117 being the elastic constant and @xmath118 the pulling velocity ) is sufficiently small , thermal fluctuations are allowed enough time to overcome the energy barrier and the unfolding force will be low . several experimental works  @xcite reported a logarithmic dependence of the unfolding force on the loading rate in the case where a single energy barrier is present along the reaction path . the analytic expression of the relation between force and loading rate was derived  @xcite within the frame of kramer s theory for a simple two - state model , by assuming that the external force reduces the height of the energy barrier , @xmath119 where @xmath120 is boltzmann constant , @xmath7 is the absolute temperature , @xmath121 is the distance between the minimum corresponding to the folded state and the activation barrier of the energy landscape , @xmath118 is the pulling speed , @xmath53 is the cantilever harmonic constant and @xmath122 is the spontaneous unfolding rate . recently  @xcite , it has been shown that the probability distribution of the force at various pulling rates does not follow the simple bell law , requiring the full kramer s theory and predicting small corrections to the logaritmic behavior of the most probable force . profile become lower and lower and the mc simulations converge to the equilibrium scenario as the number of mc steps is increased . simulation parameters : @xmath26 ; @xmath93 ; @xmath123 ; @xmath81 ; @xmath96 ; @xmath124 ; @xmath125 ; @xmath98 . notice that , using a smaller value of @xmath36 , the equilibrium profile is smoother than that in figure  [ fig : fluctuations ] , figure  [ temperature ] and preceding ones . ] we limit here to a preliminary illustration of a series of monte carlo simulations with different number of steps @xmath7 considering the pulling speed as being proportional to @xmath126 . a more detailed analysis will be presented in a forthcoming paper . figure  [ dynamics ] shows that , as the number of monte - carlo ( mc ) steps is increased , the @xmath127 plot becomes closer and closer to the profile yielded by the equilibrium simulation . in particular , a small number of mc steps , _ i.e. _ a fast pulling , corresponds to high peaks of the @xmath127 profile , whereas larger numbers of steps correspond to less and less pronounced peaks and thus , smaller unfolding forces . we performed a series of 100 independent mc runs for 5 different values of @xmath7 : @xmath128 , @xmath129 , @xmath130 , @xmath131 and @xmath132 . for each series of runs we built the histograms of the rupture spring elongations and we plotted the mean value of the histogram as a function of the logarithm of the loading rate . finally , the set of points thus obtained was linearly fitted . as shown in figure  [ statistics ] , the histograms , shift to lower values as the number of mc steps increases . figure  [ statistics ] also shows that the mean rupture forces computed from the histograms , feature an excellent linear correlation with the logarithm of the loading rate ( correlation coefficient @xmath133 ) .    if eq . [ eq : evans ] is explicitly rearranged so as to show the linear dependence of @xmath134 on @xmath135 , the coefficients of the linear equation can be equated to the corresponding parameters @xmath136 and @xmath137 of the regression line @xmath138 , so as to build the following system of equations : @xmath139    from the first equation of the system it is possible to compute the width of the activation barrier @xmath140 ; this value can then be substituted into the second equation so as to determine @xmath141 . the spontaneous unfolding rate @xmath122 , on turn , is related to the height of the activation barrier for the unfolding process , @xmath142 where @xmath143 , as explained by kramer s theory , is the reciprocal of a diffusive relaxation time . this simple computation shows how stretching experiments and simulations provide easy access to important features of the free energy landscape . in a forthcoming paper we aim at investigating the relations relating the barrier width @xmath121 and the spontaneous unfolding rate @xmath122 with molecular properties such as the folding prize @xmath77 , the threshold @xmath55 , the number @xmath36 and the length @xmath37 of the domains . let us investigate the dependence of the monte carlo simulations on the number of monte - carlo steps @xmath7 , that we can interpret as a measure of the pulling speed . in fact , in the limit of an extremely high pulling speed , we can assume that for each value of @xmath47 the polymer can adopt just a single ( or at most , a few ) conformation , so that the entropic contribution can be neglected in the computations . the mechanism outlined in section  [ typical - plot ] shows that the key event is the loss of the folding prize determining the complete extension of the protein domain and the subsequent relaxation of the spring . the problem thus arises to identify the factors affecting the probability of this crucial step . more specifically , suppose in a domain a number of contacts just below the threshold to keep the folding prize has been broken . what is the probability @xmath144 that one more contact will be cleaved within the next @xmath145 steps ?    in order to answer this question , it is necessary to compute the energy difference associated with the transition . let @xmath146 be the length of the spring when a number of contacts just below the threshold @xmath147 has been disrupted : the system retains the prize @xmath77 . if one more contact is broken , the prize is lost and the spring will be correspondingly shortened . in particular , by assuming each contact breakdown to cause a unit increase in length of the polymer , the new length of the spring will be @xmath148 . the energy difference between the two states can thus be computed as @xmath149 where a term @xmath150 has been neglected . we can now compute the probability to destroy a contact in a single step @xmath151    finally , the probability to break a contact in @xmath145 steps is given by the sum of a geometric progression @xmath152    the computations thus show that the domain unfolds with a probability depending on the prize , the cantilever stiffness , the temperature and the pulling velocity ( that in our simulations is related to the number @xmath7 of monte carlo steps ) .    ; @xmath79 ; @xmath153 ; @xmath154 ; @xmath155 ; @xmath156 ; @xmath124 ; @xmath125 ; @xmath84 . ]    in order to test our analytical treatment of the unfolding probability , we performed a series of mc simulations differing from a reference one just for a parameter ( figure  [ height ] ) . the key parameters characterizing the reference simulations are : @xmath155 , @xmath157 , @xmath125 , @xmath84 . in agreement with the analytic computation , the simulations show that an increase in the number of monte carlo steps ( @xmath158 ) , an increase in temperature ( @xmath159 ) , a decrease of the folding prize ( @xmath160 ) , and an increase of the spring constant ( @xmath161 ) , all result in an increase of the unfolding probability , thus reducing the height of the peak in the @xmath127-plot . we developed an extremely simplified model of polymer stretching in which the molecule is portrayed as a series of modules , represented as an array of contacts , and a harmonic spring ( the cantilever ) . the chemical interactions stabilizing the native conformation are simply modeled as a folding prize gained by domains where the fraction of folded monomers is above a pre - chosen threshold . our model is consistent with recent findings in refs  @xcite , showing that the unraveling of the titin immunoglobulin domain occurs very rapidly only after the breakdown of a critical number of key hydrogen - bonds . the attribution of a folding prize when a threshold value of the fraction of native contacts is crossed , also makes our approach equivalent to a g - model with force rescaling  @xcite . however , our model is thus significantly simpler than other models commonly used to study mechanical unfolding such as the wlc and fjc models and the detailed all - atom , topological and united - residue models employed for steered molecular dynamics simulations . yet , our model is detailed enough to reproduce many qualitative features of the force - extension profiles recorded in afm experiments . in particular , our model correctly reproduces the typical saw - tooth pattern with peaks characterized by a height dependent on the folding prize , the temperature and the pulling velocity .    in our study , a particular attention was paid to the relation between the heights of successive peaks in the force - extension plot . this point is quite intriguing as experimental data show that the force of unfolding tends to increase with each unfolding event , suggesting that the protein domains unfold following an increasing order of mechanical stability  @xcite . a possible explanation suggested by the literature is that the domains of the engineered modular constructs used in afm experiments are not identical but just structurally similar  @xcite . this argument may be correct in the case of real proteins where the differences in mechanical stability of the tandem repeats of the construct may arise from their different position in the tridimensional structure of the molecule . in a computer simulation , however , the domains are all perfectly identical and the explanation for the staircase pattern of unfolding peaks must be sought elsewhere . this behavior has been ascribed to an independent breaking probability of each domain . this hypothesis is however not consistent with the fact that the cantilever couples the fluctuations of all domains . our unified treatment allows to keep into consideration all contributions , and to obtain the correct equilibrium curves , that generally do not exhibit this effect . however , a detailed study of out - of - equilibrium extension curves , shows that this behavior is consistent with a finite pulling speed , i.e.  it can be ascribed to a dynamical source . our results appear to be in agreement with recent findings by cieplak  @xcite in steered molecular dynamics simulations of titin and calmodulin unfolding using a go - like model . in particular , it was found that an increase in thermal fluctuations results in a lowering of the force peaks and in their earlier occurrence during stretching . an interesting feature of afm stretching experiments is that the mean unfolding force is a logarithmic function of the pulling rate . this pattern reflects the role of thermal fluctuations : if the pulling speed is sufficiently low , then there will be enough time for thermal fluctuations to drive the polymer over the free energy barrier and the unfolding force will be accordingly low . the ability to reproduce this pattern constitutes a benchmark for any stretching model . in order to test our model , we performed a series of monte carlo simulations with different numbers @xmath7 of mc - steps and regarding the pulling velocity to be proportional to @xmath126 . our toy - model , despite its extreme simplicity , correctly reproduced the linear dependence of the mean unfolding force on the logarithm of the pulling rate .    in summary , our simple model including harmonic and long - range energy contributions , is not only capable of reproducing the force - extension saw - tooth pattern , but it also yields force spectra displaying the correct logarithmic dependence of @xmath134 on @xmath162 reported in experimental works . finally , our model provides a simple explanation of the influence of each unfolding event on the following ones , based on the role of fluctuations . our work thus shows how minimal models can be valuable tools even in the study of complex molecular systems . we gratefully acknowledge fruitful conversations with dr . l. casetti and dr . m. vassalli . we are also indebted with dr . vassalli for having provided the experimental plot shown in figure  [ figvassalli ] . this work is part of the ec project embio ( ec contract n. 012835 ) . rounsevell r.w.s . , steward a. , geierhaas c.d . , williams p.m. , paci e. and clarke j. ( 2005 ) , mechanical unfolding of tnfn3 : the unfolding pathway of a fniii domain probed by protein engineering , afm and md simulation , j. mol . biol . , 350 , 776 - 789 .    clementi c. , nymeyer h. and onuchic j.n . ( 2000 ) , topological and energetic factors : what determines the structural details of the transition state ensemble and `` en - route '' intermediates for protein folding ? an investigation for small globular proteins j. mol . biol . , 298 , 937 - 953 .
Abstract:

Answer: we present an extremely simplified model of multiple - domains polymer stretching in an atomic force microscopy experiment . we portray each module as a binary set of contacts and decompose the system energy into a harmonic term ( the cantilever ) and long - range interactions terms inside each domain . exact equilibrium computations and monte carlo simulations qualitatively reproduce the experimental saw - tooth pattern of force - extension profiles , corresponding ( in our model ) to first - order phase transitions . we study the influence of the coupling induced by the cantilever and the pulling speed on the relative heights of the force peaks . the results suggest that the increasing height of the critical force for subsequent unfolding events is an out - of - equilibrium effect due to a finite pulling speed . the dependence of the average unfolding force on the pulling speed is shown to reproduce the experimental logarithmic law .


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: after the completion of the operations for lhc run - i , the most important legacy for the whole of particle physics is the existence of a higgs boson with a mass of about 126 gev  @xcite . furthermore , as a result of combined data from various higgs search channels , it was established that the observed higgs boson has properties which are consistent with those of the standard model ( sm ) higgs state  @xcite .    however , it is natural to be left with a key question . namely , what is the actual number of scalar doublet fields present in nature ? although the higgs sector in the sm is constructed of only one isospin doublet field , that is not supported by a fundamental principle . in fact , it is well known that multi - doublet models can be considered so that the higgs boson phenomenology is consistent with current lhc data by arranging additional parameters suitably , and they are often introduced new physics models beyond the sm . for example , supersymmetry - extended models require at least two doublets in their higgs sector for gauge anomaly cancellation and to construct the yukawa lagrangian  @xcite . a second doublet field is introduced also in many radiative neutrino mass models  @xcite . the inert doublet model ( idm )  @xcite , aka a 2-higgs doublet model ( 2hdm ) with one active and one inert doublet , has been proposed to explain the existence of dark matter ( dm ) , by depriving a vacuum expectation value ( vev ) of the second doublet field and making it non - interacting with sm fields except for the ( thereby active ) higgs boson itself and the gauge states symmetry into the higgs sector , where all the sm particles [ the second doublet ] are assigned to be even [ odd ] under the @xmath5 symmetry . the particles with the odd parity are produced in pairs and the lightest one amongst these is stable , hence a dm candidate . ] . furthermore , an additional source of cp - violation is obtained in a higgs sector with multiple doublets which is required to realize a successful scenario for electro - weak ( ew ) baryogenesis  @xcite . properties of an enlarged higgs sector clearly depend on the new physics scenario embedding it , so that the determination of the true higgs boson dynamics is extremely important to extract information on the new physics model . in particular , if a multi - doublet structure is chosen for the higgs sector , what is a mass scale of the additional active and/or inert ( pseudo)scalar bosons ? how can we extract information on the total number of doublet fields ? needless to say , the direct way to probe a multi - doublet structure is discovering extra ( pseudo)scalar bosons at collider experiments . for example , even in the simplest case of 2hdms , there appear a pair of singly - charged ( @xmath6 ) , a cp - odd ( @xmath7 ) and a cp - even ( @xmath8 ) higgs boson in addition to the sm - like higgs boson ( @xmath9 ) . clearly , the most important property to find these extra scalar bosons is their masses . although , in general , these masses are taken as free parameters , we may be able to constrain their corresponding parameter space by taking into account theoretical constraints : chiefly , perturbative unitarity and vacuum stability . in 2hdms , the unitarity and vacuum stability bounds have been discussed in refs . @xcite and in ref . @xcite , respectively . by using them , dimensionless parameters in the scalar potential can be constrained , and it can be translated into limits on the mass parameters for the extra ( pseudo)scalar bosons . in particular , when there is a non - zero mixing between @xmath9 and @xmath8 in 2hdms , the exact decoupling limit of the extra scalar boson can no longer be taken , because of the unitarity constraints . therefore , one can obtain an upper limit on the mass of the extra scalar boson  @xcite . the argument can be generalized to multi - higgs doublet scenarios . the effects of a non - zero mixing can appear as deviations in the couplings of @xmath9 with gauge bosons ( @xmath10 ) and/or fermions ( @xmath11 ) from the sm predictions . hence , the very fact that the latter can be accessed experimentally enables one to extract information on the underlying ( pseudo)scalar dynamics above and beyond the mass scale directly probed at the collider concerned . take the @xmath10 couplings for example . they are presently ( i.e. , after lhc run - i ) constrained to be sm - like at the level of 10% or so . lhc run - ii ( 13 tev ) at standard integrated luminosity ( 300 fb@xmath12 ) would improve this precision by about a factor of 2 . their expected accuracy will then be about 2%  @xcite at the high - luminosity ( hl ) lhc  @xcite ( with collision energy of 14 tev and integrated luminosity of 3000 fb@xmath12 ) . furthermore , the accuracy is expected to be about 0.4%  @xcite at the international linear collider ( ilc ) with collision energy of 500 gev and integrated luminosity of 500 fb@xmath12 ( ilc500 ) . therefore , if a deviation in @xmath10 couplings is discovered at future colliders , whenever this will occur , it will at the same time suggest the possibility of an enlarged higgs sector and , if so , enable one to place an upper limit on the mass of extra ( pseudo)scalar bosons .    in this paper , we investigate constraints on dimensionless parameters in the scalar potential from the bounds imposed by perturbative unitarity and vacuum stability in 3-higgs doublet models ( 3hdms ) . the phenomenological importance of such extended higgs scenarios , amongst many possible others , has been emphasized lately in the context of the latest lhc results and dm searches . for example , 3hdms may shed light on the flavour problem , namely the unknown origin and nature of the three families of quarks and leptons , including neutrinos , and their pattern of masses , mixings and cp violation : it is possible that the three families of quarks and leptons could be described by the same symmetries that describe the three higgs doublets @xcite . such a family symmetry could be spontaneously broken along with the ew symmetry , although some remnant subgroup could survive , thereby stabilizing a possible scalar dm candidate while for certain symmetries it is possible to find a vev alignment that respects the original symmetry of the potential which will then be responsible for the stabilization of the dm candidate @xcite . although there are several types of 3hdms depending on symmetries imposed in the scalar potential @xcite@xcite , having in view the aforementioned phenomenological requirements and a need for a simple and predictive formulation , recently , @xmath5 symmetric 3hdms with one or two inert doublets have been investigated and found viable from both a dm and lhc perspective @xcite . furthermore , other than explaining current data , new regions of dm relic density and areas of parameter space yielding totally new lhc signals ( with respect to , e.g. , 2hdms ) open up @xcite .    herein , we focus on the i(1 + 2)hdm and i(2 + 1)hdm , where the former[latter ] model contains two[one ] active doublets and one[two ] inert doublets . we calculate the @xmath0-wave amplitude matrix for all possible 2to2 body scalar boson elastic scatterings in the high energy limit in both models , including contributions from the sm gauge bosons as the nambu - goldstone ( ng ) boson states due to the equivalence theorem approximation ( eta ) symmetry to obtain eigenvalues for the @xmath0 wave amplitude matrix . because the symmetry imposed in the higgs potential in our paper is different from that in ref . @xcite , derived eigenvalues are also different . we have confirmed that our formula for the eigenvalues coincides with those presented in ref . @xcite by taking appropriate replacements of the coupling constant in the potential . detailed explanations are given in the end of sec . iii - b . ] . we then examine the allowed parameter space from the unitarity and vacuum stability ( re - derived herein and found to be compliant with those in refs . @xcite ) bounds in addition to the constraints from the @xmath2 , @xmath3 and @xmath4 parameters . then , as intimated , we proceed to apply the unitarity bounds to obtain an upper limit on the masses of @xmath7 and @xmath8 in the i(1 + 2)hdm in the case with a non - zero mixing between @xmath9 and @xmath8 , and study how the limit can be changed as compared to that in 2hdms . we find that a stronger upper bound is obtained in the i(1 + 2)hdm than that in 2hdms depending on the masses of the inert ( pseudo)scalar bosons . therefore , we can exclude the i(1 + 2)hdm when the masses of @xmath7 or @xmath8 will be found to be larger than the given upper limit for a fixed amount of the mixing between @xmath9 and @xmath8 . as for the i(2 + 1)hdm , since unitarity here is perfectly realized by the sm - like higgs state @xmath9 , we study what constraints are induced on the various parameters in the potential involving inert states ) in order not to spoil the perturbative behavior of the only active scalar in this scenario . this paper is organized as follows . in sec . ii , we give the scalar potential and mass formulae for the scalar bosons . in sec . iii , constraints from vacuum stability , perturbative unitarity and the @xmath2 , @xmath3 and @xmath4 parameters are discussed . in sec . iv , we delineate the allowed parameter regions from the constraints given in the preceding section . finally , conclusions are given in sec . v.    .charge assignment under the @xmath13 symmetry in the i(1 + 2)hdm and i(2 + 1)hdm . [ cols="^,^,^,^,^,^,^ " , ] among various classes of 3hdms  @xcite , we discuss the one which satisfies the following requirements :    \(i ) it contains a sm - like higgs boson ;    \(ii ) it contains a unique candidate of dm in a given mass spectrum ;    \(iii ) it does not have flavour - changing neutral current ( fcncs ) at the tree level .    in order to satisfy the requirements of ( i ) and ( ii ) , we need at least one active and one inert doublet field . thus , the ( pseudo)scalar sector of our 3hdms must be composed of one inert plus two active doublet fields ; i.e. , i(1 + 2)hdm or two inert plus one active doublet fields ; i.e. , i(2 + 1)hdm . to guarantee a doublet field being inert and the stability of a dm candidate , we require the corresponding vev to be zero and impose an unbroken discrete @xmath5 symmetry onto the whole doublet sector . a dm candidate can then be obtained as the lightest neutral ( pseudo)scalar component of the inert doublet field with the odd parity under @xmath5 ( the sm plus the higgs / active ( pseudo)scalar states being even ) . as a consequence , we also obtain that the inert ( pseudo)scalar states can only appear in pairs in each vertex in the lagrangian , so that they can not couple to fermions directly . in the i(1 + 2)hdm , fcnc processes via neutral scalar boson mediation appear at the tree level similarly to the generic 2hdm case , because of the existence of two doublet fields with the same quantum numbers . the simplest way to avoid such fcncs is imposing another discrete @xmath14 symmetry  @xcite which can be softly - broken in general ( henceforth denoted as @xmath15 ) . here , we consider both the i(1 + 2)hm and i(2 + 1)hdm with the two discrete symmetries @xmath5 and @xmath15 symmetry . however , by imposing the @xmath15 symmetry throughout , the structure of the higgs potential can be much simpler than that in the case without it , thereby rendering possible a semi - analytical treatment of unitarity violating processes . in this connection , it turns out that we can deal with the two models uniformly as seen in the scalar potential ( given below in eq . ( [ pot ] ) ) . we thus consider the @xmath13 symmetric potential in both 3hdms . ] .    to write down the scalar potential in the 3hdms , let @xmath16 @xmath17 be isosipin doublet scalar fields , where the assignment of @xmath13 charge is listed in table  [ particle ] . the most general scalar potential is then given in the same form in both the i(1 + 2)hdm and i(2 + 1)hdm by @xmath18 \notag\\ & + \rho_1(\varphi_1^\dagger\varphi_1)(\varphi_0^\dagger\varphi_0 )   + \rho_2|\varphi_1^\dagger \varphi_0|^2 + \frac{1}{2}[\rho_3(\varphi_1^\dagger \varphi_0)^2 + \text{h.c . } ] \notag\\ & + \sigma_1(\varphi_2^\dagger\varphi_2)(\varphi_0^\dagger\varphi_0 )   + \sigma_2|\varphi_2^\dagger \varphi_0|^2 + \frac{1}{2}[\sigma_3(\varphi_2^\dagger \varphi_0)^2 + \text{h.c . } ] , \label{pot}\end{aligned}\ ] ] where @xmath19 , @xmath20 , @xmath21 and @xmath22 are complex parameters in general . throughout the paper , we take these parameters to be real for simplicity , thereby avoiding explicit cp violation . we also neglect the possibility of complex vevs , thus also removing spontaneous cp violation . hence , the ( pseudo)scalar fields can be parameterised as @xmath23,~~(i=0, ... 2 ) , \end{aligned}\ ] ] where @xmath24 are the vevs of @xmath16 with the sum rule @xmath25 gev@xmath26 . in the i(2 + 1)hdm , @xmath27 and @xmath28 while in the i(1 + 2)hdm , @xmath29 and @xmath30 . the yukawa lagrangian in the i(2 + 1)hdm is given in the same form as that in the sm , where @xmath31 gives all the sm fermion masses . in contrast , in the i(1 + 2)hdm , the structure of the yukawa lagrangian depends on the @xmath15 charge assignments for the sm fermions . in general , there are four independent choices of the assignment  @xcite and these are the so - called type - i , type - ii , type - x and type - y  @xcite . various constraints from flavour @xcite and collider physics @xcite strongly depend on this choice . however , the following discussion does not , though we will make sure to avoid phenomenologically dangerous regions of parameter space .      in the i(1 + 2)hdm , the inert sector is composed of one doublet . therefore , the cp - even , cp - odd and pair of charged scalar states in the doublet field @xmath31 correspond to the mass eigenstates . to avoid a confusion between the inert and active scalar states , we represent the inert states @xmath32 , @xmath33 and @xmath34 as @xmath35 , @xmath36 and @xmath37 , respectively . their masses are calculated as @xmath38 , \\ m_{\eta_h^{}}^2 & = \mu_0 ^ 2 + \frac{v^2}{2}\left[(\rho_1 + \rho_2 + \rho_3)\cos^2\beta+(\sigma_1 + \sigma_2 + \sigma_3)\sin^2\beta \right ] , \\ m_{\eta_a^{}}^2 & = \mu_0 ^ 2 + \frac{v^2}{2}\left[(\rho_1 + \rho_2 - \rho_3)\cos^2\beta+(\sigma_1 + \sigma_2 - \sigma_3)\sin^2\beta \right ] . \end{aligned}\ ] ]    the mass formulae for the active sector are completely the same as those in a generic 2hdm , so that we can directly apply the same mass formulae to the i(1 + 2)hdm . the mass eigenstates for the active scalar bosons are thus given as : @xmath39 where @xmath40 and @xmath41 are the ng bosons which are absorbed by a longitudinal component of the @xmath42 and @xmath43 bosons , respectively , and @xmath9 is the sm - like higgs boson . the mixing angle @xmath44 is expressed by the ratio of the vevs as @xmath45 . the squared masses of @xmath6 and @xmath7 are then calculated as @xmath46 where @xmath47 the mass matrix for the neutral cp - even ( scalar ) states is expressed in terms of @xmath44 and the usual mixing angle @xmath48 as @xmath49 where the matrix elements are @xmath50,\notag\\ m_{12}^2&=\frac{v^2}{2}\sin2\beta(-\lambda_1\cos^2\beta+\lambda_2\sin^2\beta)+\frac{v^2}{2}\sin2\beta\cos2\beta(\lambda_3+\lambda_4+\lambda_5 ) . \label{mateven}\end{aligned}\ ] ]    one of the most important variables in the active sector of the i(1 + 2)hdm , or equivalently in a 2hdm , is @xmath51 . this describes deviations in the sm - like higgs boson @xmath9 couplings to gauge bosons and fermions . namely , the ratios of the @xmath10 ( @xmath52 ) and @xmath11 couplings in the i(1 + 2)hdm to those in the sm are given by @xmath53 where the @xmath54 factors are @xmath55 or @xmath56 depending on the type of fermion ( i.e. , up - type and down - type quarks and charged leptons ) and the type of yukawa interactions ( i.e. , the aforementioned 2hdm types ) factor corresponds to @xmath57 with @xmath58 being the third component of the isospin , where @xmath59 is given in table  ii of ref . @xcite . ] . the point here is that when we take the limit of @xmath60 , both the @xmath10 and @xmath11 couplings become the same as those in the sm . we thus call this limit the sm - like configuration . furthermore , by looking at eq . ( [ mateven ] ) , we can see that only @xmath61 contains the term proportional to @xmath62 , while @xmath63 and @xmath64 are proportional to @xmath65 . thus , when we take the limit of @xmath66 , the mixing angle given in eq . ( [ tan2a ] ) becomes @xmath67 , which also gives @xmath68 and @xmath69 by choosing @xmath70 . therefore , the extra higgs bosons @xmath6 , @xmath7 and @xmath8 are decoupled due to the @xmath62 term in their mass formulae , so that we can call this limit the decoupling limit  @xcite .    from the above discussion , we can reach an important conclusion regarding the nature of the higgs sector in this model . the sm - like limit is naturally achieved when the decoupling limit is taken . in other words , the decoupling limit can not be realized when we consider the case of a deviation from the sm - like limit , because we need a strong cancellation between @xmath63 and @xmath61 in the denominator of eq . ( [ tan2a ] ) , in order to have @xmath71 . this would force the dimensionless coupling constants @xmath72 in the higgs potential to acquire quite large values which must be excluded from the view point of ensuring perturbativity of the model . this in turn implies that there is an upper limit on the masses of the extra higgs bosons when we consider the case of @xmath73 . as mentioned previously , one of the aims of this paper is indeed deriving an upper limit for these in the i(1 + 2)hdm by using bounds from vacuum stability and perturbative unitarity , which we will discuss in the next section . next , we present the mass formulae in the i(2 + 1)hdm . in this model , @xmath31 corresponds to the sm higgs doublet , so that @xmath32 and @xmath33 correspond to the ng bosons @xmath40 and @xmath41 , respectively , while @xmath34 is the sm - like higgs boson @xmath9 . the mass of @xmath9 is given just like the sm higgs boson as @xmath74 in contrast , the inert sector is composed of two doublets , so that the masses for singly - charged , cp - odd and cp - even states are , respectively , given in the @xmath75 matrix form . the mass matrix for the singly - charged scalar states is evaluated in the basis of ( @xmath76 , @xmath77 ) as @xmath78 and that for the cp - odd and cp - even states are given in the basis of ( @xmath79,@xmath80 ) and ( @xmath81,@xmath82 ) , respectively , @xmath83 the mass eigenstates can be defined by introducing three mixing angles , as follows : @xmath84 the squared mass eigenvalues for the singly - charged states ( @xmath85 , @xmath86 ) are expressed in terms of the elements of the mass matrix as @xmath87 while those for the cp - odd and cp - even states ( @xmath88 , @xmath89 ) and ( @xmath90 , @xmath91 ) , respectively , are obtained by replacing @xmath92 with @xmath93 and @xmath94 . thus , the three mixing angles are given by @xmath95    in summary , the inert sector of the i(2 + 1)hdm can be described by nine parameters : i.e. , six masses @xmath96 , @xmath97 and @xmath98 for @xmath99 and three mixing angles @xmath100 ( for @xmath101 ) . these nine physical parameters are expressed in terms of the nine parameters in the potential of eq . ( [ pot ] ) : i.e. , @xmath102 , @xmath103 , @xmath19 , @xmath104 and @xmath105 ( @xmath106 ) . in this section , we discuss constraints on the parameter space in the i(1 + 2)hdm and the i(2 + 1)hdm . as theoretical constraints , we discuss bounds from vacuum stability and perturbative unitarity . as experimental constraints , other than ensuring sampling the 3hdm parameter spaces concerned away from experimental limits extracted via direct higgs searches at lep , tevatron and lhc we also crucially take into account the so - called @xmath2 , @xmath3 and @xmath4 parameters proposed by peskin and takeuchi  @xcite in the high precision ew fits to lep data @xcite .    for the discussion of vacuum stability and unitarity , there is no difference between in the i(1 + 2)hdm and i(2 + 1)hdm . in contrast , for the @xmath2 , @xmath3 and @xmath4 parameters , different formulae are presented here in each of two models . the higgs potential should be bounded from below in any direction of the scalar boson space . the sufficient condition to guarantee such a positivity of the potential is given by the following requirements : @xmath107      constraints on the scalar self - coupling constants in the potential are extracted from the requirement of @xmath2-matrix unitarity for all the elastic scatterings of two body scalar boson states . this idea to constrain a scalar self - coupling has been proposed in ref . @xcite for the higgs sector in the sm and it has been applied to constrain the mass of the sm higgs boson . it can however be generalized to multi - doublet models , as well known @xcite . the requirement of the @xmath2-matrix unitarity is translated into a relationship for the @xmath108th partial wave amplitude @xmath109 . in the high energy limit , it can be expressed as @xmath110 this implies that @xmath109 should be on a circle with radius and center of @xmath111 and @xmath112 in the complex plane , respectively . from the above relation , we can require the following condition for the tree level amplitude of @xmath109 as  @xcite latexmath:[\ ] ]    the deviations in the @xmath2 and @xmath3 parameters from the sm predictions under @xmath224 gev and @xmath225 are given by  @xcite @xmath226 where the correlation factor between @xmath201 and @xmath202 is + 0.91 . we use the above experimental values for the numerical analysis in the next section . , but we take @xmath227 instead of @xmath228 . , title="fig : " ] , but we take @xmath227 instead of @xmath228 . , title="fig : " ]    , but we take @xmath227 instead of @xmath228 . , title="fig : " ] , but we take @xmath227 instead of @xmath228 . , title="fig : " ]    -@xmath229 plane induced by unitarity and vacuum stability in the case of @xmath230 , @xmath231 , @xmath232 gev and @xmath233 ( left panel ) and 0.996 ( right panel ) . we take @xmath234[300 ] , 500[500 ] , 600[700 ] and 650[900 ] gev for the black , blue , green and red contours in the left[right ] panel . the value of @xmath235 is fixed to be 1 in the dotted contours , while it is scanned over the range of @xmath236 in the solid contours . for all the plots , we scan the value of @xmath62 in the range of @xmath237 tev@xmath238 . outside regions from the contours are excluded by unitarity and vacuum stability . , title="fig : " ] -@xmath229 plane induced by unitarity and vacuum stability in the case of @xmath230 , @xmath231 , @xmath232 gev and @xmath233 ( left panel ) and 0.996 ( right panel ) . we take @xmath234[300 ] , 500[500 ] , 600[700 ] and 650[900 ] gev for the black , blue , green and red contours in the left[right ] panel . the value of @xmath235 is fixed to be 1 in the dotted contours , while it is scanned over the range of @xmath236 in the solid contours . for all the plots , we scan the value of @xmath62 in the range of @xmath237 tev@xmath238 . outside regions from the contours are excluded by unitarity and vacuum stability . , title="fig : " ]    in this section , we discuss how the dimensionless quartic coupling constants in the scalar potential can be constrained from perturbative unitarity , vacuum stability and the @xmath2 and @xmath3 parameters in our 3hdms . among these coupling constants , @xmath239 and @xmath240 just describe four - point interactions among the inert ( pseudo)scalar bosons in the i(1 + 2)hdm and i(2 + 1)hdm , respectively , and they do not affect to the active sector or the interaction between active and inert sectors at the tree level . we thus take these coupling constants to be zero in this section . we first show numerical results in the i(1 + 2)hdm and then in the i(2 + 1)hdm .      in the i(1 + 2)hdm , it is interesting to see how the constraints on the parameter space of the masses of active scalar bosons , @xmath6 , @xmath7 and @xmath8 , can be modified by the requirement of unitarity and vacuum stability as compared to the 2hdm case . as we can see in the sections iii - a and iii - b , the additional coupling constants which do not enter in the 2hdm , such as @xmath104 and @xmath105 , appear in the conditions of vacuum stability and unitarity . therefore , a non - zero value of these coupling constants can modify the constraints on the parameters describing the active ( 2hdm - like ) sector i.e. , @xmath241-@xmath20 , which are translated into bounds on the masses of active ( pseudo)scalar bosons .    in the following , we fix the mass of either @xmath242 or @xmath243 to be @xmath244 gev in order to satisfy the relic abundance of dm  @xcite , and take all the other masses of the inert ( pseudo)scalar bosons to be larger than @xmath245 . for definiteness , we choose @xmath36 as the dm candidate , so that we take @xmath246 . to avoid a large contribution to the @xmath3 parameter , we take @xmath228 or @xmath227 and @xmath231 . for simplicity , we take @xmath247 ( @xmath248 ) . in that case , the coupling constants @xmath249 are expressed in terms of the masses of inert scalar bosons as follows : @xmath250 we can see that the sign of @xmath21 is now positive due to the assumption of @xmath251 . the sign becomes negative when we consider the other case , of @xmath37 as the dm candidate , but the bounds from vacuum stability and unitarity do not depend on the sign of @xmath21 . from the vacuum stability condition of ( [ vs2 ] ) , the value of @xmath229 should be larger than @xmath252 , thus we take @xmath253 .    as we discussed in sec . ii - b , when @xmath73 is taken , we can not have the decoupling limit of the extra active higgs bosons and thus we can obtain an upper limit of their masses . in addition , recall that the value of @xmath51 describes the deviation in the @xmath10 couplings as shown in eq . ( [ couplings ] ) at the tree level so that we may be able to find a relationship between the upper limit on the masses and the deviation in the @xmath10 couplings . because the 1-@xmath254 error of the measurement of @xmath10 couplings is expected to be 2% at the hl - lhc and 0.4% at the ilc500  @xcite , we take @xmath233 and 0.996 , respectively , as examples .    in fig . [ const1 ] , we show the constraint on the parameter space of the i(1 + 2)hdm mapped onto the @xmath255-@xmath256 plane by unitarity , vacuum stability and the @xmath2 and @xmath3 parameters . we take @xmath233 in the left panel and 0.996 in the right panel . the value of @xmath235 is fixed to be 1 in the dotted contours , while it is scanned over the range @xmath236 in the solid contours . for all the plots , we also scan the value of @xmath62 in the range of @xmath257 . we check that these scanned ranges are enough wide to obtain the maximally allowed region on the parameter space from the unitarity and vacuum stability constraints . the blue , green and red contours show the cases of @xmath258 , 500 and 550 gev , respectively . we also add a plot for the case in the 2hdm as the black contour for comparison . the regions outside each contour are excluded by unitarity and vacuum stability . in addition , the light and dark shaded regions are excluded by the @xmath2 and @xmath3 parameters in the 2hdm and i(1 + 2)hdm , respectively . we note that the contribution to @xmath202 is given to be zero due to @xmath228 , so that the excluded region denoted by the light and dark area are from the constraint by @xmath201 . besides , we confirm that the difference in @xmath201 due to the different choices of @xmath259 is negligible .    in fig . [ const2 ] , we show the case of @xmath227 with the rest of the setup being the same as in fig . [ const1 ] . in this case , there is a non - zero contribution to @xmath202 which is proportional to @xmath260 , so that the excluded regions denoted by light and dark shading are different from those in fig . [ const1 ] . we find that the upper limits on @xmath256 and @xmath255 are related to each other and are getting stronger by taking larger values of @xmath259 . we note that , if we take the limit @xmath261 , the constraint from unitarity and vacuum stability approaches to that in the 2hdm .    in fig . [ const3 ] , we show the constraints on the parameter space on the @xmath21-@xmath229 plane by unitarity and vacuum stability in the cases of @xmath233 in the left panel and @xmath262 in the right panel , where the regions outside each contour are excluded . the values of @xmath235 and @xmath62 are taken as done in figs . [ const1 ] and [ const2 ] . the black , blue , green and red contours show the cases of @xmath234[300 , 500[500 ] , 600[700 ] and 650[900 ] gev in the left[right ] panel , where @xmath263 . we note that the parameter space shown in this plot is allowed from the @xmath2 and @xmath3 parameters . we can see that the region with @xmath264 is excluded by the vacuum stability condition in eq . ( [ vs2 ] ) . furthermore , the region above each contour is excluded by the unitarity constraints through the largest eigenvalue of the matrix @xmath186 given in eq .  ( [ 3times3 ] ) .      next , we consider an application of our bounds from unitarity and vacuum stability in the i(2 + 1)hdm . as mentioned in the beginning of this section , we take @xmath241-@xmath20 to be zero . in that case , we obtain the analytic expression of the eigenvalues for three sub - matrices of the @xmath0-wave amplitude given in eq . ( [ 3times3 ] ) as @xmath265 where @xmath266 are the orthogonal matrices to diagonalize @xmath267 . now , @xmath239 is determined by @xmath268 from eq . ( [ mh ] ) , so that the positive eigenvalues are important to get a stronger bound .    in the following calculation , @xmath269 is taken to be the dm candidate with mass @xmath245 . for simplicity , we consider a no - mixing case between two inert scalar states , which can be realized by setting @xmath270 to be exactly zero , the lightest state among @xmath89 , @xmath91 and @xmath86 is also stable . in order to have a decay mode of it , a non - zero value of @xmath19 is required . ] . in addition , to avoid the constraint from the @xmath3 parameter , we take @xmath271 and @xmath272 . we note that results do not change if we replace both[either ] @xmath273 and[or ] @xmath274 with @xmath275 and[or ] @xmath276 . in the above set up , we obtain the following relations for the quartic coupling constants : @xmath277 from the vacuum stability condition in eqs . ( [ vs2 ] ) and ( [ vs3 ] ) , we have @xmath278 by the combination of the bounds from unitarity and vacuum stability , the case with @xmath253 and @xmath279 gives the largest allowed parameter space . in this setup , we find that the eigenvalue of eq . ( [ pu1 ] ) gives the strongest bound . -@xmath229 plane induced by unitarity and vacuum stability . outside regions from the contours are excluded by unitarity and vacuum stability . ]    in fig . [ const4 ] , we show the allowed parameter space mapped onto the @xmath21-@xmath229 plane in the i(2 + 1)hdm . regions outside the contours are excluded by unitarity and vacuum stability . more precisely , the region below the contour with @xmath280 is excluded by vacuum stability and the region above the contour is excluded by unitarity . we have discussed the constraints on the parameter space of the i(1 + 2)hdm and i(2 + 1)hdm emerging from the bounds on vacuum stability , perturbative unitarity and the ew @xmath2 , @xmath3 and @xmath4 parameters . to impose the unitarity bound , we have calculated the @xmath0-wave amplitude matrix for all possible 2to2 body ( pseudo)scalar boson elastic scatterings in the high energy limit by means of the eta , where the sub - matrices with @xmath281 , @xmath282 and @xmath283 forms for the electromagnetically neutral , singly - charged and doubly - charged channels , respectively , can be in block diagonal form with @xmath1 sub - matrices as the largest component . we have first applied our formulae for the unitarity and vacuum stability bounds to constrain the parameter space in the i(1 + 2)hdm , especially mapped onto the @xmath255-@xmath256 plane . we have found that larger excluded regions on such a plane can be obtained in the case with @xmath73 as compared to those in 2hdms . we then have applied our procedure to the i(2 + 1)hdm on the @xmath21-@xmath229 plane . here , it has been clarified that the values @xmath284 , 7 , 6 and 2 are respectively excluded for any values of @xmath229 in the case of @xmath285 , @xmath286 , @xmath287 and @xmath288 . +   + @xmath289                                            f.  hartmann and w.  kilian , eur . j.  c * 74 * , 3055 ( 2014 ) . i. p. ivanov and v. keus , phys . d * 86 * , 016004 ( 2012 ) . i. p. ivanov and e. vdovin , eur . j. c * 73 * , 2309 ( 2013 ) . v.  keus , s.  f.  king and s.  moretti , jhep * 1401 * , 052 ( 2014 ) . m.  maniatis and o.  nachtmann , arxiv:1408.6833 [ hep - ph ] . i.  p.  ivanov and c.  c.  nishi , jhep * 1501 * , 021 ( 2015 ) . d.  das and u.  k.  dey , phys . d * 89 * , no . 9 , 095025 ( 2014 ) . v.  keus , s.  f.  king and s.  moretti , phys . d * 90 * ( 2014 ) 075015 . m.  aoki , r.  guedes , s.  kanemura , s.  moretti , r.  santos and k.  yagyu , phys . d * 84 * , 055028 ( 2011 ) . c.  w.  chiang and k.  yagyu , jhep * 1307 * , 160 ( 2013 ) . j.  f.  gunion and h.  e.  haber , phys . d * 67 * , 075019 ( 2003 ) . d.  lopez - val and j.  sola , eur . j.  c * 73 * , 2393 ( 2013 ) ; m.  e.  peskin and j.  d.  wells , phys . d * 64 * , 093003 ( 2001 ) ; s.  kanemura , y.  okada , h.  taniguchi and k.  tsumura , phys . b * 704 * , 303 ( 2011 ) . l.  lopez honorez , e.  nezri , j.  f.  oliver and m.  h.  g.  tytgat , jcap * 0702 * , 028 ( 2007 ) . t.  hambye , f.  -s . ling , l.  lopez honorez and j.  rocher , jhep * 0907 * , 090 ( 2009 ) [ erratum - ibid .   * 1005 * , 066 ( 2010 ) ] .
Abstract:

Answer: we calculate an @xmath0-wave amplitude matrix for all the possible 2to2 body scalar boson elastic scatterings in models with three scalar doublets , including contributions from the longitudinal component of weak gauge bosons via the equivalence theorem approximation . specifically , we concentrate on the two cases with two[one ] active plus one[two ] inert doublet fields , referred to as i(1 + 2)hdm[i(2 + 1)hdm ] , under cp conservation . we obtain three analytically irreducible sub - matrices with the @xmath1 form and eighteen eigenvalues for the amplitude matrix as an independent set , where the same formula can be applied to both models . by requiring a perturbative unitarity condition , we can constrain the magnitude of quartic coupling constants in the higgs potential . this constraint , in particular in the i(1 + 2)hdm , can be translated into a bound on masses of extra active scalar bosons . furthermore , when standard model - like higgs boson couplings with weak gauge bosons are deviated from the standard model predictions , the unitarity condition provides an upper limit on the masses . we find that stronger upper bounds on the masses of the active cp - even and cp - odd higgs bosons are obtained under the constraints from the unitarity and vacuum stability conditions , as well as the electroweak @xmath2 , @xmath3 and @xmath4 parameters , as compared to those in 2-higgs doublet models with a softly - broken @xmath5 symmetry .