individual localized spins are promising for the implementation of emerging quantum information technologies in the solid state @xcite . in this field ,
magnetic dopants in conventional semiconductors present many desirable features , such as reproducible quantum properties , stability , and potential scalability for further applications @xcite .
thanks to their point - like character , a longer spin coherence time ( compared to carriers spins ) can also be expected at low temperature , making them potentially good systems to store quantum information .
semiconductor quantum dots ( qds ) permit efficient electrical or optical manipulation of individual carrier spins @xcite .
it has been shown that the optical properties of a qd can also be used to control the spin state of individual @xcite or pairs @xcite of magnetic atoms .
the spin of a magnetic atom in a qd can be prepared by the injection of spin polarized carriers and its state can be read through the energy and polarization of the photons emitted by the qd @xcite .
the insertion of a magnetic atom in a qd where the strain or the charge states can be controlled also offers degrees of freedom to tune the properties of the localized spin such as its magnetic anisotropy responsible for the spin memory at zero magnetic field @xcite .
positively charged magnetic qds present a large magnetic anisotropy induced by the exchange interaction between the confined heavy - hole and the spin localized on a magnetic atom @xcite .
it has been shown in charge tunable mn - doped cdte / znte qds that the hybrid mn - hole spin can be prepared by optical orientation under excitation with circularly polarized light @xcite .
it was also demonstrated that a positively charged mn - doped qd forms an ensemble of optical @xmath0 systems that can be independently addressed @xcite .
this level structure suggests that the hybrid mn - hole spin states could be coherently manipulated using two resonant optical fields @xcite . for
an efficient coherent control and a practical use of this hybrid spin in a quantum device one would require however sufficiently long relaxation and coherence times for the two ground states of the @xmath0 systems .    in this article ,
resonant photoluminescence ( pl ) of the positively charged exciton ( x@xmath1 ) coupled to a single mn atom in charge tunable cdte / znte qds @xcite is exploited to analyze the spin relaxation mechanisms of the coupled hole and mn spins .
auto - correlation of the resonant pl and resonant optical pumping experiments reveal an efficient spin relaxation channel of the mn - hole system : a mn - hole flip - flop induced by the interplay of the mn - hole exchange interaction and the lattice deformation induced by acoustic phonons .
a model predicts mn - hole flip - flops in the nanosecond range , in agreement with the spin dynamics observed experimentally .
we also show that the optical @xmath0 systems are connected by inefficient forbidden spin - flips which are enhanced under a weak transverse magnetic field .
the article is organized as follows : the sample and experiments are presented in sec .
the spin structure of the positively charged mn - doped qds is shortly described in sec .
the dynamics of the resonant pl probed with auto - correlation measurements is described in sec .
iv . in sec .
v we discuss resonant optical pumping experiments allowing probing the initialization and relaxation of coupled hole and mn spins . in sec .
vi we present a model describing acoustic phonon induced mn - hole flip - flops and show that this mechanism is responsible for the strong pl observed under resonant excitation of x@xmath1-mn . finally , in sec .
vii we model the spin dynamics of positively charged mn - doped qds under resonant excitation and compare the calculated autocorrelation response and optical pumping signal with experimental results .
the sample consists of mn - doped self - assembled cdte qds grown by molecular beam epitaxy on a p - doped znte ( 001 ) substrate according to the procedure dscribed in ref .
a bias voltage is applied between a 5 nm gold schottcky gate deposited on the surface of the sample and the p - doped substrate to control the charge state of the qds . under a positive bias , as we already reported in ref .
22 , a single hole is trapped in the magnetic qds and only the emission of the positively charged exciton is observed ( figure  [ fig1](a ) ) .
individual qds containing one mn are isolated using micro - spectroscopy techniques @xcite . a high refractive index hemispherical solid immersion lens is mounted on the surface of the sample to enhance the spatial resolution and the collection efficiency of single dot emission .
the qds are excited with a tunable continuous wave ( @xmath2 ) laser tuned to an excited state of the dots or on resonance with the ground state of the positively charged exciton .
the resulting collected pl is dispersed and filtered by a one - meter double monochromator .
the temporal statistics of the emission of photons is analyzed through photon - correlation measurements using a hanbury brown and twiss ( hbt ) setup with a resolution of about 0.8 ns . under our experimental conditions , with photon counts rates of a few khz , the measured photon pair time distribution yields , after normalization , to the second order correlation function of the pl intensity , @xmath3 .    for optical pumping experiments , trains of resonant circularly polarized light
are prepared with electro - optic or acousto - optic modulators with a switching time of about 10 ns and the resonant pl is detected with a fast avalanche photodiode .
permanent magnet mounted on a translation stage are used to apply a weak magnetic field in faraday or voight configurations .    ) and positively charged ( @xmath4 ) exciton as a function of energy and bias voltage .
( b ) pl of the mn - doped qd under a positive bias voltage of v=5.5v .
inset : scheme of the energy levels of the ground ( @xmath5 ) and excited states ( @xmath4 ) in a positively charged mn - doped qd as a function of their angular momentum ( @xmath6 ) .
( c ) experimental ( left ) and calculated ( right ) color - scale plot of the linear polarization dependence of the pl of @xmath4 at b=0 t ( top ) and under a transverse magnetic field b@xmath7=0.42 t ( bottom ) , maximum transverse field used in the spin dynamics experiments .
the parameters used in the calculation are listed in table [ paraqd].,width=3 ]
a mn atom in a strained self - assembled cdte qd exhibits a fine structure dominated by a weak magnetic anisotropy with an easy axis along the qd axis @xcite . neglecting the tetrahedral crystal field of the cdte matrix @xcite , this fine structure is described by the effective spin hamiltonian    @xmath8    here , @xmath9 depicts the effect of the biaxial strain and @xmath10 describes the anisotropy of the strain in the plane of the qd .
@xmath9 varies from 0 @xmath11ev ( strain free qd @xcite ) to 12 @xmath11ev ( strained cdte layer matched on a znte substrate @xcite ) and typical values around 7 @xmath11ev are usually observed in cdte / znte qds @xcite .    when a hole is trapped in a qd containing a single mn , the spin structure is controlled by the mn - hole exchange interaction that reads    @xmath12    with i@xmath13 the exchange energy between the hole and the mn ( @xmath14 ) and @xmath15 the hole spin operator . in the presence of heavy - hole / light - hole mixing , @xmath15 , represented in the basis of the two low energy heavy - hole states ,
is related to the pauli matrices by @xmath16 and @xmath17 with @xmath18 .
@xmath19 is the coupling energy between heavy holes and light holes separated by an effective energy splitting @xmath20 .
@xmath21 is the angle relative to the [ 110 ] axis of the principal axis of the anisotropy ( shape and/or strain ) responsible for the heavy - hole / light - hole mixing @xcite . for a weak valence band mixing ,
the mn - hole energy levels are mainly controlled by @xmath22 and form a spin ladder with a quantization axis along the qds growth direction .
these states are labelled @xmath23 .
the pl of a positively charged mn - doped qd is presented in figure  [ fig1](b ) . when an exciton is injected in the qd loaded with a single hole ,
one has to consider the x@xmath1-mn complex .
the energy levels of such a complex , where the two hole spins are paired in an antiparallele configuration , is dominated by the electron - mn exchange interaction    @xmath24    with @xmath25 the electron spin and @xmath26 the exchange energy between the electron and the mn spins .
the twelve e - mn states are split into a ground state sextuplet ( total spin m=3 ) and a fivefold degenerated manifold ( total spin m=2 ) ( see inset of fig .  [ fig1](b ) ) .
these energy levels are labelled @xmath27 .
.values of the parameters used to model the positively charged mn - doped qd presented in figure [ fig1 ] .
i@xmath28 , i@xmath13 , @xmath29 , @xmath21 , @xmath30 and @xmath31 are chosen to reproduce the linear polarization intensity map of fig .
[ fig1](c ) .
the other parameters can not be extracted from the pl measurements and values for typical mn - doped qds are chosen for the calculation of the spin dynamics presented in section vii .
[ cols="<,^ , > " , ]     [ paraph ]    considering a linear dispersion of acoustic phonons @xmath32 and in spherical coordinates with @xmath33 , the explicit formula of the decay rate ( [ fermi ] ) is    @xmath34    where we used the continuum limit ( @xmath35 with @xmath36 the crystal volume ) and integrated over @xmath37 and @xmath38 .
the summation is taken over the acoustic phonon branches @xmath39 of corresponding sound velocity c@xmath40 .
the geometrical form factors for each phonon branch , @xmath41 , are given by    @xmath42     system associated with the electron - mn state @xmath43 extracted from the full level structure of a positively charged mn - doped qd ( fig .
[ fig1](b ) ) .
( b ) relaxation time @xmath44 , between the two mn - hole ground states of the @xmath0 system calculated with the material and qd parameters listed in table [ paraph ] and a temperature t=7k .
the vertical line shows the energy splitting in the studied qd of the mn - hole states involved in the @xmath0 systems considered here ( resonances ( 2 ) and ( 3 ) identified in fig .
[ fig2]).,width=316 ]    in the numerical calculation of the spin flip time @xmath44 presented in fig.[fig7 ] we use the material parameters of cdte or znte and the typical parameters for self - assembled cdte / znte qds listed in table [ paraph ] .
the calculated relaxation time strongly depend on the energy separation between the mn - hole levels @xmath45 .
this energy dependence is controlled by the size of the hole wave - function given by l@xmath46 and l@xmath47 .
the estimated flip - flop time is also strongly sensitive on the phonon induced mixing of the ground heavy - hole states with the higher energy light - hole levels @xcite . in our model
, this mixing is controlled by @xmath20 , an effective energy splitting between heavy - holes and light - holes .
this simple parameter can indeed describe more complex effects such as a coupling of the confined heavy - hole with ground state light - holes in the barriers @xcite or effective reduction of heavy - hole / light - hole splitting due to a presence of a dense manifold of heavy - hole like qd states lying between the confined heavy - hole and light - hole levels @xcite . from this modelling
we deduce that for a hole confined in small cd@xmath48zn@xmath49te alloy qds , the mn - hole flip - flop time @xmath44 can be easily bellow 2 ns for an effective heavy - hole / light - hole splitting @xmath20=15 mev and an energy separation in the mev range , typical for the mn - hole spin splitting in magnetic qds .
we will use in the following calculations @xmath44=1.5 ns for the ground states of each @xmath0 system .
using the level scheme presented in fig.[fig1](b ) for a positively charged mn - doped qd and the estimated mn - hole flip - flop rates , we can calculate the time evolution of the 24x24 density matrix @xmath50 describing the population and the coherence of the 12 electron - mn states on the excited state and the 12 mn - hole states on the ground state of a positively charged qd . in the markovian approximation , the master equation which governs the evolution of @xmath50
can be written in a general form ( lindblad form ) as :    @xmath51+l\varrho}\ ] ]    where @xmath52 is the hamiltonian of the complete system ( @xmath53-mn @xmath54 and mn - hole @xmath55 ) :    @xmath56    and    @xmath57    in ( [ lindblad ] ) , @xmath58 describes the coupling or decay channels resulting from an interaction with the environment @xcite .
the population transfers from level @xmath59 to level @xmath60 in an irreversible process associated with a coupling to a reservoir is described by a lindblad term of the form    @xmath61    where @xmath62 is the incoherent relaxation rate from level @xmath59 to level @xmath60 .
such term can describe the radiative decay of the exciton ( irreversible coupling to the photon modes ) or the relaxation of the carriers or mn spin ( irreversible coupling to the phonon modes ) .
it can also be used to describe the optical generation of an exciton in the low excitation regime where the energy shift induced by the coupling with the laser field is neglected .    a pure dephasing ( i.e. not related to an exchange of energy with a reservoir )
can also be introduced for the different spins and described by @xmath63 :    @xmath64    where @xmath65 is the pure dephasing rate of level @xmath59 .    to identify the main spin relaxation channels responsible for the observed spin fluctuations
, we first modelled the auto - correlation of the resonant pl using the full spin level structure of a p - doped magnetic qd . for a qualitative description of the observed spin dynamics , we use as an example the mn - doped qd parameters extracted from the linear polarization intensity map listed in table [ paraqd ] and reasonable order of magnitude for the spin relaxation times .    as already observed in charged mn - doped qds under pulsed resonant excitation ( ref .
23 ) , we consider that the spin dynamics in the excited state is controlled by the time evolution of @xmath54 , the generation rate of excitons @xmath66=1/@xmath67 and their radiative lifetime @xmath68=0.3 ns .
the coherence of the coupled electron - mn spins is limited by a pure dephasing term , @xmath69=0.5 ns , typical value measured in charged mn - doped qds@xcite .    for the mn - hole system in the ground state , we take into account a spin relaxation time of the mn in the exchange field of the hole , @xmath70 , describing relaxation channels involving a change of the mn spin by one unit .
this spin relaxation channel is introduced for a general description , however its characteristic time ( in the @xmath11s range ) is long compared to the time - scale of the dynamics considered in the resonant pl experiments and does not qualitatively affect the calculated time evolution .
because of the presence of valence band mixing in the qds , spin flip of the hole independently of the mn are expected to be more efficient .
a spin flip time in the 10 ns range has indeed been calculated for a hole in the exchange field of a mn @xcite .
relaxation time of the hole spin around 5 ns has also been measured at zero magnetic field in negatively charged cdte / znte qds @xcite .
we then include in the model possible spin flips of the hole by one unit with a characteristic time @xmath71=10ns .
the phonon induced mn - hole flip - flops , occurring at @xmath44 , are also introduced between the two mn - hole ground states of each @xmath0 system .    for a general qualitative description , an additional pure dephasing time @xmath72
is also included in the dynamics of the mn - hole system with a lindblad term of the form ( [ deph ] ) .
we can not extract this parameter from the experiments .
we take @xmath72= 5 ns , slightly longer than what was measured for electron - mn , as the mn - hole system is highly split and less sensitive to effective fluctuating magnetic field such as the one produced by nuclear spins for instance @xcite .     with the qd parameters listed in table [ paraqd ] and ( unless specified ) @xmath68=0.3ns , @xmath70=5 @xmath11s , @xmath73=10ns , @xmath67=0.25 ns , @xmath44=1.5 ns , @xmath72= 5 ns , @xmath69= 0.5 ns , t=7k and b@xmath74=0 .
( b ) ( c ) and ( d ) illustrate the influence on @xmath75 of @xmath44 , @xmath67 and @xmath76 respectively .
note the different vertical scale in ( b).,width=316 ]    the transition rates @xmath77 between the different mn - hole states depend on their energy separation @xmath78 .
here we use @xmath77=1/@xmath79 if @xmath80 and @xmath77=1/@xmath81 if @xmath82 @xcite .
this accounts for a thermalization among the 12 mn - hole levels with an effective spin temperature @xmath83 .
the optical excitation ( @xmath67 ) , the exciton recombination ( @xmath68 ) , the mn spin relaxation ( @xmath70 ) , the hole spin relaxation ( @xmath71 ) and the phonon induced transfer time ( @xmath44 ) produce a irreversible population transfer between level @xmath84 and @xmath85 and are described by lindblad terms ( [ inc ] ) .
detection and an excitation of @xmath86 and @xmath87 with modulated circular polarization .
the qd parameters for the calculations are those listed in table [ paraqd ] and @xmath68=0.3ns , @xmath70=5 @xmath11s , @xmath73=10ns , @xmath72= 5 ns , @xmath69= 0.5 ns , @xmath44=1.5 ns , t=10k and @xmath67=0.25ns .
( a ) influence of a variation of @xmath67 and @xmath44 .
( b ) influence of a transverse magnetic field @xmath76 .
the inset presents the transverse magnetic field dependence of the difference of population for a @xmath88 or a @xmath89 excitation.,width=316 ]    the time evolution of @xmath90 normalized by @xmath91 and with the initial condition @xmath92 ( mn - hole spin in the state @xmath93 just after the recombination of the exciton ) accounts for the auto - correlation of the emission in @xmath89 polarization of the electron - mn state @xmath94 under @xmath2 resonant excitation . @xmath90
obtained with the qd parameters listed in table [ paraqd ] is presented in fig .
[ fig8](a ) . after a fast increase ,
the calculated population presents a maximum at short delay . with reasonable spin relaxation parameters ( see caption of fig .
[ fig8 ] ) the width and the amplitude of the maximum are in good agreement with the photon bunching signals observed experimentally .
the width of the calculated bunching is controlled by all the spin - flip terms that can induce an escape out of the resonantly excited @xmath0 system . at zero transverse magnetic field ,
it is dominated by spin flips in the mn - hole system . as illustrated in fig .
[ fig8](a ) , suppressing @xmath73 gives a width of bunching only controlled by the hamiltonian evolution and the decoherence which is slightly larger than what is observed experimentally ( fig .
[ fig3 ] ) .    the dependence on the excitation intensity , @xmath67 , and transverse magnetic field , @xmath76 , are also qualitatively well reproduced by the model ( fig .
[ fig8](c ) and ( d ) respectively ) . at zero magnetic field ,
the leaks outside the excited @xmath0 systems are dominated by @xmath73 .
@xmath95 induces fluctuations in a slightly longer time scale .
the situation is different under a weak transverse magnetic field where the electron - mn states are mixed introducing new channel of escape and significantly reducing the width of the bunching .
this model is based on a large number of parameters , whose values can not all be extracted precisely from the measurements .
let us note however that suppressing the fast flip - flop process connecting the two mn - hole ground states ( @xmath96 in fig .
[ fig8](b ) ) still produces a bunching as a weak coherent transfer between the two ground states of the @xmath0 system still exist .
however , with this process only , the calculated pl intensity is about 50 time smaller and the behaviour in transverse magnetic field is not compatible with the experimental results .    with this model
, we can also calculate the population of the electron - mn states under resonant excitation with alternated circular polarization and estimate the efficiency and dynamics of the optical pumping .
figure  [ fig9 ] presents the calculated time evolution of the population of the electron - mn state @xmath97 under alternated resonant excitation of @xmath86 in @xmath88 polarization or @xmath87 in @xmath89 polarization .
this corresponds to the experimental configuration where the qd is resonantly excited with modulated circular polarization at the energy of @xmath86 and @xmath87 ( absorption ( 2 ) in fig .
[ fig2](b ) ) and the low energy resonant pl is detected in @xmath89 polarisation .
the main features of the time - resolved optical pumping experiments ( see fig .
[ fig4 ] and fig .  [ fig5 ] )
are well reproduced by the model .
the timescale of the pumping transient , in the few tens of nanosecond range , and its excitation intensity dependence are also in good agrement with the experiments ( see figure  [ fig9](a ) ) .
the influence of a transverse magnetic field , b@xmath74 , on the optical pumping transient can also be described by this model .
first , a significant reduction of the pumping time is observed for a weak magnetic field ( b@xmath74=0.2 t in fig .
[ fig9](b ) ) . as for the autocorrelation ,
this acceleration comes from the increase of the leakage out of the @xmath0 system induced by the mixing of the electron - mn states .
secondly , the transients obtained when switching the polarization from @xmath98 to @xmath99 and from @xmath99 to @xmath98 become identical for b@xmath100 t , as observed in the experiments ( fig .
[ fig5](b ) ) .    to understand this behaviour under b@xmath74 ,
let us remember that we resonantly excite @xmath86 from @xmath101 with @xmath88 light and excite @xmath87 from @xmath102 with @xmath89 photons . in both cases we detect the population of @xmath97 in @xmath89 polarization ( see the excitation / detection configuration illustrated in the inset of fig .
[ fig4](b ) ) .
if the states @xmath86 and @xmath87 are uncoupled , as it is the case at zero field , we do not detect any light during the @xmath89 excitation . with a sufficiently large mixing of @xmath86 and @xmath87 induced by the transverse magnetic field , for a @xmath89 excitation of @xmath87
, the population can be coherently transferred to @xmath86 during the charged exciton lifetime and @xmath89 light is detected after a recombination towards @xmath103 @xcite . in the optical pumping sequence ,
we can then observe , in @xmath89 polarization , a transient when the @xmath88 excitation empties the state @xmath101 but also a similar transient when the @xmath89 excitation empties the state @xmath102 .
the transverse magnetic field dependence of the difference of steady state intensity observed in @xmath98 and @xmath99 polarization ( inset of fig .
[ fig5](b ) ) is also well reproduced by the model ( inset of fig .
[ fig9](b ) ) .
this depolarization curve is controlled by the anisotropy of the electron - mn spin which is induced by @xmath30 and @xmath9 @xcite .
let us note finally that , as expected , suppressing @xmath44 from the model , a very weak average resonant pl and a fast optical pumping are obtained ( fig .
[ fig9](a ) , top curve ) .
initialized by a sequence of @xmath89/@xmath88 resonant excitation of @xmath87 and @xmath86 .
the dashed black line ( shifted for clarity ) is an exponential fit with a characteristic time @xmath104=85 ns .
( b ) corresponding calculated time evolution of the population of @xmath97 .
the parameters are those of fig .
[ fig9].,width=316 ]    including a dark time in the pumping sequence , we can also numerically evaluate the time required for the mn - hole spin to return to the ground state of the excited @xmath0 system .
the time evolution of the population of the mn - hole state @xmath105 initially prepared by a sequence of @xmath89/@xmath88 excitation resonant with @xmath94 ( and @xmath106 ) is presented in fig.[fig10 ] . when the optical excitation is switched off , after an abrupt jump due to the optical recombination of the charge exciton , the ground mn - hole state @xmath105 is repopulated in a timescale of about 100 ns , much shorter than the mn spin relaxation time used in the model ( @xmath70=5ns ) .
this relaxation is induced by the presence of valence band mixing . in the presence of valence - band mixing , @xmath107 couples two by two the different mn - hole levels .
this coupling induces a transfer of population between the different mn - hole levels .
the transfer of population becomes irreversible in the presence of dephasing and controls the observed mn - hole spin relaxation @xcite .     with @xmath108=1 ( mn - hole spin in the state @xmath109 after a @xmath89 recombination ) for a resonant @xmath88 excitation of the coupled electron - mn states @xmath110 and @xmath111 without and with a longitudinal magnetic field .
( b ) time evolution of @xmath112 under excitation with modulated circular polarization .
the parameters used in the calculations are those of fig .
[ fig9].,width=316 ]    the particular behaviour observed for a resonant excitation of the electron - mn states @xmath110 or @xmath111 ( weak photon bunching and no optical pumping at zero field , fig .
[ fig2](a ) and fig .
[ fig4](a ) respectively ) is also qualitatively explained by the model ( see figure  [ fig11 ] ) . as demonstrated in ref .
23 , the presence of a strain anisotropy term @xmath10 in the @xmath11ev range in the mn fine structure directly couples the states @xmath110 and @xmath111 which are initially degenerated .
under circularly polarized resonant excitation we either excite @xmath110 with @xmath88 photons or @xmath111 with @xmath89 photons . at zero magnetic field , the population is transferred between the two states in a time scale of a few hundreds picoseconds @xcite . under circularly polarized resonant excitation ,
the two @xmath0 systems associated with @xmath113 are simultaneously excited . for alternated circular polarization ,
a steady state is reached and no pumping transient induced by a leak outside the @xmath0 systems is expected . under a weak longitudinal magnetic field
the mn zeeman energy dominates the strain anisotropy term and the coherent transfer is blocked .
the states @xmath110 and @xmath111 are decoupled and a large amplitude of bunching and an efficient optical pumping are restored .
this behaviour observed in the experiments is qualitatively reproduced by the model .
let us finally note that in the modelling of optical pumping at zero magnetic field presented in fig .  [ fig11](b ) , fast oscillations are obtained in the first nanoseconds after the polarization switching .
these are due to the population transfer between @xmath110 and @xmath111 in the excited state ( directly coupled by e ) and between @xmath114 and @xmath115 in the ground state ( coupled by a combination of non diagonal terms of @xmath55 ) during the coherence time .
these oscillations are too fast to be observed in the experiments .
the calculated resonant pl intensity in @xmath89 polarization ( proportional to @xmath112 ) is also slightly larger for a @xmath88 excitation than for a @xmath89 excitation . the @xmath89 resonant pl probes the population of @xmath110 which is directly excited by a resonant @xmath88 laser ( see the excitation / detection configuration in the inset of fig .
[ fig4](a ) ) . on the other hand , under a @xmath89 laser ,
one excites @xmath111 and the charged exciton has a probability to recombine before being transferred to @xmath110 and detected in @xmath89 pl .
this transfer time results in a slight difference in the steady state resonant pl intensity obtained in a @xmath116 or @xmath117 configuration ( see fig .  [ fig4](a ) ) .
the model we developed here explains the main behaviour of a positively charged mn - doped qd under resonant optical excitation .
it shows in particular the importance of the effective heavy - hole / light - hole splitting @xmath20 on the mn - hole spin dynamics .
the studied cdte / znte qds have a weak valence band offset .
the resulting small value of @xmath20 is first responsible for the large influence of the qds shape or strain anisotropy on the valence band mixing .
the valence band mixing reduces the magnetic anisotropy of the mn - hole system and its spin life time .
a small @xmath20 also significantly enhances the coupling of the mn - hole spin with acoustic phonons and the phonon induced mn - hole flip - flops @xmath44 .
this fast spin dynamics limits the use of such hybrid spin system in practical quantum information devices .
the use of different qd systems with a larger valence band offset and a larger heavy - hole / light - hole splitting @xcite should significantly slow down the mn - hole spin dynamics .
using resonant pl of the positively charged exciton , we have identified an efficient spin relaxation channel for the hybrid mn - hole spin in a qd .
a modelling confirms that mn - hole flip - flops in a nanosecond timescale are induced by an interplay of mn - hole exchange interaction and the lattice deformation of acoustic phonons .
these flip - flops are responsible for the large pl intensity observed under resonant excitation of the @xmath0 systems present in positively charged mn - doped qds .
we showed that jumps out of an optically excited @xmath0 system are possible .
these leaks induce a large bunching of the resonant pl and are at the origin of the optical pumping of the mn - hole spin observed under circularly polarized resonant excitation .
escape out of the excited @xmath0 system can be enhanced by a transverse magnetic field which mixes electron - mn states in the excited state of the charged qd .
the fast mn - hole spin dynamics revealed by these experiments make difficult the practical use of this hybrid spin for quantum information devices . however , the use of a different qd systems with a better hole confinement and a larger heavy - hole / light - hole splitting would significantly reduce the influence of valence band mixing on the mn - hole spin relaxation and limit the interaction of the hybrid mn - hole spin with acoustic phonons .
this work was realized in the framework of the commissariat  lenergie atomique et aux energies alternatives ( institut nanosciences et cryognie ) / centre national de la recherche scientifique ( institut nel ) joint research team nanophysique et semi - conducteurs .
the mechanism of mn - hole flip - flop that we identified and modelled in the case of positively charged magnetic qds can also be directly observed in the resonant pl of neutral qds .
this is illustrated in fig .
[ figappendix ] which presents the pl of the exciton in a qd where all the dark exciton states ( @xmath118 ) are well separated from the bright exciton states ( @xmath119 ) and can be observed at zero magnetic field on the low energy side of the spectra .
( plain line ) and dark @xmath120 ( dotted line ) exciton energy levels calculated with @xmath26=-0.03 mev , @xmath121=0.12 mev , @xmath122 , @xmath30=10 @xmath11ev and an electron - hole exchange interaction @xmath123=-0.78 mev .
see for instance reference 32 for a description of the model.,width=316 ]    in the pl excitation experiments presented in fig .
[ figappendix ] , the detection window is set on the dark states while a circularly cross - polarized laser is scanned on the three high energy levels of the bright exciton .
three successive resonances are observed .
( 1 ) : an excitation on @xmath124 produces the dominant pl on a dark states associated with a mn spin state @xmath125 . as
a mn spin flip by two units is unlikely , we then attribute this pl to @xmath126 .
the transfer between these states involves a mn - hole flip - flop .
the second weaker contribution comes from the state @xmath127 .
it corresponds to a transfer conserving the mn spin and involving a spin flip of the electron .
( 2 ) : an excitation on @xmath128 produces the dominant pl on a dark states associated with a mn spin state @xmath129 . as discussed above , a mn spin flip by four units
is unlikely , we then attribute this pl to @xmath130 .
the transfer involve a mn - hole flip - flop .
the second contribution comes from the state @xmath131 .
it corresponds to a conservation of the mn spin and a spin flip of the electron .
( 3 ) : from an excitation on @xmath132 , the dominant pl comes from the state @xmath133 .
it corresponds to a transfer with conservation of the mn spin and a spin flip of the electron .
a pl is also observed from @xmath130 after a spin flip of the hole and from @xmath131 .
the latter transfer involves an electron - mn flip - flop .
let s note that in the state @xmath132 a mn - hole flip - flop is forbidden ( parallel hole and mn spin ) .    if we resume the observed resonances in this neutral qd , we can distinguish two types of spin transfer which can occur within the lifetime of the exciton ( typically 300ps ) : first an efficient transfer from a bright to a dark state involving a mn - hole flip - flop which produces a change of one unit of the mn spin state . secondly , a weaker transfer from a bright to a dark state involving a carrier spin - flip with conservation of the mn spin .
j. r. petta , a. c. johnson , j. m. taylor , e. a. laird , a. yacoby , m. d. lukin , c. m. marcus , m. p. hanson , a. c. gossard , science * 309 * , 2180 ( 2005 ) .
m. veldhorst , c. h. yang , j. c. c. hwang , w. huang , j. p. dehollain , j. t.muhonen , s. simmons , a .
laucht , f. e. hudson , k. m. itoh , a. morello , a. s. dzurak , nature ( london ) * 526 * , 410 ( 2015 ) .
k. saeedi , s. simmons , j. z. salvail , p. dluhy , h. riemann , n. v. abrosimov , p. becker , h .- j .
pohl , j. j. l. morton , m. l. w. thewalt , science * 342 * , 830 ( 2013 ) n. bar - gill , l.m .
pham , a. jarmola , d. budker and r.l .
walsworth , nat .
comm . * 4 * , 1743 ( 2013 ) .
p. m. koenraad , m. e. flatte , nat .
mater . * 10 * , 91 ( 2011 ) .
m. atature , j. dreiser , a. badolato , a. hogele , k. karrai , a. imamoglu , science * 312 * , 551 ( 2006 ) .
d. press , t. d. ladd , b. zhang , and y. yamamoto , nature ( london ) * 456 * , 218 ( 2008 ) .
b. d. gerardot , d. brunner , p. a. dalgarno , p. ohberg , s. seidl , m. kroner , k. karrai , n. g. stoltz , p. m. petroff , r. j. warburton , nature ( london ) * 451 * , 441 ( 2008 ) .
k. de greve , p. l. mcmahon , d press , t. d. ladd , d. bisping , c. schneider , m. kamp , l. worschech , s. hofling , a. forchel , y. yamamoto , nat .
phys . * 7 * , 872 ( 2011 ) .
l. besombes , y. leger , l. maingault , d. ferrand , h. mariette , j. cibert , phys .
lett . * 93 * , 207403 ( 2004 ) . c. le gall , a. brunetti , h. boukari , l. besombes , phys .
lett . * 107 * , 057401 ( 2011 ) .
a. kudelski , a. lemaitre , a. miard , p. voisin , t.c.m .
graham , r.j .
warburton , o. krebs , phys .
lett . * 99 * , 247209 ( 2007 ) . j. kobak , t. smolenski , m. goryca , m. papaj , k. gietka , a. bogucki , m. koperski , j .-
rousset , j. suffczynski , e. janik , m. nawrocki , a. golnik , p. kossacki , w. pacuski , nat . comm . * 5 * , 3191 ( 2014 ) .
l. besombes , c.l .
cao , s. jamet , h. boukari , j. fernandez - rossier , phys . rev .
b * 86 * , 165306 ( 2012 ) .
o. krebs , a. lemaitre , phys .
lett . * 111 * , 187401 ( 2013 ) .
l. besombes , h. boukari , c. le gall , a. brunetti , c.l .
cao , s. jamet , b. varghese , nanophotonics * 4 * , 75 ( 2015 ) .
d. e. reiter , v. m. axt , t. kuhn , phys .
b * 87 * , 115430 ( 2013 ) .
m. goryca , t. kazimierczuk , m. nawrocki , a. golnik , j. a. gaj , p. kossacki , p. wojnar , g. karczewski , phys .
* 103 * , 087401 ( 2009 ) .
j. c. oberg , m. reyes calvo , f. delgado , m. moro - lagares , d. serrate , d. jacob , j. fernandez - rossier , c. f. hirjibehedin , nat .
* 9 * , 64 ( 2014 ) .
y. leger , l. besombes , l. maingault , d. ferrand , h. mariette phys .
b * 72 * , 241309(r ) ( 2005 ) .
k. vyborny , j. e. han , r. oszwaldowski , i. zutic , and a. g.petukhov , phys .
b * 85 * , 155312 ( 2012 ) .
b. varghese , h. boukari , l. besombes , phys .
b * 90 * , 115307 ( 2014 ) .
a. lafuente - sampietro , h. boukari , l. besombes , phys .
b * 92 * , 081305(r ) ( 2015 ) . j. houel , j. h. prechtel , a. v. kuhlmann , d. brunner , c. e. kuklewicz , b. d. gerardot , n. g. stoltz , p. m. petroff , r. j. warburton , phys . rev . lett . * 112 * , 107401 ( 2014 ) .
p. wojnar , c. bougerol , e. bellet - amalric , l. besombes , h. mariette , h. boukari , journal of crystal growth * 335 * , 28 ( 2011 ) . c. le gall , r. s. kolodka , c. l. cao , h. boukari , h. mariette , j. fernandez - rossier , l. besombes , phys .
b * 81 * , 245315 ( 2010 ) .    c. le gall , l. besombes , h. boukari , r. kolodka , j. cibert , h. mariette , phys . rev .
lett . * 102 * , 127402 ( 2009 ) .
m. goryca , m. koperski , p. wojnar , t. smoleski , t. kazimierczuk , a. golnik , p. kossacki , phys .
* 113 * , 227202 ( 2014 ) .
m. qazzaz , g. yang , s.h .
xin , l. montes , h. luo , j.k .
furdyna , solid state communications * 96 * , 405 ( 1995 ) .
causa , m. tovar , s.b .
oseroff , r. calvo , w. giriat , phys . lett . *
a77 * , 473 ( 1980 ) .
a. h. trojnar , m. korkusinski , u. c. mendes , m. goryca , m. koperski , t. smolenski , p. kossacki , p. wojnar , p. hawrylak , phys .
b * 87 * , 205311 ( 2013 ) .
l. besombes , h. boukari , phys .
b * 89 * , 085315 ( 2014 ) .
s. jamet , h. boukari , l. besombes , phys .
b * 87 * , 245306 ( 2013 ) .
j. fernandez - rossier , phys .
b * 73 * , 045301 ( 2006 ) .
y. leger , l. besombes , l. maingault , h. mariette , phys .
b * 76 * , 045331 ( 2007 ) .
l. besombes , y. leger , l. maingault , d. ferrand , h. mariette , j. cibert , phys .
b * 71 * , 161307(r ) ( 2005 ) .
l. besombes , k. kheng , l. marsal , h. mariette , phys .
b * 63 * , 155307 ( 2001 ) . c. cohen - tannoudji , b. diu , f. laloe , _ mecanique quantique _ ( hermann , paris , 1973 ) .
e. tsitsishvili , r. v. baltz , h. kalt , phys .
b * 67 * , 205330 ( 2003 ) .
k. roszak , v.m .
axt , t. kuhn , p. machnikowski , phys .
b * 76 * , 195324 ( 2007 ) .
mahan , _ many - particle physics _
( plenum press , new york , 1993 ) .
woods , t.l .
reinecke , r. kotlyar , phys . rev .
b * 69 * , 125330 ( 2004 ) .
s. adachi , _ properties of group iv , iii - v and ii - vi semiconductors _ ( john wiley and sons ltd , chichester , 2005 ) .
michler , _ single quantum dots fundamentals , applications and new concepts _
( springer , berlin , 2003 ) .
luo , g. bester , a. zunger , phys .
b * 92 * , 165301 ( 2015 )
van exter , j. gudat , g. nienhuis , d. bouwmeester , phys .
a * 80 * , 023812 ( 2009 ) . c. roy and s. hughes , phys .
x * 1 * , 021009 ( 2011 ) . c. l. cao , l. besombes , j. fernandez - rossier , phys
b * 84 * , 205305 ( 2011 ) .
l. cywinski , phys .
b * 82 * , 075321 ( 2010 ) . c. le gall , a. brunetti , h. boukari , l. besombes , phys . rev .
b * 85 * , 195312 ( 2012 )
. a. o. govorov , a. v. kalameitsev , phys .
b * 71 * , 035338 ( 2005 ) .
s. moehl , f. tinjod , k. kheng , h. mariette , phys .
b * 69 * , 245318 ( 2004 ) .


the rich variety of conformations which leads to many different intrinsic properties of polymer solutions has continuously drawn considerable interest in soft matter research .
computer modeling is increasingly being used as an integral part of theoretical study , in order to both test existing theories and to trigger the development of new concepts .
furthermore , computer simulations have also become an essential tool in materials research , especially for predicting and understanding the behavior of complex systems , where a complete theory is not available .
it has been proven to be an effective and inexpensive way to study these systems . in order to observe large - scale properties , it is crucial to reduce the computational cost by coarse - graining the details of the atomic structure
. this is particularly true for polymer systems and studies of their universal static and dynamic properties @xcite . in this context , using a conventional bead - spring chain model to represent a polymer molecule in molecular dynamics ( md ) simulations is usually sufficient . in the case of dilute and semidilute polymer solutions
, a correct model also needs to take into account the effect of solvent molecules .
this effect is two - fold : on the one hand , the good solvent quality results in swelling of the random coil ; on the other , the solvent - mediated long - range dynamic correlations between different segments of the chain , known as hydrodynamic interactions ( hi ) , significantly influence the dynamical behavior @xcite . in the present methodological study , we focus on the dilute regime which is theoretically most thoroughly understood .    in order to capture hydrodynamic interactions in md simulations ,
the solvent particles need to be incorporated explicitly . typically , the number of solvent particles required for such a model is of the order of thousands even for a short chain .
although such studies are feasible @xcite , they are rather inefficient , for this reason .
therefore , a more coarse - grained description of the solvent is highly desirable .
two complementary approaches have been developed to do this .
`` mesoscopic '' methods keep the solvent degrees of freedom , but describe them in a simplified fashion .
these include dissipative particle dynamics ( dpd ) , multi  particle collision dynamics ( mpcd ) @xcite , and lattice boltzmann ( lb ) .
these approaches are typically one to two orders of magnitude faster than md @xcite .
conversely , brownian dynamics ( bd ) simulations remove the solvent degrees of freedom completely , and take their effect into account via non - trivial long - range dynamic correlations in the stochastic displacements .
this is possible due to the time scale separation between the fast solvent motion and the slow conformational polymer degrees of freedom .
since the number of degrees of freedom is reduced drastically , the method has the potential to save cpu time by additional several orders of magnitude , in particular in the dilute limit .
however , a simple implementation of the correlations @xcite leads to an algorithm which scales like @xmath3 , where @xmath4 is the number of brownian particles , and therefore becomes infeasible as soon as @xmath4 exceeds a few hundred @xcite .
it is therefore very important to treat hi by means of fixman s algorithm @xcite ( scaling roughly as @xmath5 ) , which we do in the present study .
the recently introduced method by geyer and winter @xcite would reduce the necessary cpu effort by roughly one additional order of magnitude for typical chain lengths , and exhibit a more favorable scaling ( @xmath6 ) .
however , it is based upon an inexact approximation of the hydrodynamic correlations that can not be improved systematically ( in contrast to fixman s method ) . for this reason ,
this algorithm was not implemented .
none of these approaches is sufficient to reach @xmath7 ; this latter goal is only attainable by the implementation of very recent `` superfast '' bd algorithms based upon fast fourier transforms .
these latter algorithms scale as @xmath8 , where the exponent @xmath9 depends on the details of the underlying physics , and is usually substantially smaller than unity .
these methods require the study of a confined system , and hence are not used in the present study .    while the advantages and disadvantages of the methods are well - known in general terms ( and have resulted in differing methodological preferences in different groups of researchers ) , not much is known _ quantitatively _ in terms of a clear comparison of computational efficiency .
the present paper aims at partly filling this gap .
recently , one of the present authors @xcite has proposed a new mesoscopic method for simulating polymer - solvent systems .
the solvent is represented by a fluid on a grid , simulated via the lattice boltzmann approach , while the motion of the polymer chain is governed by a continuous md model .
the two parts are coupled by a simple dissipative force .
the lattice boltzmann ( lb ) method was originally developed to simulate hydrodynamics on a grid @xcite .
it has been shown that the lb method is a fast and effective method for simulating fluid flows , which has the same speed and accuracy as other navier - stokes solvers @xcite .
ladd @xcite successfully applied the lb method to colloidal systems ( originally with a _ conservative _ coupling ) and showed that the cpu cost scales linearly with the number of particles .
moreover , he showed how fluctuations can be incorporated into the lb model , which is essential in order to investigate brownian motion @xcite .
this procedure has recently been refined and improved @xcite .
the dissipative coupling method @xcite was thoroughly tested by applying it to a single polymer chain in solution , for which the data of a previous md simulation @xcite were available , and whose parameters were used as an input for the mesoscopic model .    in this work , we study the dynamics of a single chain in a solvent to compare the predictions of the explicit solvent model via the lb method with the predictions of the implicit solvent model by bd simulations .
many recent simulation studies have investigated very similar systems , using various mesoscopic solvent models like mpcd @xcite , dpd @xcite , direct solution of the navier - stokes equation @xcite , and smoothed dpd @xcite .
these studies were mainly done in order to test and verify the validity of the chosen mesoscopic simulation approach .
meanwhile it is fair to say that a single chain in solution is a standard benchmark test bed for mesoscopic simulation methods , the necessary condition for passing being the correct reproduction of the zimm scaling laws within the limitations of finite chain length and finite solvent volume .
the present study , however , aims directly at a _ quantititave comparison _ between two different methods , and therefore it is crucial that that the underlying polymer model is exactly identical for both methods .
this is the reason why the data of refs .
@xcite , though all exhibiting essentially the same physics , are of no direct relevance for the present investigation , since all of them employ slightly different polymer models , and simulate solvents with somewhat different viscosities .
one very recent study by ladd et al .
@xcite has done a rather similar comparison between lb and bd , and arrived at similar conclusions ; the present paper should be viewed as a complementary study that puts some more emphasis on the issue of computational efficiency .    after introducing the models ,
we show how to map the input parameters of the hybrid model onto the input values of the bd model to directly compare the predicted quantities ( sec .
[ sec : models ] ) .
section  [ sec : results ] then confirms the expected physical equivalence of the two approaches in terms of comparing static and dynamic data .
furthermore , this section also presents our comparison on the numerical cost or benchmark data for both methods . finally , in sec .
[ sec : conclusions ] we summarize the results and give our conclusions .
in this work , a polymer molecule is represented by a conventional bead - spring chain model , which consists of @xmath4 beads that are connected via @xmath10 finitely extensible nonlinear elastic ( fene ) massless springs .
the lennard - jones potential , which acts between all monomers , is used to model the excluded volume ( ev ) effect .
the two potentials @xmath11 and @xmath12 are given by the expressions @xmath13 where @xmath14 is the bead - bead distance , @xmath15 is the spring constant and @xmath16 the maximum extension of the bond . @xmath17 and @xmath18 are the energy and length parameters of the lennard - jones potential , respectively .      in this method ,
the evolution of the lb variables @xmath19 is governed by the following lattice boltzmann equation @xcite : @xmath20 the variable @xmath21 is the ( partial ) fluid mass density at grid site @xmath22 at time @xmath23 , corresponding to the discrete velocity @xmath24 .
@xmath25 is the time step , and the lattice spacing is denoted by @xmath26 .
the small set of velocities @xmath24 ( @xmath27 , where the value of @xmath28 depends on the details of the model ) is chosen such that @xmath29 is a vector leading to the @xmath30th neighbor on the grid .
@xmath31 is a collision operator for dissipation due to fluid particle collisions , such that the populations always relax toward the local pseudo - equilibrium distribution @xmath32 that depends on the local hydrodynamic variables @xmath33 ( the total mass density ) and @xmath34 ( the local flow velocity ) .
the collision process is constructed in such a way that it conserves both @xmath35 and @xmath36 .
@xmath37 is the stochastic term , which is essential in order to simulate thermal fluctuations that drive brownian motion .
the local pseudo - equilibrium distribution can be represented as a second - order expansion of the maxwell - boltzmann distribution , given by @xcite @xmath38 where @xmath39 are a set of weight factors , which depend on the sublattice @xmath30 ( i.  e. the magnitude of @xmath24 ) and @xmath40 is the speed of sound . in this work ,
we have used the algorithm proposed in ref .
@xcite , however , with the modification that the original 18-velocity model ( d3q18 ) was replaced by the d3q19 19-velocity model @xcite .
the set of @xmath24 consists of the particle being at rest , the 6 nearest and 12 next - nearest neighbors on a simple cubic lattice .
the magnitudes of the velocities corresponding to these three sets of particles are @xmath41 , and @xmath42 , respectively .
the weight factors for the d3q19 model are @xmath43 , @xmath44 and @xmath45 .    the early model of ref .
@xcite only considered the thermalization of modes related to the viscous stress tensor .
it is important to note that even though this procedure is correct in the hydrodynamic limit , it provides poor thermalization on smaller length scales @xcite .
adhikari et al .
@xcite have shown that by applying thermalization to all nonconserved modes one gets a significantly improved numerical behavior at short scales ; the theoretical background is now thoroughly understood @xcite . in this work ,
we have also investigated the effects of thermalization of the kinetic modes on various dynamic properties .
the coupling to the beads is done via simple interpolation of the flow velocity from the surrounding sites , and by introducing a phenomenological stokes friction coefficient @xmath46 of the beads .
this gives rise to a friction force on the particles , plus a langevin force that balances the frictional losses .
the total momentum is conserved by subtracting the corresponding momentum transfer from the surrounding fluid .
it can be shown that this procedure satisfies the fluctuation - dissipation theorem @xcite .
for further technical details on this method and its theoretical analysis , we refer the reader to ref .
@xcite .
the configuration of a bead - spring chain is specified by the set of position vectors @xmath47 ( @xmath48 ) .
the time evolution of this configuration is governed by the it stochastic differential equation @xmath49 where summation over repeated indices is implied . here
the symbols @xmath50 , @xmath51 , @xmath52 , and @xmath53 denote the time step , boltzmann s constant , the temperature , and the diffusion tensor , respectively , where the latter decribes the hydrodynamic interactions between the beads .
the forces @xmath54 and @xmath55 are the spring and excluded - volume contributions , repectively .
@xmath56 are random variables representing a discretized wiener process , such that @xmath57 and @xmath58 , where @xmath59 is the kronecker delta and @xmath60 is the unit tensor .
finally , the tensor @xmath61 is related to the diffusion tensor such that @xmath62 @xcite .
the frictional properties of the chain and the hydrodynamic interactions between the beads are modeled via the diffusion tensor @xmath53 .
its diagonal elements contain the bead friction coefficient @xmath63 , where @xmath64 is the solvent viscosity and @xmath65 is the stokes radius of the bead .
the off - diagonal elements represent the hydrodynamic interactions via the tensor @xmath66 , for which we take the regularized rotne - prager - yamakawa ( rpy ) tensor @xcite with solvent viscosity @xmath64 and stokes radius @xmath65 . taken together ,
the diffusion tensor is given by @xmath67 further details can be found in ref .
@xcite .
the computationally most intensive part is to determine the matrix @xmath61 .
generally , cholesky decomposition of @xmath53 is used to obtain @xmath61 as an upper ( or lower ) triangular matrix , and the computational cost for this method scales as @xmath68 @xcite .
fixman made use of the fact that there are many possibilities to define @xmath61 as some square - root matrix of @xmath53 , and , noting that it is the vector @xmath69 that is required rather than the matrix @xmath61 , applied a truncated chebyshev polynomial expansion to obtain @xmath69 with a lower computational cost , scaling roughly as @xmath70 @xcite . in the present paper
, we make use of this accelerated technique as well .      for implementation on a computer , physical quantities must be represented in certain units , i.  e. in terms of suitable dimensionless ratios .
this is typically done by choosing a natural unit system where three independent elementary quantities are set to unity .    in the coupled md / lb simulation approach , this is usually done by choosing the lennard - jones parameters @xmath17 , @xmath18 , and @xmath71 ( @xmath72 , where @xmath73 is the mass of the monomer ) as the units of energy , length , and time , respectively ; this choice enables one to make direct contact with md simulations with explicit solvent @xcite .
conversely , bd simulations have traditionally @xcite used a unit system where one chooses @xmath74 as the energy unit and @xmath75 as the length unit .
this is particularly useful for the simulation of pure gaussian ( harmonic ) chains where the interaction potential has neither an energy scale nor a length scale built in .
the time unit in bd simulations is usually chosen as @xmath76 .
of course , a meaningful comparison of results requires that all data are represented in one common unit system . at this point
, one realizes that this is less straightforward than one might think at first glance . while the conversion of length and energy units is trivial , and directly facilitated by the fact that both methods use the same molecular model for the polymer chain , and perform the simulations at the same temperature , the conversion of time units is not .
this is so because of the different time scales underlying the basic updating algorithms : the md / lb method is based upon simulating the system on _ inertial _ time scales , while bd focuses directly on the larger _ diffusive _ time scales .
this is directly reflected by the occurence of the intertial parameter @xmath73 ( monomer mass ) in the md / lb time unit , which does occur in the bd model , and the diffusive parameter @xmath77 ( monomer friction constant ) in the bd time unit .
it is important to notice that @xmath77 is _ not _ an input parameter to the md / lb model .
rather , one needs to carefully distinguish between the short - time friction coefficient @xmath46 that is indeed an input parameter  it describes the stokes coupling of the monomer to the lb fluid in its immediate vicinity  and the long - time friction coefficient @xmath78 describing the particle s long - time response that is modified by solvent backflow effects .
it is this latter parameter that must be identified with the bd @xmath77 , and it is esentially an _ output _ parameter .
fortunately , the relation between @xmath46 and @xmath78 is well understood  ahlrichs and dnweg @xcite have shown that @xmath79 where @xmath80 is a numerical prefactor and @xmath26 is some measure for the range of interpolation to the surrounding lattice sites . for linear interpolation to the nearest sites
, one finds @xmath81 , if @xmath26 is the lb lattice spacing . for highly accurate results
, one should also take into account a small correction for the finite size of the simulation box @xcite ; this has however not been done in the present paper .
rather , we took the mapping determined in ref .
@xcite to calculate @xmath78 from @xmath46 and identified this with the @xmath77 parameter of the bd calculations . although the physical mapping done in this way is essentially correct , it is important to notice that this aspect introduces a certain amount of numerical inaccuracy when it comes to quantitative comparisons .
given the fact that it is intrinsically impossible to run the two simulations with the same unit system , we chose to keep the previously established systems of the two respective methods , and to map the results _ a posteriori _ by the procedure outlined above . furthermore , we chose to present all results in md / lb units .
technically , this means that for length and time unit conversions we need to set @xmath82 and @xmath83 , where the `` * '' superscript denotes bd non - dimensionalization , while `` - '' denotes a non - dimensionalization for md / lb . the corresponding factors for the conversion from bd to md / lb units are then trivially found to be @xmath84 note that @xmath85 and @xmath86 are @xmath15 and @xmath87 in non - dimensional md / lb units , respectively , while the ratio @xmath88 is just the non - dimensionalized temperature in md / lb units .
the dimensionless hydrodynamic interaction parameter @xmath89 used in the bd simulations is essentially a non - dimensionalized stokes radius .
we thus find @xmath90 or @xmath91 we therefore parameterized our simulations by first picking simulation parameters for the md / lb model ( which were then directly used for md / lb ) , then converting these to bd units using the procedure outlined above , and then running the thus - obtained equivalent bd model . for these latter simulations , a time step size @xmath92 ( in bd units )
was found to produce sufficiently accurate results .
the physical input values for the present model are chosen from the benchmark values developed in ref .
@xcite , which have been shown to reproduce the results of a typical pure md simulation @xcite . as in the comparison between lb and md simulations ,
we study a system of a single polymer chain of length @xmath93 monomers immersed in a fluid with temperature @xmath94 , density @xmath95 , and kinematic viscosity @xmath96 .
the lattice spacing @xmath97 is set to unity , which is roughly identical to the bond length ; this is necessary to resolve the hydrodynamic interactions on small length scales with sufficient accuracy .
furthermore , following ref .
@xcite , the coupling parameter @xmath98 was set to @xmath99 .
the values of the fene spring potential parameters are @xmath100 and @xmath101 .
the time step size for the polymer ( the md part of the simulations ) is set to @xmath102 .
it should be noted that such a large time step is possible since the inclusion of dissipation and noise lead to a substantial stabilization , compared to purely microcanonical md .
the value of the time step that updates the fluid should be chosen in a way such that the lb variables @xmath19 do not become negative too often . here ,
we choose @xmath103 , where such a case rarely occurred during the observation time .
it is important to mention another free input parameter which governs the time scale for the evolution of hydrodynamic interactions , known as the schmidt number @xmath104 , where @xmath105 is the diffusion constant of the single monomer .
this parameter can be set arbitrarily in the lb method by choosing @xmath106 and @xmath105 ( which can be tuned by choosing @xmath107 ) accordingly .
ideally , the value of @xmath108 should be chosen such that hydrodynamic interaction evolves much faster than the diffusion of a monomer . in our case
, we have @xmath109 , which has been shown to result in zimm - like behavior @xcite .    for the lb simulations , the polymer chain moves within a cubic box of length @xmath0 with periodic boundary conditions , while it is drifting freely in an infinite medium for the bd simulations .
in order to accurately compare various properties between the two systems , one must understand the effects of the box length @xmath0 on any observable of interest in the lb simulations .
thus it is essential to only compare quantities under identical conditions ( i.  e. independent of the box length ) .
hence various box lengths @xmath0 ranging from 10 to 35 lennard - jones units were investigated ; this allows us to extrapolate to the @xmath2 limit .
the mean square radius of gyration and the mean square end - to - end distance are given by @xmath110 with @xmath111 being the inter - particle distance .
these two quantities are both related to the number of monomers by the expression @xmath112 where @xmath106 is the flory exponent . for a self - avoiding walk ( saw ) ,
the flory exponent @xmath106 is 0.588 @xcite . in principle
, @xmath106 can be obtained from simulations , using the scaling law in eq .
[ eq : flory ] .
however , this method would require simulations for a wide range of @xmath4 values .
alternatively , one can use the static structure factor @xmath113 to obtain @xmath106 much more efficiently .    in the scaling regime @xmath114
( @xmath115 being a microscopic length of the order of the bond length ) , a power law relation between the static structure factor and the wave vector @xmath116 holds : @xmath117    figure  [ fig : skstat ] shows the static structure factor as a function of wave vector @xmath116 for the lb simulations with the presence of thermalization of all modes , and the bd simulations .
it can be clearly seen that the values of the static structure factor obtained from the lb simulations are exactly the same as those obtained from the bd simulations , indicating that they have the same static conformations . from eq .
[ eq : svsk ] , the value of @xmath106 can be extracted from the linear region of the log - log plot of @xmath118 vs @xmath116 .
as expected , the values for @xmath106 obtained via this method are the same for both the lb and the bd simulations , as reported in table  [ tab : props ] . however , they are approximately 5% higher than the asymptotically correct value , which is a consequence of the finite chain length .
the results for the mean square radius of gyration and the mean square end - to - end distance in table  [ tab : props ] further confirm this agreement with regard to static conformations between the two methods .
however , at small box length ( @xmath119 ) , the results for these static properties for the lb method deviate from their asymptotic values .
the discrepancy observed here always arises when the box length is too small compared to the chain size , where the chain is more likely to wrap over itself due to spatial restriction and hence alter its static conformations .
we also found that the two versions of lb thermalization ( `` stresses  only '' vs `` full '' thermalization ) yield identical results for the chain conformational statistics .
in general , we only quote values obtained for full thermalization , unless indicated otherwise .
the hydrodynamic radius for a single chain in an infinite medium is given by @xmath120    for a chain in a finite box , as is the case here in the lb method , it has been shown that the hydrodynamic interactions of the chain with its periodic images effectively increases @xmath121 @xcite . in order to account for this finite - size effect , a finite - size correction of order @xmath1 for most dynamic properties , resulting from the slow @xmath122 decay of hydrodynamic interactions ,
is required @xcite .
the results for the infinite - box value @xmath123 agree excellently with each other for all simulations ( see table  [ tab : props ] ) .
since the overwrapping effect is more sensitive to large inter - particle distances , it turns out that the deviation in the inverse hydrodynamic radius is too small , for the range of box lengths used , for it to be distinguishable .
as can be seen in table  [ tab : props ] , the deviation is more pronounced for the radius of gyration , and even more for the end - to - end distance .      according to dynamic scaling , the longest relaxation time @xmath124 of the chain is , by order of magnitude , identical to the time that the chain needs to move its own size , i.  e. @xmath125 , where @xmath126 is the diffusion constant of the chain s center of mass .
this leads to a dynamic scaling law @xmath127 , where @xmath128 is the dynamic scaling exponent . for a chain with hydrodynamic interactions , this relaxation time
is known as the zimm time @xmath124 . for this case , @xmath129 in the limit of long chains .
this implies that @xmath130 , which gives a dynamic exponent of @xmath131 for models with hi . for the rouse model ( i.  e. chains without hydrodynamic interactions ) , where @xmath132 , one finds a dynamic exponent of @xmath133 .
these quantities will be referred to in the discussion below .
the mean - square displacement of the chain s center of mass @xmath134 for both methods is depicted in fig .
[ fig : g3 ] . from the figure
, it can be clearly seen that @xmath135 strongly depends on the box length @xmath0 for the lb simulations .
moreover , they seem to converge to the value predicted by the bd simulations ( @xmath136 ) in the limit of large @xmath0 .
effects of thermalization of the kinetic modes in lb simulations on this property will be discussed subsequently .
the chain s center of mass diffusion constant @xmath126 can be determined by the slope of the @xmath135 vs @xmath23 curve , where the relationship @xmath137 holds . by fitting a power law to the simulation data
, we obtain the exponents and the diffusion constants shown in table  [ tab : props ] .
these exponents support the prediction of simple diffusive behavior ( @xmath138 ) .
theoretically , one would expect that two diffusive regimes exist : on the one hand , there should be a short - time diffusive regime , corresponding to time scales well below the zimm time , @xmath139 , but also well above the ballistic regime , @xmath140 ; note that @xmath141 only in the lb case , since the bd equation of motion is overdamped . on the other hand , there should be free diffusion for times @xmath142 .
both these regimes exhibit @xmath138 behavior , but with different prefactors , with a smooth crossover around the zimm time @xcite . in principle
these two different diffusion constants can be obtained via fits to the corresponding regimes . in practice , however , it turns out that the values are very close to each other , and hence the crossover is very smooth @xcite .
therefore its unambiguous identification is very difficult , i.  e. impossible within the resolution of our data .
the mean - square displacement of a single monomer @xmath30 is given by @xmath143 here , only the two innermost monomers near the center of the chain are evaluated to eliminate end effects ; the results are plotted in fig .
[ fig : g1 ] .
the values of @xmath144 behave similarly to those of @xmath135 . in the sub - diffusive time regime ,
corresponding to the short - time diffusive regime for @xmath135 , here evaluated between @xmath145 and @xmath146 , the scaling behavior @xmath147 is predicted @xcite .
the corresponding exponents obtained from a power - law fit are listed in table  [ tab : props ] and indicate a value of @xmath148 as @xmath149 . regardless of the finite - size corrections due to the box length and the effects of thermalization
, these values clearly favor the zimm model compared to the rouse model , which predicts @xmath150 .
the deviation from the asymptotic zimm value @xmath151 ( or @xmath152 ) is mainly a result of finite chain length .
figure  [ fig : g2 ] shows the mean - square displacement of a single monomer in the center of mass system ( i.  e. the two innermost monomers to eliminate end effects ) @xmath153 \\ & - & \left [ \mathbf{r}_i(t_0 )      - \mathbf{r}_\text{cm}(t_0 ) \right ] ) ^2        \rangle .
\nonumber \label{eq : g2}\end{aligned}\ ] ]    interestingly , when viewed within the center of mass system , all the results lie on top of each other , regardless of the box length @xmath0 .
this result also holds for lb simulations without full thermalization .
this shows that the global center - of - mass motion of the chain is actually the primary contribution to the deviations between lb and bd results . in the data of fig .
[ fig : g2 ] this contribution is suppressed : in terms of rouse modes , only the internal modes remain . for these modes , however , it has been shown @xcite that the hi with the periodic images is much weaker , while the leading - order @xmath122 hi cancels out .
therefore , the corresponding finite - size effect scales as @xmath154 instead of @xmath1 , and this is so small that it is invisible in fig .
[ fig : g2 ] .
theoretically , these data can also be used for estimating the zimm time as the time where the crossover to the long - time plateau occurs . however , the crossover is quite extended and smooth , making it difficult to extract .
we therefore estimated the zimm time via @xmath155 strictly speaking , this definition is only valid for a single chain in an infinite medium , where there is no finite box size effect . in the presence of finite box size
, it becomes the definition for the _ translational _ time ( @xmath156 ) rather than the zimm time : the former is subject to an @xmath1 finite - size effect , due to the strong @xmath0-dependence of @xmath157 , while the latter , being defined via the relaxation of internal modes , is only subject to an @xmath154 size effect , as discussed above .
the translational times obtained from eq .  [ eq : tauz ] ( as shown in table  [ tab : props ] ) are indeed different for different box lengths @xmath0 , as expected .
conversely , the results displayed in fig .  [ fig : g2 ] indicate that the systems with different box sizes have ( essentially ) all the same ( internal - mode ) zimm time , since their data all lie on top of each other .
next , we focus on the leading order @xmath1 finite size correction for the long - time diffusion constant of the chain s center of mass , @xmath157 . in principle
, a plot of @xmath157 vs @xmath1 should give a straight line for large @xmath0 , and an extrapolation to the limit @xmath2 should yield the same value as predicted by the bd simulations . figure  [ fig : dcmvsl ] shows the values of @xmath157 for the lb simulations with and without thermalization of the kinetic modes at various box lengths @xmath0 plotted together with the value obtained from the bd simulation at @xmath158 .
it is worth mentioning that the bd value of @xmath157 can be obtained from the mean square displacement of the chain center of mass or via fixman s expression @xcite .
the latter method has been shown to produce a much more reliable result and is easier to carry out @xcite .
the value reported here has been cross checked by both methods and the results are almost the same within error bars . for the lb simulations without thermalization of the kinetic modes , the value of @xmath157 at the asymptotic limit @xmath158 is different from that predicted by the bd simulations by about @xmath159 .
however , when all the kinetic modes in the lb simulations have been thermalized , the deviation in @xmath126 reduces to @xmath160 .
this result clearly indicates that it is very important to thermalize all the kinetic modes in order to obtain correct values for dynamic properties .
the reason for the remaining small discrepancy between lb and bd is not completely clear , since there are numerous possible sources .
firstly , it should be noted that the underlying equations of motion are quite different : lb works with inertia , while bd employs overdamped dynamics .
this results in different schmidt numbers @xmath108 and different mach numbers @xmath161 , the latter being defined as the ratio of the flow velocity to the speed of sound : both are finite in the lb method , while in the bd case they are strictly infinite ( @xmath108 ) and zero ( @xmath161 ) , respectively .
furthermore , the shape of the hi function at small interparticle distances is somewhat different for the two methods : in the bd case , we employ the rpy tensor , while the nearest - neighbor interpolation for lb results in a short - range hi that differs somewhat from the rpy tensor ( see also the discussion in ref .
finally , it should be noted that the value of the constant @xmath80 in eq .
[ eq : zeta_eff_vs_zeta_bare ] , which is crucial for the mapping between the lb friction parameter @xmath46 and the bd friction @xmath78 , is only known with some numerical inaccuracy . for highly accurate mappings ,
it is also necessary to include a finite - size correction in the definition of @xmath80 @xcite ; this was not done in the present study .    in order to examine whether the thermalization of the lb kinetic modes is also important for the internal modes of chain motion , we have performed a rouse mode analysis .
the rouse modes for a discrete chain are defined as @xcite @xmath162 \label{eq : rmode}\ ] ] for @xmath163 .
within the approximation of the zimm model , the autocorrelation function of the modes should decay exponentially @xcite @xmath164 where @xmath165 is the relaxation of the @xmath166-th mode . to validate our rouse mode analysis routine ,
we have carried out extensive simulations for a ( gaussian ) rouse chain of @xmath167 in the absence of hi and ev ; the results for @xmath165 are in excellent agreement with the analytical predictions @xcite .
figure  [ fig : rmodeall ] shows the normalized autocorrelation function for @xmath168 for the lb model with box length @xmath169 , and the bd model . for nonzero times
, there is a small deviation between lb and bd , the latter exhibiting again a slightly faster dynamics .
this deviation systematically becomes smaller upon increasing the mode index @xmath166 . since high mode index means essentially relaxation on a rather small length scale ,
it is tempting to attribute the deviation to the finite propagation of hi in the lb model , i.  e. to retardation effects , or effects of finite schmidt number , which are more important on large length scales than on small ones . nevertheless , this hypothesis is not proven .    in ref .
@xcite it was shown that the autocorrelation function is only subject to an @xmath154 finite size effect , in contrast to the usual @xmath1 behavior .
figure  [ fig : r1modevsl3 ] shows the value of the autocorrelation function of the first rouse mode @xmath170 at a fixed finite time @xmath171 , for lb simulations at various box lengths @xmath0 and bd simulations at @xmath172 . within our numerical resolution
, the data indeed confirm this @xmath154 finite size effect , both with and without thermalization of the kinetic modes .
furthermore , they demonstrate again that thermalization of all the kinetic modes in lb simulations improves the accuracy of the dynamic properties and brings them closer to the bd prediction : the deviation in the extrapolated limit @xmath2 is reduced from @xmath160 down to @xmath173 .
the reasons for the remaining discrepancies are probably of the same nature as in the case of @xmath157 .
we have also evaluated the dynamic structure factor , which is defined as @xmath174 \right ) \right > .\ ] ]    when both wave number and time are in the scaling regime ( i.  e. @xmath175 and @xmath176 ) , @xmath177 is predicted @xcite to exhibit the scaling behavior @xmath178    a plot of @xmath179 against @xmath180 should collapse to a single curve @xcite .
the results for both methods are shown in fig .
[ fig : skz ] .
the data were restricted to the scaling regime @xmath181 and @xmath182 .
these ranges were obtained from the single monomer mean - square displacement , fig .
[ fig : g2 ] , and from the static structure factor , fig .
[ fig : skstat ] , respectively . here
, we particularly focus on adjusting the exponent @xmath128 such that it would produce the best total data collapse for a chain in an infinite medium ( i.  e. in the bd model ) .
obviously , the results from the simulations show zimm - like rather than rouse - like behavior . even though we have suppressed the finite box size effect , a dynamic exponent of @xmath148 yields the best data collapse , which is somewhat smaller than the correct asymptotic one .
this result is also consistent with the value of @xmath128 obtained earlier via the exponent of @xmath144 in the sub - diffusive scaling regime ( i.  e. @xmath183 ) .
the deviation from the asymptotic value is due to the finite chain size used here , and one can expect @xmath151 only in the long chain limit @xmath184 .
more detailed comparisons of the structure factor @xmath185 are shown in fig .
[ fig : skallbd ] ( @xmath186 dependence at constant time ) and fig .
[ fig : stbyst0 ] ( time dependence for the _ normalized _ structure factor @xmath187 at constant @xmath186 ) .
figure  [ fig : skallbd ] shows the structure factor for bd simulations for a wide range of @xmath186 at three different times and the data clearly indicate that the structure factor decays rapidly with time .
the normalized structure factor @xmath187 for three different @xmath186 values are shown in fig .
[ fig : stbyst0 ] , and the data seem to indicate that the lb results approach the bd data as @xmath0 is increased , as expected .      for the ultra - dilute system considered here , the lattice boltzmann part of the hybrid lb method uses up most of the computational resources as the cpu cost for the md part for
the polymer chain is negligible .
since the dynamic properties predicted by the lb model are subject to a finite - size correction of order @xmath1 , extrapolation is required to obtain these properties in the asymptotic limit @xmath149 . to perform this extrapolation , together with checking that indeed the asymptotic @xmath1 behavior has been reached , one needs the results of at least three different box lengths . moreover
, the box length should be large enough compared to the chain size such that it does not alter the static properties .
the data displayed in table  [ tab : props ] indicate that it is safe to choose @xmath0 such that @xmath188 .
the three different box lengths @xmath0 chosen here are @xmath189 , and @xmath190 . in this work
, we set the total cpu time required for the lb simulations to be the sum of all the cpu times required to run 1000 md time steps for each of the chosen box lengths .
for bd , we take the cpu time needed to observe the system for the same time span in physical units .
each of the simulations performed for the cpu time comparison was run on an itanium 2 processor of a 1.6 ghz sgi altix server 3700 .
all the parameters used to carry out this comparison are the optimal values for both methods .
several chain sizes ranging from @xmath191 to @xmath192 have been used to obtain the cpu cost for comparison .
the results are shown in fig .
[ fig : cpu ] . for the lb method ,
it is clear that the cpu cost scales linearly with the number of particles , i.  e. the number of grid points that the solvent lives on , or @xmath193 .
since the ratio @xmath194 is kept constant , or @xmath195 , this leads to a cpu cost scaling as @xmath196 .
this is indeed found in our benchmarks , see fig .
[ fig : cpu ] .
similarly , our data also confirm the predicted @xmath70 cpu cost scaling for bd .
though the lb exponent is lower than bd , the large prefactor ensures that the total cpu cost for lb is much more expensive compared to bd for the typical chain lengths used in the literature .
it is only when the chain length is excessively large ( i.  e. @xmath4 of the order of @xmath197 or higher ) that lb will become superior to bd for a single - chain system .
the situation completely changes if one studies a semi - dilute solution instead , as has been done in ref .
for such a system , we have not done a comparison between lb and bd in terms of actual simulations ; however , by means of scaling considerations one can roughly estimate what the likely result of such a comparison would be .
a semidilute solution comprises @xmath198 chains of @xmath4 monomers each , such that the total number of monomers is @xmath199 .
therefore the bd cpu cost scales as @xmath200 , while the lb cpu cost depends on the density . within the blob picture of semidilute solutions ,
one views a chain as a sequence of `` blobs '' , each comprising @xmath201 monomers , and having size @xmath202 , which can be viewed as the typical correlation length of density fluctuations , or the typical distance from which point on chain - chain interactions become important .
since the conformation statistics within the blob is that of a self - avoiding walk , one has @xmath203 , where @xmath26 is the monomer size .
the sequence of blobs forms a random walk , hence @xmath204 .
this gives the minimum size of the simulation box , i.  e. @xmath205 , or @xmath206 .
we thus see that the cpu effort for the lb method is even slightly _ decreased _ by the factor @xmath207 compared to the single - chain case at the same @xmath4 , due to the shrinkage of the chains resulting from excluded - volume screening . in order to estimate the number of chains @xmath198
, we note that the arrangement of blobs is space - filling , i.  e. @xmath208 .
comparing this with the previous expression for @xmath193 , one finds @xmath209 .
therefore the bd effort , compared to the single - chain case , is increased by a factor of @xmath210 . taken together , this means that the ratio between lb effort and bd effort is changed by a factor of @xmath211 in favor of lb . for @xmath212 , which is needed as a minimum to resolve the gaussian statistics of the chains as a whole ,
one obtains a factor of @xmath213 , which more or less compensates the two orders of magnitude seen in fig .
[ fig : cpu ] .
taking into account that for such a system the bd simulation would have to calculate the hi with the periodic images , e.  g. , via ewald sums , which is much more complicated than the present single - chain simulation , one sees that for a semidilute solution clearly lb is more efficient , unless a `` superfast '' bd algorithm @xcite is used . for the latter case ,
the answer is not yet known .
the results of ref .
@xcite indicate that lb / md may be favorable for a rather small number of monomers ; however , this study was done ( i ) in a non - trivial geometry , which implies a more complicated bd method , and ( ii ) under complete neglect of thermal fluctuations in the lb simulations , resulting in a substantial reduction in cpu effort .
such a ( partial or complete ) neglect of thermal noise is sometimes justified in strong nonequilibrium situations such as studied in ref .
@xcite ; in that particular case , the justification was checked by additional tests @xcite .
another possible situation where lb noise is negligible is the case of strong coarse - graining , where a single lattice site can already be considered as a macroscopic thermodynamic system ( for a detailed discussion , see ref .
however , in the general case , and certainly in thermal equilibrium or weak nonequilibrium , the proper inclusion of thermal noise is necessary , as demonstrated theoretically in detail in ref .
@xcite , and also corroborated by the present numerical results .
for the general case , the estimate of the lb cpu effort given in ref .
@xcite is therefore too optimistic .
the present study has shown that brownian dynamics simulations are capable of reproducing various properties predicted by a hybrid lb / md model ( or vice versa ) .
we have demonstrated how to obtain the input values for the bd simulations from the physical input parameters of the lb model such that both models would produce the same static and dynamic properties .
for the lb model , most dynamic properties are subject to a finite - size correction of order @xmath1 . in addition to this
, it is very important to thermalize all the kinetic modes in order to obtain the correct dynamic properties .
those results that are not affected by @xmath1 finite size effects , such as the mean square displacement in the center of mass system , or the rouse mode autocorrelation function , agree very favorably with each other . for highly dilute systems where the simulation of a single chain is sufficient , bd is usually the method of choice , as it is much more efficient than the coupled lb / md approach , and finite box size effects are absent .
the situation changes however in the semidilute case , where it is easy to estimate that bd will not be able to compete , unless `` superfast '' algorithms are used .
moreover , one should take into account that the hybrid lb / md algorithm is rather easily adaptable to complicated boundary conditions , and can even be applied to flows at high reynolds numbers , where the fluid degrees of freedom become _ intrinsically _ important , and can not be handled in terms of a green s function .
this work was supported by the australian research council under the discovery projects program and the mainz max planck international research school .
computational resources were provided by the australian partnership for advanced computation ( apac ) and the victorian partnership for advanced computation ( vpac ) .
stimulating discussions with a. j. c. ladd are gratefully acknowledged .
.properties for a single chain of length @xmath214 obtained from lattice boltzmann simulations at various finite box lengths and brownian dynamics simulations in infinite medium .
@xmath215 exponent obtained by fitting a power law in the sub - diffusive scaling regime of the chain in lattice boltzmann simulations , @xmath216 $ ] .
@xmath217 exponent obtained from lattice boltzmann simulations without thermalization of all the kinetic modes . [ cols="<,^,^,^,^",options="header " , ]