the study of accretion processes is of fundamental astrophysical importance , being ubiquitous on both the largest and smallest scales , i.e. ranging from accretion by supermassive black holes in the center of galaxies at one end to star and planet formation at the other .
here we are interested in studying accretion phenomena in the presence of a large gravitational potential such as that observed in x - ray binaries ( xrbs ) and active galactic nuclei ( agn ) . in these sources ,
the fundamental timescales of interest , governing the accretion process , scale with the mass of the central object ( @xmath5 ) .
hence , studies of the galactic xrb population provide an excellent laboratory for detailed examination of the process of accretion on humanly accessible timescales .    the very - high , high - soft and low - hard states ( hereafter vhs , hss & lhs ) are the primary active accretion states observed in xrbs ( see the review by @xcite for a detailed description of accretion states in black hole binaries ) . despite being the most common mode of accretion in black hole x - ray binaries ,
the nature of the accretion flow in the low - hard state remains uncertain .
emission in the lhs is characterized by a hard power - law spectrum ( @xmath0 1.4  1.7 ) and strong x - ray variability ( 30%  40% rms ) .
correlated variability at radio / nir / optical wavelengths has also been observed from the black hole systems while in the hard state @xcite .
the hard power - law component may be the result of comptonization by a hot optically thin plasma of soft seed photons from a thermal disc or magnetic structures ( through cyclo - synchrotron processes - see , e.g. , @xcite ) .
however , the geometry of this plasma is not well understood .
a popular model for the accretion geometry in the low - hard state was given by @xcite .
in that model the standard thin accretion disc , that dominates in the spectrally soft states is radially truncated and replaced by an advection dominated accretion flow ( adaf : see @xcite ) .
the fundamental assumption of a radially recessed accretion disc may find some support in the low disc reflection fractions which are sometimes measured in the hard state ( e.g. , @xcite ) .
a number of models have been developed which do not require a recessed disc in the low - hard state .
@xcite proposed that black hole states are driven by the height and bulk velocity of magnetic flares above a disc , which remains at the innermost stable circular orbit ( isco ) .
these flares would serve to feed a mildly relativistically outflowing corona .
low disc reflection fractions do not signal a recessed disc in this model , but result from mild beaming of the hard x - ray flux away from the disc .
a similar model for the lhs , based on magnetically dominated coronae , has been proposed by @xcite . these outflowing coronae share some properties with jet based models for the hard component @xcite , in the sense that the latter also produce low reflection fractions @xcite without the need for a recessed accretion disc .
a number of recent observations call into question the assumption of a recessed disc in the lhs .
in particular gx 339 - 4 was observed during its 2004 outburst at a luminosity of @xmath6 by both _ rxte _ & _ xmm _ @xcite . the observations by _ xmm _ are critical here as they provide coverage at energies below 3 kev whereas _
rxte _ is limited to energies above this .
it was found that a cool disc blackbody ( @xmath7 0.35 kev ) consistent with an optically thick geometrically thin accretion disc extending to the innermost stable circular orbit ( isco ) existed , in contrast to theoretical expectations .
fits with reflection models revealed reflection fractions @xmath7 0.2  0.3 .
previous observations of swift j1753.5 - 0127 yielded similar results for the cool disc component , in this case a disc temperature of @xmath7 0.2 kev was required @xcite ; however , no significant disc reflection features were detected .
swift j1753.5 - 0127  was discovered by the _ swift _ burst alert telescope ( bat ) at x - ray and @xmath8-ray energies on 2005 may 30 ( @xcite ) .
subsequent observations with the x - ray telescope ( xrt ) revealed a hard power - law spectrum @xcite , while pointed _ rxte _ observations detected 0.6 hz quasi periodic oscillations ( qpo ; @xcite ) .
the system was also detected at uv @xcite and optical wavelengths @xcite .
radio observations with the _ merlin _ array also detected a variable counterpart @xcite .
observations at optical wavelengths by @xcite have also detected a significant modulation with a period of 3.2hrs , which they identify as a superhump period slightly larger than the actual orbital period .
this would make swift j1753.5 - 0127  the black hole binary with the shortest known orbital period .
in addition to the observations of @xcite above , a number of other high energy studies of this source have been published , which we summarize below ( see fig .
[ asmlc ] ) .
an analysis of _ rxte _ observations of the outburst was reported in @xcite .
the x - ray spectrum was found to be consistent with a power - law ( @xmath9 ) , while low - frequency quasi - periodic oscillations ( qpos ) were also detected ( up to 0.9 hz ) .
@xcite analysed simultaneous _ rxte _ & _ integral _ data , which were also obtained during the 2005 outburst .
the combined spectrum ( 3  400 kev ) could be fit with a model consisting of thermal comptonization modified by disc reflection ( @xmath10 ) .
qpos were also detected here , but at a lower frequency than during the outburst peak ( 0.24 hz ) .
to place meaningful constraints on the accretion flow in the lhs one requires sensitivity to both the soft x - rays ( @xmath11 2 kev ) , in order to detect the cool accretion disc , and higher energies in order to detect the most prominent disc reflection features ( 5  7 kev and 20  30 kev ) .
_ suzaku _ with its large bandpass and low background , is ideally equipped to carry out these observations . in this paper , we describe observations undertaken with the _ suzaku _ x - ray observatory in 2007 , while swift j1753.5 - 0127  was in the low - hard state . in  2 , we describe our observation and extraction of source spectra and lightcurves .
we proceed to analyze the data in  3 , where both phenomenological and more physically motivated models are considered .
the broadband spectrum ( 2  250 kev ) is consistent with a simple power - law model , although there are also contributions from the accretion disc . in
4 , these results are discussed in the context of models for the accretion flow in the low - hard state , and finally our conclusions are presented in  5 .
swift j1753.5 - 0127  was observed while still in the low - hard state by _ suzaku _ @xcite from 2007 september 19 20:36 ut until september 22 10:30 ut ( obsid:402088010 , pi : homan , see fig [ asmlc ] ) .
data were acquired over a broad spectral range ( 0.2  600 kev ) , with the x - ray imaging spectrometer ( xis : @xcite ) and the hard x - ray detector ( hxd : @xcite ) .
the source was observed at the xis nominal position for total uncorrected exposure times of @xmath7 95 ks & 82 ks respectively .
all data reduction and analysis takes place within the heasoft 6.4.1 environment , which includes ftools 6.4 , suzaku 8.0 and xspec 12.4.0x .
the latest versions of the relevant _ suzaku _ caldb files were also used .
the xis is installed at the focal plane of the four x - ray telescopes ( xrt : @xcite ) and currently consists of 3 functioning detectors xis0 , xis1 and xis3 .
xis0 & xis3 are front illuminated and provide coverage over the energy range 0.4 - 12 kev whereas xis1 is back illuminated in an effort to provide greater sensitivity at lower energies , 0.2  12 kev .
the xis has a field of view of @xmath7 18 x 18 ( 1024@xmath12 pixels ) and was operated in 5x5 , 3x3 readout mode .
in addition the data were taken in 1/4 window mode in an effort to minimize possible photon pile - up , giving a time resolution of 2s .
as the downloaded data products had been processed via the _ suzaku _ pipeline v2.1.6.15 , the data were reprocessed from the raw telemetry files as recommended in the _ suzaku _ data reduction guide ( abc guide ) .
standard screening was applied , in particular , we extracted _ asca _ event grades 0:0 , 2:4 , 6:6 with the data filtered to be taken outside of the south atlantic anomaly ( saa ) and where the earth elevation angle was greater than 5@xmath13 .
good time interval ( gti ) events were extracted using xselect , where the 3x3 and 5x5 observation mode data for each detector were extracted simultaneously .
science images , spectra and backgrounds were then extracted from these event files .
even though the observations were carried out in 1/4 window mode , we nonetheless suffered from pileup at the source position .
hence the spectra were extracted using an annular extraction region extending from 30 to 250 pixels from the source position .
this extraction region size was chosen so as to extract @xmath14 99@xmath15 of the point source flux .
as the outer radius of this annulus is larger than the window size ( @xmath7 280 , 295 , 280 pixels respectively ) , the effective extraction region is the intersection of window and the annulus .
the resulting extraction region has an area equivalent to 62% of the 250 pixels outer radius , hence , we expect to have detected a commensurate percentage of the total source flux .
background spectra were extracted from a source free region of the detector , these are automatically scaled to match the data during the spectral analysis .
response files were generated using the tasks xisrmfgen and xissimarfgen .
the background and response files were then grouped with the science spectrum for analysis in xspec .
the hxd covers the energy range from 10  600 kev , consisting of two separate detectors , ( i ) pin : silicon pin photodiodes covering the energy range 10  70 kev and , ( ii ) gso : gso / bgo phoswich scintillators covering the energy range 40 - 600 kev . due to the arrangement of the instrument , with the pin diodes residing in front of the gso scintillator in each of the 16 detector units that make up the hxd
, the raw data do not differentiate between the pin & gso , with this distinction instead being made during the extraction process .
as at the time of writing , the gso analysis procedure had not been included in the official _
data reduction pipeline , all of the hard x - ray detector data were reprocessed following the prescription in the abc guide and the 7-step guide .
gti science events were extracted using xselect with the appropriate gti files and filter options i.e. , pointing elevation @xmath14 5@xmath13 above earth and excluding data taken near the saa .
the relevant events for each detector were then extracted by requiring det_type=1:1 and det_type=0:0 for the pin and gso respectively .    the relevant background and response ( pin / gso * 20080129.rsp ) files were obtained from the _ suzaku _ website .
a gti between the data and the background file was then created using the ftool mgtime , which was then used to to extract the spectra in xselect .
standard corrections were applied to the data , i.e. deadtime for science data , pin background exposure time . as the hxd background files do not include the contribution from the cosmic x - ray background ( cxb ) , the expected cxb was simulated following the recipe detailed in the abc guide , this was then added to the background file . as a check on the accuracy of the background , a separate background estimate was made using the earth occulted data ( earth elevation = 0 ) .
this was found to be consistent with the above background .
a similar procedure is followed in the case of the gso spectral extraction with the following caveat : the background file exposure time does not need to be corrected , instead one must rebin the science spectrum to match the provided background .
the sensitivity of the gso detector is background dominated , hence the high energy detection threshold is determined not by the statistical error but by the reproducibility of the background , here we conservatively estimate the background reproducibility to be 3% ( e.g @xcite ) . in fig .
[ xispingso_back ] , we plot the extracted science & background spectra .
swift j1753.5 - 0127  is detected out to an energy of 250 kev
source and background lightcurves were extracted from the event files using xselect and the appropriate good - time - interval events , after the application of the appropriate baryocentric correction using the ftool aebarycen .
the flux from swift j1753.5 - 0127  is observed to remain constant throughout our observation .
the mean count rates for each individual detector are approximately : 16.2 , 20.8 , 17.1 , 2.6 , 2.2 counts s@xmath16(xis0 , xis1 , xis3 , pin , gso ) .
there is no evidence for any periodic modulation , the lightcurves across all detectors are characterized by rapid variability of a stochastic nature as is expected from accretion processes in the vicinity of a black hole @xcite .
this is consistent with rxte power spectra bracketing this data that revealed a power - law slope , @xmath17 .
we initially choose to fit the spectra with a number of phenomenological models in an effort to provide a broad characterization of the data .
initial spectral fits were made to the entire dataset consisting of 5 spectra in total ( xis0 , xis1 , xis3 , pin , gso ) , spanning the energy range 0.6  250 kev over which swift j1753.5 - 0127 is reliably detected .
explicitly , the spectra provide us with data in the following ranges : xis0 , xis3 0.4  10 kev ; xis1 0.2  10 kev ; pin 12  70 kev and gso 50  250 kev .
unfortunately , the low energy response of the xis detectors contains a number of uncertain residuals , hence we additionally ignored the regions below 0.6 kev and between 1.7  1.9 kev and 2.1  2.4 kev in all further modelling .
initially , the column density was fixed at a value consistent with previous detections ( @xmath18 : @xcite ) , while the normalization was allowed to vary independently for each spectrum .
the resulting normalizations were found to be consistent with those expected for the _ suzaku _ detectors .
the spectra were fit with a model consisting of an absorbed power - law ( pha*po ) , which provided an acceptable fit except at the lowest energies ( @xmath19 2 kev ) .
as the main residual was present at energies below 2 kev , this region was ignored in further fitting .
re - applying the power - law fit results in a value for the spectral index of @xmath20 ( @xmath21 ) .
it is immediately apparent that a simple phenomenological power - law model suffices to describe the spectrum from 2  250 kev .
the data was also fit with a cutoffpl model for comparison .
the fit does not indicate the presence of a cut - off in the spectrum with the cutoffpl model fit being significantly inferior to a simple power - law , i.e we find @xmath22 = 6898/6463 ( 6968/6456 ) and @xmath23 = 6962/6462 ( 7048/6455 ) for the @xmath24 models respectively .
in both cases the best fit cut - off power - law model requires the high energy cut - off to be 500 kev , i.e the intrinsic upper limit of the cutoffpl model .
ignoring the data above 150 kev does not improve the quality of the cutoffpl fit relative to the po fit .    as the gso background contains a large feature in the energy range 150  180 kev ( see fig .
[ xispingso_back ] ) , we decided to ignore the data beyond 150 kev in all further fitting .
furthermore , as there is clearly an additional component contributing at energies below 2 kev , this region is ignored while constraints are placed on the hard x - ray emission . repeating the above power - law fit
, we find @xmath25 ( @xmath26 ) unless otherwise explicitly stated ] .
inspection of the residuals reveals any contribution from disc reflection to be small , although there is evidence for curvature in the pin spectrum . to investigate the possible contribution due to disc reflection , the best fit power - law from above was convolved with the reflection model of @xcite ( reflect*po ) .
the inclination was held fixed ( cos_i _ = 0.45 ) , while the abundances of metals were frozen at the default values ( @xmath27 = 1 ) .
the best fit model reveals a slightly softer power - law ( @xmath28 ) in addition to a highly significant ( @xmath29 as determined via an f - test ) reflection fraction , @xmath30 ( @xmath31 ) . as the measured spectral parameters depend on the inclination angle , these and subsequent fits were repeated at a lower inclination angle of 30 degrees .
these fits revealed a lower reflection fraction , f @xmath7 0.14 ( see table [ specfit_params ] ) .
closer inspection of the xis spectra residuals reveals a slight excess consistent with the presence of a broad fe k line . to place a constraint on the size of any possible line
, a gaussian was added to this model .
the best energy of this line is @xmath7 6.4 @xmath32 0.1 kev , in agreement with that expected from neutral fe k@xmath33 .
this is consistent with that expected from reflection from neutral matter ( @xmath34 ) as assumed in the reflect model .
the gaussian component was then replaced with relativistic line model ( laor : @xcite ) , which more accurately represents the expected line profile in the inner disc region .
the inclination of this line is fixed at the same value as the the reflection component , the emissivity profile of the disc is fixed at @xmath35 and the outer radius of the emitting region is fixed at 400 r@xmath36 .
we find the data require the presence a broad weak iron line ( @xmath37 significant as determined via an f - test ; ew = 73@xmath38 ) .
allowing the inner disc radius to vary , we find @xmath39 .
the iron line inner radius has a strong inclination dependence with a best fit @xmath40 at 30@xmath41 .
the above inner radii are inconsistent with the isco for a schwarzschild black hole , if we extend the confidence intervals we place the following limits on the inner radius
@xmath42 ( 3@xmath43 level ) for inclinations of 63.256@xmath41 & 30@xmath41 ( cos_i _ = 0.45 , 0.8660 ) respectively .
including the lower energy flux in the fit once again , results in a chi - squared value of @xmath44 , we plot this in fig .
[ soft_excess ] .
a soft excess is present consistent with previous observations , i.e. @xcite .
a simple blackbody accretion disc component ( diskbb : @xcite ) was added to this model to account for the excess soft x - ray flux , while the power - law index , @xmath45 , and the reflection fraction , f , were frozen at their previous best fit values .
the disc component is strongly required by the data ; the best fit is achieved for a disc temperature of kt = 0.20 kev ( @xmath46 ) .
we measure the associated 0.6  10 kev ( 2  150 kev ) unabsorbed flux to be @xmath47 ( @xmath48 ) .
the column density is crucial here : although it has little impact on the spectrum above energies @xmath7 2 kev , it will have a significant impact on the shape of the spectrum and the measured flux below this value , e.g. fig .
[ soft_excess ] . in order to test the effect of different values for @xmath49
the fitting was repeated at a number of different values for the interstellar column density ranging from @xmath50 @xcite to @xmath51 @xcite .
we find the best fit as measured by _
suzaku _ was found to be @xmath52 . in table
[ specfit_params ] , we list the parameters for the best fit model and also those for the model corresponding to the value of @xmath49 measured previously by _ xmm_.    the normalization of the diskbb model is proportional to the inner disc radius , norm = ( @xmath53 [ km])/d [ 10kpc])@xmath12cos@xmath54 . in fig .
[ inner_radius ] , we plot the inner disc radius corresponding to the best fit diskbb component for various values of the column density , where we have corrected for spectral hardening and the inner disc radius ( see  [ thin_disc_isco ] ) .
we find that for reasonable values of the column density and inclination the inner disc radius may reside close to the isco , in agreement with recent work , which has provided evidence that the cool disc may reside at or near the isco , even at the low luminosities typically observed in the lhs @xcite .
we also modelled the soft excess using the ezdiskbb and diskpn models in xspec . in both cases ,
the observed blackbody component is found to be consistent with that obtained using diskbb .
we also experimented with using the kdblur kernel to relativistically blur the above continuum model ; however , this did not result in an improved fit .
we also carried out fits to the individual telemetry segments of the xis , pin & gso spectra , to check for possible variability .
the fits to each of the 3 segments were found to be consistent with each other , hence there is no evidence for spectral variability in agreement with the constant flux observed from the lightcurves .    [ cols="<,^,^,^,^,^ " , ]
we present _ suzaku _ broadband spectra ( 0.6  250 kev ) of the black hole candidate swift j1753.5 - 0127  while in the lhs .
the observed spectrum is measured to be consistent with an unbroken power - law , @xmath55 . during our observation
the flux from the source was observed to be constant , with a 0.6 - 150 kev unabsorbed flux of @xmath56 ( @xmath57 ) .
this is consistent with _
integral _ observations of gro j1655 - 40 in the lhs where unbroken power - law emission ( @xmath58 ) extending out to @xmath7 500 kev was detected at a luminosity of @xmath7 0.015 l@xmath59 ( @xcite ; although see @xcite ) .
accurate determination of the interstellar column density is crucial due to the preferential absorption of soft x - rays , which in our case are consistent with being emitted by the cool accretion disc .
measurement of @xmath49 is best achieved at soft x - ray energies . here _
xmm _ with its well studied low energy calibration should provide the most reliable determination of the interstellar column density @xmath60 @xcite .
in contrast the best fit value as measured with _
suzaku _ is @xmath61 .
however , we do note that the measured @xmath49 depends on the continuum model assumed , e.g. for the comptonizing corona model we find a best fit value @xmath62 .    in order to investigate the instrumental dependence of our measured values of the column density
, we also fit the best fit models above to a number of _ swift _ xrt spectra taken before and after our _ suzaku _ observation , with the closest occurring @xmath7 60 days beforehand ( obsid : 00030090050 ) .
the results from these fits are consistent with the results of our fits to the _ suzaku _ spectrum , with the caveat that the smaller energy coverage ( 0.6  10 kev versus 0.6  150 kev ) and exposure time ( @xmath7 2 ks ) necessarily result in larger confidence regions from fits to the _ swift _ spectra .
these measurements are in agreement with numerous independent measures of the interstellar column , sourced both from direct measurements and galactic surveys .
x - ray observations from _ swift _ measured a column density of @xmath63 @xcite , while optical measurements also require @xmath64 @xcite .
additionally an estimate of the column density may be obtained from various radio ( n@xmath65 = @xmath66 @xcite , @xmath67 @xcite ) , and far - ir ( @xmath68 @xcite ) surveys .
the radius of the cool disc component depends on the column density and inclination as illustrated in fig .
[ inner_radius ] . from above
, we see that for the directly measured values of @xmath49 towards swift j1753.5 - 0127 , the observed soft excess may originate from an accretion disc whose inner radius is consistent with the isco .
the significant soft excess observed in swift j1753.5 - 0127  ( fig .
[ soft_excess ] ) is consistent with an origin in a standard multi - colour blackbody accretion disc ( @xcite and see table [ specfit_params ] ) .
the disc normalization may be used to estimate the inner radius of the accretion disc when knowledge of the distance and inclination are available , as norm @xmath69 .
this estimate is subject to a number of corrections , ( i ) spectral hardening must be accounted for , typically this requires a multiplicative correction factor @xmath7 1.7 @xcite and , ( ii ) one must correct for the inner radius , where @xmath70 ( @xcite ) .
there are also additional errors , i.e. the zero torque inner boundary condition @xcite and radiative transfer effects @xcite , which could contribute but are difficult to quantify .    in fig .
[ inner_radius ] , we plot the inner disc radius ( corrected for ( i ) & ( ii ) above ) measured from our best fit models in  3.2 .
the best fit radius is consistent with the isco , for certain values of @xmath49 & inclination , although there are large uncertainties .
previous observations of swift j1753.5 - 0127  at a 0.5  10 kev flux of 3.9@xmath71 detected a cool accretion disc with an inner radius consistent with the isco .
the measurements presented herein , at an unabsorbed 0.6  10 kev flux of 6.8@xmath71 , are in agreement with the previous analysis of @xcite and furthermore atest to the stability of the cool disc component given that the observations were separated by almost 18 months ( see fig .
[ asmlc ] ) .
@xcite observed gx 339 - 4 at a luminosity of @xmath72 with _ xmm _ & _ rxte_. a significant excess of soft x - ray flux was detected in the _ xmm _ data alone , due to the absence of soft x - ray coverage for _
rxte_. modelling revealed this excess to be consistent with a cool accretion disc ( @xmath73 ) extending to the isco .
subsequent observations utilizing both _ rxte _ & _ swift _ at luminosities of 0.023 and 0.008 @xmath74 also detected this excess soft x - ray component @xcite .
together with the new observations presented in this paper , these results appear to confirm the presence of the inner radius of the accretion disc at the isco ( or at least at a radius far lower than expected ) for luminosities far below that at which the state transition from hss to lhs state occurs , in contrast to theoretical expectations ( e.g. @xcite ) .
@xcite have re - analysed the _ xmm _
data on swift j1753.5 - 0127  & gx 339 - 4 focusing on the timing characteristics and variability inherent to the spectra . using a new analysis technique
, they find significant evidence for the presence of the soft spectral component in agreement with @xcite .
they interpret the soft x - ray flux as originating in variability intrinsic to the accretion disc .
in particular they find that the intrinsic variability of the accretion disc is likely to be responsible for the low frequency lorentzian feature present in the power spectral density function of sources in the low - hard state .
this places a limit on the disc truncation radius of @xmath11 20 @xmath75 .
the _ xmm _
spectrum of swift j1753.5 - 0127  has also been re - analysed by @xcite .
they fit the spectrum with a relativistically blurred reflection model and find an inner disc radius of @xmath76 .
this in turn provides an measure of the spin of the putative black hole of @xmath77 .
again , inconsistent with the idea of a large truncation radius for the accretion disc in the low - hard state .
the magnitude of the reflection features is observed to be quite low in our data .
although the presence of a broad relativistic line is required by the data ( @xmath78 significant , ew @xmath79 ) , it is intrinsically weak ( norm @xmath80 ) .
the detection of a relativistic iron line is consistent with the measurements of @xcite who also measured a redshifted line at an energy of 6.2 kev during pointed _ rxte _ observations ( see fig .
[ asmlc ] ) . the inner radius inferred from the laor line fits is larger than the isco at @xmath81 10  20 gravitational radii although much less than might be expected from the standard disc picture for disc truncation , where the inner radius is expected to be an order of magnitudes larger @xcite .
the curvature present in the pin spectrum is also consistent with a low reflection fraction , f @xmath7 0.12 - 0.21 ( inclination @xmath7 30@xmath41  63.256@xmath41 ) .
we note that the reflection fraction inferred from the self consistent cdid model is at the lower end of the values inferred from the phenomenological reflect*(diskbb+po)+laor model .
it is known that the reflection features of black hole binaries are generally weaker in the low - hard state than in higher luminosity states ( i.e. @xcite ) .
the relative weakness of these features in swift j1753.5 - 0127  is primarily due to the low luminosity nature of the source ( @xmath82 ) . in comparison ,
the x - ray luminosity in gx 339 - 4 ( @xcite ) was much higher when a strong relativistically smeared iron line ( @xmath83 significant , ew @xmath7 350 ev ) was detected ( @xmath84 ) .
the value we measure for the reflection fraction in swift j1753.5 - 0127  is comparable to previous lhs observations where reflection fractions of @xmath7 0.2  0.3 were measured in cyg x-1 & gx 339 - 4 @xcite .
in contrast , observations at higher luminosities require much higher reflection fractions .
x - ray spectra of the black hole binaries xte j1650 - 500 and gx 339 - 4 in the vhs state revealed reflection fractions of approximately unity @xcite .
the low reflection fractions inferred in the low - hard state may be interpreted as evidence for a truncated disc , e.g. @xcite . alternatively the disc may remain close to the isco with the low reflection fraction resulting from beaming of the hard x - ray flux away from the disc by a mildly relativistic outflowing corona , e.g. @xcite .
a number of alternative explanations have been put forward to account for the presence of the soft excess observed in swift j1753.5 - 0127  & gx 339 - 4 .
@xcite reanalyzed the data from @xcite and find that while it is consistent with containing a soft disc component , a number of alternative continuum prescriptions that allow the disc to be truncated are also valid ; however , see @xcite who outline a number of issues with this result .
@xcite propose a scenario in which the disc is truncated at larger radii ; however , the inner edge of the disc is irradiated by the high energy comptonized photons in the corona .
this irradiation will cause the truncated inner disc edge to be hotter and hence appear to lie at smaller radii than is actually the case . in the model of @xcite , the inner edge of the cool truncated disc
is seen to overlap with the inner adaf region . here in the overlap region
the cool disc is heated to temperatures mimicking those of a cool disc at the isco .
it is also possible that what we observe is not the accretion disc but instead a ring of material that has condensed from the corona .
such an idea has been explored in detail by @xcite and @xcite .
they find that it is possible for enough material to condense from the corona to produce an inner ring of material extending from near the isco to a few tens of schwarzschild radii , while the accretion disc is truncated at larger radii .
in particular , detailed modelling was able to reproduce the soft disc components observed by @xcite in gx 339 - 4 ( see @xcite ) .
two additional black hole binaries have been observed by _
suzaku _ while in the low - hard state .
@xcite presented observations of gro j1655 - 40 , while @xcite observed cyg x-1 .
the observations took place while the systems were in the lhs at luminosities of approximately 0.007 l@xmath59 and 0.02 l@xmath59 respectively .
we also modelled the spectra of swift j1753.5 - 0127  with a combination of 2 comptonizing components , to aid comparison of this data to the _ suzaku _ observations of the black hole binaries cyg x-1 @xcite & gro j1655 - 40 @xcite . for both of these systems ,
the observed spectrum was interpreted as being due to a 2 component comptonizing corona , i.e. a population of electrons with the same temperature ( @xmath85 ) but differing optical depths ( @xmath86 ) , modified by disc reflection , diskbb+compps+compps .
the parameters of our best fit model are listed in table [ specfit_params1 ] , while the fit itself is displayed in fig .
[ spec_ii_images_2 ] . in table
[ obs_comparison ] , we display the parameters for our best fit comptonization model along with those for the best fit models for cyg x-1 and gro j1655 - 40 .    the blackbody disc component is found to be cool ( kt @xmath7 0.1 kev ) and the inner radius of the accretion disc is consistent with the outer radius of the corona within the errors ( @xmath87 ) .
the best fit radius for the iron line is also consistent with this value .
the high value for the hydrogen column density ( @xmath88 ) returned by this model is inconsistent with the available data and casts doubt the comptonization scenario as a viable model for the observed _ suzaku _ spectrum of swift j1753.5 - 0127 .
however , the presence of circumstellar matter or a local absorber could account for the excess absorption above that detected at radio wavelengths .
we find an electron temperature cooler than that measured in both cyg x-1 and gro j1655 - 40 when the spectra are interpreted in terms of this comptonizing corona model .
assuming the disc truncation model for the low - hard state to be correct , we would expect the electron temperature to decrease with increasing luminosity , e.g. @xcite .
the value of @xmath89 we measure would imply that swift j1753.5 - 0127  was at a luminosity in excess of that at which cyg x-1 was observed ( @xmath90 ) .
this is in contrast to various luminosity estimates .
indeed , we measure a luminosity of @xmath91 ; for swift j1753.5 - 0127  to have a luminosity greater than that of cyg x-1 at the time of the _ suzaku _ observations would imply a distance greater than 8.5 kpc and or a black hole mass less than 10@xmath92 .
@xcite have argued that the inner radius of the accretion disc does not extend to the isco in the lhs state based on the _ suzaku _ observations of cyg x-1 and gro j1655 - 40 ( see  4.3 ) .
they base their argument around 2 main points ; firstly the low reflection fractions measured in the lhs state show that either the accretion disc does not intrude too deeply into the corona or some form of outflow must be formed .
they consider the outflow case to be unlikely as the reflection fractions measured in both gro j1655 - 40 ( i @xmath93 ) and cyg x-1
( i @xmath94 ) are similar , whereas if the reflection was due to an outflow it would be expected to have a strong inclination dependence , which is not observed . here
, we note that the errors on the reflection fractions measured in the case of cyg x-1 ( f = 0.4@xmath95 ) , call into question the validity of this inclination based argument , given the quality of the currently available data .
secondly it is argued that the inner radius determined from fits to the iron k line in cyg x-1 ( @xmath96 ) is inconsistent with the isco .
this point is much more ambiguous as their reported value of the inner radius is @xmath97 .
this value is consistent with the isco for a schwarzschild black hole . as such
, it is clear from the currently available data that there is evidence for a cool disc component , with an inner radius that in some cases is consistent with the isco , whether or not the cool disc component actually extends to the isco is presently not clear .
finally , it is important to note the high mass x - ray binary nature of cyg x-1 . here
the accretion process , namely accretion via a stellar wind , is significantly different from that in gro j1655 - 40 ( roche lobe overflow ) and as such may preclude detailed comparison with the low mass x - ray binaries like swift j1753.5 - 0127  and gro j1655 - 40 .    lcccccccccc system & l@xmath98 & n@xmath65 & kt@xmath99 & r@xmath100 & r@xmath101 & r@xmath102 & kt@xmath103 & @xmath104 & @xmath105 & f + & [ l@xmath59 ] & [ @xmath106 ] & [ kev ] & [ km ] & [ km ] & [ km ] & [ kev ] & & & [ @xmath107 ] + swift j1753.5 - 0127 & 0.016 & 0.31 & 0.1 & @xmath7 830 & @xmath7 205 & @xmath7 360 & 53 & 0.34 & 2.57 & 0 + gro j1655 - 40 & 0.007 & 0.74 & 0.2 & @xmath7 330 & @xmath7 26 & @xmath7 65 & 135 & 0.25 & 1.2 & 0.5 + cyg x-1 & 0.02 & 0.66 & 0.3 & @xmath7 103 & @xmath7 75 & @xmath7 200 & 100 & 0.4 & 1.5 & 0.4 +    , to facilitate comparison with the seds displayed in @xcite .
the optical / nir & x - ray data are quasi - simultaneous , see  [ the_jet ] for details . ]
our knowledge of role of the relativistic radio jets in x - ray binaries has advanced rapidly in the past 10 years .
observations of numerous x - ray binaries have revealed the presence of a compact quasi - steady jet in the low - hard state , whereas in the high - soft state the jet is typically observed to disappear @xcite .
these also point towards a significant contribution from the jet towards the overall energy budget of the binary system , with the jet luminosity being possibly equal to that observed in x - rays , e.g. @xcite .
in addition , a significant correlation has been found between the x - ray and radio luminosities of accreting sources across a range of luminosity ( in eddington units ) , suggestive of a disc - jet coupling @xcite .
@xcite have shown how one could , in principle , distinguish between different jet models based on the magnitude of the reflection features present in the spectrum .
in particular , reflection fractions greater than 0.2 could not be reproduced by their jet models , reflection fractions between 0.10.2 could be due to synchrotron self - comptonization in the base of a jet , while a reflection fraction below 10% could be interpreted as evidence that the spectrum was dominated by synchrotron emission from the jet .
if the base of the jet is confined within a few gravitational radii of the black hole , relativistic light bending will become important .
this will act in the opposite direction to the beaming of the outflow leading to a slightly larger reflection fraction @xcite .
these models are also consistent with the detection of an accretion disc near the isco , as we have likely detected with _ suzaku_.    a break or cut - off in the high energy spectrum of black holes and neutron stars can be interpreted in terms of the electron temperature in a thermal comptonizing corona .
the high energy spectrum of swift j1753.5 - 0127  extends to at least 200  kev without a strong break or cut - off in our observations .
moreover , single thermal comptonization models fail to achieve acceptable fits . in the context of @xcite , the reflection fraction of @xmath7 0.12  0.21 measured in swift
j1753.5 - 0127   is consistent with a disc that is illuminated by synchrotron self - comptonization produced in the base of a jet .
a corona that is independent of a jet but which produces a non - thermal comptonization spectrum might be able to describe the data , but the base of a jet provides a natural context for the production of hard x - ray emission like that we have observed .    in section [ spec_i ] , we showed how the data may be fit with a power - law modified by reflection , this model is supported by more detailed modelling in  [ spec_ii ] . in both of these cases an incident power - law is required , the jet is a natural source of such a power - law . in fig .
[ sed_jet ] , we plot the absorption corrected spectral energy distribution for swift j1753.5 - 0127 .
the optical data is from @xcite and was taken during at approximately the same time as our _ suzaku _ spectra .
in addition , we include v - band optical and j , h & k - band nir data points obtained in april and july 2009 at the magellan and mdm observatories .
the v - band magnitude was found to be consistent with the previous measurements of @xcite , as such the nir data should be representative of the infrared portion of the sed at the time of our _ suzaku _ observations .
the radio data points represent the radio detection shortly after the 2005 outburst and likely represent upper limits to the radio flux at the time of our observations @xcite , these are also consistent with the non - detection at 5 ghz & 8 ghz by @xcite .
the solid line extending from the x - ray to the optical represents a power - law as might be expected from optically thin radio emission , i.e. @xmath108 .
this is similar to the slope of the optically thin jet in gx 339 - 4 @xcite . extrapolating this power - law to optical / nir wavelengths
shows that a spectral break is required at optical / uv frequencies .
analysis of optical data of swift j1753.5 - 0127  in its current state by @xcite favour a case where any jet emission is insignificant at optical wavelengths consistent with our extrapolated power - law , this is in agreement with an independent observation / analysis by @xcite .
@xcite has shown , based on analysis of optical data from xte j1118 + 480 , gx 339 - 4 & swift j1753.5 - 0127 , that the variability at optical wavelengths is inconsistent with being produced by reprocessing of x - rays in the accretion disc and instead favour a model where the emission in the low - hard state is dominated by the jet , e.g. @xcite .
we have presented _ suzaku _ broadband x - ray observations of the candidate black hole x - ray binary swift j1753.5 - 0127 . the broadband spectrum ( 2  250 kev )
is observed to be consistent with a simple power - law model ( @xmath109 ) .
confirming previous observations , we detect the presence of an excess at soft x - ray energies .
in addition , a weak relativistic iron line and curvature consistent with a reflection hump at 20  30 kev are detected .
these observations point towards the persistence of the accretion disc at a much lower radius than previously appreciated in the low - hard state .
estimates of the disc inner radius with both a simple diskbb+po model and a detailed reflection model reveal values consistent with the isco ( @xmath110 ) for certain values of both the column density and inclination .
in contrast modelling the spectra with a comptonization model , while revealing a truncated inner disc ( @xmath111 ) , implies a value for the hydrogen column density in disagreement with all previous estimates .
this research has made use of data obtained from the _ suzaku _ satellite , a collaborative mission between the space agencies of japan ( jaxa ) and the usa ( nasa ) .
this research has made use of data obtained from the high energy astrophysics science archive research center ( heasarc ) , provided by nasa s goddard space flight center .
this research made use of the simbad database , operated at cds , strasbourg , france and nasa s astrophysics data system .
caballero garcia m.d .
, miller j.m . , kuulkers e. , diaz trigo m.d . ,
homan j. , lewin w.h.g .
, kretschmar p. , domingo a. , mas - hesse j.m .
, wijnands r. , fabian a.c . , fender r.p . , van der klis m. , 2007 , apj , 669 , 534


collective optical ( superradiant ) effects appear in ensembles in which the distances between single emitters are much smaller than the radiation wavelength with which they interact ( dicke limit ) .
@xcite the spatial dependence of the electromagnetic field within such ensembles is negligible and thus all the systems effectively interact with common photon reservoir .
this leads to formation of rapidly decaying ( superradiant ) and optically inactive ( subradiant ) states @xcite and consequently to the appearance of a vacuum - induced coherence effect which results in occupation trapping .
@xcite    although effects resulting from collective coupling of atoms to radiative environment have been known for nearly sixty years , are very well described @xcite and have been extensively investigated experimentally , @xcite they still attract much scientific attention .
this is caused by the increasing variety of physical systems in which these phenomena may be observed , such as quantum dots ( qds ) , @xcite bose
einstein condensate , @xcite superconducting qubits , @xcite ion coulomb crystals @xcite or depolaritons.@xcite    the investigation of superradiant effects have been initiated and driven to a large extend by the promise which the short - living states show for optimization of lasers . @xcite recently the scientific interest has been focused on the concept of `` superradiant laser '' which allows to increase spectral purity of emitted light .
@xcite the technological realization of such a device presented in ref . allows to improve accuracy of atomic clocks @xcite and thus measurements of gravity @xcite and fundamental constants .
@xcite the experimental investigations of the collective effects were mostly restricted to the analysis of the superradiant states which appear spontaneously in the cascade emission and manifest themselves as a maximum in the intensity or photon emission rate .
@xcite    although the observation of the subradiance phenomena was also reported in the atomic ensembles , @xcite as the opposite of the superradiant states , the preparation of subradiant states is much more difficult , and therefore the possibilities they give were much less investigated experimentally .
recently , the preparation of an optically inactive state was reported in a system of superconducting qubits @xcite and in a diatomic molecule in an optical lattice .
@xcite the advantage of the subradiant states stems from their decoupling from the photon environment because of which they do not undergo radiative decoherence and thus may form decoherence - free subspaces .
@xcite this makes them useful for quantum information processing especially for noiseless encoding of quantum information .
@xcite the stable states also allow to construct a scalable quantum processor , @xcite quantum memories , @xcite nonlinear sign - shift gate @xcite and storing time - bin qubits @xcite for quantum cryptography .
interesting group of dicke states is formed of single - excitation combinations belonging to the class of `` w '' states which have been widely considered for quantum information processing @xcite or optimization of the quantum clock synchronization .
@xcite    systems composed of two and more coupled qds attract much scientific focus due to the richness of their properties which pave the way to new technological applications . already pairs of quantum dots allow for long - time storage of quantum information,@xcite conditional optical control of carrier states , @xcite implementation of a two - qubit quantum gate , @xcite optical writing of information on the spin state of the dopant mn atom @xcite and construct quantum nanoantennas due to the collective phenomena .
@xcite systems of three qds enable to realize two different kinds of entanglement , @xcite teleportation via superradiance @xcite , cnot gates@xcite and the control of spin blockade .
@xcite moreover , in these systems the collective transport effects ( electronic dicke or kondo - dicke effect ) may be realized @xcite and lead to the enhancement of thermoelectric efficiency .
@xcite arrays of qds allow to reduce the effect of pure dephasing on quantum information encoded in excitonic states .
@xcite    in this paper we analyze the collective optical effects in ensembles of three and four qds . compared to double qds in which only one optically inactive state may be realized , @xcite ensembles of three and more two - level systems allow to realize many stable states which occur at different exciton occupations of the single emitters . although the superradiant effects are well described in systems of identical atoms the description of such phenomena in qd ensembles requires taking into account characteristic for those system properties which distinguish them from natural ones
. therefore we include in our model the fundamental energy mismatches , different dipole moments of single dots and coupling which induces excitation transfer between single emitters , but conserves the total population of the ensemble . in our previous works concerning double qds
it has been shown that the collective optical effects are extremely sensitive to inhomogeneity of the fundamental transition energy which leads to the decay of the exciton occupation for the energy splitting much below the present technological feasibility .
this decoherence effect may be strongly reduced by sufficiently strong coupling between the dots @xcite and fully overcome in double qds with different decay rates .
@xcite in this paper we extend the results of two qds and specify conditions which allow to take advantage of the superradiant phenomena in fully inhomogeneous qd systems .
we analyze the dynamics of one electron - hole pairs and biexcitons .
we show how to adjust the system parameters in such a way that an arbitrary dark single - exciton combination may be blocked in a multiple qd , we also specify the conditions which allow to trap two excitons and prepare a system in a biexcitonic state which allows to recombine only one electron - hole pair . due to existence of two or more subradiant single - exciton eigenstates and coupling between the dots the occupation of individual dots oscillates
while the population of the whole ensemble remains stable , we show that the oscillation amplitudes may be strongly reduced if the system is initially prepared in a biexciton state .
the paper is organized as follows . in sec .
[ sec : system ] we describe the system under study , define its model and describe a method used to study the system evolution , in sec .
[ sec : results ] we present and discuss our results , we conclude the paper in sec . [ sec : conclusions ] .
the investigated system consists of @xmath0 @xmath1 quantum dots ( qds ) in which only the ground - level exciton states with fixed spin polarizations are taken into consideration . due to the strong coulomb coupling and absence of the external electric fields
we may restrict the discussion to ` spatially direct ' excitonic states , i.  e. , states with electron - hole pairs residing in the same qd which in these conditions have much lower energy than the ` dissociated ' states referring to excitons formed of carriers residing in two different dots .
@xcite these assumptions allow us to treat every qd as a two - level system which may be either empty or contain an exciton and thus describe the set of @xmath0 qds as an @xmath2-level system , with @xmath3 denoting the ground ( or `` vacuum '' ) state in which all @xmath0 dots are empty , single - exciton states @xmath4 , corresponding to one exciton localized in the @xmath5th qd , biexcitonic states @xmath6 referring to electron - hole pairs residing in the @xmath5th and @xmath7th qds , states @xmath8 with three qds , @xmath5th , @xmath7th , and @xmath9th ( @xmath10 ) , occupied with an exciton _
etc_.    present manufacturing technology does not allow to produce on demand systems of qds with identical fundamental transition energies , therefore we assume different electron - hole binding energies of each dot defined as @xmath11 where @xmath12 is the average transition energy and @xmath13 is the energy mismatch of the @xmath5th qd .
we impose @xmath14 and @xmath15 , such that @xmath16 is the mean square variation of the transition energies .
we analyze the system in a rotating frame defined by the evolution operator @xmath17,\end{aligned}\ ] ] where @xmath18 are the annihilation ( creation ) operators for the exciton in the @xmath5th qd , respectively and @xmath19 is the standard free photon hamiltonian with operators @xmath20 and @xmath21 creating and annihilating radiation modes with wave vector @xmath22 and polarization @xmath23 , while @xmath24 is the corresponding frequency .    in this frame
the hamiltonian of the system is @xmath25 the first term describes electron - hole pairs residing in the qd system .
we assume that the ground state @xmath3 corresponds to the zero energy level , so the excitonic hamiltonian is @xmath26 where @xmath27 are biexcitonic shifts due to the interaction of static dipole moments of @xmath5th and @xmath7th qd , @xmath28 is a deviation of energy caused by interaction of the dipole moments of three dots , _
etc_.    in ensembles of qds one may distinguish two types of interaction between the emitters : dipole ( frster ) coupling which in the leading order decays as @xmath29 with the qd separation @xcite and short - range coupling resulting form a combination of tunneling ( wave function overlap ) and coulomb correlations .
both types of interaction induce excitation transfer between the emitters but conserve the total exciton occupation of the ensemble .
the short - range couplings allow to rebuilt enhanced emission in the energetically inhomogeneous ensembles while in the case of dipole interaction a similar effect is achieved if the coupling is enhanced artificially by a factor of @xmath30 .
@xcite to overcome the destructive effect of the transition energy mismatch the coupling between the dots must be of the order of the energy splitting , i.  e. , @xmath31 mev for technologically feasible systems .
@xcite for a planar qd arrangement the distance between emitters is about @xmath32 nm ( the average value for the sample studied in ref . ) ,
and for such a distance the frster coupling drops to about @xmath31 @xmath33ev , which is not sufficient to stabilize collective effects .
the frster coupling reaches a fraction of mev for qds separated by only a few nanometers,@xcite which can be achieved in vertically stacked qd systems , and where indeed it can stabilize the collective effects.@xcite in the present paper we describe a planar system and therefore we prefer to consider a short - range coupling , for which we choose an exponential model , @xmath34 $ ] , where the subscripts @xmath5 and @xmath7 refer to @xmath5th and @xmath7th qd , respectively , @xmath35 is a constant amplitude , @xmath36 is the distance between the qds and @xmath37 the spatial range of the interaction .
remarkably , this simple model reproduced well recent experimental results.@xcite    the eigenstates of the hamiltonian ( [ hamiltonians ] ) do not mix quantities associated with different exciton numbers , i.  e. the eigenstates of the system are superpositions of the basis states restricted to a particular exciton number .
the second term of the hamiltonian accounts for coupling between the qd system and quantum electromagnetic field @xmath38 where @xmath39 is a coupling constant for the @xmath7th qd . here
@xmath40 is the inter - band dipole moment for the @xmath7th qd , @xmath41 is the unit polarization vector of the photon mode with polarization @xmath23 , @xmath42 is the vacuum dielectric constant , @xmath43 is the relative dielectric constant and @xmath44 is the normalization volume .
we restrict our investigations to wide - gap semiconductors with electron - hole binding energies of the order of 1 ev which allows us to describe the photon modes within the zero - temperature approximation at any reasonable temperature .    to describe the evolution of the carrier subsystem we use an equation of motion for the reduced density operator in the markov approximation . in the rotating frame it takes a form @xmath45+\mathcal{l}_{\mathrm{rad}}[\rho],\end{aligned}\ ] ] where @xmath46 is a reduced density matrix of the exciton subsystem and @xmath47 is a lindblad dissipator
@xmath48   = \sum_{i , j=1}^{n}\gamma_{ij}\left[\sigma^{(i)}_{-}\rho\sigma^{(j)}_{+ } -\frac{1}{2}\left\{\sigma^{(j)}_{+}\sigma^{(i)}_{-},\rho\right\}\right ] , \label{superlindblad}\end{aligned}\ ] ] where @xmath49 with @xmath50 being the speed of light .
since for @xmath51 the eq .
( [ gamma ] ) describes the spontaneous decay rates of single qds , @xcite the mixed ( off - diagonal ) decay rates ( [ gamma ] ) may be expressed in terms of single qd quantities , @xmath52 where @xmath53 and @xmath54 , here @xmath55 is an irrelevant phase and @xmath56 is a small angle between the dipole moments which depends of light - hole admixture .
@xcite    in the numerical simulations we assume constant energy mismatches with the parameters @xmath57 , @xmath58 , @xmath59 , @xmath60 and @xmath61 mev ( except for the fig .
[ fig:3qd_angle ] ) . for the coupling amplitudes
we take @xmath62 mev and @xmath63 nm . in sec .
[ idealsingle ] , [ coupledsingle ] and [ sec : biexciton ] we assume parallel dipole moments , since the effect of the light - hole admixture in the absence of external electric fields is negligible . @xcite
below we present an analysis of the collective effects which occur in multiple qds .
we define the dicke states for an arbitrary number of emitters and perform numerical simulations for ensembles of three and four qds . in fig .
[ fig:34qd ] we illustrate the numbering of the qds and their spatial arrangement . in sec .
[ idealsingle ] and [ coupledsingle ] we focus on single - exciton states . in sec .
[ idealsingle ] we show the evolution of uncoupled systems with identical fundamental transitions and parallel dipole moments of different amplitudes .
then , in sec
. [ coupledsingle ] we analyze the same effects in a system composed of three coupled energetically inhomogeneous dots .
the dynamics of biexciton states is presented in sec .
[ sec : biexciton ] . in sec .
[ pinjunction ] we show possibilities of controlling the exciton occupation given by a p - i - n junction .
the collective effects were described for the first time in ensembles of uncoupled identical atoms where all the emitters have the same transition energies and dipole moments .
@xcite such optical effects , resulting from the interaction in the dicke limit , are also present in uncoupled systems , even with different dipole moments , if all the emitters have identical transition energies . for the purpose of this paper
we define such systems as ideal qds .
the coupling of excitons to their radiative environment described in the dicke limit by the hamiltonian ( [ hamiltoniansrad ] ) and the fermi s golden rule , according to which the probability of releasing a photon through a transition from the initial to the final states is @xmath64 , allow to define rapidly decaying ( superradiant ) and optically inactive ( subradiant ) states , both also known as dicke states . by definition
the superradiant initial state ( @xmath65 ) corresponds to the maximum transition probability , whereas the subradiant states ( @xmath66 ) refer to the vanishing probability . due to the decoupling from the photon reservoir these are dark , optically inactive states . in the weak excitation limit , i.  e. , for a single excitation in the system from which the sample may decay only to the ground state ( @xmath3 )
, the proportionality @xmath67 allows to write the short - living states in the form @xmath68 and express the stable superpositions as @xmath69 where the coefficients @xmath70 satisfy @xmath71 .
irrespective of the number of emitters there is only one superradiant state in the system of a particular number of qds , whereas the only structure which realizes just one dark state is a double quantum dot .
the systems of three and more qds allow to realize an arbitrary number of dark states of the form ( @xmath72 ) since there are may combinations of the parameters @xmath70 for which the transition matrix element @xmath73 .     and
coherences [ ( b ) and ( d ) ] of a system of three qds prepared initially in a bright state @xmath74 for @xmath75 ns@xmath76 , @xmath77 ns@xmath76 and @xmath78 ns@xmath76 . figs .
( a ) and ( b ) show the results for uncoupled qds ( @xmath79 ) with identical transition energies ( @xmath80 ) , while figs .
( c ) and ( d ) refer to energetically inhomogeneous ( @xmath81 ) and coupled ( @xmath82 ) system arranged in a lateral array of an equidistant triangle shape with the side length @xmath83 nm .
the inset to fig .
( a ) shows the initial evolution of the exciton occupation of the qd number 2 . ]
the consequence of the co - existence of rapidly decaying and stable states is the effect of spontaneous trapping of excitation .
@xcite an arbitrary single - exciton state @xmath84 , where @xmath85 is the localization probability of the exciton on the @xmath5th qd and @xmath86 , may be expressed as a combination of the superradiant state ( [ bright ] ) and a dark state of type ( [ dark ] ) , @xmath87 where the dark state is @xmath88 }    { \sqrt{\left(\sum\limits^{n}_{i=1}\gamma_{ii}\right )    \left[\sum\limits^{n}_{i=1}\gamma_{ii}-\left(\sum\limits^{n}_{i=1}c_{i}\sqrt{\gamma_{ii}}\right)^2\right]}}.\end{aligned}\ ] ] the derivation of the eqs .
( [ single ] ) and ( [ dark2 ] ) is done in the the appendix .
the collective coupling to the radiative surrounding induces emission only from the superradiant state and thus the fraction of excitation initially spanned in the dark state , @xmath89 , remains unaffected . since the only single - exciton state which decays totally is the superradiant state , we define a state @xmath90 as being bright if its superradiant contribution does not vanish , i.  e. , if a system prepared in that state partially recombines and only a part of the initial exciton occupation remains trapped .    in fig .
[ fig:3qd_comparison ] we show the dynamics of a single - excitation induced by a common photon reservoir in an ideal system of three uncoupled ( @xmath91 ) qds with equal electron - hole binding energies ( @xmath92 ) and parallel dipole moments ( @xmath93 ) of different magnitudes . in fig .
[ fig:3qd_comparison](a ) we show the evolution of the exciton occupations of a system prepared initially in a bright state @xmath74 and in fig .
[ fig:3qd_comparison](b ) we show the corresponding coherences .
as expected , the coupling to the photon reservoir spans the excitation into the sub- and superradiant states according to eq .
( [ single ] ) and induces emission only from the short - living state @xmath94 [ fig . [ fig:3qd_comparison](a ) ] .
the excitation dynamics takes place until occupations of all dots stabilize at certain levels corresponding to the dark state defined in the eq .
( [ dark2 ] ) which confirms that the state ( [ dark2 ] ) is indeed unaffected by the photon reservoir and , after the decay of the superradiant state , neither the total exciton occupation nor occupations of single dots ( @xmath95 ) change due to radiative environment .
the emission from the above state induces decay of the total exciton occupation and excitation transfer which results in the redistribution of the occupations of single dots .
since all of the localized single - exciton states @xmath4 contribute to the superradiant state ( [ bright ] ) , the collective coupling spans the initial excitation in all of the dots even if some of them were initially empty .
therefore the population of initially unoccupied systems builds up spontaneously [ @xmath96 , green - dashed line in fig .
[ fig:3qd_comparison](a ) ] .
if the initial occupation of one of the dots is relatively small while the spontaneous decay rate from that system is sufficiently strong then the exciton occupation of that dot may vanish at same point and then restore due to the excitation transfer [ magenta - dotted line in fig . [ fig:3qd_comparison](a ) and the inset to fig . [ fig:3qd_comparison](a ) ] . during the emission process
also the evolution of the off - diagonal density matrix elements is observed , the coherences related to the initially populated dots decay , while those corresponding to initially empty systems build up spontaneously due to the increasing occupations of those dots .
when exciton dynamics in the system reaches population distribution corresponding to the optically inactive state also the off - diagonal density matrix elements stabilize at a certain non - zero level [ fig .
[ fig:3qd_comparison](b ) ] , defined by the dark contribution to the initial state .
technologically feasible qds forming multiple structures differ both in fundamental transition energies ( @xmath81 ) and dipole moments ( @xmath97 ) , and are coupled with each other ( @xmath82 ) .
as shown in the previous section and , for a double qd in ref .
, the superradiant character of the evolution of one exciton is present in the ideal systems ( @xmath80 ) , but with parallel dipole moments ( @xmath93 ) .
the collective evolution is very sensitive to the energy mismatches and is destroyed in ensembles with energy splittings of the order of the transition line width . @xcite in such systems the localized eigenstates corresponding to different energies can not form delocalized superpositions which would also be the system eigenstates .
this destructive effect may be overcome by coupling between the dots ( @xmath98 ) which delocalizes the system eigenstates and different dipole moments allowing the superradiant state to be a non - symmetric superposition of the localized states @xmath4 [ eq.([bright ] ) ] .
the single - exciton eigenstates of the system depend on the energy mismatches and coupling between the dots , while the dike states are defined by the interplay of decay rates [ eqs .
( [ bright ] ) and ( [ dark ] ) ] .
if the single - exciton decay rates [ eq . ( [ gamma ] ) for @xmath99 are adjusted in such a way that the superradiant state ( [ bright ] ) corresponds to one of the system eigenstates , then the inhomogeneous ensemble of qds interacts with its radiative environment in the `` collective regime '' i.  e. allows many effects typically present only in systems with identical electron - hole binding energies .
the amplitudes @xmath100 of a single - exciton state orthogonal to the superradiant state ( [ bright ] ) must satisfy the equation @xmath101 which implies the condition @xmath102 defining a subradiant state .
therefore , if one of the eigenstates has a superradiant character then the other eigenstates of a system are optically inactive and thus defying the `` collective regime '' requires only specifying the superradiant eigenstate .    in figs . [ fig:3qd_comparison](c ) and [ fig:3qd_comparison](d ) we show the evolution of a realistic group of three qds placed in the corners of an equilateral triangle , we assume non - equal fundamental transition energies ( @xmath61 mev ) , non - vanishing coupling between the systems ( @xmath82 ) and parallel dipole moments ( @xmath93 ) . we compare the results obtained for an inhomogeneous system coupled to the photon reservoir in the `` collective regime '' to the ideal case presented in figs .
[ fig:3qd_comparison](a ) and [ fig:3qd_comparison](b ) , where the decay rates of individual dots take the same values as in figs .
[ fig:3qd_comparison](c ) and [ fig:3qd_comparison](d ) . as can be seen in fig .
[ fig:3qd_comparison](c ) , coupled ensembles with energy mismatches of the order of mev allow to trap the same fraction of excitation as ideal dots with the same decay rates [ red - solid lines in figs .
[ fig:3qd_comparison](a ) and [ fig:3qd_comparison](c ) ] . here , as in the ideal case
, any single - exciton eigenstate may be decomposed into sub- and superradiant component according to eq .
( [ single ] ) and also in this case the superradiant state is the only state which decays totally . for an arbitrary set of decay rates , which do not correspond to the superradiant eigenstate ,
the exciton occupation of a system prepared initially in a state of the form ( [ dark ] ) is quenched and the decay of the state ( [ bright ] ) is slowed down compared to the `` collective regime '' .
although sub- and superradiant states may exist in the appropriately designed realistic systems , the dynamics of individual qd occupations differs considerably from the discussed in the previous section ideal case , when the dots interact only through the common radiative reservoir .
the situation changes when the dots communicate with each other via short - range coupling which induces excitation transfer between dots and thus oscillations in the evolution of single qd populations . in a double qd system
the oscillation amplitudes decrease and the occupations stabilize at levels corresponding to the only one dark state of the system .
@xcite in the multiple qds composed of three or more emitters oscillations of the individual dot populations never vanish [ figs .
[ fig:3qd_comparison](c ) , [ fig:3qd_comparison](d ) and fig .
[ fig:3qd_oscill ] ] .
the excitation is trapped in the system because of the existence of dark states , which in such systems may be realized by many different amplitude combinations and thus also the final population number may be realized in various ways .    )
shown by the blue line in fig .
[ fig:3qd_comparison](c ) .
the time evolution is shown for three different time intervals at picosecond resolution scale . ]    in fig .
[ fig:3qd_oscill ] we zoom the exciton occupation @xmath103 shown in fig . [ fig:3qd_comparison](c ) ( blue line ) . as can be seen in fig .
[ fig:3qd_oscill](a ) the initial oscillation pattern is relatively complicated which is caused by the existence of three energy gaps defined by the differences between the superradiant eigenstate and the subradiant ones , @xmath104 , and by the energy splitting between the two subradiant eigenstates , @xmath105 .
the period of the envelope oscillations is a multiplication of the corresponding three periods , i.  e. @xmath106 and @xmath107 .
the period @xmath108 itself defines the fine oscillations of the occupation .
due to the emission process , the superradiant contribution attenuates and interference pattern simplifies [ fig .
[ fig:3qd_oscill](b ) ] .
finally , when the short - living state decays and the total exciton occupation becomes trapped , the evolution of the single dot occupation shows single - mode pattern which repeats with the time @xmath108 and with amplitude depending of the initial occupation of individual dots [ fig . [ fig:3qd_oscill](c ) ] .
defined in fig .
[ fig:34qd ] ( a ) , on the energy mismatch ( b ) and the corresponding steady state occupations for the initial state @xmath109 ( c ) and ( d ) , respectively . for all figs .
we assume @xmath110 nm and @xmath111 nm and a constant decay rate of the qd number 2 ( @xmath112 ns@xmath76 ) . ]    the system of three and more qds allows many different planar arrangements of the emitters . since the coupling amplitudes ( @xmath98 ) depend on the distances between emitters , the eigenstates of the system , and thus the decay rates for which the ensemble interacts collectively with its radiative environment , also depend on the arrangement of the dots . using the geometry design defined in fig .
[ fig:34qd](a ) we calculate the dependence of the decay rates @xmath113 and @xmath114 on the spatial arrangements of the system .
we assume constant distances @xmath115 and @xmath116 and thus constant values of the coupling amplitudes @xmath117 and @xmath118 and change the angle @xmath119 from @xmath120 degrees to linear design , i.  e. , we increase the distance between dots @xmath31 and @xmath121 . as seen in fig .
[ fig:3qd_angle](a ) the values of the decay rates necessary to form the `` collective regime '' slightly decrease with increasing angle @xmath119 .
the two similarly decreasing decay rates ( while the third one is constant ) , according to eq .
( [ single ] ) lead to increasing steady state ( final ) exciton occupation which is shown in fig . [ fig:3qd_angle](c ) .
if the coupling between two out of three qds is much stronger than coupling of that dots with the third one ( e.g. @xmath122 , @xmath123 ) then one of the system eigenstates has large contribution from the localized state associated with the weakly coupled dot ( number @xmath124 ) and to achieve `` collective regime '' the decay rate @xmath125 must be much smaller from @xmath113 and @xmath114 . in the limiting case of vanishing couplings @xmath126 and
@xmath118 the localized state @xmath127 becomes the system eigenstate and the corresponding decay rate ( @xmath125 ) vanishes .
consequently the pair of coupled dots acts as a double qd while the third dot does not contribute to the evolution .    in fig .
[ fig:3qd_angle](b ) we show the dependence of the decay rates forming the `` collective regime '' on the energy mismatch .
both calculated decay rates decrease with increasing energy separation , but one changes slowly while the second decays fast , this leads to decreasing ( eq . [ single ] ) steady state occupation shown in fig .
[ fig:3qd_angle](d ) . for appropriately arranged ensemble of dots the decay rates may decrease in such a way that for smaller energy mismatches the final exciton occupation decreases with increasing energy mismatch , but after exceeding a critical point it increases .
the similar effects in the dynamics of single excitons are observed for electron - hole pairs confined in ensembles of four and more dots .
( a ) superradiant initial state .
[ ( b ) and ( c ) ] initial biexciton state which allows for a recombination of only one exciton .
( @xmath128 , @xmath129 , @xmath130 , @xmath131 , @xmath132 , and @xmath133 ) .
the values of the decay rates of individual dot are the same for all figs . and take the values : @xmath134 ns@xmath76 , @xmath135 ns@xmath76 , @xmath136 ns@xmath76 , and @xmath137 ns@xmath76 .
the biexction shifts are @xmath138 mev , @xmath139 mev , @xmath140 mev , @xmath141 mev , and @xmath142 mev with @xmath143 being an arbitrary parameter .
the line types used in panel ( b ) are also valid for panel ( a ) . ]    in multiple qds built out of three and four units more than one exciton may be delocalized and thus these systems allow for more complex collective effect .
below we focus on biexciton states which in general are described by a vector @xmath144 , where @xmath145 . the realistic ensembles of qds ( @xmath81 , @xmath146 , @xmath82 and @xmath147 ) permit superradiance phenomena in the two - exciton subspace if the biexcitonic , as well as single - exciton , eigenstates correspond to the dicke states . as in the single - exciton case , the biexcitonic state is considered to be superradiant if the exciton occupation of a system prepared in this state decays totally . although eigenstates of the hamiltonian ( [ hamiltonians ] ) do not mix localized basis states associated with different exciton numbers , the biexcitonic superradiant state may be formed provided one of the single - exciton eigenstates has a superradiant character .
thus the biexcitonic superradiant states occur only if collective effects are present in the single - exciton subspace . due to the equal number of single - exciton eigenstates and qds forming the ensemble the `` collective regime '' in the single - exciton domain
may be achieved by adjusting only the decay rates ( dipole moments ) .
as will be explained in detail below , the same rates @xmath148 ( [ gamma ] ) define the biexcitonic dicke states .
thus , in order to achieve the collective effects in the two - exciton subspace one has to also appropriately adjust the spatial arrangement of the dots or energies .
the superradiant two exciton superpositions may be spanned as well in ensembles of four emitters as in triple qds and take a form @xmath149 similarly to the single - exciton case , the biexcitonic superradiant states are defined by the maximum value of the transition probability ( @xmath150 ) , but this time from the initial biexciton state to the final single - exciton superradiant one ( [ bright ] ) .
the form of the condition is governed by the coupling to the radiative environment [ hamiltonian ( [ hamiltoniansrad ] ) ] which induces decay of only one exciton at a time .
thus total quenching of two excitons must occur through formation of single - exciton superradiant states . as can be seen in fig .
[ fig : biex_fig](a ) , due to the decay of the biexciton superradiant state a part of initial excitation is initially transferred to the single - exciton state which reaches a maximum population and then is totally quenched , together with the biexciton state .
the condition for the existence of the superradiant state ( [ brightbi ] ) implies that the transition to subradiant states ( [ dark ] ) vanishes .
two excitons can be blocked in a system if a transition from the biexciton state to any single - exciton state is forbidden , i.  e. the transition matrix element @xmath151 vanishes .
due to infinite number of possible single - exciton states this condition reduces to the requirement @xmath152 . to simplify the description we define non - normalized amplitudes @xmath153 in such a way that the amplitudes of biexciton states @xmath154 . for a triple
qd the condition leads to a system of three equations of the form : @xmath155 where every subscript @xmath5 , @xmath7 , and @xmath9 takes the values @xmath31 , @xmath124 , and @xmath121 , respectively .
the system is satisfied only in the case of vanishing amplitudes @xmath156 which means that it is impossible to block two excitons in a triple qd .
the coefficients of a stable biexcitonic state spanned in a system of four qds must satisfy the system of four equations of the form : @xmath157 where every subscript @xmath5 , @xmath7 , @xmath9 , and @xmath158 takes the values @xmath31 , @xmath124 , @xmath121 , and @xmath159 , respectively .
the above equations lead to the condition for the @xmath153 numbers in a form @xmath160 although in a system of four qds many dark states may be spanned , in realistic systems the parameters of biexciton states may be adjusted in such a way that only one , particular ( for specified values @xmath161 and @xmath162 ) state is blocked . thus , as in the ideal systems and similarly to the double qds ,
@xcite the contribution of any pair of dots to the total biexcitonic population is constant in time ( @xmath163 ) .    in realistic qds
one may realize either superradiant or subradiant biexciton eigenstate but never both simultaneously as in the ideal system . in both cases the basis is supplemented by a third kind of state which allow for a recombination of one electron - hole pair and trapping of the second exciton .
the evolution of a four qd system prepared in this state is shown in figs .
[ fig : biex_fig](b ) and [ fig : biex_fig](c ) . while the population of the biexciton state decreases , the occupation of single - exciton state increases until it stabilize at the level corresponding to total trapping of one electron - hole [ green - dashed line [ fig : biex_fig](b ) ] .
the third basis state must be orthogonal to the superradiant state ( [ brightbi ] ) , which equals to the requirement of vanishing transition matrix element between the biexciton state and the single - exciton superradiant state ( @xmath164 ) and means that the state @xmath165 has a subradiant character .
the orthogonality to the subradiant biexcitonic state ( [ subbiex ] ) excludes contributions from the two - exciton dark states and thus the population trapping occurs only due to formation of a single - exciton subradiant state ( [ dark ] ) .
since the transition to the state ( [ bright ] ) and thus the decay of single - exciton states is forbidden , one exciton is blocked in the ensemble .
because the biexcitonic subradiant states can not be formed in triple qds , the requirement defining biexcitonic states which allow to block one exciton reduces to the orthogonality to the superradiant state , and the states take a form @xmath166 where @xmath167 . whereas , due to the orthogonality to the superradiant state the six amplitudes of the biexcitonic states spanned in ensembles of four qds take a form : @xmath168 where the coefficients must satisfy equation @xmath169 and , using further the orthogonality to the subradiant state relations , @xmath170    due to the existence of single and biexciton optically inactive states in ensembles of four qds the initial biexcitonic states allow for many combinations of final occupation
. an arbitrary fraction of exciton occupation ( @xmath171 ) may be trapped by an appropriate combination of blocked single - excitons and biexcitons due to contribution to the initial state from dark states and basis states which allow to recombine only one electron - hole pair .
irrespective of the initial number of excitons , ensembles of three qds allow to span only single - exciton subradiant states and thus block only single - excitons .
it is important to emphasis that if the system was prepared in a bright biexcitonic state which leads to trapping of single - exciton occupation [ figs .
[ fig : biex_fig](b ) and [ fig : biex_fig](c ) ] , then the pronounced oscillations due to coupling between the dots appear only in occupation of localized biexciton states , while the amplitudes of oscillations in the populations of single dots are negligible [ [ fig : biex_fig](c ) ] compared to single - exciton initial state [ fig .
[ fig:3qd_comparison](d ) ] .
this means that the biexciton initial state allows to achieve well defined stable single - exciton subradiant states , which is important for the application of quantum computation .
the presence of optical collective effects in realistic qd systems requires high accuracy of a system parameters which may be controlled on the manufacturing stage or by external fields by , e.g. implementing the dots into the intrinsic region of a p - i - n junction .
this structure provides a possibility of a separated injection of electrons and holes into qds from both sides of a sample and control of exciton dynamics and qd parameters through application of contacts on n- and p- type regions .
it has been shown that controlled with a bias voltage carrier tunneling into a single qd in a p - i - n structure incorporated into a microcavity leads to regulated emission of single photons and pairs of photons.@xcite the ideas were followed by a technological realization of an electrically driven single - photon emitter with a layer of self - organized inas qds.@xcite the gate voltages constructed over dots allow to control energies of the dots and dipole moments , but the magnitudes of decay rates of single qds ( [ gamma ] ) depend on the average energy of the ensemble ( @xmath12 ) and thus operations on one dot change the decay rates of all qds in the system .
nm and prepared initially in a single - exciton subradiant state . in regions i and iii the decay rates are the same as in fig .
[ fig : biex_fig ] i.  e.  selected such that the systems interact collectively with the radiative surrounding .
the initial state is @xmath172 .
the red solid line corresponds to changes of the decay rates of individual qds , whereas the blue dotted line to the same changes plus changes of the relative angular orientation of the dipoles . in regions ii the decay rates and angles
have been modified such that the ensemble does not interact collectively with its photon environment ( see text ) . ]    in order to simulate this possibility we have calculated the time evolution of the system of four qds and we have changed the decay rates at some selected time points . the results are shown in in fig . [
fig : pin ] where the system is initially in the `` collective regime '' and has been prepared in a single - exciton subradiant state .
we assume that the control electric fields are weak enough to exclude `` dissociated '' exciton states .
as expected , initially the excitation is blocked in the system .    at time 2.5
ns a change of parameters occurs which destabilizes the system and induces quenching of the excitation .
the effect may be produced by a variation of the internal electric field in the p - i - n junction , which is simulated here by a change of the decay rates of each single dot , and shown by the red solid lines in fig .
[ fig : pin ] .
the decay rates have been changed as follows : @xmath173 in this case the quenching is relatively weak because of a small change in the ratio of the decay rates but the effect is visible . at time 7.5 ns the parameters of each dot are changed back and the system population is again stable .
the decay may be enhanced by changing the orientation of the dipoles as shown by the blue dotted lines in fig .
[ fig : pin ] .
all dipoles are parallel in the regions i and iii , i.  e. all angles @xmath93 .
but now , in addition to the previous variations of @xmath148 the angles are also modified in the regions ii : @xmath174 radians .
the misalignment of the dipole moments creates thus a stronger decay .
the initial conditions may be restored at any time which may result in the trapping of a desired fraction of the initial occupation .
we have studied the optical collective effects due to interaction of multiple quantum dots built of three and four emitters with the radiative surrounding .
ensembles of three and more emitters allow to span many subradiant states which facilitate preparation of the system in an optically inactive single - exciton states and for ensemble of four emitters also in the biexciton subspace .
we specified the conditions which allow the superradiance phenomena to occur in coupled inhomogeneous systems with different fundamental transition energies and dipole moments ( and thus decay rates ) .
we discussed the dynamics of single electron - hole pairs and biexcitons .
although many features typical for identical atoms , such as spontaneous trapping of excitation , may also occur in inhomogeneous qds there are differences in the dynamics of these systems . in principle , coupling between the dots induces excitation transfer between the dots which together with a possibility to define many dark states in ensembles of three and more dots lead to oscillations in the occupation of single dots .
the amplitudes of these oscillates may be considerably reduced if the system is prepared initially in an appropriate biexction state which allows for trapping of one electron - hole pair .
we envision that the presented collective effects may be controlled if the ensemble of dots is placed in the intrinsic region of a p - i - n junction with contacts constructed over the dots which due to sufficiently weak electric fields allow to control the dynamics of excitons and thus collective effects .
this work was funded by the icelandic research fund ( rannis ) and by a polish ncn grant no .
dec-2011/01/b / st3/02415 .
a. s. acknowledges support within a scholarship for outstanding young scientists granted by the polish mnisw .
the authors are grateful to pawe machnikowski for fruitful discussions .
to express an arbitrary single - exciton state in terms of the superradiant state ( [ bright ] ) and a subradiant state ( [ dark ] ) we begin with a derivation of a formula for a localized state @xmath4 .
the subradiant state which allows to cancel out all different from the state @xmath4 localized contributions to the superradiant state has a form @xmath175 and is orthogonal to the superradiant state for an arbitrarily chosen state @xmath4 .
therefore the localized single - exciton states may be decomposed into a superposition of the superradiant ( [ bright ] ) and subradiant state defined in the formula ( [ subi ] ) according to the equation @xmath176 the above formula allows to define a sub- and superradiant component in every single - exciton state @xmath177 since there is only one superradiant state in any system of @xmath0 qds , the short - living contribution @xmath178 is proportional to the state ( [ bright ] ) with a weight factor @xmath179 , while the stable part @xmath180 is a combination of @xmath0 subradiant states ( [ subi ] ) , which as a sum of dark states remains optically inactive irrespective of a number of emitters .
the component @xmath181 is proportional to the subradiant state defined by the formula ( [ dark2 ] ) with an amplitude @xmath182 .
r. h. dicke , phys . rev . * 93 * , 99 ( 1954 ) . h. m. nussenzweig , _ introduction to quantum optics _ , ( gordon and breach , nowy jork , 1973 ) . g. s. agarwal , _ quantum statistical theories of spontaneous emission and their relation to other approaches _ , springer tracts in modern physics vol . *
70 * , edited by g. hhler ( springer , berlin , 1974 ) .
m. gross and s. haroche , phys
. rep . * 93 * , 301 ( 1982 ) .
m. r. singh , w. lau , phys .
a * 231 * , 115 ( 1997 ) .
n. skribanowitz , i. p. herman , j. c. macgilvray , m. s. feld , phys .
lett . * 30 * , 309 ( 1973 ) .
d. pavolini , a. crubellier , p. pillet , l. cabaret , s. liberman , phys .
lett . * 54 * , 1917 ( 1985 ) .
m. scheibner , t. schmidt , l. worschech , a. forchel , g. bacher , t. passow , and d. hommel , nat . phys . * 3 * , 106 ( 2007 ) .
k. baumann , ch .
guerlin , f. brennecke , and tilman esslinger , nature * 464 * , 13011306 ( 2010 ) .
j. m. fink , r. bianchetti , m. baur , m. gppl , l. steffen , s. filipp , p. j. leek , a. blais , and a. wallraff , phys .
103 * , 083601 ( 2009 ) .
j. a. mlynek , a. a. abdumalikov , jr .
, j. m. fink , l. steffen , m. baur , c. lang , a. f. van loo , and a. wallraff , phys .
a * 86 * , 053838 ( 2012 ) .
p. f. herskind , a. dantan , j. p. marler , m. albert , and m. drewsen , nature physics * 5 * , 494 - 498 ( 2009 ) .
o. kyriienko , a. v. kavokin , and i. a. shelykh , arxiv:1211.0688v2 ( 2013 ) .
a. a. belyanin , v. v. kocharovsky , vl .
v. kocharovsky , and d. s. pestov , laser physics * 13 * , 61167 ( 2003 ) .
d. meiser and m. j. holland , phys . rev .
a * 81 * , 033847 ( 2010 ) .
d. meiser and m. j. holland , phys .
a * 81 * , 063827 ( 2010 ) .
j. g. bohnet , z. chen , j. m. weiner , d. meiser , m. j. holland , and james k. thompson , nature * 484 * , 78 ( 2012 ) .
a. d. ludlow , t. zelevinsky , g. k. campbell , s. blatt , m. m. boyd , m. h. g. de miranda , m. j. martin , j. w. thomsen , s. m. foreman , j. ye , t. m. fortier , j. e. stalnaker , s. a. diddams , y. l. coq , z. w. barber , n. poli , n. d. lemke , k. m. beck , and c. w. oates , science * 319 * , 1805 ( 2008 ) .
y. y. jiang , a. d. ludlow , n. d. lemke , r. w. fox , j. a. sherman , l .- s .
ma , and c. w. oates , nature photonics * 5 * , 158 ( 2011 ) . c. w. chou , d. b. hume , t. rosenband , and d. j. wineland , science * 329 * , 1630 ( 2010 ) .
t. m. fortier , n. ashby , j. c. bergquist , m. j. delaney , s. a. diddams , t. p. heavner , l. hollberg , w. m. itano , s. r. jefferts , k. kim , f. levi , l. lorini , w. h. oskay , t. e. parker , j. shirley , and j. e. stalnaker , phys .
98 * , 070801 ( 2007 ) .
s. blatt , a. d. ludlow , g. k. campbell , j. w. thomsen , t. zelevinsky , m. m. boyd , j. ye , x. baillard , m. fouch r. l. targat , a. brusch , p. lemonde , m. takamoto , f .-
hong , h. katori , and v. v. flambaum , phys .
. lett . * 100 * , 140801 ( 2008 ) .
m. gross , c. fabre , p. pillet , and s. haroche , phys .
lett . * 36 * , 10351038 ( 1976 ) .
s. filipp , a. f. van loo , m. baur , l. steffen , and a. wallraff , phys .
a * 84 * , 061805(r ) ( 2011 ) .
y. takasu , y. saito , y. takahashi , m. borkowski , r. ciuryo , and p. s. julienne , phys .
lett . * 108 * , 173002 ( 2012 ) . n. yamamoto , arxiv:1210.2632 ( 2012 ) .
x . gong , x .- b .
niu , j. li , y .- f .
huang , and g .- c .
guo , phys .
a * 77 * , 042317 ( 2008 ) p. zanardi , m. rasetti , rev
. lett . * 79 * , 3306 ( 1997 ) .
p. zanardi , f. rossi , phys .
rev . lett .
* 81 * , 4752 ( 1998 ) .
d. petrosyan and g. kurizki , phys .
lett . * 89 * , 207902 ( 2002 ) .
a. kalachev , phys .
a * 76 * , 043812 ( 2007 ) .
a. kalachev , v.v .
samartsev , quantum electronics * 35*(8 ) , 679 ( 2005 ) .
a. kalachev and s. krll , phys .
a * 74 * , 023814 ( 2006 ) s. b. zheng , phys .
a * 74 * , 054303 ( 2006 ) .
p. agrawal and a. pati , phys .
a * 74 * , 062320 ( 2006 ) .
l. li and d. qiu , j. phys .
a : math . theor . * 40 * , 10871 ( 2007 ) . h. situ and d. qiu , j. phys .
theor . * 43 * , 055301 ( 2010 ) .
radel ben - av and iaakov exman , phys .
rev a * 84 * , 014301 ( 2011 ) .
e. pazy , i. damico , p. zanardi , and f. rossi , phys .
b * 64 * , 195320 ( 2001 ) .
t. unold , k. mueller , c. lienau , t. elsaesser , and a.  d. wieck , phys .
lett . * 94 * , 137404 ( 2005 ) .
e. biolatti , r.  c. iotti , p. zanardi , and f. rossi , phys .
lett . * 85 * , 5647 ( 2000 ) .
m. goryca , t. kazimierczuk , m. nawrocki , a. golnik , j.  a. gaj , p. kossacki , p. wojnar , and g. karczewski , phys .
lett . * 103 * , 087401 ( 2009 ) .
s. mokhlespour , j. e. m. haverkort , g. slepyan , s. maksimenko , and a. hoffmann , phys .
b * 86 * , 245322 ( 2012 ) .
w. dr , g. vidal , and j. i. cirac , phys . rev . a * 62 * , 062314 ( 2000 )
. y. n. chen , c. m. li , d. s. chuu , and t. brandes , new j. of phys .
* 7 * , 172 ( 2005 ) .
j. p. kestner and s. das sarma , phys .
a * 84 * , 012315 ( 2011 ) .
m. busl , r. snchez , and g. platero , phys .
b * 81 * , 121306(r ) ( 2010 ) .
p. trocha , and jzef barna , phys .
b * 78 * , 075424 ( 2008 ) .
e. vernek , p. a. orellana , and s. e. ulloa , phys .
b * 82 * , 165304 ( 2010 ) .
q. wang , h. xie , y .- h.nie , and w. ren , phys .
b * 87 * , 075102 ( 2013 ) .
a. grodecka and p. machnikowski , phys .
b * 73 * , 125306 ( 2006 ) .
a. sitek and p. machnikowski , phys .
b * 75 * , 035328 ( 2007 ) .
b. szafran , t. chwiej , f.  m. peeters , s. bednarek , j. adamowski , and b. partoens , phys .
b * 71 * , 205316 ( 2005 ) .
b. szafran , acta phys .
a * 114 * , 1013 ( 2008 ) .
b. w. lovett , j. h. reina , a. nazir , and g. a. d. briggs , phys .
b * 68 * , 205319 ( 2003 ) .
m. kozub ,  .
pawicki , p. machnikowski , phys . rev . * 86 * , 121305 ( 2012 ) .
marlan o. scully , m. suhail zubairy , quantum optics ( cambridge university press , cambridge , 1997 ) j. danckwerts , k. j. ahn , j. frstner , and a. knorr , phys .
b * 73 * , 165318 ( 2006 )
. m. r. singh , phys .
rev . a * 75 * , 043809 ( 2007 ) .
m. r. singh , ch .
racknor , and d. schindel , appl .
phys , lett . * 101 * , 051115 ( 2012 ) .
o. benson , c. santori , m. pelton , and y. yamamoto , phys .
* 84 * , 2513 ( 2000 ) .
z. yuan , b. e. kardynal , r. m. stevenson , a. j. shields , c. j. lobo , k. cooper , n. s. beattie , d. a. ritchie , m. pepper , science * 295 * , 102 ( 2002 ) .