metal - organic frameworks ( mofs ) are organic - inorganic hybrid solids in which a usually low - dimensional inorganic subnetwork is connected with adequately functionalized organic molecules to form a 3d framework .
the fascinating structures and properties of these materials have been the focus of intense research efforts over the last 15 years or so.@xcite at present , a wide variety of new compounds continues to be discovered , because of the huge combinatorial possibilities offered by the amalgamation of organic and inorganic chemistry .
this new class of compounds has even proved able to produce interesting properties for industrial applications.@xcite one of the main reasons for this success is their porous character.@xcite various applications have been envisaged for mofs , principally in the field of gas storage and separation,@xcite but also for catalysis,@xcite luminescence,@xcite or magnetorefrigeration.@xcite    mofs containing magnetic ions furthermore represent a goldmine for low - dimensional and frustrated magnetism.@xcite first , when the organic moiety of these materials is composed of magnetically inert molecules such as sufficiently long alkane chains , it can act as a simple spacer that prevents the establishment of electronically - driven magnetic interactions between e.g. spin chains or spin layers .
in other words , such systems have a strongly low - dimensional character regarding through - bond magnetic interactions and , consequently , they only become magnetically ordered at low or very - low temperatures , when the one - dimensional ( 1d ) or two - dimensional ( 2d ) correlation length diverges and promotes sizable classical dipolar interactions between correlated spin blocks .
second , in addition to the intrinsic low - dimensionality , magnetic mofs also show a large diversity in the type of topologies formed by the inter - connection of the magnetic ions ( m ) through the possible superexchange pathways ( m - o - m pathways in most cases ) . despite of the clear interest of these materials in the field of magnetism , hybrid frameworks are relatively little studied by neutron scattering techniques,@xcite although these methods are commonly used to investigate the magnetic structures and excitations of purely inorganic compounds .
this is partly due to the difficulty in synthesizing deuterated compounds when complex organic ligands are involved , and partly due to the weak magnetic contributions to the total neutron scattering from these materials .
recently , we reported a new cobalt^ii^ hydroxy - sebacate framework of formula co^ii^@xmath0(oh)@xmath1(c@xmath2@xmath3h@xmath2@xmath4o@xmath0)@xmath5.@xcite the compound contains co^ii^-based magnetic chains which are very well separated from each other ( 1 to 2 nm ) .
we revealed the coexistence of magnetic long - range order ( lro ) and slow dynamics of the magnetization below @xmath7 k. this new compound is an excellent candidate for a detailed magnetic study for at least two reasons .
first , the magnetic chains are much better isolated from each other than in most of the spin - chain compounds .
second , the chains themselves present a rather complex topology of magnetic interactions between spin carriers which are well known for their large single - ion anisotropy .
spin - chains attract much attention because they provide genuine examples for testing models that are less tractable in higher dimensions.@xcite for example , spin chains with nearest - neighbor antiferromagnetic interactions display a large variety of ground states with non - classical magnetic phenomena.@xcite although ideal 1d spin systems do not undergo lro at finite temperature due to strong quantum fluctuations,@xcite slight perturbations such as , for instance , weak interchain interactions can make the quantum critical state unstable to three - dimensional ( 3d ) lro.@xcite another striking feature of 1d magnetic systems is the so - called single - chain magnet ( scm ) behavior observed in a relatively large number of metallo - organic compounds.@xcite these systems , though not long - range ordered , present bulk magnetic behavior at low temperature that results from slow intrachain dynamics arising from the combination of a large magnetic anisotropy with ferromagnetic interactions .
this was predicted in the pioneering work of r. j. glauber in 1963.@xcite    in this article we provide deeper insight on the magnetism in the ordered phase of co^ii^@xmath0(oh)@xmath1(c@xmath2@xmath3h@xmath2@xmath4o@xmath0)@xmath5 . in section
ii we recall the necessary information concerning the crystal structure and macroscopic magnetic behavior of the system at @xmath8 k. experimental details are given in section iii .
our new results are presented and discussed in section iv .
firstly , the zero - field magnetic structure is determined from neutron powder diffraction data .
the magnetic field - dependence of this spin arrangement is qualitatively evaluated .
secondly , we present new magnetic measurements for temperatures down to 70 mk of both magnetization and ac susceptibility , and deduce the characteristic energies of the observed dynamics . finally , in section v , we discuss the magnetic structure in relation with the magnetic couplings and anisotropy in the compound .
we further propose a scenario to describe the dynamics .
finally , in section vi we summarize the present work .
( oh)@xmath1(c@xmath2@xmath3h@xmath2@xmath4o@xmath0)@xmath5.@xcite co^ii^1o@xmath4 and co^ii^2o@xmath9 polyhedra are represented in yellow and red , respectively .
( i ) projection of the 3d metal - organic framework perpendicular to the magnetic chain axis .
( ii ) atomic structure of the 1d inorganic subnetwork.,width=264 ]    .
( i ) j1 , j3 and j4 correspond to corner - sharing pathways having the angle values given on the figure , while j2 corresponds to an edge - sharing pathway with co^ii^1-o - co^ii^2 angles of about 91^o^ and 96^o^. ( ii ) spin topology viewed along the [ 101 ] crystallographic direction , i.e. by applying a 45^o^ rotation around the chain axis from the orientation of figure [ fig.2](i).,width=302 ]    co^ii^@xmath0(oh)@xmath1(c@xmath2@xmath3h@xmath2@xmath4o@xmath0)@xmath5 adopts a _
_ p__2@xmath2__/c _ _ monoclinic structure in which magnetic chains develop along the * b*-axis .
isoreticular compounds , i.e. having a similar framework organization though eventually based on different organic ligands , have been reported with manganese^ii^ and iron^ii^ ions @xcite .
figure [ fig.1 ] sheds light on the main structural characteristic of this compound : the strongly 1d character due to large distances between the inorganic chains .
the distance between the centers of neighboring chains varies from 11 to 21  , with the shortest path ( i.e. in the [ 001 ] direction , see figure [ fig.1](i ) ) not involving chemical connections . from the point of view of the magnetism ,
the combination of these large interchain distances with either the absence of interchain chemical connection or its presence through magnetically inert molecules ( alkane chain with 9 consecutive c - c single bonds ) precludes any role of interchain interactions other than classical dipolar forces .
the 1d inorganic subnetwork of this material contains two crystallographically distinct magnetic ions , co^ii^1 and co^ii^2 , that have distorted octahedral and distorted trigonal bipyramidal oxygen environments , respectively .
both sites have intrinsically asymmetric crystal field environments as the oxygen atoms originate from carboxylate functions and hydroxide groups .
the chains can be initially viewed as being constituted of a central ` zig - zag ' chain of corner - sharing co^ii^1o@xmath4 octahedra .
in addition , co^ii^2o@xmath9 trigonal bipyramids share oxygen atoms with three consecutive co^ii^1o@xmath4 octahedra of the main chain .
both magnetic sites are in the general 4__e _ _ position which , within a single inorganic chain , splits into two equivalent positions related by the 2-fold screw axis and labeled as ` a ' and ` b ' on figure [ fig.1](ii ) .
four intrachain superexchange pathways can be distinguished and are emphasized on figure [ fig.2](i ) .
first , yellow co^ii^1o@xmath4 octahedra are connected to each other along the chain axis by sharing corners through a large superexchange angle of about 136^o^ ( j1 ) .
second , each co^ii^2o@xmath9 bipyramid ( red ) connects with two co^ii^1o@xmath4 octahedra by corner - sharing through smaller superexchange angles ( @xmath10102^o^ for j3 and @xmath10107^o^ for j4 ) ; and with one co^ii^1o@xmath4 octahedron by edge - sharing ( j2 ) with the nearly orthogonal angles ( @xmath1091^o^ and @xmath1097^o^ ) .
the resulting spin topology is represented on figure [ fig.2](ii ) .
in addition to these four m - o - m pathways , weaker intrachain magnetic ` super - superexchange ' pathways exist through the carboxylic functions of the sebacate ligand ( m - o - c - o - m pathways ) , leading to further terms expected to contribute to the full hamiltonian .
the first magnetic characterizations of co^ii^@xmath0(oh)@xmath1(c@xmath2@xmath3h@xmath2@xmath4o@xmath0)@xmath5 have been previously reported.@xcite the compound undergoes magnetic long - range ordering below t@xmath12=5.4 k , as evidenced by heat capacity ( c@xmath13 ) measurements and by the appearance of magnetic scattering below t@xmath12 as observed by neutron diffraction .
the presence of a broad bump in the zero - field c@xmath13(t ) curve , between the @xmath14-peak at t@xmath12 and 2 to 3 times t@xmath12 , was attributed to the effects of magnetic correlations within the chains .
a complex behavior of the magnetic susceptibility was observed above the lro transition .
the non - monotonic shape of the product @xmath15 could be described by introducing two effective energy parameters ( 20 and -30 k ) which account for both anistropy and ( ferromagnetic and antiferromagnetic ) exchange interactions .
ac susceptibility measurements showed the existence of dynamical processes at t @xmath16 t@xmath12 .
ac susceptibility data were taken on a physical properties measurement system ( quantum design , usa ) for frequencies of the ac magnetic field between 10 and 6500 hz .
the analysis revealed the existence of a thermally activated behavior within the 3 - 4.5 k temperature range , with an energy barrier of about 67 k and a characteristic time of the order of 1.4@xmath17 s. the observed slow dynamics was attributed to the motion of thermally activated domain walls within the magnetic chains .
finally , ac susceptibility measurements performed under various bias dc magnetic fields revealed a strong displacement of the main peak towards higher temperatures , for both in - phase and out - of - phase components .
this latter observation was suggested to evidence the existence of field - induced magnetic order for temperatures above the zero - field transition temperature t@xmath12 .
the two powder samples of co^ii^@xmath0(oh)@xmath1(c@xmath2@xmath3h@xmath2@xmath4o@xmath0)@xmath5 used in this work were synthesized hydrothermally according to the method reported previously.@xcite chemicals ( sebacic acid and naoh ) were not deuterated , obtained from commercial suppliers , and used as received .
high - purity distilled water was used for the reaction medium .
long - exposure laboratory x - ray powder diffraction revealed the samples to be single phase .
neutron scattering experiments were conducted using the high - flux d20 two - axis diffractometer at the ill ( institut laue langevin , grenoble , france ) .
the neutron wavelength was @xmath18  .
all experiments were performed with a 6 mm diameter vanadium sample container and a cryomagnet providing temperatures down to 1.8 k and magnetic fields up to 6 t. zero - field data were first recorded by simply loading the 2 g powder into the sample holder .
two long duration neutron diffraction patterns were recorded at 1.8 k and 20 k ( 8 hours for each ) .    for measurements with an applied magnetic field ,
the sample was immersed in deuterated isopropanol in order to avoid reorientation of the powder particles .
isopropanol freezes around t @xmath10 115 k in a glassy state , giving rise to a weak and broad bump in the data around @xmath19 .
after initial rapid cooling of the sample in the remnant field of the magnet , the temperature was never raised above the glass transition temperature of isopropanol for the remainder of the experiment .
it was therefore guaranteed that the powder grains remained fixed and that the same powder average was measured in each run .
a diffraction pattern was recorded under these conditions , in zero - field at 10 k , and used for subtraction of the nuclear scattering .
data were then recorded under zero - field and under two different applied magnetic fields ( 0.5 and 2 t ) . for each of these three field values ,
diffraction data were collected for different temperatures across the transition , from 1.8 k to higher temperatures .
after these measurements , another zero - field 10 k pattern was recorded in order to confirm the absence of reorientation of the powder particles .
the neutron patterns were analyzed by rietveld refinements using the fullprof program .
@xcite the shape of the diffraction peaks was fitted using a pseudo - voigt function . besides the six peak profile parameters , the fit of the paramagnetic pattern recorded at 20 k comprised zero shift , cell parameters , preferred orientation of the powder particles , and scale factor .
the agreement indices r@xmath20 and r@xmath21 are the factors defined in ref . .
the difference resulting from the two long - duration patterns recorded at 1.8 and 20 k was used for the refinement of the magnetic structure .
all parameters were fixed to their previously refined values ( 20 k ) , except the lattice parameters which were taken from a le bail decomposition@xcite of the 1.8 k pattern .
the agreement factor noted r@xmath22 is the bragg agreement factor@xcite for the fit of the magnetic scattering observed on the difference pattern .
the representational analysis techniques described by bertaut@xcite were used to determine the possible magnetic structures compatible with both the magnetic wavevector and the symmetry of the crystal structure .
the irreducible representations of the propagation vector group were determined with the program basireps@xcite .
the possible solutions were first tested by refining the mixing coefficients of the basis vectors , as implemented in the fullprof program.@xcite the final rietveld refinement of the magnetic scattering was performed in the ` spherical coordinates procedure ' of the software .
the three parameters ( r , @xmath23 , @xmath24 ) per co^ii^ characterizing their atomic magnetic moments correspond to the value of the magnetic moment @xmath25 ( in bohr magnetons ) , the polar @xmath23 angle with z axis and the azimuth @xmath24 angle between the x axis and the projection of r onto the ( x , y ) plane . since this mode works only if z is perpendicular to the ( x , y ) plane , the present monoclinic space group was described in the laue class 1 1 2__/m _
_ for the magnetic refinement .
two powder samples ( sample 1 - mass @xmath26 mg and sample 2 - mass @xmath27 mg ) were measured using low temperature superconducting quantum interference device ( squid ) magnetometers equipped with a miniature dilution refrigerator developed at the institut nel - cnrs grenoble @xcite . sample 1 was obtained from the 2 g sample used for neutron scattering experiments .
the samples were prepared in a copper pouch and mixed with apiezon grease ( for thermalization and to avoid reorientation of the powder under magnetic field ) .
they were attached to a copper tress suspended from the mixing chamber of the dilution fridge .
@xmath28 hysteresis loops were measured for both samples in the ` high ' field magnetometer which reaches 8 t. ac and relaxation measurements were performed in the ` low ' field magnetometer , which reaches 3500 oe , for sample 1 .
the magnitude of the alternating field was @xmath29 oe .
the rietveld refinement of the 20 k neutron pattern was obtained by using the atomic coordinates determined from x - ray single - crystal diffraction at 100 k ( see figure [ rietveld_nuclear])@xcite . the good agreement between the calculated and recorded neutron patterns ( r@xmath20 = 0.115 and r@xmath21 = 0.127 ) suggests the absence of significant modifications of the crystal structure below 100 k. the refined unit cell parameters are : @xmath30  , @xmath31  , @xmath32  , and @xmath33 ^o^.    the magnetic scattering was extracted from the data recorded at 1.8 k by subtracting the 20 k pattern .
the resulting ` difference pattern ' shows weak magnetic bragg peaks both on top of some nuclear peaks of the ( _ h_,0,0 ) type and on the ( 0,0,_l _ ) and ( _ h_,0,_l _ ) peaks with odd _
l _ indices .
these two latter families of reflections , absent in the high - temperature pattern , are forbidden by the _ c_-glide of the nuclear space group .
these indications for such a non - bravais lattice ( breaking of the nuclear point - group symmetries ) imply a commensurate magnetic structure with a propagation vector * * _ k_**=(0,0,0 ) .    to simplify the analysis , group theory was used in order to predict the magnetic arrangements compatible with the nuclear symmetry , thus reducing the number of independent parameters in the refinement . in the landau theory , magnetic fluctuations in the paramagnetic state
while approaching a second - order phase transition reflect the crystal symmetries . at the transition ,
one magnetic mode is stabilized .
the magnetic modes compatible with the nuclear space group _ g _ (= _
_ p__2@xmath2__/c _ _ ) are given by the irreducible representations of the little group _
g@xmath34 _ that leave * _ k _ * invariant , i.e. the irreducible representations of _ g _ for * * _ k_**=(0,0,0 ) .
the application of the group theory , using the irreducible representations tabulated by kovalev@xcite and available in the basireps program,@xcite shows that the group representation can be reduced to a sum of four irreducible representations of order 1 : @xmath35 .        [ cols="<,^,^,^,^,^,^,^,^,^,^,^,^,^,^,^,^,^,^ " , ]     the four possible magnetic arrangements resulting from this analysis are summarized in table [ irreps ] . in principle ,
@xmath36 and @xmath37 can be excluded because they do not contain a ferromagnetic coupling while magnetic scattering is detected on the top of allowed nuclear peaks . moreover , testing the four models by refining the mixing coefficients of the basis vectors of the two magnetic moments indicates @xmath38 to most likely describe the magnetic structure .
the three other models provide unreasonable magnetic moment values and/or much larger agreement factors . imposing the irreducible representation @xmath38 , the refinement of the six mixing coefficients converges and gives reasonable magnetic moment values .
in particular , the magnetic moment of co^ii^1 is oriented along the chain axis ( * b*-axis ) , implying that the co^ii^1 sublattice forms a simple antiferromagnetic chain . in a second step ,
the refinement was realized by describing the magnetic moments by their spherical coordinates ( see section [ exp ] ) , in order to improve the stability of the refinement procedure .
this procedure confirmed the orientation of the co^ii^1 moment to lie along the * b*-axis .
thus , we constrained the magnetic moment of co^ii^1 with @xmath39^o^ and @xmath40 = 0@xmath41 and 180@xmath41 for the positions related by the screw axis , in accordance with @xmath38 .
the magnetic moment of co^ii^2 , for its part , was refined with its three spherical coordinates by imposing the couplings relevant to @xmath38 : ( @xmath42 , @xmath43 , @xmath44 ) for co@xmath45 ions in positions 1 and 3 , and ( @xmath42 , @xmath43 , @xmath46 ) for co@xmath45 in positions 2 and 4 . using these four parameters ( @xmath47 for co@xmath48 and @xmath42 , @xmath43 , @xmath44 for co@xmath45 )
, the refinement procedure converges satisfactorily ( r@xmath49 , @xmath50=1.49 , see figure [ rietveld_magnet ] ) .
starting from this model , lifting the additional constraints applied to the magnetic moment of co^ii^1 ( @xmath51 ) does not lead to a significantly different magnetic structure . in the final model
, the refined @xmath43 value of about @xmath52 corresponds to the monoclinic angle of the unit cell , indicating that the co@xmath45 magnetic moment is roughly in the ( * b*,*c * ) plane .
consequently , @xmath44 ( @xmath53 ) defines both the rotation of the co@xmath45 magnetic moments in the ( * b*,*c * ) plane relative to the chain axis , and their canting angle between the positions related by the screw axis ( @xmath54 ) .
model ( see text for details).,width=283 ]     correspond to a setting with x = * a * , y = * c * , z = * b * ( laue class 1 1 2__/m _ _ , see section [ exp]).,width=321 ]    this refinement gives magnetic moment values of 2.19(7 ) and 2.68(8 ) for co@xmath48 and co@xmath45 respectively ( see figure [ magn_struct ] and table [ table_magmom ] ) .
the co@xmath48 sublattice is fully compensated ( @xmath51 ) , while co@xmath45 has both compensated ( @xmath55 ) and uncompensated ( @xmath56 ) components .
the global ferromagnetic component for this material is thus of about @xmath57.co@xmath58mol@xmath59 .
this value , after averaging over the three directions ( @xmath60.co@xmath58mol@xmath59 ) , is in good agreement with the size of the magnetization jump at very low field ( @xmath61.co@xmath58mol@xmath59 , see section [ mh ] and ref . ) . the magnetic structure is presented on figure [ magn_struct ] , and is further discussed in section [ discussion ] .    *
7c & _ _ m__@xmath62 & _ _ m__@xmath63 & _ _ m__@xmath64 & * m & @xmath23 & @xmath24 + co@xmath581 - 1 & 0 & 2.19(7 ) & 0 & 2.19(7 ) & 0 & 0 + co@xmath581 - 2 & 0 & -2.19(7 ) & 0 & 2.19(7 ) & 180 & 0 + co@xmath581 - 3 & 0 & 2.19(7 ) & 0 & 2.19(7 ) & 0 & 0 + co@xmath581 - 4 & 0 & -2.19(7 ) & 0 & 2.19(7 ) & 180 & 0 + co@xmath582 - 1 & 0 & -2.13(5 ) & -1.64(8 ) & 2.68(8 ) & -36(2 ) & 101(3 ) + co@xmath582 - 2 & 0 & 2.13(5 ) & -1.64(8 ) & 2.68(8 ) & 216(2 ) & 101(3 ) + co@xmath582 - 3 & 0 & -2.13(5 ) & -1.64(8 ) & 2.68(8 ) & -36(2 ) & 101(3 ) + co@xmath582 - 4 & 0 & 2.13(5 ) & -1.64(8 ) & 2.68(8 ) & 216(2 ) & 101(3 ) + *      the stability of the magnetic structure against the application of an external magnetic field was verified by measuring the magnetic scattering at 1.8 k under fields of 0.5 and 2 t ( see section [ exp ] ) . despite the complex experimental setup ( cryomagnet and dispersion of the compound in deuterated isopropanol ) ,
the high hydrogen content of the compound , and the diluted character of the magnetism , weak magnetic elastic scattering could be observed ( see figure [ patterns_3h ] ) . although the bad statistics preclude refinement of these data , the absence of drastic changes in the 1.8 k magnetic structure is evidenced up to 2 t. this result is in agreement with the shape of the isothermal magnetization at 1.8 k ( see section [ mh ] and ref . ) , i.e. a fast increase at very low field due to the ferromagnetic component of the magnetization , followed by a mild increase due to competition between the applied field and the local anisotropy of co@xmath58 ions .
data were also recorded as a function of temperature for 0 , 0.5 and 2 t. figure [ 3ref_h ] summarizes the obtained data through the thermal evolution of three magnetic bragg reflections .
we can first look at the ( 10@xmath65 ) reflection , which is characteristic of the interchain correlations ( anti-*c * reflection , see figure [ magn_struct ] ) . in zero - field ( see figure [ 3ref_h ] ) , the ( 10@xmath65 ) intensity increases below @xmath7 k , the transition temperature determined from macroscopic measurements @xcite . when a magnetic field is aplied , the intensity of this reflection remains finite above @xmath66 .
this indicates the persistence of interchain correlations in presence of an applied field well above @xmath66 : up to 10 and 13 k for applied field values of 0.5 and 2 t respectively .
this is not surprising since the interchain magnetic arrangement is ferromagnetic - like , so that a magnetic field is expected to stabilize these correlations .
it is worth noting that these results are in good agreement i ) with the previously reported phase diagram established from the bias dc field dependence of the ac susceptibility maxima,@xcite and ii ) with the position of the maxima in the derivative of the temperature dependence of the magnetization : 9.8  k for 0.5 t and 13.9 k for 2 t ( data not shown).@xcite    second , we consider the ( 001 ) and ( 100 ) reflections .
these reflections are relevant for observing the evolution of the canting within the chains .
for instance , considering @xmath67 , ( 001 ) is maximal for @xmath68 while ( 100 ) is completely absent , and vice versa for @xmath69 .
the temperature dependence of these reflections depends strongly upon the magnetic field , especially at the onset of the interchain correlations when field - induced magnetic scattering is detected .
this indicates significant rearrangements of the spin structure within the chains at temperatures for which magnetic field induces long - range ordering , and points toward a rich @xmath70 phase diagram .
interestingly , at low temperature , the respective intensities of these reflections are little affected by the magnetic field up to 2 t. this suggests that the applied magnetic field only provokes a small canting with respect to the zero field magnetic structure .        ) reflections recorded upon warming under zero - field ( left ) , 0.5 t ( middle ) , and 2 t ( right ) .
the plotted magnetic neutron scattering is obtained after subtraction of a pattern recorded in zero - field at 10 k.,width=321 ]      @xmath28 hysteresis loops start to open below 2 k as shown in figures [ fig_hyst1 ] and [ fig_hyst2 ] .
the maximum coercive field is about 0.5 t at low temperature .
these hysteresis loops are strongly affected by the dynamics existing in the chains , resulting in relaxation at a given field during the extraction measurements .
this can be seen clearly in the data at 800  mk or 1 k , for example .
vs. @xmath71 for sample 1 at several temperatures between 100 mk and 1.5 k , measured with field increments of 500 oe every 110 s for @xmath72 k , and of 50 oe every 80 s for @xmath73 k. , width=302 ]    below 1 k , magnetic avalanches ( that is to say , abrupt reversals of the magnetization ) were measured in sample  1 , at fields between 0.4 and 0.5 t ( see figure [ fig_hyst1 ] ) .
these avalanches are typical of samples with slow relaxation processes at very low temperature ( see for example ref . ) and have been observed in other ordered spin - chains @xcite .
the avalanche mechanism can be described as follows : due to slow dynamics ( large energy barrier , pinned domain - walls ,  ) , when the sample is driven out of thermodynamic equilibrium , by an applied field , the magnetization does not follow the field . at a given field ( the resonant tunneling field or the nucleation field for example ) ,
the reversal of the magnetization becomes easier , and spins can start to flip .
then , because the system is at very low temperature , the specific heat is very low and the energy @xmath74 released in this reversal heats the neighboring spins , which in turn can flip .
this process results in an increase of the internal temperature of the sample , which itself induces a faster relaxation , and provokes a global reversal of the magnetization termed a magnetic avalanche . in some cases avalanches
can be suppressed by slowing the field sweeping rate in the hysteresis loop , or by improving the thermal contact between the sample and the mixing chamber of the dilution fridge which provides the cooling power .
the latter was done in the present compound with sample 2 , which was better thermalized than sample 1 ( smaller size and more apiezon grease ) .
indeed , magnetic avalanches could be suppressed in this sample ( see figure [ fig_hyst2 ] ) .
interestingly , around @xmath75 t , there is however a change of slope in the magnetization curve , showing that the magnetization starts to respond to the applied magnetic field .
this indicates that @xmath76 corresponds to an energy scale of the system , and that the previously observed avalanches trigger at a field which is characteristic of a faster relaxation of the system .
however , because the measurements were made on a powder and the magnetic structure is very complex , the origin of this energy scale is not readily determinable .     vs. @xmath71 for sample 2 at several temperatures between 100 mk and 1.5 k , measured with field increments of 1000 oe every 110 s for @xmath72 k , and of 500 oe every 110 s for @xmath73 k. , width=302 ]      as shown in ref . , and observed here in the @xmath28 hysteresis loops , co^ii^@xmath0(oh)@xmath1(c@xmath2@xmath3h@xmath2@xmath4o@xmath0)@xmath5  presents slow dynamics in the ordered regime , characterized by an energy barrier of about 60 k which is attributed to domain - wall dynamics . in order to get a better understanding of the dynamics we have conducted complementary experiments at temperatures lower than 2 k.    first , we have measured the ac susceptibility from 2.11 hz down to the lowest frequency reachable in our magnetometer , 1  mhz . in this frequency range , the ac susceptibility still shows a frequency dependence , which we characterize by tracking the maximum in the imaginary part of the susceptibility @xmath77 ( see inset of figure [ fig_tau ] ) . expecting that , at the maxima , the measurement frequency @xmath78 matches the relaxation time @xmath79 of the system , we have plotted @xmath80 as a function of the inverse of the temperature of the maximum @xmath81 , together with the higher frequency results from ref .
( see figure [ fig_tau ] ) .
these data can be fitted with an arrhenius law @xmath82 with @xmath83  s and @xmath84  k in the whole frequency range , that is to say , over 7 decades .
this shows that the mechanism behind the dynamics is robust in the ordered regime , and remains thermally activated down to at least 2.1 k. the existence of a distribution around the central barrier energy e@xmath85 was previously suggested by the shape of the cole - cole plots ( @xmath77 vs. @xmath86).@xcite the parameter extracted from the fit of these plots to a generalized debye model@xcite ( @xmath87 for t = 4 k , see ref . ) , can be used to estimate the full width at half maximum of this distribution ( @xmath88 k ) , by following the phenomelogical formalism developped by mydosh for the study of spin glasses.@xcite     vs. @xmath89 in a semi - logarithmic scale from ac susceptibility measurements .
red points are from mpms ( qd ) measurements ( 10 hz@xmath90 6.5 khz ) @xcite , and blue squares from measurements in dilution fridge magnetometers ( lt ) ( 1 mhz @xmath91 hz ) .
the line is a fit to the arrhenius law with @xmath83 s and @xmath84 k. inset : imaginary part @xmath77 of the susceptibility at low temperature .
, width=302 ]    to probe the dynamics at lower temperatures , wherein we expect a slower response , we have performed measurements of the relaxation of the magnetization as a function of time . the protocol is the following : ( i ) the sample is heated to about 2 k ( temperature at which there is no hysteresis in the magnetization ) @xcite , ( ii ) a field @xmath92  oe is applied , ( iii ) the sample is cooled towards the target temperature , ( iv ) the field is removed and the magnetization is measured as a function of time .
we note that , according to the neutron scattering experiments under applied magnetic field ( see section [ mag_struc ] ) , a field of 3500 oe does not seem to destabilize the magnetic structure .
this implies that the induced magnetization is due to the reversal of magnetic domains or to canting phenomena .
the relaxation curves obtained in zero field are shown in figure [ relax ] for temperatures between 0.9 and 1.8 k. the relaxation is not purely exponential , and can be fitted with stretched exponential functions .
this confirms the existence of a distribution of relaxation times . to investigate more precisely the barrier distribution
, we followed the approach of ref . by plotting the measured magnetization as a function of a single variable , @xmath93 .
all curves can be scaled using the parameter @xmath94 s ( see figure [ fig_mtlnt ] ) , thus confirming the picture of a distribution of barriers crossed by thermal activation .
a barrier distribution @xmath95 can then be deduced from the derivative of the magnetization with respect to @xmath96 , @xmath97 , which is plotted on the inset of figure [ fig_mtlnt ] .
vs time @xmath98 at several temperatures between 0.9 and 1.8 k. the sample is heated at about 2 k and cooled to the target temperature in an applied field of -3500 oe . at @xmath99 , the field is removed .
, width=283 ]     vs the energy @xmath100 , for 0.9 k@xmath101 k. @xmath102 s is chosen so that all the magnetization curves collapse .
inset : derivative @xmath97 vs @xmath96 showing the barrier distribution ( see text ) .
the line is a gaussian law centered on @xmath103 k with @xmath104=16.5 .
, width=302 ]    the obtained barrier distribution is almost symmetric , but rather large : the full width at half maximum @xmath105 evaluated by the fit to a gaussian function is about 10  k. surprisingly , the barrier distribution is centered on @xmath106 k , that is to say about half of the energy barrier @xmath107 obtained from ac susceptibility measurements ( in zero field ) .
in addition , the value of @xmath108 is three orders of magnitude larger than the one deduced from ac measurements .
these results indicate the existence of different relaxation processes for the two types of measurement .
this is in strong contrast with other spin - chains on which the same analysis was made @xcite , and where the obtained barrier energy was similar from ac susceptibility and relaxation of the magnetization .
vs the energy @xmath100 with @xmath109 s , obtained for relaxation curves measured in applied fields of 500 , 1000 , 1500 , 2000 and 2500 oe between 0.8 and 1.2 k. the curve at @xmath110 ( inset of figure [ fig_mtlnt ] ) is shown for comparison .
the inset gives the central barrier energy @xmath107 vs applied field @xmath71 .
the line is a fit for @xmath111 to the equation : @xmath112.,width=298 ]    we have also performed measurements of the magnetization relaxation in non - zero applied fields , from 500 to 2500 oe .
the protocol described above was used for temperatures between 0.8 and 1.2 k. note that for fields larger than 2000 oe , and at high enough temperatures , the relaxation is so fast that avalanches occur .
this phenomenon , which involves self - heating of the sample , prohibits access to the full barrier distribution at these fields .
nevertheless , the same scaling procedure was applied and led to the same @xmath108 ( @xmath113 s ) .
furthermore , the barrier distribution could be obtained ( see figure [ fig_dmdevsh ] ) .
two main features are clear from this figure . first , the distribution width is much smaller in applied field than in zero field . in the presence of an applied magnetic field
, the full width at half - height @xmath105 varies between 3.5 and 4.5 k , a value to be compared with the 10 k obtained in zero applied field .
second , the mean energy barrier decreases when the field increases .
the dependence is linear for non - zero applied fields ( see inset of figure [ fig_dmdevsh ] ) .
this is consistent with the two - level system picture used for the analysis , in which the energy barrier would be reduced by the zeeman energy : @xmath114 , where @xmath115 is the energy barrier in zero field , and @xmath116 is the magnetic moment parallel to the applied field .
the obtained @xmath115 value is smaller than the experimental @xmath117 ( @xmath118 k against @xmath106 k - see inset of figure [ fig_dmdevsh ] ) , but within the barrier distribution measured at zero field .
the slope @xmath116 is equal to 30 @xmath119 .
this might indicate that the energy barrier does not concern single spin flips but rather blocks of spins of effective moment @xmath116 . in that sense ,
the larger distribution width , as well as the difference between @xmath117 and @xmath115 , might be explained by a more complex domain structure in zero field than in finite applied field for which a domain selection occurs ( as shown by magnetization measurements ) .
the complex magnetic structure and the diversity of the exchange paths in this cobaltous chain make the interpretation of the above results challenging . nonetheless , the observed slow dynamics has a large signal , and is very clean and reproducible , suggesting that there should be a well - defined origin of this phenomenon . in the following ,
we examine the exchange paths and the relevant energy scales in the system , and present a possible scenario explaning the observations described in the previous section .      the magnetic structure results from the magnetic interactions and the local anisotropy of the co@xmath120 ions .
it is interesting to compare the anisotropy expected in the two types of environments adopted by the co@xmath120 ions . as mentioned in section [ struct ] ,
the co@xmath1201 coordination consists of distorted octahedra , whereas the co@xmath1202 ion is 5-coordinated in a trigonal bipyramid . in the former ,
the magnetic moment is almost in the basal plane , like already observed in other magnetic structures involving co@xmath120 ions in similar environments.@xcite .
the j1 exchange coupling ( see below ) could be responsible for the proposed orientation , thus implying that this coupling is stronger than the in - plane anisotropy of the co@xmath121 ions .
nevertheless , our neutron data can not exclude that the magnetic moments lie exactly in the basal plane . in the latter , co@xmath1202 ,
the anisotropy can be stronger @xcite than in the octahedral coordination : it has been observed to be even larger when the co@xmath120 ion is out of the basal plane , which is the case here ( see figure [ fig_coordination ] ) , and to be of order 25  k @xcite .    1 site , right : co@xmath1202 site .
both cations contain one hydroxide group among their ligand atoms.,width=283 ]    we examine now the exchange paths that appear to be relevant for the understanding of the intrachain magnetic structure .
their estimation is deduced from both the goodenough - kanamori rules@xcite , and the values reported in cobalt - based magnetic systems with similar exchange networks @xcite .
they are summarized in both table [ tablej ] and figure [ fig.2 ] .
the strongest magnetic coupling j1 is antiferromagnetic and induces the collinear antiferromagnetic arrangement of the co@xmath48 sublattice .
the j2 ( weak ferromagnetic ) and j3 ( weak antiferromagnetic ) interactions determine the relative orientation between the co@xmath1201 and co@xmath1202 sublattices .
the strong canting of the co@xmath122 ion with respect to the co@xmath1201 ion demonstrates the importance of the single - ion anisotropy of this ion . in that sense ,
given the anisotropy of the ions , the experimentally deduced magnetic structure satisfies the three interactions j1 , j2 and j3
. the j4 coupling , which is expected to be weakly antiferromagnetic , is frustrated in the obtained magnetic structure , thus suggesting a smaller coupling compared with j2 and j3 .    *
5c exchange & exchange & sharing & exchange & estimated + constant & path&type&type & value + j1 & co1a - oh - co1b & corner & af & [ -20,-40 ] k + j2 & co1-o - co2 & edge & f & [ 10 , 30 ] k + j3 & co1-oh - co2 & corner & af & [ -15,-25 ] k + j4 & co1-o - co2 & corner & af & [ -15,-25 ] k +    given the structural characteristics of co^ii^@xmath0(oh)@xmath1(c@xmath2@xmath3h@xmath2@xmath4o@xmath0)@xmath5 presented in section [ struct ] , the 3d magnetic lro in this material can only be explained by interchain dipolar interactions .
the energy of the dipolar interaction between two spins @xmath123 and @xmath124 separated by a distance @xmath125 is    @xmath126 \label{dipolar}.\ ] ]    in the present case , the interchain dipolar interaction is typically a few millikelvins , and implies a ratio of about 10@xmath127 between the intrachain and interchain interactions .
in such a system with a strong 1d character , the critical temperature describing lro originates from the combination of the strong intrachain correlations , resulting in spin blocks exhibiting a super - spin ( c.f . uncompensated component of the magnetic structure ) , and the interchain interactions between these super - spins @xcite . in the present material
, this picture is most readily applied along the crystallographic * c*-axis which corresponds to the shortest interchain distance . in that direction
the dipolar interaction stabilizes an interchain ferromagnetic arrangement between the intrachain super - spins .
however , the complexity of the magnetic structure makes this description much more complicated for the @xmath128 direction .
the reproduction of the magnetic motif observed in this direction is not favourable when considering the orientation of the super - spins emerging from the ferromagnetic component .
the dipolar interactions between individual spins have thus to be considered along * a*-axis .
preliminary calculations of the dipolar energy suggest that the obtained orientation is more energetically favourable with respect to the dipolar interaction than its 180@xmath129 counterpart .      in lro spin chains ,
very strong correlations exist within the chains , while the spin chains are weakly coupled to each other .
the result is that when the temperature is larger than the interchain interactions ( of a few millikelvins in our case ) , the dynamics are governed by intrachain dynamics , which can be described in terms of domain - walls motion .
this was observed in several ordered spin chains @xcite : dependent on the ratio between the local anisotropy and exchange interactions , the dynamics can be described in terms of single - chain magnets or motion of extended domain - walls .    in co^ii^@xmath0(oh)@xmath1(c@xmath2@xmath3h@xmath2@xmath4o@xmath0)@xmath5 , the presence of two co@xmath120 sites with a different single - ion anisotropy , combined with the existence of several exchange couplings , makes the description in terms of domain - walls complicated . to get a better understanding , we return to the striking result of the analysis of the dynamics : the ac susceptibility ( i.e zero - field ) dynamics is accounted for by an energy barrier of @xmath84 k over a large range of temperature , whereas the relaxation of the magnetization ( i.e. from a finite magnetization state )
can be described by a barrier of half value @xmath106 k. to go further , we will have to keep in mind that the magnetic structure is only canted by the field applied during measurements of the magnetization relaxation ( see section [ neut_field ] ) .
the zero - field dynamics measured by ac - susceptibility can be ascribed to the motion of domain - walls along the chains . due to the strong anisotropy
, we can expect that these domain - walls are relatively narrow .
then the observed energy barrier would be due to the cost in energy ( exchange and anisotropy ) to move the wall at the domain boundary .
the ranges of the exchange constants and anisotropy ( see section [ discussion1 ] ) is consistent with the obtained energy barrier .
we recall that ( i ) magnetization relaxation measurements were done after a field of -3500 oe was applied , and that ( ii ) the sample being a powder , the field acts on the three directions of the crystal .
moreover , as already mentioned , a fast increase of the magnetization is observed when a field is applied .
this can be attributed to a selection of domains with respect to the resulting spin component in the ferromagnetic - like direction , the @xmath130 direction . in the other directions , or at larger fields ,
the magnetization increases smoothly by spin canting . because of the strong anisotropy of the co@xmath1202 ions ,
this induced magnetization can be assigned to the co@xmath1201 ions .
this canting is effective only along the directions * a * and * c * , since the co@xmath1201 ions are oriented along the * b*-axis .
when the field is removed , the return to equilibrium , that is to say , the relaxation of the magnetization , will consist of two processes : the reorientation of the co@xmath1201 ions and the creation of domains to recover the zero magnetization .
it is only the second process which needs to overcome a barrier by nucleation of domain - walls inside the chains .
the obtained @xmath117 value should be associated with this nucleation .
the question that immediately arises is why @xmath117 is half of @xmath107 .
a possible reason could be that the nucleation occurs at the boundary of a chain ( or near a defect ) , where it is more energetically favourable , thus creating only one domain - wall .
on the contrary , the dynamics probed by ac susceptibility concern existing domains where the two boundaries of the domains are moving in response to the ac - field ( because the magnetization has to remain zero ) . the obtained characteristic relaxation time @xmath108 is three orders of magnitude larger in magnetization measurements than in ac susceptibility measurements , suggesting that the nucleation process is much slower .    to get a deeper insight in these dynamics , and
especially the field induced dynamics , further measurements on single crystals are needed .
these would allow us to determine if an anisotropy exists in the dynamics ( which should be the case in our scenario ) , and to provide a more quantitative analysis of these observations .
in conclusion , we have determined the magnetic structure of co^ii^@xmath0(oh)@xmath1(c@xmath2@xmath3h@xmath2@xmath4o@xmath0)@xmath5 .
it is made of spin - chains involving two sublattices : one is purely antiferromagnetic , the second gives rise to a ferromagnetic component along the * c*-axis . due to the magnetic propagation vector being * * _ k_**=(0,0,0 ) , the magnetic motif is reproduced in the neighboring unit cells forming the three dimensional lattice . the magnetic structure at low temperature
is not modified upon the application of a magnetic field , at least up to 2  t. only canting within the chains is observed .
furthermore , the applied field reinforces the interchain correlations resulting in a higher ordering temperature which increases up to 13  k at 2  t ( to be compared to the transition temperature @xmath7  k ) . at very low temperature , in the ordered phase , slow dynamics
are observed in a large temperature / frequency range which are ascribed to dynamics of domains inside the chains .
interestingly , these dynamics are extremely sensitive to the presence of a magnetic field .
finally , we notice that the nature of the present material , comprising inorganic chains interconnected through organic linkers , is reminiscent of other metal - organic frameworks that have shown high intrinsic flexibility in response to high pressures @xcite or to the absorption of guest molecules @xcite .
this point is particularly interesting as it might be possible to tune the organization of the chains relative to each other and thus the magnetic ground state . for this purpose ,
powder x - ray diffraction experiments under back pressure of co@xmath1 ( up to 100 bars ) and h@xmath1 ( up to 200 bars ) were performed . however , the crystal structure of co^ii^@xmath0(oh)@xmath1(c@xmath2@xmath3h@xmath2@xmath4o@xmath0)@xmath5did not seem to be affected .
experiments under higher pressure , using a diamond anvil cell , will be the scope of future experiments .    c. paulsen
is warmly acknowledged for the use of his magnetometers .
the authors thank the institut laue langevin , grenoble , france for provision of the research facilities .
t. hansen and a. daramsy are acknowledged for their assistance during the d20 experiments .
thanks mesr for his ph.d . funding and j. s. white for reading of the article .
we acknowledge the support of the european community - research infrastructures under the fp7 capacities specific programme , microkelvin project number 228464 .
a. k. cheetham , c. n. r. rao , and r. k. feller , chem
( cambridge ) * 2006 * , 4780 . j. r. long , and o. m. yaghi , chem .
soc . rev . *
38 * , 1213 ( 2009 ) .
g. frey , chem .
rev . * 37 * , 191 ( 2008 ) .
a. u. czaja , n. trukhan , and u. mller , chem .
soc . rev . *
38 * , 1284 ( 2009 ) .
h. li , m. eddaoudi , m. okeeffe , and o. m. yaghi , nature , bf  402 , 276 ( 1999 ) .
g. frey , c. mellot - draznieks , c. serre , f. milange , j. dutour , s. surbl , and i. margiolaki , science * 5743 * , 2040 ( 2005 ) . l. j. murray , m. dinc , and j. r. long , chem .
soc . rev . *
38 * , 1294 ( 2009 ) .
li , r. j. kuppler , and h - c .
zhou , chem .
soc . rev . *
38 * , 1477 ( 2009 ) .
j. y. lee , o. k. farha , j. roberts , k. a. scheidt , s. t. nguyen , and j. t. hupp , chem .
soc . rev . *
38 * , 1450 ( 2009 ) .
m. d. allendorf , c. a. bauer , r. k. bhakta , and r. j. t. houk , chem .
soc . rev . *
38 * , 1330 ( 2009 ) .
r. sibille , t. mazet , b. malaman , and m. franois , chem .
j. * 18 * , 12970 ( 2012 ) .
m. kurmoo , chem .
soc . rev . *
38 * , 1353 ( 2009 ) .
e. ressouche , j. x. boucherle , b. gillon , p. rey , and j. schweizer , j. am .
soc . * 115 * , 3610 ( 1993 ) .
r. feyerherm , a. loose , m. a. lawandy , and j. li , j. phys .
solids * 63 * , 71 ( 2002 ) .
r. feyerherm , a. loose , m. a. lawandy , and j. li , appl .
a : mater .
sci . process .
* 74 * , s778 ( 2002 ) .
j. l. manson , h. n. bordallo , j. w. lynn , q. huang , r. feyerherm , a. loose , l. chapon , and d. n. argyriou , appl .
a : mater .
sci . process .
* 74 * , s722 ( 2002 ) .
r. feyerherm , a. loose , p. rabu , and m. drillon , solid state sci . * 5 * , 321 ( 2003 ) .
a. zheludev , b. grenier , e. ressouche , l. p. regnault , z. honda , and k. katsumata , phys .
b * 71 * , 104418 ( 2005 ) .
e. lhotel , v. simonet , e. ressouche , b. canals , d. b. amabilino , c. sporer , d. luneau , j. veciana , and c. paulsen , phys .
b * 75 * , 104429 ( 2007 ) .
r. a. mole , j. a. stride , t. unruh , and p. t. wood , j. phys .
matter * 21 * , 076003 ( 2009 ) .
j. saines , j. r. hester , and a. k. cheetham , phys .
b * 82 * , 144435 ( 2010 ) .
a. mesbah , r. sibille , t. mazet , b. malaman , s. lebegue , and m. franois , j. mater . chem . * 20 * , 9386 ( 2010 ) .
a. mesbah , b. malaman , t. mazet , r. sibille , and m. franois , crystengcomm .
* 12 * , 3126 ( 2010 ) .
r. a. mole , j. a. stride , p. f. henry , m. hoelzel , a. senyshyn , a. alberola , c. j. gmez garcia , p. r. raithby , and p. t. wood , inorg .
50 * , 2246 ( 2011 )
. p. j. saines , h. h. m. yeung , j. r. hester , a. r. lennie , and a. k. cheetham , dalton trans .
* 40 * , 6401 ( 2011 ) .
o. fabelo , l. candillas - delgado , i. puente orench , j. a. rodriguez - velamaza , j. campo , and j. rodriguez - carvajal , inorg . chem . * 50 * , 7129 ( 2011 ) .
o. fabelo , l. candillas - delgado , j. pasan , p. diz - gallifa , c. ruiz - perez , f. lloret , m. julve , i. puente orench , j. campo , and j. rodriguez - carvajal , inorg . chem . * 52 * , 12818 ( 2013 ) .
g. nnert , h - j .
koo , c. v. colin , e. m. bauer , c. bellitto , c. ritter , g. righini , and m - h .
whangbo , inorg . chem .
* 52 * , 753 ( 2013 ) .
r. sibille , a. mesbah , t. mazet , b. malaman , s. capelli , and m. franois , j. sol . state chem .
* 186 * , 134 ( 2012 ) .
r. sibille , t. mazet , b. malaman , t. gaudisson , and m. franois , inorg .
chem . * 51 * , 2885 ( 2012 ) .
r. georges , j. j. borrs - almenar , e. coronado , curly , and m. drillon , _ one - dimensional magnetism : an overview of the models _ , in magnetism : molecules to materials i : models and experiments ( eds j. s. miller and m. drillon ) , wiley - vch verlag gmbh & co. kgaa , weinheim , frg , 2002 , p. 1 .
_ order and disorder in quantum spin systems _ , edited by t. tonegawa , s. miyashita , and h. takayama , suppl .
145 ( progress of theoretical physics , japan , 2002 ) .
f. d. m. haldane , phys .
* 45 * , 1358 ( 1980 ) .
f. d. m. haldane , phys .
b * 25 * , 4925 ( 1982 ) .
h. a. bethe , z. phys .
* 71 * , 205 ( 1931 ) .
g. gomez - santos , phys .
b * 41 * , 6788 ( 1990 ) .
y. uchiyama , y. sasago , i. tsukada , k. uchinokura , a. zheludev , t. hayashi , n. miura , and p. bni , phys .
lett . * 83 * , 632 ( 1999 ) .
a. caneschi , d. gatteschi , n. lalioti , c. sangregorio , r. sessoli , g. venturi , a. vindigni , a. rettori , m. g. pini , and m. a. novak , angew .
, int . ed .
* 40 * , 1760 ( 2001 ) . c. coulon , h. miyasaka , and r. clrac , struct .
bonding ( berlin ) * 122 * , 163 ( 2006 ) .
d. gatteschi , and a. vindingni , _ single - chain magnets _ , in molecular magnets : physics and applications ( eds j. bartolom , l. fernando , and j. f. fernandez ) , springer , 2014 , p. 191 .
j. saines , p. jain , and a. k. cheetham , chem .
* 2 * , 1929 ( 2011 ) .
r. sibille .
_ organic - inorganic hybrid solids based on dicarboxylate molecules and 3d or 4f elements : relations between structure and magnetic properties , dimensionality effects _ , thse .
vandoeuvre - ls - nancy : universit de lorraine , * 2012*. http://docnum.univ-lorraine.fr/public/ddoc_t_2012_0163_sibille.pdf r. j. glauber , j. math .
phys . * 4 * , 294 ( 1963 ) .
j. rodriguez - carvajal , physica b * 192 * , 55 ( 1993 ) .
l. b. mccusker , r. b. von dreele , d. e. cox , d. loer , and p. scardi , j. appl . crystallogr . * 32 * , 36 ( 1999 ) .
a. le bail , h. duroy , and j. l. fourquet , mater .
* 23 * , 447 ( 1988 ) .
e. f. bertaut , in _ magnetism _ , edited by g. t. rado , and h. a. shul , academic : newyork , 1963 , vol .
iii , chap .
j. rodriguez - carvajal , _ basireps : a program for calculating irreducible representation of little groups and basis functions of polar and axial vector properties .
_ part of the fullprof suite of programs , 2007 available at http://www.ill.eu/sites/fullprof ; c. ritter , solid state phenom . * 170 * , 263 ( 2011 ) . c. paulsen , in _ introduction to physical techniques in molecular magnetism : structural and macroscopic techniques - yesa 1999 _ , edited by f. palacio , e. ressouche , and j. schweizer , servicio de publicaciones de la universidad de zaragoza , zaragoza , 2001 , p. 1 . o. v. kovalev , _ representation of the crystallographic space groups , irreducible representations , induced representations , and corepresentations _ , 2nd ed .
( gordon and breach science , amsterdam , 1993 ) .
j .- g . park and c. paulsen , physics * 6 * , 55 ( 2013 ) .
e. lhotel , d. b. amabilino , c. sporer , d. luneau , j. veciana and c. paulsen , phys .
b * 77 * , 064416 ( 2008 ) .
e. lhotel , e. n. khatsko and c. paulsen , phys .
b * 74 * , 020402 ( r ) ( 2006 ) .
m. hagiwara , j. magn .
magn . mater . * 177 - 181 * , 89 ( 1998 ) .
d. huser , a. j. van duyneveldt , g. j. nieuwenhuys , and j. a. mydosh , j. phys .
c * 19 * , 3697 ( 1986 ) .
e. lhotel , v. simonet , j. ortloff , b. canals , c. paulsen , e. suard , t. hansen , d. j. price , p. t. wood , a. k. powell , and r. ballou , eur .
j. b * 86 * , 248 ( 2013 ) .
ideally , the sample should be heated over @xmath7 k , but the working conditions of the dilution fridge prevents from heating to such a high temperature . j. j. prjean and j. souletie , j. physique * 41 * , 1335 ( 1980 ) .
j. b. goodenough , _ magnetism and chemical bond _
( john wiley and sons publ . , new york,1963 ) .
e. j. kinast , c. a. dos santos , d. schmitt , o. isnard , m. a. gusmao and j. b. m. da cunha , j. alloys comp .
* 491 * , 41 ( 2010 ) .
m. c. viola , m. j. martinez - lope , j. a. alonso , j. l. martinez , j. m. de paoli , s. pagola , j. c. pedregosa , m. t. fernandez - diaz , and r. e. carbonio , chem . mat . * 15 * , 1655 ( 2003 ) .
m. e. foglio , m. c. dos santos , g. e. barberis , j. m. rojo , j. l. mesa , l. lezama , and t. rojo , j. phys
. condens .
matter * 14 * , 2025 ( 2002 ) .
l. robertson , m. duttine , m. gaudon , and a. demourgues , chem . mat . * 23 * , 2419 ( 2011 ) .
t. jurca , a. farghal , p. h. lin , i. korobkov , m. murugesu , and d. s. richeson , j. am .
soc . * 133 * , 15814 ( 2011 ) . m. drillon , and p. panissod , j. magn .
mater . * 188 * , 93 ( 1998 ) .
s. ostrovsky , w. haase , m. drillon , and p. panissod , phys .
b * 64 * , 134418 ( 2001 ) . c. coulon , r. clrac , w. wernsdorfer , t. colin , and h. miyasaka , phys . rev .
lett . * 102 * , 167204 ( 2009 ) .
a. u. ortiz , a. boutin , a. h. fuchs , and f. x. coudert , phys .
lett . * 109 * , 195502 ( 2012 ) .
k. j. gagnon , c. m. beavers , and a. clearfield , j. am .
soc . * 135 * , 1252 ( 2013 ) .
s. h. lapidus , g. j. halder , p. j. chupas , and k. w. chapman , j. am .
soc . * 135 * , 7621 ( 2013 ) . c. serre , f. millange , c. thouvenot , m. nogues , g. marsolier , d. lour , and g. frey , j. am .
soc . * 124 * , 13519 ( 2002 ) .


massive data from online social media and online social networking ( osn ) enables accurate analysis of the human behavior and of the interaction times between individuals through technology such as in word - of - mouth marketing , opportunistic networking and viral spreading .
empirical measurements contradict the commonly assumed exponential distribution of poissonian inter - event times @xcite in spreading processes or epidemic models .
multiple publications @xcite report that _ the interactivity time _ , defined as the duration between two consecutive tasks like sending emails , accessing web pages , instant messaging and phone calls , follow power - law distributions .
these findings recently led to non - markovian analyses , addressed in the work of cator _ et al . _
@xcite , iribarren and moro @xcite , van mieghem and van de bovenkamp @xcite and schweizer _ et al .
given that the inter - activity distributions are heavy tailed , word - of - mouth spreading , viral infections or dynamics of memes are expected to endure or survive longer , compared to basic markovian models @xcite .
barabasi @xcite infers that heavy - tailed distributions may arise from a priority queue , where individuals execute tasks of which the majority can be completed in short time , but some tasks wait long due to a perceived priority .
barabasi s priority queue model fits the distribution of durations between events quite well , leading to a power - law distribution with an exponent @xmath0 around 1 .
a power - law random variable @xmath1 has the probability density function @xmath2 where @xmath3 and @xmath4 is the lower bound for @xmath5 .
the probability density function ( pdf ) of a lognormal random variable @xmath5 for @xmath6 is @xmath7   } { \sigma t\sqrt{2\pi } }   \label{pdf_lognormal}\ ] ] where @xmath8 are called the parameters of the lognormal pdf , that are the mean and variance of @xmath9 as shown in appendix [ sec_properties_lognormal ] .
when we assume that inter - event durations are power - law distributed , we encounter the following issues :    1 .   in many cases , only a part of the data ( the tail larger than @xmath10 ) is modeled by a power - law .
the lower bound @xmath10 in ( [ pdf_powerlaw ] ) does not correspond to the physical minimum of the random variable @xmath5 , but @xmath10 is fitted from the data by ignoring smaller values that do not obey the power - law .
often , these smaller values may have a large probability to occur , so that their neglect is difficult to justify .
in other words , only a part of the process ( above @xmath10 ) is modeled by a power - law ( [ pdf_powerlaw ] ) , while the other part ( below @xmath10 ) is not .
most processes or measurements possess both a lower as well as an upper bound .
apart from the lower bound @xmath10 , an upper bound @xmath11 is often invoked , at which a cut - off is observed : the power - law behavior is confined to the range @xmath12 , although @xmath13 and @xmath14 . however , it is often unclear whether the process in the deep tail still obeys a power - law distribution or some other , much faster decreasing distribution .
the upper bound @xmath11 is usually empirically determined , rather than based on the physical maximum of @xmath5 .
as long as @xmath15=c\int_{\kappa}^{\infty}t^{-\gamma}dt=\left (   \frac{\kappa}{\tau}\right )   ^{1-\gamma}\ ] ] is small ( with respect to the measurement precision ) , the upper bound @xmath11 is justified , else other validation arguments are needed .
3 .   we demonstrate here that the binning of data ( either by the data - provider or by the researcher ) alters the shape of a lognormal distribution into an apparent power - law .    particularly in relation to human activities or behaviors
, we question in this paper the widely assumed power - law distribution .
we present measurements of inter - event durations from digg.com and reddit.com , two online social news aggregators @xcite , and from the enron data set , a collection of emails sent by employees of the company enron , and argue that a lognormal distribution is a valid candidate for the distribution of human inter - event durations in section [ sec_measurements ] .
the problem of fitting a lognormal is explained in section [ sec_fittinglognormal ] , followed by previously reported lognormally distributed data sets in section [ sec_reported_parameters ] .
existing theoretical models are compared in section [ sec_debate ] , a plausible interpretation for lognormal human behavior is proposed in section [ sec_plausible_explanation_lognormal_interactivity_times ] and section [ sec_conclusion ] concludes .
mathematical properties of the lognormal distribution are deferred to the appendix [ sec_properties_lognormal ] .
appendix [ sec_likelihood_testing ] presents results of likelihood tests for the observed distributions .
all events in an osn are based on users activities , such as posts , friend - requests and comments . a complete data set including all activities of users from digg.com for a duration of 4 years , described in tang _
_ @xcite and doerr _ et al . _
@xcite , allows us to analyze the time frame in which users of digg.com add their friends .
_ @xcite found that reaction times in a retweet network from twitter and digg are close to a lognormal distribution with parameters @xmath16 and @xmath17 .
the process of adding friends shows bursts of activity also observed and analyzed in email communication @xcite . in these publications the question arises , whether the observed distribution of inter - activity times is described by a power - law , lognormal or a cascading poisson process .
@xcite describe that circadian and weekly activity cycles in human behavior are the factors that lead to heavy tails .
they proposed a cascading poisson process , consisting of a _ nonhomogeneous _ poisson process that reflects periodicity and a _ homogeneous _ poisson process describing active intervals , which is shown to model the interactivity distributions of e - mail communication .    because the network of digg.com is directed ( like in twitter.com or other osns ) ,
a user can be followed by other users to become their friend , while a user can not add followers .
this means that the process of adding friends is solely based on the user himself , whereas obtaining followers depends on the activities of other users .
the random variable @xmath18 denotes the time between the addition of two friends and , similarly , @xmath19 is the time between receiving two followers .        figures [ fig_interfriendaddtime_seconds ] and [ fig_binsize1s_followers ] depict the distribution of durations between adding friends @xmath18 and of receiving followers @xmath19 for 7.4 million friendship relations in digg . fitting the data , i.e. all realizations of @xmath18 , by the state of the art technique by clauset _
_ @xcite to a power - law ( [ pdf_powerlaw ] ) results in an exponent @xmath20 for @xmath21 , whereas all realizations of @xmath19 are fitted best by a lognormal ( [ pdf_lognormal ] ) with parameters @xmath22 and @xmath23 . for @xmath19
an estimated p - value of 0.0 indicates that a power - law fit provides not the best solution , whereas the p - value for @xmath18 of @xmath24 indicates a reasonable fit for a power - law distribution .
the reason that the distribution of @xmath19 does not fit the lognormal distribution over the whole range lies in the nature of digg.com .
@xcite showed that only a few users were active over a long period , while just a fraction of them was actually submitting stories .
only 2% of all registered users succeeded to have their submissions promoted  to the frontpage @xcite .
since the username of a submitter appears next to the story , these users receive a lot of followers during the relatively small period that their story was listed on the front page .
therefore , @xmath25 is large for small @xmath26 in fig .
[ fig_binsize1s_followers ] . in summary , even ignoring the small time values , the random variable @xmath18 and @xmath19 possess different distributions . the process that generates @xmath18 is more likely described by a variant of barabasi s priority queue model , leading to power - law behavior .
since the random variable @xmath19 is generated by the collective dynamics of different individuals ( that are likely weakly dependent ) , central limit law arguments may point towards the lognormal distribution @xcite .
similar properties occur in sending and receiving emails .
obviously , a user can only send an email when he is online , whereas emails arrive at a user s inbox at his email server at all times .
we analyzed the durations between receiving and sending emails in the enron data set , which contains emails of all employees of enron during roughly 6 years , starting in january 1998 until february 2004 .
similar distributions arise as shown later in figs .
[ fig_enron_sent_diff ] and [ fig_enron_received_diff ] .    a third data set from reddit.com , an osn in which users mainly submit , comment or vote on bookmarks .
the pdf of the duration @xmath27 between consecutive comments and submissions in reddit are shown in fig .
[ fig_reddit_comment_and_submission ] .
the tails of the distributions in fig .
[ fig_reddit_comment_and_submission ] seem nicely fitted by a power - law distribution with exponents @xmath28 .
however , as mentioned in the introduction , the increasing regime in @xmath29 for small values of @xmath30 , the peak of @xmath31 nor the concave form of the pdf can be modeled by a power - law distribution .
the two main approaches are based on the pdf and the edf ( empirical distribution function ) , after a logarithmic transformation of the data .
figure  [ fig_ccdf_seconds_normal_followers ] depicts the data fitted to the distribution function of a normal distribution ( [ distribution_lognormal ] ) and shows that the parameters of the edf are about the same as those in fig .
[ fig_binsize1s_followers ] ( pdf approach ) .
( same data as in fig . [ fig_binsize1s_followers ] ) . ]
we generate @xcite @xmath32 realizations of a lognormal random variable @xmath5 , which we denote by the set @xmath33 @xmath34 .
all @xmath32 realizations lie in the interval @xmath35   $ ] , where the minimum value is @xmath36 and the maximum value is @xmath37 .
the binning operation @xcite divides the entire data interval @xmath38   $ ] into @xmath39 sub - intervals of length @xmath40 and the @xmath41-th subinterval @xmath42   $ ] for @xmath43 is associated with a bin @xmath44 , that contains the number of realizations of @xmath5 or data points of set @xmath45 lying within the @xmath41-th subinterval,@xmath46 clearly , for a given set @xmath45 , increasing the binsize @xmath47 decreases the number of bins @xmath39 .
the effect of binning on the pdf @xmath48 is depicted in fig .
[ fig_different_binning_m10_s2_lognormnoise ] , where @xmath49 realizations from a lognormal random variable with parameters @xmath50 and @xmath51 are drawn . as demonstrated in appendix [ sec_properties_lognormal ] , the scaled random variable @xmath52 has parameter @xmath53 , but @xmath54 , implying that binning  only changes the parameter @xmath55 , but leaves the parameter @xmath56 invariant ! by binning ( scaling ) the distribution with different binsizes , the up - going regime , where @xmath48 increases with @xmath30 , disappears and the observable part of the distribution evolves  towards a straight line on a log - log plot . binning the data with larger binsizes
decreases the parameter @xmath55 , even to the extent that @xmath55 may become negative .
if @xmath57 decreases , the maximum of the lognormal at @xmath58 tends to zero ( infinitely far to the left on a log - log scale ) .
therefore , just the decreasing part of the quadratic shape of ( [ pdf_lognormal ] ) will be visible in a log - log plot .
binned per minute . ]
binning the data in different time units , say per minute instead of per second , scales the distribution by a factor of 60 ( since 60 seconds equals 1 minute ) , which will shift the parameter @xmath59 of the cdf towards the left to @xmath60 .
however , a remarkable property of the lognormal distribution is that the parameter @xmath56 will not change after linear scaling as shown in the appendix [ sec_properties_lognormal ] .
binned per hour . ]
an important consequence of a binning operation is illustrated in the pdfs in fig .
[ fig_pdf_minutes_lognormal_followers ] and fig .  [ fig_pdf_hours_lognormal_followers ] , that show the pdf of @xmath19 binned per minute and per hour , respectively .
conforming to the theory , the parameter @xmath55 decreases by a factor of about @xmath61 , while the parameter @xmath56 keeps its value , but the shape of the distribution changes . the distribution , binned per hour , is shown in fig .
[ fig_pdf_hours_lognormal_followers ] , where the distribution can be confounded with a power - law distribution ( with exponential cut - off ) .
moreover , this misinterpretation may be justified , because the exponent @xmath62 of a power - law fitted to the hourly binned data is close to the exponent @xmath20 found for the @xmath18 distribution in fig .
[ fig_interfriendaddtime_seconds ] .
[ fig_pdf_hours_lognormal_followers_exponential ] shows the complementary cumulative distribution function of @xmath19 binned per hour , fitted using the method of clauset _ et al .
_ @xcite . the resulting power - law exponent in fig .
[ fig_pdf_hours_lognormal_followers_exponential ] is @xmath63 .
the visualization in fig .
[ fig_pdf_hours_lognormal_followers_exponential ] is misleading because the original distribution is , as shown earlier , a lognormal distribution .     where data is binned per hour . ]    for completeness we like to refer to table [ tab : estimatedparametersusingthemethodofclausetetal ] in appendix [ sec_likelihood_testing ] , where we list fitted parameters , p values and the result of log - likelihood tests , conducted using the method of virkar and clauset @xcite .
as shown in the table , the larger the used binsizes , the higher the p values indicating that a power - law fits the data . for small binsizes ( low p values )
the results are not conclusive or rather in favor of a lognormal distribution .    to further demonstrate the effect , figs .
[ fig_pdf_minutes_lognormal_followers ] , [ fig_pdf_hours_lognormal_followers ] and [ fig_pdf_hours_lognormal_followers_exponential ] depict distributions in which the data was artificially binned by minutes and hours . in reality
, such binning operation might occur if the data provider ( the osn or web service ) returns values per hour or even larger scales .
in other words , data measured in larger time steps over a certain total time range thus erases the possibility to distinguish between a lognormal and power - law distribution .
the sending time of an email in the enron data set until may 2001 was stored per minute , whereas afterwards ( june 2001 - february 2004 ) the sending time is stored with an accuracy of a second .
the resulting pdf of time differences between sent emails is plotted in fig .
[ fig_enron_sent_diff ] , where the black and red dots represent the measurements per second and per minute , respectively .
the parameters of the fitted lognormal ( @xmath64 ) and power - law distribution @xmath65 ) correspond to those of digg.com in fig .
[ fig_binsize1s_followers ] and fig . [
fig_interfriendaddtime_seconds ] .
the pdf of inter - arrival times of received emails , split in observations per minute ( red ) and per second ( black ) in fig .
[ fig_enron_received_diff ] , can be compared to the durations of added followers , where the parameters are again in the same range as shown before .
we found that the edf of the data ( not shown ) agrees with the pdfs in fig .
[ fig_enron_sent_diff ] and fig .
[ fig_enron_received_diff ] , underlining the likeliness that a lognormal random variable models the data better than a power - law .
table [ tab : powerlawexponentsfoundinhumandynamics ] lists inter - arrival times and the fitted exponents found in publications in which the aspect of binning was ignored .
rather small power - law exponents @xmath0 between @xmath66 and @xmath67 are found in multiple data sets of human activity .    [ c]lll data - set & exponent ( @xmath0 ) & reference + time between sending emails & 1 & barabasi ( 2005 ) @xcite + time between sending emails & 1.2 & eckmann _ et al . _
( 2004 ) @xcite + time between clicks in website usage & 1 , 1.25 & gonalves and ramasco ( 2008 ) @xcite + time between similar actions in web data & 1.1 , 1.2 , 1.8 & radicchi ( 2008 ) @xcite + time between messages in instant - messaging & 1.53 & leskovec and horvitz ( 2008 ) @xcite + time between phone calls & 0.9 & candia _ et al . _
( 2008 ) @xcite + time between phone calls & 0.7 & karsai _ et al . _
( 2012 ) @xcite + time between sending emails & 1 & karsai _ et al . _
( 2012 ) @xcite + time between sending short messages & 0.7 & karsai _ et al . _
( 2012 ) @xcite    lognormal distributions were reported for similar human activity as shown in table [ tab : addlabel ] , which extends the collection of limpert _ et al .
the oldest analysis was conducted by boag @xcite in 1949 .
based on the scaling invariance of @xmath56 as demonstrated in appendix [ sec_properties_lognormal ] , the parameter @xmath56 can be compared over different measurements , whereas the parameter @xmath55 can not , since @xmath55 depends on the units in which the lognormal random variable is measured .
interestingly , all @xmath56 s in table [ tab : addlabel ] lie within an amazingly small range of @xmath68 , which shows that the parameter @xmath56 only varies over about one order of magnitude in different measured phenomena . if @xmath56 is rather small as in our measurements where @xmath69 , the first term in ( [ loglogform ] ) dominates and a quadratic function appears in a log - log plot . consequently , these rather small values of @xmath56 and those in table [ tab : addlabel ] contradict the common deductions made from ( [ loglogform ] ) in appendix [ sec_properties_lognormal ] , namely that only for large values of @xmath56 power - law and lognormal distributions are indistinguishable . equations ( [ pdf_lognormal_power_law_like_form ] ) and ( [ power_law_exponent_lognormal ] ) in appendix [ sec_properties_lognormal ] explain why the lognormal distribution may seem linear in a certain regime .
[ c]llll @xmath55 & @xmath56 & process & reference + 5.547 & 2.126 & email forwarding & iribarren and moro @xcite + @xmath70 & @xmath71 & email forwarding & stouffer _ et al . _
@xcite + @xmath72 hour & & & + @xmath73 days & & email forwarding & stouffer _ et al .
_ @xcite + 2.47 & 0.38 & infection times & nishiura @xcite + 2.47 & 0.36 & latency periods of diseases & limpert @xcite + 14 days & 1.14 & latency periods of diseases & sartwell @xcite + 100 days & 1.24 & latency periods of diseases & sartwell @xcite + 2.3 hours & 1.48 & latency periods of diseases & sartwell @xcite + 2.4 days & 1.47 & latency periods of diseases & sartwell @xcite + 12.6 days & 1.50 & latency periods of diseases & sartwell @xcite + 21.4 days & 2.11 & latency periods of diseases & sartwell @xcite + 9.6 months & 2.5 & survival times after cancer diagnosis & boag @xcite + 15.9 monts & 2.8 & survival times after cancer diagnosis & feinleib and macmahon @xcite + 17.2 months & 3.21 & survival times after cancer diagnosis & feinleib and macmahon @xcite + 14.5 months & 3.02 & survival times after cancer diagnosis & boag @xcite + 60 years & 1.16 & age of onset of alzheimer & horner @xcite + 4 days & & incubation periods ( viral infections ) & lessler _ et al . _
@xcite + 3 to 5 & @xmath71 & task completion & linden @xcite + 0.5 & 1.4 & strike duration & lawrence @xcite + & & time of individual activities & mohana _ et al .
_ @xcite + 0.43 & 1.634 & call duration & spedalieri _ et al .
_ @xcite + 3.5 & 0.70 & message holding time & barcelo and jordn @xcite + 7.439 & 0.846 & transmission holding time & barcelo and jordn @xcite + 3.29 & 0.890 & channel holding time & barcelo and jordn @xcite + 3.3 & 0.89 & channel holding time & barcelo and jordn @xcite + @xmath74 & @xmath75 & & + @xmath76 & @xmath77 & call holding time & bolotin @xcite + & & citations & eom and fortunato @xcite + & & citations & redner @xcite + 1 to 2 & @xmath78 & citations & stringer _ _ et al.__@xcite + & 1.095 & citations & radicchi _ _ et al.__@xcite + @xmath79 & @xmath80 & & + @xmath81 & @xmath82 & retweeting behavior & doerr _ _ et al.__@xcite + 5.29 & 0.42 & distribution of votes on pages of digg.com & van mieghem _
@xcite + & & &
mitchell @xcite discusses the fascinating and abundant appearance of power - law distributions and mentions processes such as self - criticality , that are related to phase transitions , as the main producers of power - laws , but she remarks that the precise nature as well as the deviations from a power - law are still open to debate .
table [ tab : addlabel ] illustrates that the research on lognormal distributions has a long history .
unfortunately , also the processes that produce ( almost ) lognormal behavior are not well understood .
mitzenmacher @xcite overviews processes leading to power - law and lognormal distributions , emphasizing that minor changes in the process will lead to either a power - law or a lognormal distribution .
mitzenmacher mentions the work of gabaix @xcite , who analyzed the size distribution of cities in the united states .
interestingly , gabaix @xcite found that the city size distribution follows a zipf distribution , which is a power - law distribution ( [ pdf_powerlaw ] ) with an exponent @xmath83 .
gabaix argues that cities can not become infinitely small , a fact that imposes a lower bound to their size .
when modeling the size of cities as a markov chain with a fixed number of cities , which grow stochastically as proposed by gibrat @xcite , then the steady state of the markov chain will follow zipf s distribution with an exponent @xmath83 .
if there is no lower bound on the city size , so that cities can be arbitrarily small , then the distribution degenerates to a lognormal distribution .
gibrat @xcite , whose work is often associated with the law of proportionate effect @xcite , argues in his model that transition probabilities or the variance of transition probabilities in a growth process are independent from the size .
gibrat , who estimated the distribution of city sizes , concluded that the law of proportionate growth leads to a distribution that is lognormal .
however , simon @xcite showed that gibrat s law of the proportionate effect may lead to other heavy - tailed distributions as well .
champernowne @xcite and cordoba @xcite analyzed the income distribution of england and wales with markov theory .
again , their crucial assumption is that incomes have a lower bound .
eeckhout @xcite analyzed the distribution of city sizes by using accurate data from the us census in 2000 and found that the heavy tail obeys zipf s law .
but , the entire distribution is better described by a lognormal than a pareto distribution . by comparing the census data from 1990 and 2000
, eeckhout shows that the growth of a city is independent of its size .
the parameters of the lognormal distribution found by eeckhout are @xmath84 and @xmath85 .
the difference with the above mentioned markov chain approach lies in the fact that eeckhout @xcite modeled the process by a multiplicative process , which leads to a lognormal distribution .
this multiplicative process , proposed by kapteyn @xcite in 1903 for the first time and later coined the law of proportionate effect  by gibrat @xcite , is based on the central limit theorem applied to a multiplicative process , which leads to lognormal distributed sizes @xcite .
gibrat s growth process is defined as @xmath86 where the size @xmath87 of an item at state @xmath30 depends upon the previous size @xmath88 times a positive , random factor @xmath89 . by taking the logarithm of both sides in ( [ gibrat_governing_eq ] ) and denoting @xmath90
, we obtain , after iteration , @xmath91 by the central limit theorem @xcite , @xmath92 , we arrive at , for large @xmath30,@xmath93\rightarrow\frac{1}{\sqrt{2\pi}}\int_{-\infty}^{\frac{y - t\mu}{\sigma\sqrt{t}}}e^{-u^{2}/2}\;du\ ] ] from which approximately @xmath94 where @xmath55 denotes the mean and @xmath95 the variance of the sequence @xmath96 .    for large @xmath30 , a power - law distribution ( [ pdf_powerlaw ] ) tends to zero as @xmath97 , while a lognormal distribution ( [ pdf_lognormal ] ) tends considerably faster as @xmath98   \right )   $ ] to zero , illustrating that the deep tails are significantly different .
malevergne @xcite addressed this fact by the concept of a slowly varying function , in particular , for @xmath99 and @xmath100 , the power - law distribution ( [ pdf_powerlaw ] ) features @xmath101 the lognormal distribution ( [ pdf_lognormal ] ) on the other hand is not slowly varying . in the limit @xmath99
, a lognormal distribution will , for @xmath102 , always tend to zero : @xmath103 this different tail behavior questions whether exponential bin sizes should be used in a log - log plot , because these may hide the rapid decrease in the tail .
perhaps , a plausible explanation of the appearance of lognormal human interactivity times is that the logarithm @xmath104 of a human interactivity time @xmath27 is gaussian or normally distributed ( see appendix [ sec_properties_lognormal ] ) , most likely as a consequence of the central limit theorem . roughly speaking
, the central limit theorem applies for a large number of weakly dependent random variables , where none of them is dominant .
thus , rather than concentrating on the time @xmath27 , it seems more natural to focus on the random variable @xmath104 as the decisive quantity .
the logarithm of a human related measure often occurs : for example , in our hearing system , the intensity of sound is logarithmically experienced .
but more importantly , the large variability in human performance also seems logarithmically distributed @xcite and the measured interactivity times are strongly related to the large differences in human performance and/or behavior .    on the other hand , if @xmath105 is an exponential random variable with mean @xmath106 for @xmath107 ( else @xmath108 ) and @xmath109 , then @xmath110   = \pr\left [   y\leq\log t\right ]   = 1-e^{-\left ( \gamma-1\right )   \left (   \log t-\log\tau\right )   } = 1-\left (   \frac{t}{\tau } \right )   ^{1-\gamma}\ ] ] so that @xmath111   = \frac{\gamma-1 } { \tau^{1-\gamma}}t^{-\gamma}$ ] is a perfect power - law probability density function as in ( [ pdf_powerlaw ] ) .
however , if @xmath104 @xmath112 were exponentially distributed with mean @xmath106 , the memoryless property of the exponential distribution @xcite would indicate that @xmath113   = \pr\left [   \log t\geq t\right]\ ] ] in other words , given that the logarithm of a human interactivity time is larger than @xmath114 time units , the probability that @xmath104 exceeds @xmath115 time units is equal to the probability that @xmath104 exceeds @xmath30 time units , for any @xmath114 and @xmath30 larger than @xmath116 , precisely as if the log - interactivity time @xmath114 never has been spent or waited , which is quite counterintuitive for a duration between consecutive activities .
alternatively , with @xmath117 and @xmath118 , the memoryless property of @xmath104 implies scale - freeness in @xmath27 : @xmath119   = \pr\left [   \log t\geq t\right]\ ] ] which shows , ignoring the condition @xmath120 , the independence on the scale @xmath121 .
but , we have shown in section [ sec_effect_binning ] that rescaling the human interactivity time ( by different bin sizes ) definitely alters the distribution .    from these arguments
, we can infer that a power - law time @xmath27 is less defendable than a lognormally distributed @xmath27 .
in this paper , the interactivity durations of individuals , between creating friendship relations , writing emails , commenting and voting on online content are analyzed .
we found that the distribution of durations to add friends follows a power - law with an exponent of @xmath122 , whereas the durations to acquire followers are well described by a lognormal with @xmath123 and @xmath124 . due to the small probability of executing two tasks in a small time interval ( typically ignored in fitting a power law )
, we claim that a lognormal distribution covers the entire activity time range better than a power distribution .
in addition , we show that binning of lognormally distributed data can seriously affect the perception : the parameter @xmath55 shifts towards smaller values , but the parameter @xmath56 of a lognormal distribution does not change after a binning or scaling operation . in the extreme case
, only the heavy tail of the lognormal distribution ( [ pdf_lognormal ] ) remains , which follows a power - law distribution ( [ pdf_powerlaw ] ) with an exponent of @xmath0 close to 1 .
as explained in the appendix , there exists an interval in which the lognormal distribution is indistinguishable from a power - law distributions with power - law exponent @xmath125 , for small @xmath126 .
similar observations and concerns , discussed for a long time in the literature of city size and income distributions and reviewed in section [ sec_debate ] , supports that a lognormal distribution is modeling the whole data range better than power - laws .
finally , section [ sec_plausible_explanation_lognormal_interactivity_times ] argues that the logarithm of a quantity associated with human activities rather than the quantity itself is the better descriptor , because human performance seems to fit a lognormal .
99 p.  van  mieghem . .
cambridge university press , 2014 .
barabsi . . , 435:207211 ,
may 2005 .
eckmann , e.  moses , and d.  sergi .
entropy of dialogues creates coherent structures in e - mail traffic .
, 101(40):1433314337 , 2004 .
b.  gonalves and j.  j. ramasco . .
, 78.2026123 , 2008 .
f.  radicchi . . , 80.026118 ,
aug 2009 .
j.  leskovec and e.  horvitz . .
in _ proceeding of the 17th international conference on world wide web _ , pages 915924 .
acm , 2008 .
j.  candia , m.  c.  gonzlez , p.  wang , t.  schoenharl , g.  madey , and a .-
barabsi . uncovering individual and collective human dynamics from mobile phone records .
, 41(22):224015 , 2008 .
m.  karsai , k.  kaski , a .-
barabsi , and j.  kertsz . . , 2:397 ,
may 2012 .
e.  cator , r.  van  de bovenkamp , and p.  van  mieghem .
susceptible - infected - susceptible epidemics on networks with general infection and cure times . ,
87:062816 , jun 2013 .
j.  l.  iribarren and e.  moro . .
, 103(3):038702 , jul 2009 .
p.  van  mieghem and r.  van  de bovenkamp . .
, 110:108701 , mar 2013 .
i.  scholtes , n.  wider , r.  pfitzner , a.  garas , c.  j.  tessone , and f.  schweitzer .
causality - driven slow - down and speed - up of diffusion in non - markovian temporal networks .
september 2014 .
m.  karsai , m.  kivel , r.  k.  pan , k.  kaski , j.  kertsz , a.  l.  barabsi , and j.  saramki . . , 83(2):025102 ,
s.  tang , n.  blenn , c.  doerr , and p.  van  mieghem . .
, 13(5):11631175 , 2011 .    c.  doerr , n.  blenn , s.  tang , and p.  van  mieghem . .
in _ computer communications _
, volume  35 , pages 796809 , 2012 .
p.  van  mieghem .
human psychology of common appraisal : the reddit score .
, 13(6):14041406 , 2011 .
r.  d.  malmgren , d.  b.  stouffer , a.  e.  motter and l.  a.  n.  amaral .
, 105(47):18153 , 2008 .
r.  d.  malmgren , d.  b.  stouffer , a.  s.  l.  o.  campanharo and l.  a.  n.  amaral on universality in human correspondence activity .
, 325(5948):16961700 , 2009 .    c.  doerr , n.  blenn , and p.  van  mieghem . . ,
8(5):e64349 , 05 2013 .
a.  clauset , c.  r.  shalizi , and m.  e.  j. newman . . , 51(4):661703 ,
e.  limpert , w.  a.  stahel , and m.  abbt .
log - normal distributions across the sciences : keys and clues . , 51(5):341352 , 2001 .
j.  w.  boag . . ,
11(1 ) , 1949 .
d.  b.  stouffer , r.  d.  malmgren , and l.  a.  n.  amaral . .
arxiv : physics/0510216 , 2005 .
d.  b.  stouffer , r.  d.  malmgren , and l.  a.  n.  amaral .
log - normal statistics in e - mail communication patterns , arxiv : physics/0605027 2006 .
h.  nishiura . . , 2007
p.  e.  sartwell . .
, 51:310318 , 1950 .
m.  feinleib and b.  macmahon .
variation in the duration of survival of patients with the chronic leukemias .
, 15 , 1960 .
r.  d.  horner .
, 126(3):40914 , 1987 .
j.  lessler , n.  g.  reich , r.  brookmeyer , t.  m.  perl , k.  e.  nelson , and d.  a.  cummings .
, 9(5):291300 , may 2009 .
w.  j.  van  der  linden . .
, 31(2 ) , 2006 .
r.  j.  lawrence . . ,
147(3 ) , 1984 .
s.  mohan , m.  gopalakrishnan , h.  balasubramanian , and a.  chandrashekar . .
, 58(6):827831 , 2006 .
a.  spedalieri , i.  martn - escalona , and f.  barcel . in _ wireless communications and networking conference _ , pages 23812386 .
ieee , 2005 .
f.  barcel and j.  jordn . .
electronics letters , vol.34 no.2 _ , pages 146147 , 1998 .
v.  a.  bolotin . .
, 12:433438 , 1994 .
eom and s.  fortunato . . , 6(9):e24926 ,
s.  redner . .
, 58(6):4954 , 2005 .
m.  stringer , m.  sales - pardo , and l.  a.  n.  amaral .
effectiveness of journal ranking schemes as a tool for locating information .
, 3(2):e1683 , 2008 .
p.  van  mieghem , n.  blenn , and c.  doerr . .
, 83:251261 , 2011 .
10.1140/epjb / e2011 - 20124 - 0 .
m.  mitchell , _ complexity : a guidede tour _ , oxford university press , 2009 .
m.  mitzenmacher . .
, 1(2 ) , 2003 .
x.  gabaix . .
, 114:114739 , 1999 .
r.  gibrat . .
librairie du recueil sirey , 1931 , french .
h.  a.  simon . .
, 42(34):425440 , 1955 .
d.  g.  champernowne . . ,
318351 , 1953 .    j.  c.  cordoba . . , 49(4):14631468 , 2008 .
j.  eeckhout . .
, 94(5):14291451 , september 2004 .
j.  c.  kapteyn . .
, 19(3):205206 , 1918 .
y.  malevergne , v.  pisarenko , and d.  sornette .
. swiss finance institute research paper series 09 - 40 , swiss finance institute , 2009 .
w.  shockley , on the statistics of individual variations of productivity in research laboratories , proceedings of the institute of radio engineers ( ire ) , vol .
279 - 290 , 1957 .
y.  virkar and a.  clauset .
power - law distributions in binned empirical data .
, 8(1):89119 , 03 2014 .
f.  black and m.  s. scholes . .
, 81(3):63754 , may - june 1973 .
d.  garcia , p.  mavrodiev , and f.  schweitzer .
social resilience in online communities : the autopsy of friendster . in _ proceedings of the first acm conference on online social networks _ , acm , february 2013 .
h.  kesten . . ,
cxxxi:207248 , 1973 .
f.  radicchi , s.  fortunato , and c.  castellano . . ,
105(45 ) , 2008 .
a lognormal random variable @xcite is defined as @xmath109 where @xmath127 is a gaussian or normal random variable .
hence , @xmath128 .
the distribution function @xmath129=\pr[y\leq\log t]$ ] is @xmath130   du=\frac{1}{2}\left (   1+\text{erf}% \left (   \frac{t-\mu}{\sigma\sqrt{2}}\right )   \right ) \label{distribution_lognormal}%\ ] ] where erf(@xmath26 ) is the error function .
the probability density function ( pdf ) of a lognormal random variable @xmath5 follows from the definition @xmath131 $ ] , for @xmath6 , as ( [ pdf_lognormal ] ) , where @xmath8 are called the parameters of the lognormal pdf , while the mean and variance are @xcite @xmath132   = e^{\mu}e^{\frac{\sigma^{2}}{2}}\ ] ] and @xmath133   = e^{2\mu}e^{\sigma^{2}}\left (   e^{\sigma^{2 } } -1\right)\ ] ] the limit @xmath134 reduces to a dirac delta function at @xmath135 , thus @xmath136 .
given the mean and variance , the parameters of the lognormal are found as @xmath137   } { \left (   e\left [ x\right ]   \right )   ^{2}}\right )   \label{parameter_sigma^2}%\ ] ] and @xmath138   -\frac{\sigma^{2}}{2 } \label{parameter_mu}%\ ] ] although @xmath139   \geq0 $ ] , we remark that the parameter @xmath55 can be negative .
moreover , ( [ parameter_sigma^2 ] ) and ( [ parameter_mu ] ) show that the scaled lognormal random variable @xmath52 , where @xmath140 is a positive real number , has mean @xmath141 and @xmath142 .
hence , _ scaling _ by a factor @xmath140 does not change the parameter @xmath56 , which has interesting consequences for binning and measured data : the unit ( e.g. second versus hours ) in which the random variable is measured does not alter the parameter @xmath56 , only the parameter @xmath55 .
the change of the argument @xmath143 in @xmath48 leads to @xmath144   } { \sigma\sqrt{2\pi } } \label{pdf_lognormal_e^u}%\ ] ] illustrating that the scaled lognormal pdf @xmath145 is a gaussian pdf @xmath146 with mean @xmath147 .
the maximum of @xmath48 occurs at @xmath148 and equals @xmath149 , which follows directly from ( [ pdf_lognormal_e^u ] ) .
moreover , we find easier from ( [ pdf_lognormal_e^u ] ) than from ( [ pdf_lognormal ] ) that @xmath150 and that @xmath151 .
this means that any lognormal starts at @xmath152 from zero , increases up to the maximum at @xmath153 after which it decreases towards zero at @xmath154 .
thus , the lognormal is bell - shaped , but , in contrast to the gaussian , the lognormal pdf is not symmetric around its maximum at @xmath148 and can be seriously skewed .    the expression for the lognormal pdf in ( [ pdf_lognormal ] ) can be rewritten @xcite in a power - law - like form as @xmath155 where the exponent @xmath156 equals@xmath157 which illustrates that a lognormal random variable behaves as a power - law random variable , provided the last fraction in ( [ power_law_exponent_lognormal ] ) is negligibly small , say @xmath158 .
the latter happens when @xmath159 .
thus , when @xmath160 $ ] or in terms of the lower bound @xmath161 and upper bound @xmath162 defined in section [ sec_introduction ] , the pdf ( [ pdf_lognormal ] ) of a lognormal random variable is almost indistinguishable from the pdf ( [ pdf_powerlaw ] ) of a power - law random variable with power exponent @xmath125 . the @xmath30-interval @xmath163 $ ] exceeds the maximum @xmath58 of the lognormal pdf and is clearly longer when @xmath56 is larger ( as well as the tolerated accuracy @xmath158 increases ) .
_ @xcite mention that , if @xmath164 , then the exponent @xmath156 in ( [ power_law_exponent_lognormal ] ) varies less than 0.3 units over a range of three orders of magnitude .
taking the logarithm on both sides of ( [ pdf_lognormal ] ) results in @xmath165 if @xmath56 is large , then the second term @xmath166 in ( [ loglogform ] ) dominates , which leads to a straight line with in a log - log plot resembling a power - law with exponent @xmath83 . on the other hand ,
if @xmath56 is small or @xmath167 , the first , quadratic term in ( [ loglogform ] ) dominates .
we apply the method of clauset _ et al . _ @xcite and virkar and clauset @xcite to our data to distinguish between power - law and lognormal distributions .
@xcite approached the problem of estimating the exponent @xmath0 in ( [ pdf_powerlaw ] ) by testing different distributions . in their data , lognormals and power - laws
were not clearly distinguishable either : for some tested data - sets , a lognormal distribution actually achieved a higher p - values ( goodness of fit ) than power - laws , but log ratio tests suggested that other tested distributions are closer to power - laws . by using the technique in @xcite to fit the distribution of @xmath18 , a power - law with exponent @xmath20 was found with a reasonable @xmath168-value of @xmath24 .
the distribution of @xmath19 is most likely not a power - law because the @xmath168-value of @xmath169 .
table [ tab : estimatedparametersusingthemethodofclausetetal ] lists the parameters for all used datasets using clauset and virkar s fitting technique @xcite with the according log - likelihood ( log - lh ) and p values .
a positive log - likelihood ratio indicates that the power - law is favored over the lognormal .
[ c]l|cll|ll|ll|r distribution & & power - law & & lognormal & & log - lh & & # of data - points + ( binned in ) & exponent ( @xmath0 ) & p & xmin & @xmath55 & @xmath56 & p & lr & + & & & & & & & & + adding friends @xmath18 ( seconds ) & 1.53 & 0.23 & 59 & 0.54 & 2.26 & 0.0015 & 418.67 & 7,156,722 + adding follower @xmath19 ( seconds ) & 2.03 & 0.0 & 12 & 10.45 & 2.75 & 0 & -1327.5 & 6,734,405 + adding follower @xmath19 ( hours ) & 1.81 & 0.37 & 116 & 6.43 & 3.1 & 0 & -297.66 & 41,184 + sending emails ( seconds ) & 2.33 & 0.0 & 86 & 9.53 & 3.24 & 0 & -373.7441 & 68,346,901 + sending emails ( minutes ) & 2.03 & 0.21 & 1 & 6.9 & 3.32 & 0 & -391.13 & 62,031,701 + receiving emails ( seconds ) & 2.01 & 0.0 & 1 & 10.31 & 2.73 & 0 & -958.67 & 69,528,701 + receiving emails ( minutes ) & 2.21 & 0.05 & 182 & 6.23 & 3.54 & 0 & -136.31 & 15,365,001 + commenting on reddit.com ( seconds ) & 1.08 & 0.0 & 112 & 10.2 & 3.13 & 0 & -0.0024 & 180,156,539 + commenting on reddit.com ( hours ) & 1.91 & 0.01 & 64 & 0.88 & 2.41 & 0 & -20.73 & 50,044 + & & & & & & & &