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: 2d semiconductors of group v elements , including phosphorene and arsenene , have been rapidly attracting interest due to their significant fundamental band gap , large density of states near the fermi level , and high and anisotropic carrier mobility . combination of these properties places these systems very favorably in the group of contenders for 2d electronics applications beyond graphene and transition metal dichalcogenides@xcite . keeping in mind that the scope of group iv semiconductors such as si has been broadened significantly by introducing isoelectronic iii - v compounds , it is intriguing to see , whether the same can be achieved in iv - vi compounds that are isoelectronic to group v elemental semiconductors . even though this specific point of view has not yet received attention , there has been interest in specific iv - vi compounds , such as sns , snse and gete , for thermoelectric and photovoltaic applications . it appears likely that specific search for isoelectronic counterparts of layered semiconductors such as phosphorene and arsenene may guide us to yet unexplored 2d semiconducting iv - vi compounds that are stable , flexible , and display a tunable band gap .        as a yet unexplored iv - vi compound , we study the layered structure of silicon mono - sulfide ( sis ) . we use _ ab initio _ density functional theory ( dft ) to identify stable allotropes , determine their equilibrium geometry and electronic structure . we have identified two nearly equally stable allotropes , namely the black - phosphorus - like @xmath0-sis and the blue - phosphorus - like @xmath1-sis , and show their structure in fig . [ fig1](a ) and ( d ) . both @xmath0-sis and @xmath1-sis monolayers display a significant , indirect band gap that depends sensitively on the in - layer strain . unlike 2d semiconductors of group v elements with the corresponding nonplanar structure , different sis allotropes show a strong polarization either within or normal to the layers . we find that sis may form both lateral and vertical heterostructures with phosphorene at a very small energy penalty , offering an unprecedented tunability in structural and electronic properties of sis - p compounds . associated with states in the energy range between the fermi level @xmath2 and 0.2  ev below the top of the valence band , shown by the shaded region in fig . [ fig1 ] , in a monolayer of ( a ) @xmath0-sis and ( b ) @xmath1-sis . @xmath3  e / @xmath4 contours are superposed with ball - and - stick models of related structures . [ fig2 ] ]    since all atoms in @xmath5 layered structures of group v elements are threefold coordinated , the different allotropes can all be topologically mapped onto the honeycomb lattice of graphene with 2 sites per unit cell.@xcite an easy way to generate iv - vi compounds that are isoelectronic to group v monolayers is to occupy one of these sites by a group iv and the other by a group vi element . in this way , we have generated the orthorhombic @xmath0-sis monolayer structure , shown in fig . [ fig1](a ) , from a monolayer of black phosphorene ( or @xmath0-p ) . the hexagonal @xmath1-sis monolayer , shown in fig . [ fig1](d ) , has been generated in the same way from the blue phosphorene ( or @xmath1-p ) structure . the monolayer structures have been optimized using dft with the perdew - burke - ernzerhof ( pbe)@xcite exchange - correlation functional , as discussed in the methods section . we found that presence of two elements with different local bonding preferences increases the thickness of the sis monolayers when compared to the phosphorene counterparts . the 2d lattice of @xmath0-sis is spanned by the orthogonal bravais lattice vectors @xmath6   and @xmath7   , which are about 3% longer than the lattice vectors of black phosphorene . the 2d hexagonal lattice of @xmath1-sis is spanned by two bravais lattice vectors @xmath8   , which are @xmath9% longer than those of blue phosphorene . unlike in the counterpart structures of the phosphorene monolayer , we find @xmath1-sis to be more stable by about 12  mev / atom than @xmath0-sis . we checked the phonon spectrum of free - standing @xmath0-sis and @xmath1-sis monolayers and found no soft modes that would cause a spontaneous collapse of the structure . an apparently different layered sis structure had been synthesized@xcite by exposing the layered casi@xmath10 structure to s@xmath10cl@xmath10 . it appears that the s@xmath10cl@xmath10 reagent reacted with the ca atoms intercalated between silicene layers in casi@xmath10 by forming cacl@xmath10 and saturating the valencies of the silicene layers by @xmath11s@xmath11s bridges . even though no structural information has been provided in that study , we identified a likely candidate structure and found it to be less stable than the postulated sis allotropes , on which we focus next . since the electronegativity of s is higher than that of si , we expect an electron transfer from si to s atoms . based on a mulliken population analysis , we estimate a net transfer of 0.3 electrons in @xmath0-sis and 0.2 electrons in @xmath1-sis from silicon to sulfur atoms . this charge redistribution , combined with the nonplanarity of the structure , causes a net polarization . we find an in - plane polarization for @xmath0-sis and and an out - of - plane polarization for @xmath1-sis .    whereas dft generally provides an accurate description of the total charge density and equilibrium geometry , interpretation of kohn - sham energy eigenvalues as quasiparticle energies is more problematic . still , we present our dft - pbe results for the electronic structure of @xmath0-sis and @xmath1-sis monolayers in fig . even though the fundamental band gaps are typically underestimated in this approach , the prediction that both systems are indirect - gap semiconductors is likely correct . our calculated band structure and the corresponding density of states for @xmath0-sis , presented in fig .  [ fig1](b , c ) , suggest that the fundamental band gap value @xmath12  ev should be significantly larger than in the isoelectronic black phosphorene counterpart . the band structure near the top of the valence band shows a significant anisotropy when comparing the @xmath13 and @xmath14 directions . in analogy to phosphorene , @xmath0-sis should exhibit a higher hole mobility along the @xmath15-direction than along the @xmath16-direction the dft - based fundamental band gap @xmath17  ev in the monolayer of @xmath1-sis is even larger . the band structure in the symmetric honeycomb lattice of @xmath1-sis is rather isotropic , as seen in fig . [ fig1](e ) . the top of the valence band is very flat , resulting in a heavy hole mass and a large density of states ( dos ) in that region , as seen in fig . [ fig1](f ) . the character of frontier states is of interest not only for a microscopic understanding of the conduction channels , but is also crucial for the design of optimum contacts.@xcite whereas dft - based band gaps are typically underestimated as mentioned above , the electronic structure of the valence and the conduction band region in dft is believed to closely correspond to experimental results . in fig . [ fig2 ] , we show the charge density associated with frontier states near the top of the valence band . these states , which correspond to the energy range highlighted by the green shading in the band structure of @xmath0-sis in fig . [ fig1](b ) and that of @xmath1-sis in fig . [ fig1](e ) , cover the energy range between the mid - gap and 0.2  ev below the top of the valence band . the valence frontier states of @xmath0-sis in fig .  [ fig2](a ) and @xmath1-sis in fig .  [ fig2](b ) are similar in spite of the notable charge density differences caused by the larger dos of @xmath1-sis in this energy range . the difference between s and si atoms is also reflected in the character of the valence frontier states at these sites . whereas @xmath18 states contribute most at sulfur sites , silicon sites contribute a mixture of @xmath19 and @xmath18 states . these frontier states differ from those of phosphorene , which are related to lone pair electron states . similar to phosphorene , the fundamental band gap value of sis also depends sensitively on the in - layer strain , as seen in fig . [ fig3 ] . due to their nonplanarity , accordion - like in - layer stretching or compression of sis structures may be achieved at little energy cost as shown in the supporting information . the energy cost is particularly low for a deformation along the soft @xmath15-direction , requiring @xmath20  mev / atom to induce a @xmath21% in - layer strain . we believe that in view of the softness of the structure , similar strain values may be achieved during epitaxial growth on particular incommensurate substrates . we also note that tensile strain values such as these have been achieved experimentally in suspended graphene membranes that are much more resilient to stretching due to their planar geometry and stronger bonds . consequently , we believe that strain engineering is a viable way to effectively tune the fundamental band gap in these systems .    our results for @xmath0-sis in fig . [ fig3](a ) indicate that the band gap decreases when the structure is compressed and increases slightly when it is stretched . the largest change in the band gap , namely its reduction to 0.3  ev , may be achieved during a @xmath22% compression . as seen in fig . [ fig3](b ) , we expect the fundamental band gap of @xmath1-sis to be reduced during both stretching and compression . within the @xmath21% range , we find that the band gap may be tuned in the range from @xmath23  ev . this high degree of band gap tunability in sis appears particularly attractive for potential applications in flexible electronics . since the geometry and lattice constants of sis and phosphorene are very similar , it is likely that the two could interface naturally in lateral and vertical heterostructures , thus further advancing the tunability of their electronic properties . in fig  [ fig4 ] , we present electronic structure results for bilayers consisting of an sis and a phosphorene monolayer in both @xmath0- and @xmath1-phases as the simplest examples of vertical heterostructures . we have optimized the bilayer structures assuming commensurability , i.e. setting the primitive unit cells of each monolayer to be the same . the optimum geometry of the @xmath0-(sis / p ) bilayer is shown in fig  [ fig4](a ) and that of the @xmath1-(sis / p ) bilayer in fig  [ fig4](b ) . we find the inter - layer interaction in the two bilayer systems to be rather weak , amounting to @xmath20  mev / atom based on our dft - pbe calculations . taking better account of van der waals interactions would likely increase the calculated interaction energy and reduce the large interlayer separation of @xmath24   . whereas the precise interlayer interaction and separation are not of primary concern here , our most important finding is that the weak interaction is not purely dispersive in nature , since we find a substantial rehybridization of states between the adjacent sis and phosphorene layers . consequently , as we show more explicitly in the supporting information , the bilayer band structure is not a mere superposition of the two monolayer band structures in the same assumed geometry . as shown in fig . [ fig4](b ) , the @xmath0-(sis / p ) bilayer is a direct - gap semiconductor with @xmath25  ev , smaller than the band gap in either isolated monolayer . since the gap is indirect in the @xmath0-sis monolayer , whereas it is direct in @xmath0-p , the cause for the direct band gap in @xmath0-sis / p may be the dominance of p states near the gap . we find this conjecture confirmed by examining the frontier orbitals of the bilayer . as we show in the supporting information , the frontier states in the conduction band region of the bilayer are essentially purely p - based , whereas those in the valence band region have only a small contribution from the sis layer . the electronic band structure of the aa - stacked @xmath1-(sis / p ) bilayer is presented in fig . [ fig4](d ) . the fundamental band gap is indirect in this system and its value @xmath26  ev is again smaller than in isolated blue phosphorene and @xmath1-sis monolayers . as seen in the supporting information , also the frontier states of @xmath1-(sis / p ) are dominated by p orbitals both in the valence and the conduction band region . since both sis and phosphorene are rather flexible , they may adjust to each other and form also in - layer heterostructures at little or no energy penalty . we constructed two types of sisp@xmath10 heterostructures for both @xmath0 and @xmath1 allotropes and show their geometry and electronic structure in fig . one type of heterostructures , shown in fig . [ fig5](a ) and [ fig5](e ) , contains p - p and si - s atom pairs completely separated from like atom pairs . the other type of heterostructures , presented in fig .  [ fig5](c ) and [ fig5](g ) , contains alternating , contiguous sis and phosphorus chains . these structures maintain a rectangular lattice with four atoms per unit cell . generally , we found the in - layer heterostructures to be less stable than pure phosphorene and sis monolayers . the least stable heterostructures among these are those with isolated p - p or si - s atom pairs , shown in fig . [ fig5](a ) and [ fig5](e ) , which are @xmath27  ev / atom less stable than sis and phosphorene monolayers due to their highly frustrated geometries . the heterostructures with contiguous phosphorus and sis chains , presented in fig . [ fig5](c ) and [ fig5](g ) , may better optimize the nearest - neighbor environment . this causes less frustration , making these systems only @xmath28  ev / atom less stable than isolated sis and phosphorene monolayers . we found all four in - layer heterostructures to be indirect - gap semiconductors . as in the vertical heterostructures , we found the fundamental band gaps to be substantially smaller than in isolated sis and phosphorene monolayers . as seen in fig . [ fig5](a , b ) and [ fig5](e , f ) , the fundamental band gap @xmath29 is close to 0.5  ev in the less stable heterostructures with isolated si - s and p - p pairs . we found larger band gap values in the more stable heterostructures with contiguous sis and p chains , namely @xmath30  ev in @xmath0-sisp@xmath10 shown in fig . [ fig5](c , d ) and @xmath31  ev in @xmath1-sisp@xmath10 shown in fig . [ fig5](g , h ) , both nearly half the value of the smaller band gap in either sis or phosphorene monolayers . these findings indicate an intriguing possibility of isoelectronic doping as an effective way to tune the electronic properties of sisp@xmath32 systems .    in view of the different successful synthesis approaches used to form layered iv - vi compounds such as sis@xcite or sns@xcite with a structure similar to @xmath0-sis , we believe that a suitable synthesis path will be found to form also the allotropes introduced in this study . the weak inter - layer interaction in layered pure sis compounds or sis / p heterostructures should allow for a mechanical exfoliation yielding monolayers and bilayers , same as in graphene and phosphorene . chemical vapor deposition ( cvd ) , which had been used successfully in the past to grow graphene and silicene@xcite , may become ultimately the most common approach to grow few - layer sis on specific substrates . as suggested above , compounds that are isoelectronic to group v layered semiconductors are not limited to pure iv - vi systems , but may contain other group v elements in the same layer or in a vertical few - layer heterostructure . the initial interest in black phosphorene could thus be significantly expanded by considering the whole range of group iv elements including si , ge , sn , and pb , group vi elements such as s , se , and te , and group v elements such as p , as , sb and bi , leading to a virtually limitless number of compounds and structural phases . the large family of iv - vi systems that are isoelectronic to group v elemental semiconductors will include 2d semiconductors with a sizeable fundamental band gap , a high carrier mobility and chemical stability . expanding this family to heterostructures containing also group v elements should provide a way to effectively tune the electronic properties of the pristine structures , which will likely bring unlimited richness to the field of 2d semiconductors .    in conclusion , we have proposed iv - vi compounds as isoelectronic counterparts to layered group v semiconductors in analogy to iii - v compounds , which have significantly broadened the scope of group iv semiconductors . using _ ab initio _ density functional theory , we have identified yet unrealized structural phases of silicon mono - sulfide ( sis ) including the black - phosphorus - like @xmath0-sis and the almost equally stable , blue - phosphorus - like @xmath1-sis . we found that both @xmath0-sis and @xmath1-sis monolayers display a significant , indirect band gap that depends sensitively on the in - layer strain . unlike 2d semiconductors of group v elements with the corresponding nonplanar structure , different sis allotropes are polar . we found that sis may form both lateral and vertical heterostructures with phosphorene at a very small energy penalty , offering an unprecedented tunability in structural and electronic properties of sis - p compounds . combining other group iv and group vi elements with group v elements is expected to lead to a large family of layered semiconductors with an unprecedented richness in structural and electronic properties . our computational approach to study the equilibrium structure , stability and electronic properties of sis structures is based on _ ab initio _ density functional theory with the perdew - burke - ernzerhof ( pbe)@xcite exchange - correlation functional as implemented in the siesta@xcite and vasp@xcite codes . we use periodic boundary conditions throughout the study , with few - layer structures represented by a periodic array of slabs separated by a vacuum region of @xmath33   . in our siesta calculations we used norm - conserving troullier - martins pseudopotentials@xcite and a double-@xmath34 basis including polarization orbitals . the reciprocal space has been sampled by a fine grid@xcite of @xmath35  @xmath36-points in the brillouin zone of the primitive unit cell of the 2d structures . we used a mesh cutoff energy of @xmath37  ry to determine the self - consistent charge density , which provided us with a precision in total energy of @xmath38  mev / atom . in vasp calculations , we used an energy cut - off of 500  ev and the same k - point sampling as mentioned above . all geometries have been optimized using the conjugate gradient method@xcite , until none of the residual hellmann - feynman forces exceeded @xmath39  ev /  . equilibrium structures and energies based on siesta have been checked against values based on the vasp code . this study has been supported by the national science foundation cooperative agreement # eec-0832785 , titled `` nsec : center for high - rate nanomanufacturing '' . computational resources have been provided by the michigan state university high performance computing center . [ 1]#1 @ifundefined      liu ,  h. ; neal ,  a.  t. ; zhu ,  z. ; luo ,  z. ; xu ,  x. ; tomanek ,  d. ; ye ,  p.  d. phosphorene : an unexplored 2d semiconductor with a high hole mobility . _ acs nano _ * 2014 * , _ 8 _ , 40334041 zhu ,  z. ; tomnek ,  d. semiconducting layered blue phosphorus : a computational study . _ phys . lett . _ * 2014 * , _ 112 _ , 176802 guan ,  j. ; zhu ,  z. ; tomnek ,  d. phase coexistence and metal - insulator transition in few - layer phosphorene : a computational study . _ phys . rev . lett . _ * 2014 * , _ 113 _ , 046804 li ,  l. ; yu ,  y. ; ye ,  g.  j. ; ge ,  q. ; ou ,  x. ; wu ,  h. ; feng ,  d. ; chen ,  x.  h. ; zhang ,  y. black phosphorus field - effect transistors . _ nature nanotech . _ * 2014 * , _ 9 _ , 373377 han ,  m.  y. ; zyilmaz ,  b. ; zhang ,  y. ; kim ,  p. energy band - gap engineering of graphene nanoribbons . lett . _ * 2007 * , _ 98 _ , 206805 elias ,  d.  c. ; nair ,  r.  r. ; mohiuddin ,  t. m.  g. ; morozov ,  s.  v. ; blake ,  p. ; halsall ,  m.  p. ; ferrari ,  a.  c. ; boukhvalov ,  d.  w. ; katsnelson ,  m.  i. ; geim ,  a.  k. control of graphene s properties by reversible hydrogenation : evidence for graphane . _ science _ * 2009 * , _ 323 _ , 610613 radisavljevic ,  b. ; radenovic ,  a. ; brivio ,  j. ; giacometti ,  v. ; kis ,  a. single - layer mos2 transistors . _ nature nanotech . _ * 2011 * , _ 6 _ , 147150 zhao ,  l .- d . ; lo ,  s .- h . ; zhang ,  y. ; sun ,  h. ; tan ,  g. ; uher ,  c. ; wolverton ,  c. ; dravid ,  v.  p. ; kanatzidis ,  m.  g. ultralow thermal conductivity and high thermoelectric figure of merit in snse crystals . _ nature _ * 2014 * , _ 508 _ , 373377 mathews ,  n. electrodeposited tin selenide thin films for photovoltaic applications . _ solar energy _ * 2012 * , _ 86 _ , 10101016 sinsermsuksakul ,  p. ; heo ,  j. ; noh ,  w. ; hock ,  a.  s. ; gordon ,  r.  g. atomic layer deposition of tin monosulfide thin films . energy mater . _ * 2011 * , _ 1 _ , 11161125 guan ,  j. ; zhu ,  z. ; tomnek ,  d. tiling phosphorene . _ acs nano _ * 2014 * , _ 8 _ , 1276312768 perdew ,  j.  p. ; burke ,  k. ; ernzerhof ,  m. generalized gradient approximation made simple . _ phys . rev . lett . _ * 1996 * , _ 77 _ , 38653868 hengge ,  e. ; olbrich ,  g. polymere siliciumsulfide mit schichtstruktur . _ z. anorg . . chemie _ * 1969 * , _ 365 _ , 321324 tomnek ,  d. interfacing graphene and related 2d materials with the 3d world . _ j. phys . : condens . matter _ * 2015 * , _ 27 _ , 133203 lee ,  c. ; wei ,  x. ; kysar ,  j.  w. ; hone ,  j. measurement of the elastic properties and intrinsic strength of monolayer graphene . _ science _ * 2008 * , _ 321 _ , 385388 frank ,  i.  w. ; tanenbaum ,  d.  m. ; van  der zande ,  a.  m. ; mceuen ,  p.  l. mechanical properties of suspended graphene sheets . _ b _ * 2007 * , _ 25 _ , 25582561 huang ,  m. ; yan ,  h. ; heinz ,  t.  f. ; hone ,  j. probing strain - induced electronic structure change in graphene by raman spectroscopy . _ nano lett . _ * 2010 * , _ 10 _ , 40744079 su ,  h. ; xie ,  y. ; xiong ,  y. ; gao ,  p. ; qian ,  y. preparation and morphology control of rod - like nanocrystalline tin sulfides via a simple ethanol thermal route . _ j. solid state chem . _ * 2001 * , _ 161 _ , 190  196 kim ,  k.  s. ; zhao ,  y. ; jang ,  h. ; lee ,  s.  y. ; kim ,  j.  m. ; kim ,  k.  s. ; ahn ,  j .- h . ; kim ,  p. ; choi ,  j .- y . ; hong ,  b.  h. large - scale pattern growth of graphene films for stretchable transparent electrodes . _ nano lett . _ * 2009 * , _ 457 _ , 706 reina ,  a. ; jia ,  x. ; ho ,  j. ; nezich ,  d. ; son ,  h. ; bulovic ,  v. ; dresselhaus ,  m.  s. ; kong ,  j. large area , few - layer graphene films on arbitrary substrates by chemical vapor deposition . _ nano lett . _ * 2009 * , _ 9 _ , 30 vogt ,  p. ; de  padova ,  p. ; quaresima ,  c. ; avila ,  j. ; frantzeskakis ,  e. ; asensio ,  m.  c. ; resta ,  a. ; ealet ,  b. ; le  lay ,  g. silicene : compelling experimental evidence for graphenelike two - dimensional silicon . lett . _ * 2012 * , _ 108 _ , 155501 artacho ,  e. ; anglada ,  e. ; dieguez ,  o. ; gale ,  j.  d. ; garcia ,  a. ; junquera ,  j. ; martin ,  r.  m. ; ordejon ,  p. ; pruneda ,  j.  m. ; sanchez - portal ,  d. the siesta method ; developments and applicability . _ j. phys . cond . mat . _ * 2008 * , _ 20 _ , 064208 kresse ,  g. ; furthmller ,  j. efficient iterative schemes for _ ab initio _ total - energy calculations using a plane - wave basis set . b _ * 1996 * , _ 54 _ , 1116911186 kresse ,  g. ; hafner ,  j. _ ab initio _ molecular dynamics for liquid metals . _ b _ * 1993 * , _ 47 _ , 558561 kresse ,  g. ; hafner ,  j. _ ab initio _ molecular - dynamics simulation of the liquid - metal  amorphous - semiconductor transition in germanium . _ b _ * 1994 * , _ 49 _ , 1425114269 kresse ,  g. ; joubert ,  d. from ultrasoft pseudopotentials to the projector augmented - wave method . _ b _ * 1999 * , _ 59 _ , 17581775 troullier ,  n. ; martins ,  j.  l. efficient pseudopotentials for plane - wave calculations . _ b _ * 1991 * , _ 43 _ , 1993 monkhorst ,  h.  j. ; pack ,  j.  d. special points for brillouin - zone integrations . _ b _ * 1976 * , _ 13 _ , 51885192 hestenes ,  m.  r. ; stiefel ,  e. methods of conjugate gradients for solving linear systems . _ * 1952 * , _ 49 _ , 409436
Abstract:

Answer: in analogy to iii - v compounds , which have significantly broadened the scope of group iv semiconductors , we propose iv - vi compounds as isoelectronic counterparts to layered group v semiconductors . using _ ab initio _ density functional theory , we study yet unrealized structural phases of silicon mono - sulfide ( sis ) . we find the black - phosphorus - like @xmath0-sis to be almost equally stable as the blue - phosphorus - like @xmath1-sis . both @xmath0-sis and @xmath1-sis monolayers display a significant , indirect band gap that depends sensitively on the in - layer strain . unlike 2d semiconductors of group v elements with the corresponding nonplanar structure , different sis allotropes show a strong polarization either within or normal to the layers . we find that sis may form both lateral and vertical heterostructures with phosphorene at a very small energy penalty , offering an unprecedented tunability in structural and electronic properties of sis - p compounds .


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: strong sensitivity to initial conditions is a signature of chaotic behavior in dynamical systems @xcite . in closed quantum systems , however , such a behavior is excluded by the linearity of the unitarity of the time evolution . nevertheless , despite this characteristic linearity of closed quantum systems , chaotic features can arise in open quantum systems @xcite . in particular , non - linear quantum state transformations can arise when identically prepared quantum systems are subjected to an entangling unitary transformation and a subsequent selective measurement performed on parts of the system @xcite . iterating such non - linear quantum state transformations may result in strong dependence on the initial conditions and in complex chaos @xcite . quantum state purification protocols involving identically prepared weakly entangled two - qubit systems , as introduced by bennett and deutsch @xcite , are typical examples involving repeated sequences of unitary entangling transformations and selective conditional quantum measurements @xcite . in these protocols the parameters are designed in such a way that quantum state purification can be achieved and a complicated dependence on initial conditions is avoided @xcite . nevertheless , it has been demonstrated that in general an intricate dependence of the resulting entanglement on the initially prepared states will occur @xcite . sensitivity to initial conditions can also be used for amplifying small initial differences of quantum states ( schrdinger microscrope @xcite ) . however , it has been demonstrated that such an amplification requires necessarily a large number of identically prepared systems which have to be discarded during the process @xcite . recent advances of quantum technology have enabled the control of complicated many - body quantum systems to an unprecedented degree of accuracy so that even subtle effects of quantum measurement - induced nonlinearities may become amenable to experimental observation @xcite . in particular , the area of cavity quantum electrodynamics offers interesting perspectives for future experiments in these directions in view of its possibilities to measure and control the interaction of individual atoms with a single mode of the quantized radiation field with high precision @xcite .    in this paper we explore the intricate nonlinear properties of measurement - induced chaos in a cavity quantum electrodynamical scenario . this scenario is a generalization of the dynamics of a one - atom - micro - maser @xcite which may be viewed as an early example of an iterative unitary quantum dynamics conditioned on state selective quantum measurements . motivated by recent experimental advances which realize the tavis - cummings model @xcite we consider an ensemble of identically prepared two - level atoms ( qubits ) which interact pairwise with a single mode of the radiation field . afterwards , one member of each pair and the corresponding cavity field are measured . conditioned on the measurement results , the unmeasured atoms are kept or discarded . in practice , this may be implemented with the help of a single cavity and a pair of optical conveyor belts @xcite . the atoms are moved through the cavity by the conveyor belts in such a way that only one pair of atoms interacts with the cavity mode at a time and then the cavity is re - initialized after each interaction . the remaining atoms form a new identically prepared ensemble of smaller size . similar as in entanglement distillation protocols , the state changes of the remaining two - level atoms are described by an iterated nonlinear quantum transformation . it is demonstrated that this nonlinear quantum transformation exhibits measurement - induced chaos and the resulting julia sets can be connected or disconnected . furthermore , it is demonstrated that this nonlinear quantum map can implement a schrdinger microscope capable of amplifying the distinguishability of non - orthogonal quantum states . the paper is organized as follows . in sec . [ model ] we solve the dynamical equation of the two - atom tavis - cummings model . we present exact and approximate analytical solutions that facilitate the analysis of the dynamics . furthermore , we briefly explain the atomic postselection scheme assisted by projection of the field . this is used in sec . [ protocol ] to introduce a protocol that implements a nonlinear mapping of atomic probability amplitudes . in sec . [ analysis ] we analyze the fractal structure in the complex plane after iteration of the ideal map and compare with the numerical implementation using the exact - numerical quantum operation . in sec . v we present a possible application , where the non - linear dynamics transforms two , initially close quantum states into almost orthogonal ones in a few steps of the protocol . finally , in sec . vi we propose a setup based on currently available experimental technology capable of implementing this few - step scenario . the two atoms tavis - cummings model describes the resonant interaction between two atoms , say @xmath0 and @xmath1 , and a single mode of the radiation field @xcite . the atoms have ground states @xmath2 and excited states @xmath3 ( @xmath4 ) separated by an energy difference of @xmath5 which matches the energy of a photon inside the empty cavity . in the interaction picture the hamiltonian can be expressed in the following form @xmath6 where @xmath7 and @xmath8 are the atomic raising and lowering operators ( @xmath4 ) , and @xmath9 ( @xmath10 ) is the annihilation ( creation ) operator of the single mode field . the interaction picture is taken with respect to the reference hamiltonian @xmath11 which is a constant of motion as it commutes with the interaction hamiltonian @xmath12 . for this reason both operators , @xmath12 and @xmath13 , can be diagonalized simultaneously . in fact , there is a set of common eigenvectors with eigenvalue zero , namely @xmath14 . these states are written in terms of the fock states @xmath15 of the field and the atomic states @xmath16 ( @xmath17 ) together with the atomic bell states @xmath18 the evaluation of the rest of the eigenvectors can be simplified by realizing that @xmath12 has a block - diagonal form in the basis @xmath19 with blocks given by the following matrices @xmath20 the eigenvalues of these matrices are given by @xmath21 for @xmath22 , @xmath23 for @xmath24 , and @xmath25 for @xmath26 , with @xmath27 the transformations that diagonalize each of the blocks @xmath28 are given by @xmath29 these matrices are the blocks of the orthogonal transformation @xmath30 that diagonalizes the hamiltonian @xmath12 as @xmath31 . having solved the eigenvalue problem for @xmath12 , it is now possible to evaluate the time - dependent state vector @xmath32 for any given initial pure state @xmath33 . in this work we consider as initial condition a normalized product state of the two atoms and the single - mode field that can be expressed as @xmath34 we have considered a general pure state @xmath35 of the atoms with probability amplitudes @xmath36 , @xmath37 and @xmath38 . for the single mode of the radiation field we have chosen a coherent state @xmath39 with mean photon number @xmath40 . using the eigenbasis of @xmath12 , the exact solution of the time - dependent state vector can be written as @xmath41 with the relevant photonic states @xmath42 and with the aid of the following abbreviations @xmath43      the time - dependent solution of the state vector can be significantly simplified in the case of high values of the mean photon number , i.e. , @xmath44 . in this limit @xmath45 , i.e. , the mean of the poisson distribution @xmath40 is much larger than the standard deviation @xmath46 . therefore we approximate @xmath47 and @xmath48 by @xmath49 and we also use the approximations @xmath50 the last line is obtained from the first - order taylor expansion in @xmath51 of the frequencies around @xmath52 . this is valid whenever the product between the second - order contribution times the interaction time @xmath53 remains small , a condition that is satisfied when @xmath54 @xcite . with these considerations and by introducing the abbreviations @xmath55 the photonic states can be simplified to @xmath56 where we have introduced the field states @xmath57 which are coherent states up to an additional phase . the description in terms of coherent states allows a simpler analysis of the dynamics . our aim is to prepare the atoms in an atomic postselection scenario where the atoms are prepared conditioned to a successful projection of the field onto the initial coherent state @xmath58 in a simplified and ideal implementation . in such a case , one would have to consider the following overlaps @xmath59 this result can be obtained by noting that the overlap between coherent states is given by @xmath60}\right|    \approx&e^{- g^2t^2 } ,    \nonumber    \label{}\end{aligned}\ ] ] which can be neglected if @xmath61 . therefore , after interaction with the resonator and projection onto state @xmath58 , both atoms are left in the state @xmath62 with @xmath63 success probability . the final state is actually a superposition of two states with probability amplitudes proportional to the intial ones . therefore , the atomic postselection can be understood as a projection of the atomic state with the following rank two projector @xmath64 where we have introduced the state @xmath65 . the operation @xmath66 represents the effective description of the interaction of the atoms with the resonator and the postselection via measurement of the field .      considering the projection onto a coherent state is an idealization that provides a convenient simplified picture . in practice , however , it is sufficient to project onto a state with vanishing overlap with the time - dependent field components @xmath67 and with finite overlap with @xmath58 . a typical experimental setting able to achieve this goal is a balanced homodyne measurement @xcite . the basic idea is to use a @xmath68 beam splitter to combine the field to be measured with a reference coherent field parametrized by its phase @xmath69 . photons from the two outputs of the beam splitter are collected using photodetectors . in the strong limit of the reference field @xcite , the probability of measuring a photocurrent difference between the detectors is proportional to the projection of the field onto the eigenstate @xmath70 of a field quadrature @xmath71 . this probability density for a coherent state @xmath58 is given by @xmath72 ^ 2\right\ }    \label{probcoh1}\ ] ] with @xmath73 . this can attain values close to unity by choosing the phase in such a way that @xmath74 and restricting values of @xmath75 close to zero . the square of the overlap with the other field components , that are also coherent states , can be evaluated as @xmath76 ^ 2\right\ } ,    \label{probcoh2}\ ] ] with @xmath77 . by choosing an appropriate interaction time @xmath53 , these overlaps can be made exponentially small . [ onestep ]     two two - level atoms in the same state @xmath78 interact with an optical resonator prepared in a coherent state @xmath58 . before the interaction , the gate @xmath79 is applied to one of the atoms and after the interaction this same atom is projected onto its ground state . the field is projected onto the initial coherent state . finally one of the atoms is left in the state @xmath80 . , scaledwidth=47.0% ]    in this section we use the atomic postselection scheme of the two - atom tavis - cummings model in order to implement an entangling quantum operation which by iteration leads to a nonlinear mapping of atomic probability amplitudes . the protocol is depicted schematically in fig . [ scheme ] . we consider the two two - level atoms initially prepared in a product state of the form ( @xmath81 ) @xmath82 for later convenience we have included the phase @xmath83 of the coherent state . before interacting with the optical resonator , a unitary gate @xmath84 is applied to atom @xmath1 . we choose the following gate @xmath85 which can be implemented by driving the atom with a resonant electromagnetic field and properly controlling the coupling and duration of the interaction @xcite . after the application of @xmath86 and before entering the resonator we get the following atomic probability amplitudes @xmath87 the probability amplitude @xmath88 does not need to be specified , as the resulting quantum operation projects the atoms onto a subspace orthogonal to @xmath89 as can be noted from eq . . with these initial conditions , both atoms interact with the electromagnetic field inside an optical cavity prepared in a coherent state @xmath58 . after the interaction a projection @xmath90 of the field onto the initial coherent state @xmath58 is performed and the atoms are left in the state @xmath91 the success probability of this projection is @xmath92 afterwards a projection @xmath93 onto the ground state of atom @xmath1 is implemented leaving atom @xmath0 in the state @xmath94 this event occurs with success probability @xmath95 . the overall success probability of the postselections is then given by @xmath96 the last inequality follows from analyzing eq . and noting that @xmath97 attains its minimum value when @xmath98 and @xmath99=-1 $ ] . up to normalization the final state is given by @xmath100 by iterating this procedure we attain a scheme implementing the following quantum map for the @xmath101th step @xmath102 with the complex functions @xmath103 the map is independent of the parameter @xmath83 , as one can note that the phase factor @xmath104 appears in the probability amplitude of state @xmath105 in the same manner as in the initial state @xmath106 of eq . . we note that the iteration of the map involves repeated action of the protocol on an ensemble of atoms . the protocol acts on a pair of identically prepared atoms from the ensemble and prepares one atom probabilistically . the other atom becomes useless from the point of view of the protocol , as a result of the projective measurement on it . after acting on all the atoms of the ensemble , one arrives at a smaller ensemble of less than one half in size . rapid downscaling of the ensemble size is a unavoidable condition for any quantum dynamics truly sensitive to initial conditions @xcite . in practice , realizing many steps of the protocol would require an exponentially large initial ensemble which would not be realistic . another practical aspect is that employing more than one cavity would be challenging with today s experimental possibilities . on the other hand , as we will demonstrate in the next sections , already a few steps can be enough to make highly overlapping initial quantum states almost orthogonal . furthermore we will outline an experimental proposal in sec . [ section - experimental ] with currently available technology by applying an optical conveyor belt and a single cavity . the dynamics within the approximations we have made is fully described by the iterative complex function in eq . ( [ map ] ) . this is a quadratic rational map @xcite , similar to the maps occurring in the measurement induced nonlinear quantum dynamical schemes first described in @xcite . in the following , we first carry out an analysis of the general properties of the iterated map @xmath107 of eq . ( [ map ] ) by using concepts from the theory of complex dynamical maps @xcite . then we compare its behavior to the numerical solution of the complete iterated dynamics , based on the hamiltonian of eq . ( [ hamilton ] ) and the subsequent selective measurements .       as a function the parameter @xmath108 . blue corresponds to the one - cycles @xmath109 , red corresponds to @xmath110 . dots , lines , and circles represent superattractive , attractive , and neutral cycles , respectively . numerical investigation of the enlarged regions between the neutral one - cycles show two different attractive 4-cycles ( orange lines ) and a single attractive 6-cycle ( green lines ) close to the two ends of the region . the central part of the enlarged regions contain `` islands '' of attractive @xmath111 cycles . the dotted circles indicate that is hard to identify the border of different regions . ] the periodic orbits or fixed cycles of the map @xmath107 can be determined from the relation @xmath112 . the one - cycles or fixed points as well as the 2-cycles can be determined analytically . for @xmath24 we find @xmath113 for @xmath114 , in addition to the above one - cycles , one can find two more points which are transformed into each other by @xmath107 . these form the single nontrivial two - cycle @xmath115    the stability of the fixed cycles can be checked by calculating the multiplier @xmath116 . a fixed cycle is repelling , neutral , attractive , or superattractive if @xmath117 , @xmath118 , @xmath119 , or @xmath120 , respectively . such an analysis can be carried out analytically for the one- and two - cycles , however , for @xmath121 it is a nontrivial task . the analysis of the multipliers shows that for each of the one - cycles there are certain parameter regions where they are attractive . on the other hand , the two - cycle given by eq . ( [ 2cyc ] ) , is repelling for any value of @xmath108 . for the determination of the longer ( @xmath1223 ) attractive cycles we can use the method based on the iteration of the critical points of the map . the critical points of @xmath107 are those which solve the equation @xmath123 . in this case , there are two critical points : @xmath124 a general theorem on iterated rational polynomial maps states that a rational map of degree @xmath125 can have at most @xmath126 attractive cycles . following the orbits of the critical points one can find all stable cycles of the iterated map ( in this case at most @xmath127 ) . [ fig2 ] shows where , according to the analytical calculations , the one - cycles are superattractive ( dots ) , attractive ( lines ) , and neutral ( circles ) as a function of the parameter @xmath108 . the numerical iteration of the critical points in the regions between the neutral one - cycles shows that there are two different attractive 4-cycles ( orange lines ) and a single 6-cycle ( green lines ) close to the two ends of the regions . the actual @xmath128 values belonging to the attractive 4- and 6-cycles depend on the parameter @xmath108 . in between these regions , it is numerically hard to rule out the existence of very long stable periodic orbits . the precision of our numerical simulation made it possible to identify a few `` islands '' of attractive fixed cycles of @xmath111 . the remaining part of this region may belong to maps without any stable periodic orbit , which means that all initial states belong to the julia set . the dotted circles indicate that the border between different regions is hard to determine numerically , which is an indication of the fractal nature of the regions . let us note that for @xmath129 and @xmath130 the map is actually not a genuine complex map since @xmath131 in these cases . the fractal nature of the map is more apparent when one determines the julia set of @xmath107 , i.e. the set of points which do not converge to an attractive cycle for a given @xmath108 . one way of numerically finding the points belonging to the julia set is backwards iterating the map starting from a point which is an element of a repelling cycle of the map . we show in fig . [ fig3 ] the julia set of @xmath107 for @xmath132 . in this case , the julia set is a totally disconnected set , all other initial points converge to the single attractive cycle @xmath133 , or physically speaking to the state @xmath134 . the analysis of the orbits of the critical points reveals important properties of the julia set . in this case both critical points converge to the same attractive fixed point , consequently the julia set is totally disconnected , similarly to the well - known cantor set @xcite .     for @xmath132 . ] another important case is when the two critical points converge to two distinct fixed points , then the julia set is connected . this case is illustrated by the map at parameter value @xmath135 shown in fig . [ cplane1 ] . for quadratic rational maps a general theorem ensures that the julia set is either totally disconnected , or connected @xcite . complex plane after @xmath136 iterations of the map in eq . for @xmath135 . two amplification levels are shown , confirming the fractal structure of the julia set separating the regions whose points converge to the attractive fixed points @xmath137 ( grey ) and @xmath138 ( black ) . the region indicated by the square in the middle of the left figure is magnified in the right figure . , title="fig : " ]   complex plane after @xmath136 iterations of the map in eq . for @xmath135 . two amplification levels are shown , confirming the fractal structure of the julia set separating the regions whose points converge to the attractive fixed points @xmath137 ( grey ) and @xmath138 ( black ) . the region indicated by the square in the middle of the left figure is magnified in the right figure . , title="fig : " ]      in order to investigate the real performance of the two - atom tavis - cummings model without the approximations of sec . [ model ] , we compute a numerically exact version of the operator @xmath66 in eq . . the matrix elements are evaluated as @xmath139 where we considered the atomic basis @xmath140 . the interaction time @xmath53 and coupling strength @xmath141 satisfy the relation @xmath142 . each iteration of the map is then evaluated by renormalizing the following outcome @xmath143 for qubit @xmath0 . in fig . [ cplane2 ] we plotted the real part of the @xmath136th iteration for two different values of the mean photon number @xmath40 , namely @xmath144 , and @xmath145 . with precision of two ( one ) decimal places the two fixed points also converge to + 1 and -1 in the case of @xmath146 ( @xmath147 ) . both figures reveal a fractal structure which resembles more the ideal case for larger values of @xmath40 .     the same as left part of fig . [ cplane1 ] for the the numerically exact quantum map and two values of the mean photon number @xmath40 : @xmath144 ( left ) and @xmath145 ( right ) . , title="fig : " ]   the same as left part of fig . [ cplane1 ] for the the numerically exact quantum map and two values of the mean photon number @xmath40 : @xmath144 ( left ) and @xmath145 ( right ) . , title="fig : " ] the number of atoms needed grows exponentially with the number of iterations even in an ideal case , which follows from the quantum magnification bound @xcite . in a realistic experiment , one can expect that only a few steps of the iteration can be carried out . we propose a test here , which would demonstrate the strong nonlinear behavior of the evolution even after 2 or 3 steps . a useful aspect of nonlinear quantum state transformations is that small initial differences between two similar quantum states can be amplified , enabling to distinguish them , realizing a schrdinger microscope , a term originally suggested by lloyd and slotine @xcite . let us consider the simple case when @xmath148 . in this case the unitary of eq . ( [ gate ] ) will become the well - known z gate . we choose two initial quantum states close to each other in the form of eq . ( [ atom - state ] ) with @xmath149 , and @xmath150 . their overlap is close to unity , @xmath151 . the two states will become almost orthogonal ( with a scalar product of @xmath152 ) after three steps of the iteration as can be seen in fig . [ numit ] where we show the plane of initial states colored according to the number of iterations needed for a complex number @xmath128 to reach either of the fix points @xmath137 or @xmath138 with a precision of @xmath153 . the julia set of the map is indicated by the yellow region on the figure , coinciding with the imaginary axis .     to reach either of the fix points @xmath137 or @xmath138 with a precision of @xmath153 . ]    in fig . [ scalprod ] we show the evolution of the overlap of the above two initial states using the complete solution with mean photon numbers @xmath154 and @xmath155 and by using the idealized map ( [ map ] ) . interestingly , the low - photon - number case leads to a faster decrease in the overlap during the first few steps of the iteration , but then converges to a larger value compared to the ideal map .     and @xmath156 after @xmath51 iterations of the ideal map ( blue squares ) , and the complete map with mean photon number 10 ( red circles ) and 100 ( black crosses ) . ] our basic protocol involves atomic and photonic postselection and therefore there is always a finite probability of failure . this means that in order to implement several iterations of the map , one requires several copies of the initial qubit pair . the procedure explained in sec . [ onestep ] has to be applied to every single copy of the ensemble . the number @xmath157 of qubit pairs required to achieve @xmath51 iterations , can be bounded from below by taking into account the success probability @xmath158 in . in addition , one has to take into account that half of the atoms in the ensemble are lost after being measured . therefore , the number of pairs scales exponentially as @xmath159 .    on a first thought one would naively consider the use of @xmath157 optical cavities for @xmath157 atomic pairs . however , there is another simpler solution motivated by current experimental implementations @xcite , where a standing - wave dipole trap or `` optical conveyor belt '' is used to coherently transport neutral atoms into an optical resonator .    using this setting , two conveyor belts are required to transport atoms into the cavity . in an initial stage , @xmath157 atoms are prepared in the minima of the two optical traps and are aligned as depicted in fig . [ schemebelt ] . for convenience , we number the atoms from left to right . the unitary gate @xmath160 is applied at this preparatory stage to atoms labeled with an even ( odd ) number in the upper ( lower ) conveyor belt , we call them marked atoms . the two conveyor belts are moved forward into the direction of the cavity until the first pair reaches the other side of the cavity . then , the conveyor belts stop in order to allow the measurement of the first marked atom and the field inside the cavity . afterwards , the cavity is reset to the state @xmath58 and the conveyor belts move again repeating the process . after all atoms have interacted with the cavity , the marked atoms are blacklisted as they are no longer useful . they are depicted in gray in fig . [ schemebelt ] . in order to pair only the useful atoms , the lower conveyor belt is shifted one period to the left , leaving the first marked atom without a partner . in this way , the potentially successfully prepared atoms are aligned . the process is repeated with both conveyor belts moving to the opposite side to start the second iteration . in the aforementioned implementation of the second iteration we have ignored the possibility of failure in the postselection . in order to overcome this problem , one has to keep track of successfully prepared atoms and then shift the conveyor belts in order to align useful pairs before transporting them into the cavity . possible implementation of the protocol using neutral atoms coherently transported using optical conveyor belts . , scaledwidth=47.0% ] we have proposed a non - linear map of qubit states in a cavity quantum electrodynamical scenario where the qubits are encoded in two - level atoms . the core step requires the interaction of two equally prepared atoms with the field inside an optical resonator according to the tavis - cummings model . by subsequent field detection and selective measurement of one of the atoms , the unmeasured atom is postselected into a state nonlinearly depending on its initial state .    from a mathematical point of view , we have studied the complex function describing this mapping of pure qubit states , where we have exploited the fact that any pure state of a qubit can be described by a complex parameter . we have performed an analysis of stable cycles under the iteration of the function and studied the behaviour in the complex plane . in particular , we have numerically investigated the julia set which changes from connected to disconnected for different parameters of the system . thus , our study offers a demonstration of chaotic behaviour in a quantum mechanical setting involving sequences of unitary transformations and postselective measurements .    from a physical perspective , we have proposed the realization of this scheme using an ensemble of equally prepared atoms in two optical conveyor belts that are coherently transported and interact in pairs with a single optical resonator . we have estimated the number of atoms required for each iteration of the protocol taking into account the success probability of the measurements involved . as an application we have discussed a schrdinger microscope implementation in which two initially close pure quantum states can be discriminated almost perfectly by a few iterations of our scheme .    although possible realizations of this nonlinear qubit map require cutting edge quantum technological developments , such as optical conveyor belts and controlled two - qubit interactions with a single mode radiation field , in view of the rapid experimental advances in cavity quantum electrodynamics its realization is within reach of nowadays technology . this work was supported by the hungarian academy of sciences ( lendlet program , lp2011 - 016 ) and the national research , development and innovation office ( k115624 , nn109651 ) and by the deutscher akademischer austauschdienst ( mb - daad project no . . o. k. acknowledges support from the j ' anos bolyai research scholarship of the hungarian academy of sciences . m. tavis and f. w. cummings , phys . rev . * 170 * , 279 ( 1968 ) . torres , j.z . bernd and g. alber , phys . a * 90 * , 012304 ( 2014 ) . torres , j.z . bernd and g. alber , appl . b * 122 * , 117 ( 2016 ) . a. i. lvovsky and m. g. raymer , rev . * 81 * , 299 ( 2009 ) . m. a. nielsen and i. l. chuang , _ quantum computation and quantum information _ ( cambridge university press , cambridge , uk , 2000 ) . d. meschede and a. rauschenbeutel , adv . at . . phys . * 53 * , 75 ( 2006 ) . m. khudaverdyan , w. alt , i. dotsenko , t. kampschulte , k. lenhard , a. rauschenbeutel , s. reick , k. schrner , a. widera and d. meschede , new j. phys . , * 10 * , 073023 ( 2008 ) . s. brakhane , w. alt , t. kampschulte , m. martinez - dorantes , r. reimann , s. yoon , a. widera and d. meschede phys . , * 109 * , 173601 ( 2012 ) . j. milnor and tan lei , experiment . math . * 2 * , 37 ( 1993 ) .
Abstract:

Answer: we propose a cavity qed scheme based on the tavis - cummings model to realize a conditional , non - linear dynamics exhibiting complex chaos in an ensemble of two - level atoms . the atoms interact pairwise with the coherent cavity field , which is subjected to homodyne measurement after the interaction . one member of the pair of atoms is then measured and upon condition that it is found in its ground state the other member of the pair is kept , otherwise discarded . the quantum state of the resulting ensemble of atoms is related to the state of the original ensemble via a non - linear transformation . by iterating this process with the help of an optical conveyor belt , for example , one arrives at complex deterministic chaos , where the evolution of pure quantum states is sensitive to the initial conditions . we uncover basic features of this dynamics by exploring stability regions in the space of initial states and analyze the occurrence of stable periodic orbits in dependence of the parameters of the model . it is demonstrated that sensitivity to initial conditions can be applied to amplify small differences of quantum states .


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: because of the potential for temperature gradients or other thermal stresses to cause electronic devices to fail , knowledge of the coefficients of thermal expansion ( cte ) of materials is important . either active materials with desired thermal properties or composites of active and compensating materials can mitigate the ill effects of thermal expansion in real devices.@xcite a compensating material must be chemically and electrically compatible with the functional material and will typically contract with increasing temperature ; i.e.  it will exhibit negative thermal expansion ( nte ) . while the cte for many materials have been cataloged , the set of known nte materials is small , and materials appropriate for specific applications may not yet be available . a means of rapidly predicting the thermal expansion properties of yet - to - be investigated materials will be extremely useful . the displacement of the equilibrium positions of atoms in a material with temperature is the source of thermal expansion . while anharmonicity in the vibrations of bonded pairs of atoms will drive the atoms farther apart with increasing temperature , crystal structure also plays an essential role in the specific thermal expansion characteristics of a material , providing behaviors over the range from positive thermal expansion ( pte ) to nte . accurate theoretical modeling of thermal expansion in materials necessitates quantum mechanical calculations and dynamical calculations or , minimally , quasi - harmonic modeling . as a result , typical calculations of thermal expansion address the microscopic causes of experimentally - observed cte on a case - by - case basis . especially prevalent recently are experimental and theoretical investigations of materials that exhibit nte , such as @xmath0 , @xmath1 ( with @xmath2 ) , @xcite @xmath3,@xcite and @xmath4.@xcite except for materials with very similar electronic and structural properties , a global picture capable of guiding searches for new materials with desired thermal expansion characteristics is slow to emerge . comprehensive reviews of thermal expansion in materials , emphasizing nte in both theory and experiment , can be found in refs .   and .    for microscopic atomic displacements to create nte , the motions of the atoms must carry them into spaces already existing within the lattice , while simultaneously drawing neighboring ions closer together . it is well - known@xcite that the definition of cte , @xmath5 can , with the help of a maxwell relation , be rewritten as @xmath6 which shows that for a material to exhibit nte , entropy must increase with pressure , contrary to typical expectations . decreasing the volume available to a particle , e.g.  by increasing pressure , is associated with decreasing entropy . the definition in eq . [ eq : cteentropy ] is isothermal which effectively means that the momentum space contribution to entropy changes negligibly relative to the real space contribution . from this we deduce that applying pressure in an nte material increases the volume available to its constituent atoms . on the other hand , applying pressure to pte materials ( at constant temperature ) results in decreasing entropy , implying a corresponding decrease in the volume available to the atoms follows from the same line of argument used above . in order for there to be more volume available to the atoms within a crystalline material under increasing pressure , atomic motions must be directed more significantly into the open spaces within the lattice . in other words , the atomic motions possess significant components in directions perpendicular to the bonds with nearest neighbor atoms . the difference between nte and pte is therefore quantifiable by the degree to which ionic thermal displacements are longitudinal or transverse with respect to the bonds . experimental evidence supporting this view is found , for example , in a study of nte in @xmath4 , in which inelastic neutron scattering shows that the sc - f bonds lengthen with increasing temperature and that the material contracts over a wide temperature range as a result of large transverse motions of the f ions.@xcite in our recent molecular dynamics investigation we demonstrated that thermal expansion in a single structure , with expanding bonds modeled with first- and second - neighbor interactions via morse potentials , can be varied from nte to pte by increasing the second - neighbor interaction strength relative to the first - neighbor interactions . @xcite by adjusting second - neighbor interactions , we reduced the transverse motions of the light ions , resulting in the emergence of pte in the model . we propose that an approach to predicting thermal expansion can be found by scanning the literature for the structures of crystals , focusing on quantities that may have a relationship to the entropy of the material and its potential to increase or decrease with respect to pressure , such as the space occupied by ions , their coordinations , and the distribution of atomic masses . we employ the wealth of structural information in databases such as the inorganic crystal structure database ( icsd ) @xcite and the crystallographic open database ( cod ) in this work . by correlating materials with known cte to their structures , we can begin to determine useful descriptors for thermal expansion . in section  [ section : proposeddescriptors ] , we discuss the physical underpinnings for the quantities that we propose for descriptors of cte . in section  [ section : results ] , we present the choice of descriptors from the original list of quantities of interest along with the result of performing a fit using the descriptors developed . we conclude with a discussion of the implications this correlation will have in the search for materials with desired thermal expansion properties . for each atom in a unit cell of each material , the number of nearest neighbors and formal oxidation state are determined . from this information we obtain ionic radii @xcite @xmath7 for each ion and an estimate of the volume each ion occupies ( assuming it occupies a sphere of volume @xmath8 ) . the fractional volume occupied is the total volume estimated for all the atoms in the unit cell divided by the volume of the unit cell @xmath9 : @xmath10 a complete list of the data used for this investigation is provided in the appendix[section : appendix ] . as seen in fig . [ fig : alphavolumeplot ] , thermal expansion data for the materials we sampled shows there is a relationship between thermal expansion and the volume occupied by atoms in the lattice ; more open space corresponds to a greater likelihood for nte .    the cluster of points with @xmath11/k in fig . [ fig : alphavolumeplot ] consists of binary materials that have the rock salt structure , with the corresponding high coordination providing an explanation for the large positive thermal expansion in these materials . we note that coulomb interactions between second neighbors will be repulsive and strong . the data point with the highest volume occupied and very low pte belongs to bn . because of its zincblende structure , the b and n atoms in bn have tetrahedral coordination . we also note that bn has a high bulk modulus , which implies that the tetrahedral bonds are strong . however , the low cte of bn suggests that ionic thermal motions are mostly directed into open spaces which are available because of low coordination . we deduce from this that bond coordination in a material is a possible descriptor to predict thermal expansion coefficients . finally , for some materials that exhibit negative thermal expansion , it has been noted that motion of light ions in the directions transverse to their bonds is at the root of their atypical thermal expansion . one underlying reason for this is that lighter ions undergo greater amplitude of motion than their heavier counterparts , which will assist the emergence of nte provided there exist open spaces and low coordinations around the light ions . therefore , we include atomic mass in the list of potential descriptors . ( color online ) the experimental coefficient of thermal expansion ( cte ) is plotted as a function of fractional volume occupied as defined in eq . [ eq : fractionalvolume ] for the materials gleaned from the literature . apparently , nte is possible only when @xmath12 is less than @xmath13 . the trend in the data is highlighted with a linear fit . the data with expansion coefficients greater than @xmath14 are ionic solids in the rock salt structure . the rms deviation between the data and the fit is @xmath15 and the correlation coefficient is 0.53.,width=326 ] quantities extracted in the search for descriptors , in addition to @xmath12 , are the mean bond coordination @xmath16 , its standard deviation @xmath17 , the mean atomic mass @xmath18 , the minimum mass @xmath19 , the maximum mass @xmath20 , and the standard deviation of mass @xmath21 . in fig .  [ fig : correlation ] , we display the linear correlations between descriptor candidates and the experimental cte . in addition , we include unitless ratios @xmath22 , @xmath23 , and the ratio of minimum to maximum mass @xmath24 . the quantities that have the strongest correlations with the cte are @xmath12 and @xmath16 , with correlations 0.53 and 0.57 , respectively . these two candidates correlate strongly with each other ; little is to be gained from using them together . the next - most - highly correlated quantity is @xmath22 , with a ( negative ) correlation of @xmath25 to the cte . furthermore , this quantity has a low degree of correlation with the volume ratio @xmath12 , making it a better quantity for using jointly with @xmath12 to predict cte . the candidates for mass have lower correlations with the cte . the ratio of minimum to maximum atomic mass correlation with cte is @xmath26 and cross - correlates only moderately with @xmath12 and @xmath22 . from our search for moderate - to - strong correlations with cte and weak - to - moderate cross - correlations , we determine @xmath12 , @xmath22 , and @xmath24 to be a good set of descriptors . these quantities have the additional advantages of being unitless and possessing values between zero and one for all materials in our list , simplifying analysis of their fitting strengths . ( color online ) the correlations ( pearson-@xmath27 ) between candidate descriptors and the coefficient of thermal expansion cte for 69 materials taken from the literature . moderate linear correlation with cte is found for both the mean coordination and occupied volume . weak - to - moderate correlation is exhibited by both the mean mass and the ratio of minimum to maximum atomic mass . cross correlations between descriptors permit identification of more nearly independent descriptors : the occupied volume , the ratio of the standard deviation to the mean coordination , and the ratio of minimum to maximum mass.,width=336 ]    we performed a second - order polynomial fit using these three descriptors and the 69 materials in our data sample , using the relation @xmath28 the correlation plot for this fit is displayed in fig . [ fig : multi - fit ] and the coefficients are listed in table  [ table : coefficients ] . the correlation coefficient is 0.80 and the rms - deviation between the fit and the data is @xmath29 , a noticeable improvement over the the fit displayed in fig . [ fig : alphavolumeplot ] , with outlying points more closely fit . points to the right of the diagonal are ones for which the model predicts thermal expansion greater than the experimentally measured value . the previously noted deviation in bn has been reduced by a factor of @xmath30 . the predicted values for the ionic solids have also been increased toward their experimental values . the material with largest positive deviation ( @xmath31 ) is @xmath32 . the largest negative deviation arises for @xmath33 . .[table : coefficients]the parameters used in eq . [ eq : fit ] to evaluate @xmath34 . [ cols="^,^,^,>",options="header " , ]     in order to visualize the model and to isolate ranges of the descriptors that characterize the cte , we present contours of constant @xmath35 across the ranges of @xmath36 and @xmath24 at fixed volume ratios @xmath37 in fig . [ fig : layers ] . superimposed on the contours are points indicating the descriptor values of data used in the fit , including volume ratios within @xmath38 of the value of @xmath12 displayed in the corresponding panel . as noted previously , materials with the lowest volume ratios are most - likely to have nte . in panels ( a ) and ( b ) in fig . [ fig : layers ] we see that the nte materials generally have mass ratios below 0.4 and coordination descriptors @xmath22 of between 0.2 and 0.5 . in addition to open space , important contributors to nte are atoms with a distribution of masses and a range of coordinations . the majority of materials in this work have volume ratios around 0.5  0.6 . these are displayed in panels ( c ) and ( d ) in fig . [ fig : layers ] . these materials also generally possess cte that would be considered typical . as seen for the low - volume - ratio materials , the lower values of @xmath35 fall in the middle of the ranges of the coordination and mass descriptors , suggesting nonlinear dependence . many of the materials here have zero coordination variance as a result of simple structures , such as rocksalt for which the atoms have identical coordination environments . a zero value in the coordination descriptor indicates the material has a high pte value . greater variance in the coordination indicates lower values of pte . the trend linking large pte with low coordination descriptor values and lower cte with larger values of the coordination descriptor continues into the most volume - crowded materials , which are seen in panels ( e ) and ( f ) of fig . [ fig : layers ] . many of these materials have zero values of the coordination descriptor , indicating simple , high - coordination structures . the materials with lower cte values either have mass ratios closer to 0.5 or have more complex bonding environments , as quantified by large coordination descriptor values . the trends here indicate that materials that can be considered typical , with cte @xmath39 , have volume ratios @xmath40 , coordination descriptors @xmath41 , and mass ratios @xmath42 . a coordination descriptor @xmath43 generally means large pte . nte requires low volume ratio and mid - range values for the other descriptors . the cte of bn and @xmath32 are both overestimated by the present description . in the case of bn , its zincblende structure means all atoms have tetrahedral coordination with zero variance , and its large ionic volume gives the material higher model cte . we note that a number of materials with zincblende or diamond lattice structures , such as si , show nte for a range of low temperatures.@xcite a relatively low resistance for b - n - b ( or si - si - si ) bond angles to change may account for the discrepancy between the model and experiment . the present model does not have a parameter that can access this particular characteristic . the reappearance of a zone of nte in the middle - right in panels ( e ) and ( f ) of fig . [ fig : layers ] is a result of including bn in the fit . removing bn from the fit eliminates this nte region and improves some values for the cte at the high end of the volume ratio . @xmath32 is in the mid - range of volume and coordination descriptors but its mass ratio is very low , @xmath44 , which in this model drives the predicted cte upward .    the cte of ki is underestimated . while it has the rocksalt structure , with the corresponding @xmath45 , its mass ratio is 0.3 and volume ratio is 0.59 , which cause the model to predict the lower value . the essential dynamical contribution to thermal expansion can not be fully captured by a static model , but we demonstrate that focusing on local properties that are linked to the dynamics , such as open space , bond coordinations , and atomic masses , provides a glimpse into these thermal motions for a range of different materials . the present model reveals trends in cte values without specific descriptors to characterize the flexibility of _ a - b - a _ bond angles to bend nor the rigidity of polyhedral atom groups , which are frequently cited as important contributors to thermal expansion . ( color online ) correlation plot of experimental and calculated cte generated by a two - parameter second - order polynomial fit to experimental data using eq . [ eq : fit ] with the parameters in table  [ table : coefficients ] . the root - mean - square deviation of the calculated coefficients of expansion , indicated with blue dashed lines , is @xmath46 . the pearson-@xmath27 for this fit is @xmath47.,width=326 ]     ( color online ) contour plots showing the dependence of the calculated @xmath35 upon @xmath36 and @xmath24 at fixed volume ratios @xmath12 , as indicated . the circles plotted display the descriptor values for the materials used in the fit . the volume ratios for the data points displayed in each panel fall within @xmath38 of the value of @xmath12 indicated in each panel.,width=355 ] the thermal expansion properties of solids emerge from the dynamics of ions and the bonds that join them . as a result , accurate theoretical prediction of thermal expansion is only possible with careful inclusion of the appropriate quantum mechanical considerations . these time - consuming calculations are a barrier to high - throughput searches of materials for candidates that possess particular thermal expansion characteristics . by encoding the relevant physics into simple descriptors , we have demonstrated that it is possible to approximate the thermal expansion coefficient of a variety of materials , including ionic binary compounds , perovskites , silicates , and other oxides by considering only the crystal structure of the material , the atomic masses , and the ionic radii of the atoms that compose it . these descriptors provide guidance for the search for new materials capable of mitigating thermal expansion in devices . additional features that may prove important to this investigation include properties of the electronic structure of materials in , for example , magnetically driven nte . table  [ table : alldata ] lists the sources of structures and cte used in this work . @xmath4 & & 261068 & @xmath48@xcite & @xmath49 + @xmath50 ( high @xmath51 ) & & 418846 & @xmath52@xcite & @xmath53 + @xmath54 ( high @xmath51 ) & 1008795 & & @xmath55@xcite & @xmath56 + @xmath57 ( high @xmath51 ) & & 80448 & @xmath58@xcite & @xmath59 + @xmath60 & & 50941 & @xmath61@xcite & @xmath62 + @xmath63 & & 391467 & @xmath64@xcite & @xmath65 + @xmath66-quartz & & 89289 & @xmath67@xcite & @xmath68 + @xmath69 & 9011469 & & @xmath70@xcite & @xmath71 + @xmath72 & 9014631 & & @xmath73@xcite & @xmath74 + @xmath75 & 9014166 & & @xmath76@xcite & @xmath77 + @xmath78 & 9008834 & & @xmath79@xcite & @xmath80 + @xmath81 & 1522201 & & @xmath82@xcite & @xmath83 + @xmath84 & & 54697 & @xmath85@xcite & @xmath86 + @xmath87 & 9007325 & & @xmath88@xcite & @xmath89 + @xmath90 & 1010168 & & @xmath91@xcite & @xmath92 + @xmath93 & 1011265 & & @xmath94@xcite & @xmath95 + @xmath96 & 4119149 & & @xmath97@xcite & @xmath98 + @xmath99 & 9006315 & & @xmath100@xcite & @xmath101 + @xmath102 & 9006111 & & @xmath103@xcite & @xmath104 + @xmath32 & 1528988 & & @xmath103@xcite & @xmath105 + @xmath106 & 9007310 & & @xmath107@xcite & @xmath108 + @xmath109 & 9005857 & & @xmath110@xcite & @xmath111 + @xmath112 & 9007701 & & @xmath113@xcite & @xmath114 + @xmath115 & 9001685 & & @xmath116@xcite & @xmath117 + @xmath118 & 1531195 & & @xmath119@xcite & @xmath120 + @xmath121 & 9003592 & & @xmath122@xcite & @xmath123 + @xmath124 & 9010061 & & @xmath125@xcite & @xmath126 + @xmath127 & 9009768 & & @xmath125@xcite & @xmath128 + @xmath129 & 9007485 & & @xmath130@xcite & @xmath126 + @xmath131 & 9015790 & & @xmath132@xcite & @xmath133 + @xmath134 & 9004744 & & @xmath135@xcite & @xmath111 + @xmath136 & 9000958 & & @xmath137@xcite & @xmath138 + @xmath139 & 9003990 & & @xmath140@xcite & @xmath141 + @xmath142 & 9015662 & & @xmath143@xcite & @xmath144 + @xmath145 & 9005861 & & @xmath146@xcite & @xmath147 + @xmath148 & 9000344 & & @xmath149@xcite & @xmath150 + @xmath151 & 9002686 & & @xmath152@xcite & @xmath153 + @xmath154 & 1528398 & & @xmath155@xcite & @xmath156 + @xmath157 & 9000526 & & @xmath158@xcite & @xmath159 + @xmath160 & 9008279 & & @xmath161@xcite & @xmath162 + @xmath163 & 1011033 & & @xmath164@xcite & @xmath165 + @xmath166 & 9006864 & & @xmath167@xcite & @xmath168 + @xmath169 & 1010299 & & @xmath167@xcite & @xmath170 + @xmath171 & 1524203 & & @xmath172@xcite & @xmath173 + @xmath174 & 9000904 & & @xmath175@xcite & @xmath176 + @xmath177 & 9010486 & & @xmath178@xcite & @xmath179 + @xmath180 & 9007121 & & @xmath181@xcite & @xmath182 + @xmath183 & 9007534 & & @xmath184@xcite & @xmath185 + @xmath186 & 1011064 & & @xmath184@xcite & @xmath187 + @xmath188 & 9009766 & & @xmath189@xcite & @xmath126 + @xmath190 & 1514118 & & @xmath191@xcite & @xmath192 + @xmath193 & 1530736 & & @xmath194@xcite & @xmath195 + @xmath196 & 7200689 & & @xmath197@xcite & @xmath198 + @xmath199 & 9006717 & & @xmath200@xcite & @xmath201 + @xmath202 & 9006758 & & @xmath203@xcite & @xmath117 + @xmath204 & 9013383 & & @xmath205@xcite & @xmath206 + @xmath207 & 9011666 & & @xmath208@xcite & @xmath209 + @xmath210 & & 181799 & @xmath211@xcite & @xmath162 + @xmath212 & 9010943 & & @xmath211@xcite & @xmath213 + @xmath214 & 9011682 & & @xmath215@xcite & @xmath216 + @xmath217 & 9007457 & & @xmath218@xcite & @xmath219 + @xmath220 & 8103714 & & @xmath218@xcite & @xmath221 + @xmath222 & 9008652 & & @xmath223@xcite & @xmath224 + @xmath225 & & 165593 & @xmath226@xcite & @xmath221 + @xmath227 & & 165592 & @xmath228@xcite & @xmath229 + @xmath230 & 9008665 & & @xmath231@xcite & @xmath232 + @xmath33 & 9009735 & & @xmath233@xcite & @xmath234 + @xmath235 & 9008681 & & @xmath236@xcite & @xmath237 + @xmath238 & 9008669 & & @xmath239@xcite & @xmath240 +    this work has been supported by the department of energy office of basic energy sciences , under grant number de - fg02 - 07er46431 . computational support was provided by the national energy research scientific computing center ( nersc ) . 51ifxundefined [ 1 ] ifx#1 ifnum [ 1 ] # 1firstoftwo secondoftwo ifx [ 1 ] # 1firstoftwo secondoftwo `` `` # 1''''@noop [ 0]secondoftwosanitize@url [ 0 ]  + 12$12  & 12#1212_12%12@startlink[1]@endlink[0]@bib@innerbibempty link:\doibase 10.3390/ma5061125 [ * * ,   ( ) ] link:\doibase    10.1039/c4cs00461b [ * * ,   ( ) ] link:\doibase    10.1126/science.272.5258.90 [ * * ,   ( ) ] ,   @noop * * ,   ( ) link:\doibase 10.1103/physrevlett.80.4903 [ * * ,   ( ) ] http://dx.doi.org/10.1038/24115 [ * * ,   ( ) ] link:\doibase 10.1103/physrevb.60.7234 [ * * , ( ) ] link:\doibase 10.1103/physrevlett.86.4692 [ * * ,   ( ) ] link:\doibase 10.1103/physrevlett.89.215902 [ * * ,   ( ) ] link:\doibase 10.1103/physrevb.68.014303 [ * * ,   ( ) ] link:\doibase    10.1103/physrevlett.93.225501 [ * * ,   ( ) ] link:\doibase    10.1103/physrevlett.95.255501 [ * * ,   ( ) ] link:\doibase 10.1103/physrevlett.109.195503 [ * * ,   ( ) ] link:\doibase    10.1103/physrevlett.112.045505 [ * * ,   ( ) ] link:\doibase 10.1021/cm501107w [ * * ,   ( ) ] ,   \doibase    http://dx.doi.org/10.1063/1.4867437 [ * * , ( ) ] link:\doibase    10.1103/physrevb.86.094110 [ * * ,   ( ) ] link:\doibase 10.1063/1.3155191 [ * * ,   ( ) ] link:\doibase    10.1103/physrevlett.107.195504 [ * * ,   ( ) ] link:\doibase    10.1103/physrevb.92.224302 [ * * ,   ( ) ] http://stacks.iop.org/0953-8984/17/i=4/a=r03 [ * * ,   ( ) ] link : www.scopus.com [ * * ,   ( ) ] link:\doibase 10.1103/physrevb.93.214304 [ * * , ( ) ] @noop _ _ ,  edited by  , ,  and  ( , ) link:\doibase 10.1107/s0108768102006948 [ * * ,   ( ) ] link:\doibase    10.1107/s1600576715022396 [ * * , ( ) ] link:\doibase 10.1107/s1600576714025904 [ * * ,   ( ) ] @noop * * , ( ) link:\doibase    10.1107/s0021889809016690 [ * * , ( ) ] @noop * * ,   ( ) link:\doibase 10.1107/s0567739476001551 [ * * ,   ( ) ] link:\doibase 10.1103/physrev.163.779 [ * * , ( ) ] link:\doibase    10.1021/ja106711v [ * * ,   ( ) ] \doibase http://dx.doi.org/10.1016/s0925-8388(01)02003-5 [ * * ,   ( ) ] link:\doibase 10.1021/acs.chemmater.5b00662 [ * * ,   ( ) ] ,   @noop * * ,   ( ) @noop * * ,   ( ) link:\doibase 10.1007/s002690050161 [ * * ,   ( ) ] `` , ''  ( ,  )  chap . \doibase    http://dx.doi.org/10.1016/j.jct.2003.12.010 [ * * ,   ( ) ] \doibase http://dx.doi.org/10.1016/0025-5408(71)90080-8 [ * * ,   ( ) ] \doibase http://dx.doi.org/10.1063/1.321373 [ * * ,   ( ) ] @noop * * , ( ) @noop * * ,   ( ) link:\doibase 10.1103/physrevb.53.3013 [ * * , ( ) ] link:\doibase 10.1016/j.jssc.2004.07.036 [ * * ,   ( ) ] link:\doibase    10.1007/s10973 - 013 - 3421 - 8 [ * * ,   ( ) ] \doibase    http://dx.doi.org/10.1063/1.3672198 [ * * , ( ) ] @noop * * ,   ( ) link:\doibase 10.1007/s10765 - 015 - 1881 - 8 [ * * , ( ) ] @noop * * , ( )
Abstract:

Answer: thermal expansion in materials can be accurately modeled with careful phonon calculations within density functional theory . however , because of interest in controlling thermal expansion and the time consumed evaluating thermal expansion properties of candidate materials , either theoretically or experimentally , an approach to rapidly identifying materials with desirable thermal expansion properties would be of great utility . we show that the fraction of crystal volume occupied by ions , based upon ionic radii , the deviation of bond coordination divided by mean bond coordination , and the ratio of least to greatest atomic mass are descriptors that correlate with the coefficient of thermal expansion for a variety of materials found in widely accessible databases . correlation is greatly improved by combining these descriptors in a multi - dimensional quadratic fit . open space combined with a range of atomic masses and a variety of local bond coordinations indicates materials with lower coefficients of thermal expansion . materials with single - valued local coordination and less open space have the highest coefficients of thermal expansion .