ca@xmath0ruo@xmath1 ( cro ) is the mott insulating@xcite end member of the series ca@xmath2sr@xmath3ruo@xmath1 , which possess a rich diversity of structural , magnetic and transport properties@xcite .
sr@xmath0ruo@xmath1 , the other end member , is a spin - triplet superconductor with broken time - reversal symmetry @xcite .
the metal - insulator transition in cro at @xmath4  k @xcite goes along with severe structural distortions @xcite , the ruo@xmath5 octahedrons flattens , which is also apparent in jumps of the lattice constants in the order of 0.1   .
these structural distortions are enhanced till the onset of antiferromagnetic order @xcite .
the nowadays widely used picture assumes that an orbital ordering is associated with the structural changes @xcite .
the 4@xmath6 orbitals become doubly occupied and the 4@xmath7 singly occupied resulting in flattened octahedrons and a @xmath8 state . in the past the nature of the mott transition of this multiband system with four @xmath9 electrons on the ru site
was intensively discussed @xcite .
recently , it was proposed , that spin - orbit coupling ( soc ) in this 4@xmath9 system is strong enough to change the multiplet structure and couples @xmath10 and @xmath11 to @xmath12 resulting in a nonmagetic @xmath13 ground state .
the occurring of a magnetic order was proposed to be of a singlet magnetism type ( see e.g. sec .
5.5 . in ref .
@xcite ) , which was called excitonic magnetism in ref .
@xcite , and a special type of magnon dispersion was predicted in this theory @xcite .
the main branches of the dispersion obtained by inelastic neutron scattering ( ins ) experiments , however , could be successfully described with a conventional heisenberg model @xcite and strongly disagree with @xmath13 calculations @xcite . nevertheless , neutron scattering experiments revealed features and additional magnetic scattering not describable with a simple square - lattice antiferromagnet .
the main in - plane transverse modes exhibit a large anisotropy gap , underlining the impact of broken tetragonal symmetry in combination with soc .
furthermore , an additional signal was detected in the neutron scattering experiments at low and at high energies @xcite , that can not be explained by the two transverse magnon branches .
here we present a study of the magnon dispersion using ins with and without polarization analysis on crystals containing 0 , 1 , and 10  % ti combined with density functional theory ( dft ) calculations .
the full dispersion of magnetic excitations in cro reveals a very strong impact of soc , and oxygen moments that lead to an additional branch .
( color online ) crystal and magnetic structure of cro and the associated reciprocal space . in ( a ) one layer of the crystal and magnetic structure
the ru ( blue balls ) sit in octahedrons of oxygen ( red balls ) .
the ab - plane is indicated in light blue . in black and green the orthorhombic and tetragonal cells are shown .
the pink arrows denote the magnetic moments on the ru and o sites . at the tip of the low right ru moment
, the directions of the different polarizations of the magnon modes are indicated by black arrows and some labellings of atoms and atom distances are given . in ( b ,
c ) the ( hkl ) planes in reciprocal space are illustrated for l=0 and l=1 , respectively .
the orthorhombic and tetragonal cells are drawn in black and red , respectively .
green circles denote tetragonal zone centers and pink circles zone boundaries .
black squares ( stars ) mark antiferromagnetic zone centers of the @xmath14- ( @xmath15)-centered antiferromagnetic order . ]
crystal growth of insulating ca@xmath2sr@xmath3ruo@xmath1 is severely hampered by the metal - insulator transition occurring in pure cro slightly above room temperature , @xmath4  k. the space group does not change at this first - order phase transition , but there are sizable jumps in the lattice - parameters , in particular for @xmath16 , with @xmath17  @xcite . therefore , upon cooling the crystals after the growth process , they tend to crack and only small pieces of mm@xmath18 size can be recovered .
because inelastic neutron scattering experiments require samples of @xmath19cm@xmath18 size , we substituted a small amount of ru by ti in order to broaden the phase transition and to be able to cool the crystals through the metal - insulator transition without breaking them into small pieces .
crystals with a mass of 0.5  g to 1  g were obtained with only 1  % ti substitution .
we were also able to obtain two larger pieces of pure cro , one with a mass of 0.45  g and one with 0.3  g , but these crystals posses a bad mosaic spread of 5@xmath21 , which most likely is the reason why they were not destroyed upon cooling .
nevertheless , these crystals are suitable for some inelastic neutron scattering studies .    the ca@xmath0ru@xmath3ti@xmath22o@xmath1 crystals with @xmath23% used in this neutron scattering study ( labeled 0ti , 1ti , 10ti , respectively )
were grown by the floating - zone method in a canon machinery inc .
sc1-mdh11020-ce furnace equipped with two 2000  w halogen lamps and a coldtrap following the procedure described in references @xcite .
caco@xmath24 , ruo@xmath0 and tio@xmath0 were mixed in stoichiometric ratios and a ru excess of 32  % was added .
the powder was mixed and reacted for 24  h with an intermediate grinding .
then a rod was pressed and sintered at 1350  @xmath21c . a growth speed of 17  mm / h and a feed - rod speed of 20  mm / h were used and the atmosphere contained 90  % ar and 10  % o at a total pressure of 9  bar .
phase purity was checked by x - ray powder diffraction from crushed parts of the single crystals , which indicates no impurity phases .
the lattice constants at room temperature were obtained by lebail fits of these powder pattern using the fullprof suite@xcite .
single - crystal x - ray diffraction experiments were performed on a bruker x8 apex diffractometer equipped with a charge - coupled - device detector .
details of this structure determination and a detailed list of the refined parameters can be found in the appendix a. magnetization measurements were performed with a superconducting quantum interference device ( squid ) magnetometer from quantum design , and the resistivity was measured with a standard four - point method .
spin wave calculations were performed using spinw @xcite .
elastic and inelastic neutron scattering experiments were performed with the cold triple - axis spectrometer ( tas ) 4f1 at the llb , with the thermal tas in3 at the ill and with the polarized thermal tas in22 at the ill .
the polarized neutron scattering experiments were performed using heusler ( 111 ) monochromator and analyzer crystals .
a pg filter was inserted on the scattered beam before the neutron - spin flipper and the monitor was put on the incoming beam between monochromator and sample . a set of helmholtz coils was used to produce the guide field and the sample was zero field cooled ( less than 2  g ) in an orange - type cryostat . a mounting of 3 crystals with a total mass of 2.5  g containing 1  % ti ( 1ti )
was used for the polarized neutron scattering experiment .
the sample was twinned with non - equal twinning - fractions of 2.5:1 determined by scanning the orthorhombic ( 200 ) and ( 020 ) reflections .
the @xmath25 plane [ see fig .  [ explanations](a ) ] was chosen as the scattering plane in order to efficiently integrate the inelastic signal along the vertical direction , where the resolution is poor .
we used the standard coordinate system in polarized neutron scattering @xcite .
* x * is parallel to the scattering vector ( * q * ) , * y * and * z * are perpendicular to * q*. while * y * is in the scattering plane , * z * is perpendicular to it .
therefore , * y * lies in the @xmath25  plane of the crystal and * z * parallel to the @xmath16  axis . in neutron scattering only magnetic moments perpendicular to * q * intervene .
therefore , the magnetic scattering intensities sense a geometry factor sin@xmath26(@xmath27 ) with @xmath27 being the angle between * q * and the magnetic moment , which corresponds to either the static ordering moment in diffraction or to the oscillating moment in a magnon . with longitudinal neutron polarization six different channels
can be analyzed : three spin - flip ( sf ) and three non - spin - flip ( nsf ) channels .
while phonons always contribute to the nsf channel there is an additional selection rule for magnetic scattering .
the magnetic component parallel to the direction of the neutron polarization contributes to the nsf channel , while the components perpendicular to the neutron polarization generate sf scattering . by combining this polarization rule with the geometry factor one can distinguish the different magnetic components . in the geometry we use for the polarized neutron scattering experiment
we see the @xmath16 polarized modes in the sfy channel and they do not loose intensity due to the geometry factor , because the scattering vector is always parallel to the @xmath25 plane and thus perpendicular to @xmath16 .
in contrast the sfz channel contains the in - plane modes , the transverse and longitudinal ones , but weighted with the geometry factor .
for example at @xmath28 the geometry factor for the transverse mode is sin@xmath29 and that for the longitudinal mode sin@xmath30 . at @xmath31
this ratio is inverted [ see fig .
[ explanations](b , c ) ] . respecting the twinning ratio of 2.5:1
, we expect that the intensities at @xmath28 and @xmath31 have a ratio of 0.6 for the transverse mode . for the longitudinal mode
this ratio is inverted .
thus , it is possible to distinguish between transverse modes and longitudinal modes by comparing scattering at properly chosen * q*. because the polarization of the neutron beam is not perfect , one has to correct the intensities for the finite flipping - ratio ( fr ) @xcite .
thereby , the magnetic signals are obtained from the intensities of the different sf channels corrected for the fr : @xmath32\ ] ] the fr of our experiment on in22 amounts to 12 , which is determined by comparing the signals of rocking scans on the ( 200 ) bragg - reflection in the sf and nsf channels .    on the cold tas
4f1 crystals with different ti content were used , 0ti and 10ti .
0ti has a mass of 0.45  g and a twinning ratio of 1:1 and 10ti a mass of 0.93  g and a twinning ratio of 9:1 .
the scattering plane for both samples was chosen to be the [ 010]/[001]-plane .
so the scattering plane for the other twin domains was the [ 100]/[001]-plane .
for all 4f1 scans a pyrolytic graphite monochromator and analyzer were used , a cooled be - filter was put on * * k**@xmath33 to suppress higher order contaminations .
for the cooling of the samples an orange - type cryostat was used and @xmath34 was set to 1.55  @xmath35 for all scans .    with the in3 tas we analyzed the magnetic order of the crystal containing 1% ti used in the previous neutron scattering study@xcite ( 1tib ) .
this crystal was found to exhibit a majority twin domain of 95% and the experiment was performed in the [ 010]/[001 ] orientation used in the inelastic neutron scattering experiment .
( color online ) temperature dependence of physical properties of cro with different ti content . ( a )
field - cooled magnetization curves in 0.1  t measured on heating .
( b ) resistivity on heating in dark color and on cooling in bright color .
only the 10  % ti curves do not coincide . ]
the crystal structure of cro @xcite is similar to the structure of sr@xmath0ruo@xmath1 of k@xmath0nif@xmath1 type @xcite , where the ru atoms sit in oxygen - octahedron cages , which are corner shared in the @xmath25  plane . in cro
the octahedrons are rotated and tilted and they also become severely distorted . as a consequence , the tetragonal symmetry of sr@xmath0ruo@xmath1
is reduced to the orthorhombic space group pbca@xcite .
one layer of the crystal structure is drawn in fig .
[ explanations ] , which also depicts the magnetic order with antiferromagnetic moments on ru pointing along @xmath36 .
the orthorhombic unit cell is rotated by @xmath37 with respect to the tetragonal one and enhanced to : @xmath38 , @xmath39 , @xmath40 . unless otherwise specified all notations refer to the orthorhombic lattice of the majority twin orientation .
our crystals usually posses two twin orientations , which are obtained by interchanging the @xmath41 and @xmath42 directions .
all crystals were examined by x - ray diffraction , magnetization and resistivity measurements in addition to the neutron scattering experiments described below .
[ magnetization ] shows the magnetization and electric resistance data of the three ti concentrations .
the temperature of the metal - insulator transition is 4  k lower for the sample containing 1  % ti than for the pure compound .
the sample containing 10  % ti does not show indications for a phase transition upon cooling down to 80  k where the experimental limit of the increasing resistivity is reached .
the absolute values of the resistivity curves have a large uncertainty because of the first - order structural phase transition .
there the crystals tend to crack , which prohibits the current to flow through the hole sample .
this effect has been frequently observed during a single measurement cycle .
the resistivity is enhanced by a multiplication factor after passing the structural transition . as the samples have already passed this transition after the crystal growth
, there are some cracks inside the sample , which is evident from measuring several pieces of the same crystal without obtaining the same room - temperature specific resistivity .
the magnetic and insulating properties in ca@xmath2sr@xmath3ruo@xmath1 are closely related to the crystal structure .
cro is heavily distorted with respect to sr@xmath0ruo@xmath1 , which possesses the ideal structure of k@xmath0nif@xmath1 type without structural distortions , but which already is close to such a structural instability@xcite . in the layered
ruthenates the distortions can be described as octahedron rotation around the @xmath16  direction and tilting around an in - plane axis .
these distortions and the associated structural phase transitions result from bond - length mismatch , and the distortions in cro are induced by the chemical substitution of sr by the isovalent , but smaller ca , which is not able to fill the space between the oxygen octahedrons like sr does .
so the octahedrons start to rotate and tilt in order to reduce the coordination volume around ca . by chemically substituting the smaller ti for ru , one expects a small decrease of these deformations , which is indeed realized .
the tilt and rotation angles of 10ti are significantly smaller than the values determined in the samples containing only 0 and 1  % ti , see appendix a , tab .
[ tab : char ] .
temperature and ti doping dependence of the ruo@xmath5 octahedron distortion obtained by single crystal x - ray diffraction analyses ; ( a , b ) presents the ratios of in - plane to out - of - plane ruo bond distances ( flattening ) and of the two octahedron edges parallel to @xmath36 and @xmath43 , respectively . both entities are enhanced in the insulating phase at low temperature due to the lifting of orbital degeneracy . ]
most interesting is the impact of the ti doping on the deformation of the ruo@xmath5 octahedron , which is a fingerprint of the lifting of orbital degeneracy@xcite .
the structural change at the metal - insulator transition in pure cro is characterized by a jump of the @xmath16  lattice - constant and by a flattening of the octahedron as it is visible in the ru - o bond distance ratio@xcite . upon cooling into the insulating phase this
octahedron flattening continues until it saturates at the antiferromagnetic phase transition .
slightly below room temperature the octahedron shape passes from elongated to flattened @xcite .
in all samples studied here , this crossover occurs below room temperature , thus at 100  k the octahedron shape is flattened for all ti substitutions , see fig .
[ ruooo ] . due to flattening of the lattice and octahedron the @xmath6 orbital shifts down in energy compared to the @xmath44,@xmath45 orbitals .
a similarly strong and anomalous temperature dependence is also observed for the ratio of the two o - o octahedron edge lengths parallel @xmath43 and @xmath36 which is related to the orthorhombic splitting [ @xmath46 @xcite . at room temperature
the octahedron is longer along @xmath43 , while it is elongated along @xmath36 at low temperature .
all these effects can be attributed to a temperature dependent orbital ordering @xcite .
the considerable elongation along @xmath36 at low temperature agrees with soc and the alignment of the magnetic moment mainly along the @xmath36 direction @xcite .
figure [ ruooo ] illustrates that these two distortions of the ruo@xmath5 octahedron become suppressed by the 10  % ti substitution .
this underlines the orbital ordering character of these distortions , which is obviously suppressed by ti@xmath47 with an empty @xmath48 shell .    while the @xmath25  plane ( or the average in - plane parameter ) increases with high ti content , the @xmath16  lattice - constant decreases at room temperature .
the reduction of the @xmath16  lattice - constant , in hand with an increase of the @xmath25  plane , is the structural signature of the insulating state of cro . with 1  % of ti substitution
there is only a small decrease in @xmath49 ( fig .
[ magnetization ] ) , with 10  % of ti substitution there are drastic effects . in the resistivity of 10ti
there are no indications for a sharp metal - insulator transition on cooling from 700 down to 80  k , where the upper experimental limit of the resistance experiment is reached .
note that the absolute values of the resistivity curves are prone to a large uncertainty because of the first order structural phase transition which causes cracks in the crystals .
but the orders of magnitude larger resistivity of the 10ti sample and the absence of the metal - insulator transition is unambiguous .
this remarkable stabilization of the insulating state by only small amounts of ti substitution is also seen in ca@xmath24ru@xmath0o@xmath50 @xcite .
since the ionic radii of ti and ru are very similar , this suppression of the metallic state seems to originate from the very effective suppression of the hopping .
also in sr@xmath0ruo@xmath1 ti substitution has a strong impact : it stabilizes spin - density wave ordering associated with the fermi - surface nesting of the pure material for only 2.5%  ti @xcite .
also in sr@xmath24ru@xmath0o@xmath50 ti substitution stabilizes a spin - density wave magnetic instability @xcite .
( color online ) elastic neutron scattering scans with polarization analysis .
the orthorhombic splitting is clearly visible in this elastic scans .
( a , b ) present the values of the nsfx ( sfx ) channel in magenta ( blue ) of longitudinal scans across the strong nuclear reflections ( 200 ) and ( 020 ) , respectively .
the values of the nsfx channel are divided by the fr .
these elastic longitudinal scans show that both twins are well separated .
( c - f ) present data of the sfx ( blue ) , sfy ( red ) and sfz ( black ) channels of scans across the magnetic ( c ) ( 100 ) , ( d ) ( 010 ) , ( e ) ( 210 ) , and ( f ) ( 2 - 10 ) reflections . ]    the magnetic structure of 1ti was studied on in22 with polarization analysis .
the analysis of the magnetic structure of 1tib on in3 is described in appendix b. scans along the ( 200 ) and ( 020 ) reflections ( fig .
[ polel ] ) reveal that crystal 1ti exhibits both twin - domain orientations in the ratio 2.5:1 .
[ polel](c - f ) show the sf intensities in the three channels for scans across the ( 100 ) , ( 010 ) , ( 120 ) and ( -120 ) magnetic reflections for crystal 1ti .
the magnetic structure of cro has been previously determined by neutron powder diffraction @xcite . as illustrated in fig .
[ explanations](a ) the magnetic moments are essentially aligned along the @xmath36  direction ( corresponding to the tilt axis ) with antiferromagnetic alignment between nearest neighbors . due to the low symmetry of the crystal structure ( pbca )
the dzyaloshinski - moriya interaction is finite and results in a canting of the moments along the @xmath43 direction .
the canting of the moments along @xmath43 is parallel for all ru sites in a single layer resulting in a net ferromagnetic moment per layer .
there are two different magnetic structures reported for cro @xcite , which differ by the stacking of the single layer arrangement shown in fig .
[ explanations](a ) .
the magnetic structure of the main antiferromagnetic @xmath36 component is either @xmath14 or @xmath15-centered . in the @xmath14-centered phase the two magnetic moments at ( 0,0,0 ) and ( 0,0.5,0.5 )
are parallel , in the @xmath15-centered one the two moments at ( 0,0,0 ) and ( 0.5,0,0.5 ) are parallel .
the magnetic space groups are @xmath51 ( @xmath14-centered ) and @xmath52 ( @xmath15-centered ) .
while in @xmath51 the net ferromagnetic canted moments per layer cancel due to an antiferromagnetic stacking , the @xmath15-centered @xmath52 structure results in a total ferromagnetic moment , that can be measured with a magnetometer .
there is also canting along @xmath16 direction in both magnetic structures .
this canting corresponds to an antiferromagnetic @xmath16 component that is @xmath15-centered in @xmath51 and @xmath14-centered in @xmath52 .
this moment should , however , be small because the main part of the ordered moment points along the tilt axis , and only microscopic methods can detect such an antiferromagnetic @xmath16 component .
the analysis of the magnetic bragg peaks with polarization analysis first confirms that moments point along the @xmath36 direction and show a dominating @xmath15-centered scheme for 1ti and an almost exclusive @xmath15-centered scheme for 1tib , see appendix b. in contrast , small pure cro samples show only the @xmath14-centered scheme @xcite .    in fig .
[ magnetization](a ) the magnetization on heating of samples with different ti content is presented .
the magnetic field of 0.1  t is applied in the @xmath53 plane .
the weak ferromagnetic component parallel to the @xmath43 direction , resulting from the canting of the spins , dominates the magnetization below @xmath54 .
the fact that the magnetization of the lowest temperature is the highest , points to the @xmath15-centered phase , in which the ferromagnetic component does not cancel among two neighboring layers @xcite .
thus all samples contain a dominating @xmath15-centered phase in agreement with the neutron diffraction studies .
the small crystal used for the magnetization measurement with 1  % ti content is essentially untwinned , while the sample containing no ti is partially twinned with nonidentical twin fractions , in contrast to the much larger pure sample used for neutron scattering , for which the twinning ratio amounts to 1:1 .
if the field is applied along [ 1,1,0 ] , the observable ferromagnetic component is reduced by a factor @xmath55 but both twin domain orientations contribute ( sample containing 10% ti ) .
the highest total ferromagnetic component is observed in the pure sample , and the sample 1ti already exhibits a reduced ferromagnetic component possibly stemming from a slightly reduced fraction of the @xmath15-centered phase . we can conclude that cro exhibits sizable moment canting resulting in an ordered magnetic component along the @xmath43 direction of @xmath56 ( canting angle @xmath57  deg ) .
the canting of the magnetic moment arises from a strong dzyaloshinski - moriya interaction that in turn stems from the strong soc .
minimizing the static energy of a pair sensing only the heisenberg interaction @xmath58 and the dzyaloshinski - moriya interaction , @xmath59 , yields the condition tan@xmath60 , and thus a rough estimate of the dzyaloshinski - moriya interaction in cro : @xmath61 .
the magnetization curves further show a decrease of the nel temperature with increasing ti content , see tab .
[ tab : char ] in appendix a , which can be explained by the dilution of the magnetic lattice by non - magnetic ti .
( color online ) energy scans at various * q * values with polarization analysis . raw data of different ( a ) nsf , and ( b , c , e - g ) sf channels with the polarization parallel to * x * ( * y*),(*z * ) in blue ( red ) , ( black ) at different scattering vectors .
( d , h , i ) show the results of the polarization analysis with eq .
( 1 ) , in red m@xmath62 , and m@xmath63 in blue ( green ) at * * q**=(2,1,0 ) ( * * q**=(1,2,0 ) ) . in ( d )
the blue line depicts the modeled asymmetric shape resulting from the folding of the magnon dispersion with the resolution function , see ref .
@xcite , which is also used for the modeling in ( i ) , but scaled with the expected correction factors for geometry and twinning fractions .
the inset in ( i ) shows a zoom into the data . ]    with the alignment of the antiferromagnetic moment along @xmath36 one expects non - degenerate transverse magnon modes corresponding to polarization along @xmath43 or @xmath16 @xcite . with our previous unpolarized inelastic neutron scattering experiments @xcite only the @xmath43 polarized mode could be clearly identified .
the gap of the @xmath16 polarized transverse mode could either be hidden in the shoulder of the in - plane signal , which inevitably arises from the folding of the instrument resolution with the steep spin - wave dispersion , or appear at higher energies where weak signals were detected @xcite . with the new polarized experiment the in - plane and @xmath16 polarized transverse magnons can be easily separated in the polarization and it was the first aim to search for the @xmath16 polarized mode .
figure  [ modes](a - c ) shows energy scans through the anisotropy gap at the antiferromagnetic zone center .
the nsf channels shown in fig .  [ modes](a ) have a larger background than the sf channels [ fig .
[ modes ]  ( b - c ) ] and are thus less informative but confirm the main conclusions .
the @xmath16 polarized mode must entirely contribute to the sfy and nsfz channels with a complete geometry factor .
spin - wave calculations show that the two non - degenerate transverse modes exhibit an intrinsic signal strength inversely proportional to their zone - center energy @xcite .
therefore , the @xmath16 polarized mode can be excluded in the asymmetric peak of the in - plane transverse mode .
these scans confirm the in - plane character of the signal in the range 14 to 20  mev , see fig .
[ modes](d ) . a special effort
was laid on the analysis of the signal maximum at the in - plane magnon gap , e=14  mev , see fig .
[ modes](e , f ) . with the partially detwinned crystal one expects this signal from the transverse magnon to be reduced by a factor 0.6 for @xmath28 , which agrees with the measured ratio of the intensities at 14  mev between @xmath28 and @xmath31 of 0.55(10 ) .
we also measured the backfolded mode at the ferromagnetic zone - center , which according to the spin - wave calculations , see below , exhibits a @xmath16 polarization .
indeed such a @xmath16 polarized signal can be determined at @xmath64 and e=14  mev [ fig .
[ modes](g ) ] by counting for a very long time ( 45  min for each data point ) .    in view of the recent observation of longitudinal magnetic excitations in cro
@xcite it seemed worthwhile to search for such modes first at the antiferromagnetic zone centers ( 2,1,0 ) and ( 1,2,0 ) . in fig .
[ modes](i ) the fitted background from the in - plane polarized transverse mode from fig .  [ modes](d ) is used to separate a weak signal peaking at 29  mev .
the signals for the two different * q * values are assumed to posses the same shape , but are scaled with the expected factors for geometry and twinning . around 30  mev the stronger signal for @xmath28 compared to @xmath31 thus points to a longitudinal mode in the sense that it is polarized parallel to the static moment .
the energy of this mode corresponds fairly well to the double of the gap of the in - plane transverse branch and the strength is below 10  % of that signal .
therefore it seems most likely that this longitudinal signal completely stems from the two - magnon excitation , which is expected to appear in the longitudinal polarization channel @xcite .
a similar discussion has also been initiated for the parent material of feas - based superconductors where again the two - magnon explanation seems most likely @xcite .
the 30  mev longitudinal signal is thus not anomalous .
the @xmath16 polarized antiferromagnetic zone - center mode is detected at higher energy , where the unpolarized experiment found some evidence for additional scattering @xcite , but where the signal strength is expected to be rather small @xcite . in order to cover the energies of the order of 40
mev we needed to use a larger value of @xmath65  @xmath35 , which also allows one to avoid the contamination appearing at e=44  mev for the standard value of @xmath66  @xmath35 .
the magnetic signals polarized along @xmath67 and along * y * at large energies are shown in fig .
[ modes](h ) and show that the zone - center @xmath16-polarized magnon possesses a large energy of 45.5(1.5 )  mev .
its signal strength is in rough agreement with the signal strength of the in - plane transverse mode with the energy being enhanced by a factor three .
our conclusion of @xmath16 polarized modes at the ferromagnetic and antiferromagnetic zone centers agrees with the interpretation of similar polarized neutron scattering experiments performed with twinned crystals in a different scattering geometry @xcite .
this large @xmath16 anisotropy is remarkable , as it strongly deforms the magnon dispersion in cro with respect to a simple isotropic model .
the splitting of the two antiferromagnetic zone - center magnons , 14 and 45.5  mev , is larger than the dispersion of the in - plane branch to the zone - boundary . as will be discussed below
, the entire dispersion of transverse branches is , nevertheless , well described with the @xmath8 spin - wave model using the holstein - primakoff transformation .
the in - plane and @xmath16-polarized branches exchange their character : while the @xmath16-polarized mode is the high - energy mode at the antiferromagnetic zone center , it appears at the lower energy of the in - plane transverse magnon at the ferromagnetic zone center , in accordance to the data shown in fig .
[ modes](g ) .
the in - plane polarized branch thus starts at 14  mev at the antiferromagnetic zone center , exceeds to the zone boundary at 41  mev and then continues to stiffen till 45.5  mev at the ferromagnetic zone center .
the out - of - plane polarized branch just exhibits the opposite dispersion . by comparing the intensities at the ferromagnetic zone centers @xmath68 and @xmath69 with an energy transfer of 45  mev
we may confirm the in - plane transverse character of this high - energy mode appearing at the ferromagnetic zone center .
the expected ratio of the signals taking into account the different twinning fractions and geometry factors amounts to 1:2.5 for a transverse and to 2.5:1 for a longitudinal magnon .
we obtain a ratio of 0.1(2 ) indicating the transverse in - plane character .
evidence for a longitudinal high energy branch could not be obtained in our experiment @xcite .
magnon dispersion calculated with the spinw program using the values given in the text .
the blue circles denote the fit values of the dispersion obtained in this and the previous @xcite neutron scattering studies .
the black lines are the calculated dispersion and the color code denote the calculated neutron scattering intensity of the convoluted spectra . in each panel
black stands for maximum intensity and white for none . in ( b ,
c ) only the in - plane- and out - of - plane polarized transverse modes are shown , respectively , illustrating the opposite dispersion of these branches . ]    the anisotropy gap of the @xmath16-polarized transverse magnon of 45.5  mev leads to the uncommon feature in the magnon dispersion , that the transverse in - plane branch continues to increase in energy between the zone boundary and a ferromagnetic zone center . this peculiarity can , however , be well described with a rather conventional model .
the hamiltonian , which is used for the description of the magnon dispersion is given by :    @xmath70    the sum runs over pairs of magnetic ions , so that each pair or bond appears twice and @xmath10 is set to 0.67 following the experimental results in ref .
the spin - wave calculations using the spinw program and the parameters @xmath71  mev , @xmath72  mev , @xmath73  mev , @xmath74  mev and @xmath75  mev , give a good description of the magnon dispersion obtained in this and the previous @xcite neutron scattering studies [ fig .
[ spinw](a ) ] .
panels fig .
[ spinw](b , c ) show the in - plane and c - polarized branches , respectively , starting at an antiferromagnetic zone center and proceeding to a ferromagnetic one . in this figure
the color code denotes maximum neutron scattering intensity of the convoluted spectra with black color , and zero intensity with white .
the energy resolution is set to 1  mev .
the spin - wave calculation thus perfectly describes not only the energy dispersion but also the polarization of the magnon branches .
note that this rather uncommon dispersion is described with strong single - ion parameters arising from soc .
these strong anisotropy terms interfere with the heisenberg interaction parameters in contrast to models invoking only weak anisotropy .
therefore , the @xmath58 parameter differs from that obtained by fitting only the lower part of the dispersion with a small uniaxial anisotropy @xcite .
figure  [ spinw](a ) shows that the non - sinusoidal parts of the dispersion are perfectly described , clearly better than without the strong anisotropy terms .
it is , however , important to note that the model is not unique .
it is possible to obtain similar fitting by partially reducing the single - ion anisotropy and by inducing an anisotropic nearest - neighbor interaction .
the coupling parameter j@xmath76 acting between neighboring layers splits the magnon modes into two by introducing a finite dispersion perpendicular to the planes .
whether the lower or upper mode is seen at @xmath77 is determined by the sign of @xmath78 . with the partially detwinned crystals ( in particular with the experiment on 1tib )
it is possible to determine the sign of this interaction .
here it is chosen to couple spins at ( 0,0,0 ) and ( 0,0.5,0.5 ) .
the dominant @xmath36 components of these spins are parallel for the @xmath14-centering and antiparallel for the @xmath15-centering .
if @xmath78 is positive ( negative ) , corresponding to an antiferromagnetic ( ferromagnetic ) coupling , the upper ( lower ) mode is seen at @xmath77 .
the experimental data reveal that the upper mode is seen at @xmath77 , so the coupling is ferromagnetic , stabilizing the @xmath14-centered phase .
this coupling contradicts the observation of a @xmath15-centered structure as the main magnetic scheme , see subsection iii.b .
cro thus exhibits the uncommon situation that the minimum magnon energy does not occur at the magnetic bragg peaks .
this observation can be attributed to an anisotropic @xmath78 which differs for the @xmath43 and @xmath36 spin components .
the linearized augmented plane wave ( lapw ) method as implemented in the wien2k package@xcite was used for the dft calculations .
the exchange correlation potential for the generalized gradient approximation ( gga ) was chosen to be in the form proposed in ref .
the soc was treated in a second variational way .
the parameter of the plane - wave expansion was chosen to be @xmath79 , where @xmath80 is the smallest atomic sphere radii ( @xmath81  a.u .
, @xmath82  a.u .
, @xmath83  a.u . ) and @xmath84 the plane - wave cutoff .
we used a mesh consisting of 800  @xmath85points .
the on - site hubbard repulsion and intra - atomic hund s exchange were taken to be @xmath86  ev and @xmath87  ev@xcite in the gga+u@xcite and gga+u+soc calculations .
the exchange parameters were recalculated via total energies of three different collinear magnetic configurations ( a ferromagnetic and two , in the @xmath16 direction differently stacked , antiferromagnetic configurations ) using the gga+u approach .
we found that @xmath88  mev , which agrees with the results of the spin - wave calculations .    in order to estimate the single - ion anisotropy we add soc to the calculation scheme and computed the energies of two configurations , where all spins are either directed along @xmath16 or lie in the @xmath25  plane .
the antiferromagnetic order was assumed in these calculations .
the lowest total energy corresponds to the configuration , where all spins are directed predominantly along the @xmath36  axis .
the single - ion anisotropy is found to be very large , @xmath89  mev , again in good agreement with experiment .
very recently microscopic magnetic parameters were calculated for similar @xmath90 and @xmath91 values finding a nearest - neighbor interaction of 3 - 6mev and a large single ion anisotropy in agreement with our results @xcite .
the spin moment on the ru ions was calculated to be @xmath92 ( i.e. @xmath93 ) , while the orbital moment is @xmath94 .
we also find a sizable magnetic moment on the apical o @xmath95 , while planar ligands stay nonmagnetic .
the total moment , @xmath96 , agrees well with the powder neutron diffraction result , and the small value of the orbital moment indicates that cro is not close to a @xmath13 state driven by soc as it is also deduced from the calculations in reference @xcite .
( color online ) influence of ti substitution on the anisotropy gap .
( a ) shows energy scans through the anisotropy gap at the antiferromagnetic zone center using the crystal 10ti , l=2 ( blue ) and l=3 ( black ) and 0ti , l=2 ( red ) . in ( b )
the energy of the lower zone - center mode of crystals with different ti content is shown , the data for 1  % ti are from ref .  @xcite . ]    in fig .
[ ti ] the influence of different ti substitutions on the anisotropy gap is addressed by comparing constant * q * scans across the in - plane gap at the antiferromagnetic zone center .
the modes are slightly split due to finite inter - layer interaction , see subsection iii.d .
even @xmath97 corresponds to the higher and odd to the lower modes , respectively .
the lower zone - center magnon energy is displayed as function of ti doping in fig .
[ ti](b ) .
the gap clearly diminishes with increasing ti content .
since the non - magnetic ti dilutes the magnetic lattice a general softening can be expected , as it is visible in the anisotropy gap , which in first view corresponds to the square root of exchange and anisotropy energies .
in addition , ti also perturbs the lifting of orbital degeneracy , as it is shown in fig .
[ ruooo ] and subsection iii .
a. in consequence also the single - ion anisotropy will be reduced with increasing ti content .
the impact of a small ti content of the order of one percent , which is used to favor the synthesis of large crystals , can however , safely be neglected .
10ti , which is essentially untwinned , shows a @xmath97 dependence of the anisotropy gap [ fig :  [ ti](a ) ] consistent with the previous study @xcite and the discussion given above .
( color online ) additional magnon branch near 5  mev .
( a , b ) show scattering maps of energy versus * q*. ( c ) presents some characteristic scans included in ( b ) for an energy transfer of 3  mev ( 5mev ) , ( 8  mev ) in blue ( red ) , ( green ) .
the data in ( a - c ) are taken from ref .
( d ) shows * q * scans with an energy transfer of 5  mev for 10ti ( blue ) and 0ti ( red ) .
the intensities are scaled with the maximum intensities of the scans in fig .
[ ti ] . scaled in the same way
are the energy scans in ( e ) at @xmath98 using the crystals 10ti , @xmath99 ( blue ) , @xmath100 ( red ) and 0ti , @xmath100 ( black ) .
( f ) presents the energy of this additional mode at * * q**=(0,1,0 ) for crystals containing different amounts of ti . ]
( color online ) polarization analysis of the additional mode .
sf intensities of energy scans at ( a , d ) @xmath101 , and ( b , e ) @xmath77 for the sfx ( sfy ) , ( sfz ) chanel in blue ( red ) , ( black ) .
( c ) shows the results of the polarization analysis with eq .
( 1 ) , in red m@xmath62 and in blue m@xmath63 .
( d , e ) present only the 5  mev data of the two * q * positions , which are measured with higher statistics . ]
the previous unpolarized experiments gave evidence for an additional scattering at the lower energy of 5mev that could not be explained by the magnon dispersion expected for a simple square - lattice antiferromagnet .
the data , which are presented in the mappings shown in fig .
[ omode2](a , b ) , were obtained on the thermal tas in8 @xcite .
panel ( c ) shows some characteristic scans included in ( b ) .
the 5  mev signal is much weaker than the in - plane transverse mode in particular when considering the by a factor three reduced energy .
furthermore , it is not possible to connect this low - energy feature with a dispersing branch with a slope comparable to that of the transverse branch .
instead the intensity of this feature seems to be rather localized in * q * space and to exhibit a flat dispersion .
the signal could also be followed along the vertical direction without any indication for finite dispersion .    in ref .
@xcite it was speculated that this extra mode could arise from the orbital disorder induced by the ti substitution , a scenario which can now be ruled out by the comparison of ti - containing and ti - free cro .
[ omode2](d - e ) show * q * and energy - scans across the 5  mev signal on samples with 0 and 10% of ti .
most interestingly there seems to be no difference in the strength of the low - energy signal when normalized to that of the in - plane gap mode at @xmath1914  mev .
the additional signal is thus intrinsic to the magnetic order in cro , and it is not induced by disorder , which must be significantly enhanced in ti10 .
the ti content furthermore does not change the energy of the additional scattering , see the fitted energy maxima displayed in fig .
[ omode2](f ) . the minor variation between measurements on the same concentration on different spectrometers can be attributed to the energy calibration .
the polarization of the additional mode was determined with polarized neutron scattering following the same procedures as those used to determine the character of the dispersion at higher energies .
the polarization analysis directly shows that the additional excitation possesses an in - plane character , because it is not observed in the sfy channel ( or in @xmath102 ) , see fig .
[ omodepol ] . comparing the intensities at @xmath101 and
@xmath77 one may furthermore deduce that the additional mode is @xmath43  polarized , which is transverse with respect to the static moment .
the expected ratio of the intensities at 5  mev between @xmath77 and @xmath101 based on the twinning fractions and the geometry factor is 1:2.5 for an in - plane transverse mode and 2.5:1 for an in - plane longitudinal mode , respectively .
the measured ratio is 0.25(21 ) characterizing this mode as an in - plane transverse mode .
the simple antiferromagnetic square lattice exhibits a magnon dispersion consisting of two transverse branches , which were both observed at larger energy in cro , see above . in order to explain the additional mode the model needs to be extended , and the weakness and localization of the signal point to a smaller moment with smaller coupling .
the gga+u calculations on cro reveal a small magnetic moment situated at the two apical oxygens of 0.11@xmath103 each .
such a polarization of the oxygen results from the strong hybridization between ru @xmath104 and o @xmath105 orbitals .
dft calculations on several ferromagnetic or nearly ferromagnetic ruthenates find a sizable amount of magnetic moment on o , which in total sum up to about 30  % of the entire magnetization @xcite .
experimentally this sizable magnetization of the o states has been observed by polarized neutron diffraction on ca@xmath106sr@xmath107ruo@xmath1 where indeed 33  % of the magnetization stems from o states @xcite .
cro , however , exhibits antiferromagnetic order , in which the o bridging two antiparallel ru moments can not be polarized , see fig .
[ explanations](a ) .
therefore , there is no moment on the in - plane o sites , which is also confirmed by the dft calculations of cro .
however , the apical o is connected only to a single ru site and can be polarized .
there are two principal contributions to the exchange between ru and o moments due to hopping between o 2@xmath108 and ru 4@xmath9 orbitals in the scenario where the ru @xmath6 orbital is doubly occupied and the @xmath44 and @xmath45 singly occupied , see e.g. ref .
one is a hopping from the @xmath109 and @xmath110 to the @xmath44 and @xmath45 orbitals , respectively .
there , because of pauli exclusion principle , only antiparallel spins can hop , leaving parallel spins in the now singly occupied o @xmath108-orbitals . in total
this yields a ferromagnetic coupling : @xmath111 where @xmath112 is the hopping integral between @xmath44 and @xmath109 , and @xmath45 and @xmath110 orbitals of ru and oxygen . @xmath113 and @xmath114 are the centers of the ru @xmath115 and apical o @xmath108 bands .
the factor two in eq .
( [ t2u ] ) reflects the fact that there are two processes ( virtual hoppings from @xmath109 to @xmath116 and from @xmath110 to @xmath45 orbitals ) , contributing to this ferromagnetic ru - o exchange .    the second contribution to
the ru - o exchange comes from the hopping from the @xmath117 orbital of o to the empty @xmath118 orbital of ru . because of hund s rule , predominantly parallel spins hop , leaving anitparallel spins on the o @xmath117 orbital .
this yields an antiferromagnetic coupling between ru and o moments : @xmath119 here @xmath120 is the hopping matrix element between a @xmath117 orbital of o and an empty @xmath121 orbital of ru , cf .
@xcite .
the @xmath122 hopping is larger than the @xmath123 one @xmath124 ) @xcite , but the antiferromagnetic exchange ( eq . [ t2jd2 ] ) occurs through the @xmath125 band , which lies much higher in energy than @xmath115 .
our gga+u calculations show that @xmath126 ev , while @xmath1274.7  ev .
there is only one such antiferromagnetic exchange path , in contrast to two for ferromagnetic exchanges
. [ t2u ] ) and hund s exchange on the ru , @xmath128 ev , also slightly reduces this exchange ( eq .
[ t2jd2 ] ) .
therefore , one should expect ferromagnetic coupling to dominate , so that the moment of oxygen would be parallel to the moment of neighboring ru .
our dft calculations support this conclusion .
( color online ) magnon dispersion calculated with the spinw program including the additional mode .
the labeling corresponds to the one in fig .
[ spinw ] . ]
we have extended the spin - wave model by a small moment on the apical oxygens , 0.11  @xmath103 , coupled ferromagnetically to the next ru moment , see fig .
[ explanations](a ) .
this extension necessitates to reduce the single - ion anisotropy parameter @xmath129 in order to keep the anisotropy gap of the in - plane transverse mode at the correct energy .
the additional mode can not be obtained at an energy near 5  mev at the antiferromagnetic zone center without introducing an anisotropy in terms of a single ion anisotropy for the magnetic moment on the o sites , or an anisotropic coupling parameter .
because the single ion anisotropy stems from soc which is small in o , the first scenario seems to be less likely .
a better choice is an anisotropic coupling , because the ru with its larger soc is also involved in this process .
therefore we used @xmath130  mev to fit the experimental results . with this anisotropic coupling
, @xmath129 has to be reduced to 0.5  mev to keep the in - plane transverse mode at the experimental value .
with this extension we can perfectly describe the additional feature , see fig .
[ omodedisp ] .
the flat dispersion of the additional branch results from the lack of coupling between two o moments and perfectly describes the experimental finding .
also the localization of the signal strength at the antiferromagnetic zone center is well reproduced by this model .
we can therefore conclude that the additional moment , situated on the apical o site and ferromagnetically coupled to the ru , can explain the additional low - energy magnetic signal in cro .
we also tried to reproduce the data with the o moments being antiferromagnetically coupled to ru , but such a scenario can not explain the location of the scattering in q - space .
the significant polarization of the o site is another consequence of the strong hybridization of the ru @xmath104 and o @xmath108 orbitals resulting in large hopping parameters , which can already be deduced from the closeness of the metallic phase of cro .
it furthermore yields a convincing explanation of the reduced moment , because the in - plane oxygen ion , which should carry a similar amount of charge carriers , can not be magnetically polarized in the antiferromagnetic structure reducing the ordered moment from the simple @xmath8 expectation .
the ordered moment may thus also oscillate in length by transferring moment between ru and in - plane o , which can give rise to a longitudinal mode , as an alternative explanation for the observation in reference @xcite .
various neutron scattering experiments on partially untwinned crystals of ca@xmath0ru@xmath3ti@xmath22o@xmath1 were performed to determine the magnetic structure and magnon dispersion in the insulating state of cro .
most remarkable is the large magnetic anisotropy , which results in a splitting of the two transverse zone - center modes that exceeds the full dispersion . the entire dispersion of transverse magnon branches , however , is perfectly described by a conventional spin - wave model , in which the strong impact of the soc is reflected by large anisotropy parameters . dft calculations within the gga+u approximation yield an orbital moment of only 0.13@xmath103 , which indicates that cro is not close to a @xmath13 state , which can be seen in the discrepancy of the experimental magnon dispersion with the @xmath13 calculation @xcite .
the calculated magnetic interaction and anisotropy parameters reasonably well agree with those obtained by fitting the dispersion , and also the ordered moment is correctly reproduced in the dft calculation .
the dft study finds a sizable amount of magnetic moment on the apical oxygen reminiscent of previous reports on ferromagnetic ruthenates .
this additional moment explains an additional signal appearing in the neutron scattering experiments at lower energy and limited to the antiferromagnetic zone center .
extending the spin - wave model to the weakly ferromagnetic coupled oxygen moments describes the flat dispersion and the limited appearance in * q * space of this signal .    the spin - wave dispersion in cro is thus dominated by the impact of strong soc and by the presence of magnetic moments on the oxygen sites .
this work was supported by the deutsche forschungsgemeinschaft through crc 1238 projects no .
a02 and no .
b04 and by the russian foundation of the basic research via program 16 - 02 - 00451 , russian president council on science through md-916.2017.2 , faso ( theme `` electron '' 01201463326 ) and mon ( project 236 ) .
complete crystal structure analyses were performed with an x8-apex by bruker axs single - crystal diffractometer with a goniometer in kappa - geometry and x - ray radiation from a molybdenum anode with a wavelength of @xmath131   .
the distance between the sample and the detector was set to 50  mm .
structure refinements were carried out using jana2006 @xcite .
a type i extinction correction was applied during the refinements and the data were corrected for absorption . the thermal parameters for ru1 and ti1 were constrained to be equal , and the total occupancy of this site was fixed to 1 .
the results of the structural refinements are given in table [ tab : ref ] .
table [ tab : char ] presents further characteristics of the crystal structure as well as the metal - insulator- and antiferromagnetic transition temperatures .
@xmath54 is determined from the magnetization curves presented in fig .
[ magnetization](a ) by finding the zero point of the second derivation .
@xmath49 is determined in the same way from the resistance curves presented in this figure .
[ cols= " < , < , < , < , < , < , < " , ]     [ tab : char ]
( color online ) elastic neutron scattering scans of 1tib across ( a ) @xmath69 , ( b ) @xmath77 , ( c ) @xmath134 , ( d ) @xmath135 , ( e ) @xmath135 .
blue color denotes scans taken in the magnetic phase at 10  k and red color scans taken in the non magnetic phase at 150  k. ]    the crystal 1tib used in the previous neutron scattering study @xcite was further investigated on the thermal tas in3 at the ill .
the crystal was mounted in the [ 010]/[001 ] orientation into a ill orange cryostat , @xmath34 was set to 2.662  @xmath35 for all scans and a pyrolitic graphite filter was used to suppress higher order wave length . rocking scans revealed the mosaic spread to be less than 0.5@xmath21 . with the good resolution of the in3
it is easily possible to resolve the bragg peaks from the two twins present in the crystal .
the twinning ratio amounts to 20:1 , which is deduced from rocking scans of strong nuclear bragg reflections at the 2@xmath136 scattering angles of the ( 020 ) and ( 200 ) reflections .
figure  [ magn - str](a ) shows a longitudinal scan over the ( 020 ) reflection , the smaller twin is barely seen in the tail .
scans across the positions , where magnetic scattering can be expected , are shown in ( b - e ) . there
several scans of the kind ( 01l ) are compared at 10  k , in the antiferromagnetic phase , and at 150  k , in the non magnetic phase .
it is clearly seen , that the minority twin gains intensity for an odd @xmath97 and not for an even @xmath97 .
the majority twin gains intensity for even l and much less , but also clearly detectable , for odd l. the ( 012 ) and ( 014 ) reflections are magnetic bragg peaks for the @xmath15-centered phase but not for the @xmath14-centered one .
so the strong gain of intensity of these reflections in the antiferromagnetic phase reveals an almost exclusive @xmath15-centered phase .
in opposite the ( 011 ) and ( 013 ) peaks are bragg peaks in the @xmath14-centered phase , but not in the @xmath15-centered one .
so the gain of intensity at low temperatures of those peaks points to the @xmath14-centered phase , which is however much weaker than the gain of the bragg peak intensities of the @xmath15-centered phase .
so the main stacking scheme of this crystal is the @xmath15-centered one .
53ifxundefined [ 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.1143/jpsj.66.1868 [ * * ,   ( ) ] link:\doibase 10.1103/physrevlett.84.2666 [ * * ,   ( ) ] link:\doibase    10.1103/physrevb.63.174432 [ * * ,   ( ) ] @noop * * , ( ) link:\doibase 10.1103/physrevlett.93.146401 [ * * ,   ( ) ] @noop * * , ( ) @noop * * ,   ( ) @noop * * ,   ( ) @noop * * ,   ( ) link:\doibase    10.1103/physrevb.60.r8422 [ * * ,   ( ) ] link:\doibase 10.1103/physrevb.58.847 [ * * ,   ( ) ] link:\doibase    10.1103/physrevlett.87.077202 [ * * ,   ( ) ] link:\doibase    10.1103/physrevlett.104.226401 [ * * ,   ( ) ] @noop ( ) ,   link:\doibase 10.1103/physrevlett.91.056403 [ * * ,   ( ) ] link:\doibase 10.1103/physrevlett.88.017201 [ * * ,   ( ) ] link:\doibase    10.1103/physrevlett.81.2747 [ * * ,   ( ) ] link:\doibase 10.1103/physrevlett.98.216403 [ * * ,   ( ) ] http://stacks.iop.org/0295-5075/63/i=1/a=097 [ * * ,   ( ) ] link:\doibase 10.1140/epjb / e20020021 [ * * ,   ( ) ] @noop _ _  ( ,  ,  ) link:\doibase 10.1103/physrevlett.111.197201 [ * * ,   ( ) ] link:\doibase 10.1103/physrevb.90.035137 [ * * ,   ( ) ] link:\doibase    10.1103/physrevlett.115.247201 [ * * ,   ( ) ] @noop ( ) ,   \doibase http://dx.doi.org/10.1016/s0921-4526(99)00989-8 [ * * ,   ( ) ] link:\doibase 10.1002/crat.201600020 [ * * , ( ) ] \doibase    http://dx.doi.org/10.1016/0921-4526(93)90108-i [ * * ,   ( ) ] http://stacks.iop.org/0953-8984/27/i=16/a=166002 [ * * ,   ( ) ] ,  ed . ,
@noop _ _  ( ,  ) link:\doibase    10.1103/physrevb.86.060410 [ * * ,   ( ) ] @noop * * ,   ( ) link:\doibase    10.1103/physrevb.57.1236 [ * * ,   ( ) ] link:\doibase    10.1103/physrevb.76.014505 [ * * ,   ( ) ] \doibase    http://dx.doi.org/10.1016/s0921-4534(96)00637-5 [ * * ,   ( ) ] link:\doibase 10.1103/physrevb.84.201102 [ * * ,   ( ) ] link:\doibase    10.1103/physrevlett.88.197002 [ * * ,   ( ) ] link:\doibase 10.1103/physrevb.79.054422 [ * * ,   ( ) ] link:\doibase 10.1103/physrevb.72.094104 [ * * ,   ( ) ] link:\doibase    10.1103/physrevb.24.3939 [ * * ,   ( ) ] link:\doibase 10.1103/physrevx.3.041036 [ * * ,   ( ) ] link:\doibase 10.1140/epjb / e2016 - 60588 - 6 [ * * ,   ( ) ] @noop ( ) ,   @noop @noop _ _  ( ,  ,  ) http://www.ncbi.nlm.nih.gov/pubmed/10062328 [ * * ,   ( ) ] link:\doibase 10.1038/nmat1605 [ * * ,   ( ) ] link:\doibase 10.1103/physrevb.86.064429 [ * * ,   ( ) ] http://iopscience.iop.org/0953-8984/9/4/002 [ * * , ( ) ] g. zhang and e. pavarini , phys . rev .
b * 95 * , 075145 ( 2017 ) .
link:\doibase 10.1103/physrevb.56.2556 [ * * , ( ) ] link:\doibase    10.1103/physrevlett.89.087202 [ * * ,   ( ) ] @noop _ _  ( ,  ,  )  p.   @noop * * ,   ( )