the quantum mechanical properties of the vacuum state in quantum field theories are presently attracting the attention of researchers from very different communities , from astrophysics and gravitation , to quantum optics and condensed matter physics .
a central idea of all these studies is the possibility of converting zero - point quantum fluctuations into observable quantum vacuum radiation by some spatial and/or temporal dependence of the background over which the quantum field is propagating  @xcite .
a number of fundamental physical effects belong to this category , from the dynamical casimir effect when the boundary conditions of the quantum field are rapidly varied in time  @xcite , to the hawking radiation when it propagates on a curved space - time showing a black hole horizon  @xcite .    in the wake of the pioneering work by unruh  @xcite
, researchers have started investigating condensed - matter systems where the propagation of some low - energy excitation field follows an effectively curved space - time geometry .
in suitable configurations showing a black hole horizon for this low - energy excitation , a quantum vacuum emission is expected to appear via a mechanism analogous to hawking radiation . among the many systems that have been investigated in this perspective
@xcite , low - temperature superfluids and nonlinear optical systems are nowadays considered the most promising ones .
dilute superfluids such as bose - einstein condensates of ultracold atoms join a very simple spectrum of elementary excitations with sonic dispersion , with the possibility of pushing the sensitivity of measurements down to the quantum limit .
their potential as analog models was first proposed in refs .   and later confirmed by _
ab initio _ numerical simulations of the condensate dynamics  @xcite .
even if acoustic black hole configurations have been experimentally realized  @xcite , no evidence of hawking radiation has been reported yet .
analog models based on nonlinear optical systems were pioneered in ref .  : a strong pulse of light propagating in a kerr nonlinear medium can be used to generate a moving spatial interface separating regions of different light velocity . for suitably chosen parameters
, the interface then behaves as a horizon that can be crossed by light in one direction only , which should result in hawking radiation being emitted .
some experimental claim in this direction was recently reported  @xcite and has raised a number of interesting questions concerning the interpretation of the observed radiation in terms of hawking processes  @xcite .
a completely new perspective to analog models based on nonlinear optical systems was opened by marino  @xcite , who first proposed the use of quantum fluids of light to generate acoustic black hole configurations and , then , to look for the hawking radiation of bogoliubov phonons on top of the photon fluid . in the following years
, this idea has been pushed forward by a number of authors , who have considered different geometrical configurations  @xcite and material systems , and very recently exciton - polariton fluids in semiconductor microcavities  @xcite .    in the present work ,
we report a comprehensive study of analog hawking radiation effects in quantum fluid of light , specifically exploiting the quantum fluid properties of an exciton - polaritons condensate .
our theoretical model fully includes their intrinsically non - equilibrium nature  @xcite and describes quantum fluctuations of the polariton field within the so - called truncated - wigner formalism of degenerate quantum gases  @xcite . as a result ,
our numerical calculations are able to provide quantitative predictions for the observable quantities , and to point out clear and accessible signatures of analog hawking radiation in the emitted light from the cavity .
even if our discussion is mostly focused on the specific case of semiconductor microcavity devices  @xcite where polariton superfluidity has been first demonstrated  @xcite , all our conclusions straightforwardly extend to generic planar cavity devices filled with or made of a kerr nonlinear medium .
the article is organized as follows . in sec .
[ theory ] we introduce the physical system under consideration and we summarize the theoretical tools that are used to describe it . in sec .
[ horizon ] we describe the original laser beam configuration that we propose to generate an analog black hole horizon with a large surface gravity .
the observable consequences of hawking mode - conversion processes are illustrated in the following sections . in sec .
[ scattering ] we discuss the scattering of a coherent phonon wavepacket on the horizon : the signature of classical hawking processes is visible as an additional wavepacket emerging from the horizon on a negative - norm branch . in sec .  [ hawking ] , we present the numerical evidence of the hawking radiation due to the conversion of zero - point quantum fluctuations into observable radiation , and we compare it with the theoretical expectations .
conclusions are finally drawn in sec .
[ concl ] .
we consider a polariton wire device where cavity photons propagating along one dimension are resonantly coupled to the fundamental exciton transition in one or more quantum wells embedded in the cavity layer ( in - plane polarized heavy - hole to conduction band transition ) . starting from a planar microcavity where light is confined along the growth axis @xmath0 by a pair of distributed bragg reflectors ,
a suitable etching procedure is used to laterally pattern the device in the @xmath1 direction , and create the ridge structure that is schematically shown in fig .
[ fig1](a ) . provided the losses are weak enough , the bosonic excitations that result from the strong coupling of the cavity photon with the quantum well exciton have the typical mixed light - matter nature of one - dimensional exciton - polaritons , free to propagate along the wire axis according to the dispersion law that is plotted in fig .
[ fig1](b ) . a high quality experimental realization of this polariton wire concept was reported in ref .
, from which we shall take the state - of - the art system parameters to be used in the numerical simulations .
along the wire axis .
the thin dashed and dash - dotted lines indicate the bare cavity photon and exciton dispersions , respectively . the thick dashed line indicates the parabolic approximation of the lower polariton dispersion .
system parameters are inspired from ref .
: @xmath2  mev , @xmath3  mev , @xmath4 ( with @xmath5 the free electron mass ) , @xmath6 . ,
title="fig:",scaledwidth=45.0% ] +   along the wire axis .
the thin dashed and dash - dotted lines indicate the bare cavity photon and exciton dispersions , respectively .
the thick dashed line indicates the parabolic approximation of the lower polariton dispersion .
system parameters are inspired from ref .
: @xmath2  mev , @xmath3  mev , @xmath4 ( with @xmath5 the free electron mass ) , @xmath6 .
, title="fig:",scaledwidth=48.0% ] +      a theoretical description of the dynamics of this system can be developed in terms of the standard hamiltonian for the two coupled bosonic fields describing the quantum well exciton and the cavity photon  @xcite @xmath7 where @xmath8 is the longitudinal coordinate along the wire axis and the two - component operator vector @xmath9 summarizes the quantum fields describing the exciton @xmath10 and photon @xmath11 fields .
each of them satisfies one - dimensional bosonic commutation rules , e.g. @xmath12=\delta_{i , j}\delta(x - x')$ ] , with @xmath13 . throughout the whole paper ,
we restrict ourselves to the case where the system is pumped on a single spin state , so that the spin degrees of freedom can be neglected in the theoretical model .
the single - particle hamiltonian @xmath14 describing the evolution of the non - interacting exciton and cavity photon fields has the simple representation @xmath15 in terms of the bare cavity photon and exciton dispersion law @xmath16 in a spatially homogeneous one - dimensional system .
the rest energy @xmath17 and the effective mass @xmath18 of the cavity photon are determined by the spatial confinement along the @xmath0 and @xmath1 axis by the dbr mirrors and by the total internal reflection at the etching interfaces , respectively . in practical calculations
the exciton dispersion can be safely neglected , @xmath19 , as the exciton mass @xmath20 is orders of magnitude larger than the photon mass .
the strength of the light - matter coupling is quantified by the rabi frequency @xmath21 , proportional to the square root of the oscillator strength per unit area of the excitonic transition  @xcite .    in spatially homogeneous systems , for which the external potentials acting on the cavity photon and the quantum well exciton vanish ( @xmath22 ) , diagonalization of the single particle hamiltonian @xmath14 as a function of wave vector @xmath23 leads to the well - known upper and lower polariton branches of dispersions @xmath24 and @xmath25 . as it is shown in fig .
[ fig1](b ) the bottom part of the lower polariton dispersion is accurately approximated by a parabolic form @xmath26 where the @xmath27 and @xmath28 parameters have the clear physical meaning of a rest energy and an effective mass for the ( lower ) polariton . in the resonant case @xmath29
, one has the simple relation @xmath30 .
several methods to generate external potentials @xmath31 with almost arbitrary spatial shapes along the @xmath32 direction have been experimentally demonstrated in the last years , from the in - plane patterning of the microcavity layer  @xcite , to the all - optical potential created by a strong pump laser with opposite polarization  @xcite , to the application of a mechanical stress to the device  @xcite .
exciton - exciton interaction , due to the coulomb interactions between the electron and holes forming the excitons , are the main source of optical nonlinearity in this problem .
experimental measurements  @xcite report a value on the order of @xmath33 , which translates into a reduced one - dimensional interaction parameter to be used in the hamiltonian ( [ eq : model ] ) @xmath34 , where the integral is over the transverse wire width ( @xmath35 ) , and the forth power of the exciton wave function comes from the kerr - type nature of the nonlinearity . assuming the transverse confinement of the exciton envelope function to be described by @xmath36 , a good estimate can be given to be @xmath37 .
taking into account the exciton weight in the polariton , this translates in the resonant @xmath38 case into a one - dimensional polariton - polariton interaction constant @xmath39 .
polaritons are coherently injected into the system by an external laser drive coupled to the cavity photon via the non - perfect reflectivity of the dbr mirrors . in the theoretical model ,
such processes are described by the terms in the last row of eq .
( [ eq : model ] ) , where @xmath40 is the spatio - temporal profile of the coherent pump beam including the coupling coefficient proportional to the dbr transmission amplitude .
correspondingly , photons ( excitons ) are subjected to radiative losses at a rate @xmath41 ( @xmath42 ) that have to be described at the level of the master equation for the density matrix , @xmath43 , for the coupled quantum fields @xmath44+\mathcal{l}[\rho].\ ] ] performing the standard born - markov approximation  @xcite , the dissipation superoperator can be cast into the lindblad form @xmath45   \ , .\end{gathered}\ ] ] the resulting decay rate of polaritons is a weighted average of the exciton and photon ones @xmath46 : in the vicinity of the bottom of the lower polariton branch of fig .
[ fig1](b ) , it has an almost momentum - independent value @xmath47 .      in a regime of weak exciton - exciton interactions ,
the dynamics of the system can be accurately captured by the mean - field approximation : the quantum fields @xmath48 are replaced by the classical fields corresponding to their expectation values @xmath49 and evolving according to the pair of nonlinear partial differential equations @xmath50 \left (   \begin{array}{c } \phi_{\mathrm{x } } ( x , t ) \\    \phi_{\mathrm{c } } ( x , t ) \end{array }    \right )   \ , , \end{gathered}\ ] ] which generalize to the non - equilibrium context of polaritons the well - known gross - pitaevskii equation ( gpe ) of dilute bose condensed gases  @xcite .
[ sec : homog ] in the simplest case of a spatially homogeneous system under a coherent pump in a plane wave form @xmath51\ ] ] of frequency @xmath52 , wavevector @xmath53 and amplitude @xmath54 , the analog of the equation of state can be obtained by simply injecting the same plane wave ansatz @xmath55 $ ] into the generalized gpe ( [ eq : mf_gpe ] ) .
@xmath56m@xmath57 , and frequency @xmath58 mev tuned slightly above the lower polariton branch at @xmath59  mev .
equal loss rate are assumed for the exciton and the photon fields , @xmath60 mev . the value @xmath61 of the nonlinearity corresponds to a @xmath62 m wide polariton wire ( see text ) .
lower panel : spatial intensity profile of the coherent pump along the wire axis , with excitation parameters : @xmath63 , @xmath64 , @xmath65 , @xmath66 .
, scaledwidth=50.0% ]    this leads to the simple form @xmath67 for the steady - state excitonic density @xmath68 as a function of the pump parameters , and to an analogous expression for the cavity - photon density  @xcite .
the pump detuning from the cavity photon and exciton modes at wavevector @xmath69 are indicated here as @xmath70 and @xmath71 .    in the quantum fluid language ,
this equation plays the role of the equation of state for the polariton fluid under a coherent pump in a plane wave form : while the polariton flow speed is controlled by the incident wavevector @xmath69 , the density has a more complex dependence on the pump intensity @xmath72 and frequency @xmath52 .
a most remarkable example is illustrated in fig .
[ fig2]a : the choice of a pump frequency slightly above the lower polariton branch leads to an optical bistable behavior .    for reasons that will soon become clear , we shall concentrate our attention onto the case where the pump intensity is tuned in the close vicinity of the end - point of the upper branch of the hysteresis loop ( indicated as @xmath73 in fig .
[ fig2]a ) . in order to reach this working point in a given spatial region ,
a pump with a non - trivial spatial envelope , @xmath74 , multiplying the plane wave can be used , as shown in the lower panel of fig .
[ fig2](b ) .
upstream of the region of interest ( i.e. for @xmath75 the pump amplitude is tuned to a value @xmath76 higher than the switch - on point of the hysteresis loop ( i.e. the end point of the lower branch ) , so to switch the system to the upper branch of the hysteresis loop  @xcite . in the long region of interest ( @xmath77 ) ,
the pump amplitude is maintained at a spatially constant value @xmath73 , so to keep the system sits at the desired working point with the desired in - plane flow wavevector @xmath69 .
the behavior of the density @xmath78 and the wavevector @xmath79 in the downstream region ( @xmath80 ) is determined by a complex interplay of propagation and losses .
first numerical and analytical studies of this regime have appeared in refs .   and
: the spatial drop of the polariton density away from the pump results in an increase of the polariton speed under the effect of the interactions .
the lower the losses , the slower this acceleration effect .
the dispersion of collective excitations on top of the steady state of the coherently pumped condensate is obtained as usual by linearizing the gpe around the steady state solution  @xcite .
a complete discussion of the many different behaviors that are possible depending on the pump frequency and intensity can be found in refs .  . here
, we shall restrict our attention to the two cases of present interest , namely ( i ) a coherent pump whose intensity is tuned right at the turning point of the upper branch of the hysteresis loop , and ( ii ) the ballistic flow regime in the downstream region .
if the polariton population is restricted to the lower polariton branch and we perform a parabolic approximation of the branch bottom with effective mass @xmath28 , the dispersion of the collective excitations on top of a condensate of polariton density @xmath81 in uniform flow at speed @xmath82 has , in both cases , the usual bogoliubov form @xmath83 in the first ( i ) case , the local wavevector @xmath84 is fixed to the in - plane component @xmath53 of the pump wavevector . in the second
( ii ) case , hydrodynamic effects  @xcite make the local wavevector @xmath85 to slowly vary in space . in both cases ,
the local in - plane flow speed @xmath82 is related to the local wavevector by @xmath86  , and the polariton speed is defined as @xmath87 in terms of the polariton current @xmath88 \simeq   j_{\rm c}=\frac{1}{2i m_{\rm c}}[\hat{\psi}^*_{\rm c}\,\nabla { \psi}_{\rm c}- \psi_{\rm c}\,\nabla{\psi}^*_{\rm c}]$ ] . given its very heavy mass , @xmath6 , the exciton fraction does not significantly contribute to the polariton current . ] .    for small wave vectors , @xmath89
, the bogoliubov dispersion tends to a ( lossy ) sonic dispersion @xmath90 with a sound speed @xmath91 of fluid excitations such that @xmath92 , and a doppler shift term due to the background fluid flow at @xmath82 . for large wavevectors
, it recovers a ( lossy ) parabolic single - particle dispersion @xmath93 with a hartree energy shift @xmath94 due to the interactions with the condensate .
plots of the bogoliubov dispersion at different spatial positions are shown in the small panels of figs .
[ fig : nobh_scheme ] and [ fig : bh_scheme ] . in the following of the paper ,
the will be used to interpret the numerical results for the evolution of excitation wavepackets and for the density - density correlations .    before proceeding ,
it is important to emphasize the conditions underlying the dispersion ( [ eq : sonic ] ) of the collective excitations . in the coherent pump case ( i )
, this form only holds at the end point of the upper branch of the hysteresis loop , while gaps and/or band sticking regions appear in all other cases because of the coherent pump locking the condensate phase  @xcite .    on the other hand ,
the presence of a soft excitation branch of sonic nature is guaranteed in the ballistic case ( ii ) by the fact that the condensate phase is completely free to evolve both in space and time  @xcite . however , as the density is spatially varying , the very concept of dispersion relation holds only locally and is limited to excitations whose wavelength is small as compared to the characteristic length scale of the density profile .      a standard technique to go beyond the mean - field approximation and include the fluctuations of the polariton field around its mean value is based on the wigner representation of quantum fields . within the _ truncated wigner approximation _
, the dynamics of the quantum field problem can be described by a stochastic partial differential equation that can be numerically solved .
this technique was first introduced in the context of quantum fluids in ref .   and its application to atomic condensates was fully developed and characterized in ref .  .
its extension to the polariton case was developed in ref .
: remarkably , the presence of loss and pump terms in the stochastic equations suppresses some of the difficulties of the truncated wigner method and guarantees its accuracy as long as interactions between individual polaritons are weak enough . if a spatial grid of real - space spacing @xmath95 is used to numerically solve the stochastic partial differential equation , the weak interaction condition can be formulated as @xmath96 .
while this condition rules out the possibility of using the truncated wigner method to study strongly correlated polariton states , e.g. in the polariton blockade regime  @xcite , it does not hinder the study of quantum hydrodynamics effects such as the analog hawking radiation , which originate from the collective dynamics of a large number of polaritons .
the stochastic partial differential equations of the truncated wigner method for the coupled exciton and polariton fields has the explicit form  @xcite    @xmath97 \left (   \begin{array}{c } \phi_{\mathrm{x } } \\
\phi_{\mathrm{c } } \end{array }    \right)\right\}\,dt +   \frac{\hbar}{\sqrt{4\,\delta x } }   \left (   \begin{array}{c } \sqrt{\gamma_{\mathrm{x}}}\,dw_{\mathrm{x } } \\
\sqrt{\gamma_{\mathrm{c}}}\,dw_{\mathrm{c}}\end{array }    \right ) \ , , \ ] ]    where @xmath95 is the spacing of the real - space grid and @xmath98 and @xmath99 are complex valued , zero - mean , independent gaussian noise terms , with white noise correlation in both space and time @xmath100 with @xmath13 . in the numerical simulations , we reconstruct the equilibrium wigner distribution by first letting the system evolve to its steady state under the monochromatic pump , and then taking a large number of independent configurations by sampling the stochastic evolution at different times spaced by @xmath101 . in order to ensure statistical independence of the different realizations , a large enough @xmath101 has to be taken such that @xmath102 . in practice ,
a number of realization on the order of @xmath103 is used , taken at time intervals on the order of @xmath104 . as usual in wigner approaches  @xcite , the stochastic averages over the configurations of different functions of the fields provide the expectation value of the corresponding symmetrically ordered operator .
inspired from previous numerical studies of acoustic black holes in atomic condensates  @xcite , we shall look for the signature of analog hawking radiation in the normalized , zero - delay correlation of the cavity field intensity defined as @xmath105 within the wigner formalism , the different contributions to ( [ eq : g2def ] ) have the form  @xcite @xmath106 and @xmath107 where the @xmath108 averages indicate the classical averages over the different stochastic configurations of the field sampled at interval times @xmath101 .
after having reviewed the technical tools to study the dynamics of the polariton quantum fluid , we can now proceed to discuss realistic configurations that can be used to generate an analog acoustic black hole with a large surface gravity .
all simulations are performed by solving the generalized gross - pitaevskii equation ( [ eq : mf_gpe ] ) using experimental parameters taken from the recent work in ref .  .
the most straightforward configuration to generate an acoustic horizon was proposed by solnyshkov et al .
@xcite and is illustrated in fig .
[ fig : nobh_scheme](b ) : the presence of losses is responsible for a spatially decreasing density profile of a ballistically flowing polariton condensate and a correspondingly increasing velocity profile .
eventually , the polariton density tends to zero and the flow becomes necessarily super - sonic .
the polariton density and speed in the pumped region is fixed by the pump beam parameters : if these are suitably chosen to have an initially sub - sonic flow , a horizon must necessarily appear at some point .    in spite of the simplicity of this configuration
, some care has to be paid to the value of hawking temperature that can be expected . for smooth horizons within the hydrodynamic limit  @xcite ,
the hawking temperature @xmath109 is determined by the surface gravity @xmath110   \right|_{x_{\rm hor } }   { \label{eq : surf_grav}}\ ] ] according to @xmath111 , with @xmath112 the boltzmann constant . according to our numerical simulations and analytical approximations  @xcite ,
the characteristic length scale of the variation of the flow parameters ( density and flow speed ) is inversely proportional to the loss rate @xmath113 .
direct combination of these two general facts shows that the achievable values of the hawking temperature is limited from above by the loss rate : for the case in the figure , one indeed has @xmath114 mev .     and
the sound @xmath115 speeds at steady state .
the parameters and the spatial shape of the pump laser are the same as in fig .
( c ) dispersion of the collective excitations at two spatial positions : point @xmath116 ( left panel ) is chosen to be in the sub - sonic flat - top region ; point @xmath117 ( right panel ) is located downstream of the horizon in the super - sonic region .
red solid ( blue dashed ) lines refer to positive ( negative ) norm bogoliubov modes .
( d ) spatial profile of the exciton density at a time @xmath118  ps after the arrival of the probe pulse .
the probe laser parameters are chosen so to resonantly excite the bogoliubov branch indicated as @xmath119 in panel ( c ) : frequency @xmath120 mev , wavevector @xmath121 .
the probe spot has a gaussian shape centered at @xmath122 m with a waist @xmath123 m and a peak amplitude @xmath124 mev@xmath125m@xmath126 .
its temporal shape is also gaussian of duration @xmath127 ps .
all numerical calculations are performed in a 800-@xmath56 m simulation box with 2048 lattice points and absorbing boundary conditions at the edges of the box .
, scaledwidth=46.0% ]    a possible solution to overcome this difficulty , is to insert a narrow repulsive potential in the ballistic flow region , close to the edge of the pump spot .
the efficiency of this technique is illustrated by the numerical calculations shown in fig .
[ fig : bh_scheme](b ) .
flow across the defect mostly occurs via tunneling processes through the potential barrier : this produces a sudden drop of the condensate density and a corresponding sudden increase of the flow speed . as a result
, a horizon appears in the vicinity of the defect center , separating a sub - sonic upstream region from a super - sonic downstream one , with a surface gravity as high as @xmath128 mev , i.e. about a factor @xmath129 larger than the purely ballistic flow case of fig .
[ fig : nobh_scheme ] .    .
the parameters and the spatial shape of the pump laser are again the same as in fig .
the defect consists of a square potential barrier of thickness @xmath130 m and is located at a position @xmath131 m close to the flat - top edge .
the height of the defect potential is @xmath132 mev and acts on both exciton and photon components of the polariton field .
the horizon appears in the close vicinity of the defect with a large surface gravity .
the probe laser parameters are chosen so to resonantly excite the bogoliubov branch indicated as @xmath119 in panel ( c ) : frequency @xmath133 mev , wavevector @xmath134 .
the shape and the amplitude of the pump spot are the same as in fig .
[ fig : nobh_scheme](d ) . , scaledwidth=47.0% ]    another difficulty of the purely ballistic flow configuration of fig .
[ fig : nobh_scheme](b ) was the significant spatial variation of the polariton density on both sides of the horizon , which may strongly distort the geometrical structure hawking signal .
the defect configuration of fig .
[ fig : bh_scheme](b ) appears to be favorable also in this respect : as the flow in the upstream region is sub - sonic , no bogoliubov - cerenkov emission can take place from the defect in this direction  @xcite and the polariton density remains almost flat across the whole flat - top region .
the slow spatial variation of the density in the downstream super - sonic region is hardly avoided unless more complex laser configurations are used , but does not seem to prevent identification of the hawking effect .    on the other hand
, the use of a coherent pump raises the crucial issue of avoiding all those branch - sticking and gap effects typical of non - equilibrium condensates  @xcite that would spoil the analogy with gravitational systems . as we have reviewed in sec .
[ sec : bogoliubov ] , choosing the value of the laser amplitude in the flat - top region in the vicinity of the end - point of the upper branch of the hysteresis loop is enough to have a linear , sonic dispersion of the low-@xmath135 dispersion of the form ( [ eq : sonic ] ) at all spatial positions around the horizon  , and hence has no influence on the final results , which further supports the robustness of our scheme . ] .    as it was illustrated in fig .
[ fig2 ] , keeping the system stably on the upper branch of the hysteresis loop requires including a region of higher pump intensity upstream of the flat - top .
exception made for its role in switching on the optical bistability , the high pump intensity region play no role in the hawking physics . as we have discussed in sec .
[ sec : bogoliubov ] , the bogoliubov dispersion in the downstream region is guaranteed to be sonic at low wavevectors by the ballistic nature of the flow in this region .
the flow speed is guaranteed to be super - sonic by the suppressed value of the density beyond the defect .
ps , ( b ) @xmath136 ps , and ( c ) @xmath137 ps after the arrival of the probe pulse .
each panel displays the excitonic density modulation with respect to the time - independent steady state .
the pump and probe parameters are the same as for the purely ballistic flow shown in fig .
[ fig : nobh_scheme](d ) .
the labels on the wavepackets refer to the bogoliubov dispersions of fig .
[ fig : nobh_scheme](c ) , while the arrows indicate the propagation direction of each wave packet .
, scaledwidth=50.0% ]     ps , ( b ) @xmath136 ps , and ( c ) @xmath137 ps after the arrival of the probe pulse .
each panel displays the excitonic density modulation with respect to the time - independent steady state .
the pump and probe parameters correspond to the defect configuration shown in fig .
[ fig : bh_scheme](d ) .
the labels on the wavepackets refer to the bogoliubov dispersions of fig .
[ fig : bh_scheme](c ) , while the arrows indicate the propagation direction of each wave packet .
, scaledwidth=50.0% ]
the numerical simulations that we have presented in the previous section demonstrate the possibility to create configurations of a polariton flow through a horizon separating regions of sub- and super - sonic regimes with a sizable surface gravity . this configuration will be a powerful workbench for studies of hawking phenomena .
inspired by related theoretical work on atomic condensates  @xcite and experiments in classical hydrodynamics  @xcite , in this section we shall first consider the scattering of a phonon wavepacket off the horizon , a process where the classical counterpart of the hawking effect manifests itself as an additional reflected wavepacket .
this discussion will be the starting point for the next section , where we shall investigate the very hawking effect , namely the conversion of zero - point quantum fluctuations into observable radiation by the horizon .
the scattering dynamics of bogoliubov phonons hitting the horizon can be studied in terms of the generalized gpe ( [ eq : mf_gpe ] ) : far from the horizon , the propagation of weak wavepacket perturbations follows the bogoliubov dispersion discussed in sec .
[ sec : bogoliubov ] albeit with spatial dependent flow and sound speeds . figs .
[ fig : nobh_scheme ] to [ fig : bh_snap ] illustrate the different aspects of this physics for the two previously mentioned cases , namely a smooth flow with a low surface gravity a defect configuration with a sizable surface gravity .
for each configuration , the structure of the flow is illustrated in figs .
[ fig : nobh_scheme](b ) and [ fig : bh_scheme](b ) , while the smaller ( c , d ) panels show the bogoliubov dispersion at different spatial positions .
once the system is in its steady state , a low - frequency wavepacket perturbation is generated in the downstream super - sonic region by means of an extra coherent laser pulse with a spatio - temporal shape of the form @xmath138 the probe wave vector ( @xmath139 ) and frequency ( @xmath140 ) are chosen to be resonant with the bogoliubov branch labeled as @xmath119 in the ( c ) panels .
the probe beam is centered at @xmath117 ; its waist @xmath141 determines the spatial size of the generated wavepacket ( and its inverse width in @xmath135-space ) .
the temporal duration @xmath142 has to be long enough not to excite the upper polariton branch nor the negative - norm bogoliubov branch .
its amplitude @xmath143 is chosen weak enough to remain within the validity domain of a linearized description of excitations .    the generated wavepacket shown in figs .
[ fig : nobh_scheme](d ) and [ fig : bh_scheme](d ) then propagates in the leftward direction against the horizon .
a series of snapshots of this evolution are presented in figs .
[ fig : nobh_snap ] and [ fig : bh_snap ] .
the difference between the behaviors of the weak and of a large surface gravity cases is apparent : in the former case , only two wavepackets visibly emerge from the horizon at long times , located on the @xmath144 and @xmath145 bogoliubov branches . on the other hand , in the latter case , an extra wavepacket appears on the @xmath146 branch .
remarkably , the bogoliubov norm of this wavepacket is negative , which signals the occurrence of hawking conversion processes .
the identification of the different wavepackets in terms of their bogoliubov branch is confirmed by the numerical measurement of their group velocity and of their wavevector .    of course , the physical reason why we did not observe a @xmath146 wavepacket in fig .  [ fig : nobh_snap ] is that the low value of the surface gravity strongly suppresses the amplitude of the hawking process .
as it was expected on general grounds  @xcite and then explicitly verified  @xcite , the hawking scattering amplitude for a bogoliubov wavepacket of given carrier frequency @xmath147 scales down exponentially with the surface gravity according to the boltzmann factor @xmath148 . on the other hand , excitation of the @xmath146 branch
would be totally forbidden if the flow was everywhere sub - sonic and no horizon was present : in this case , in fact , energy could not be conserved in the hawking scattering process .
this quantity directly reflects in the intensity correlations of a near - field image of the emitted light from the microcavity . the upper ( a ) panel refers to the case of a purely ballistic flow considered in figs .
[ fig : nobh_scheme ] and [ fig : nobh_snap ] .
the lower panel ( b ) refers to the case with a defect in the ballistic flow considered in figs .
[ fig : bh_scheme ] and [ fig : bh_snap ] .
the colored dashed lines in ( b ) indicate the different hawking features as discussed in the text.,scaledwidth=50.0% ]    the process studied in the previous section provides a solid evidence of the classical counterpart of the hawking effect , namely the possibility of inter - converting positive and negative norm waves at the horizon .
now we shall proceed with the study of the hawking effect _
stricto sensu _ ,
i.e. the emission of radiation by the horizon via the inter - conversion of zero - point quantum fluctuations into observable radiation .    to this purpose
, one has to fully include in the model the quantum fluctuations of the fields around the mean - field value predicted by the gross - pitaevskii equation . in this work ,
we adopt the truncated wigner approach described in sec .
[ sec : montecarlo ] , where expectation values of quantum operators are calculated as classical averages over a suitably chosen stochastic partial differential equation .
inspired from previous theoretical  @xcite and numerical  @xcite work , we shall consider the spatial correlation function of the quantum fluctuations of the fluid density : the smoking gun of hawking radiation consist of a negative long - range correlation signal between points located on opposite sides of the horizon .
numerical plots of the normalized density - density correlation function @xmath149 are shown in the two panels of fig .
[ fig : g2 ] for the two cases of a weak and strong surface gravity @xmath150 , respectively .
while in the upper panel , fig .
[ fig : g2](a ) ( weak surface gravity , @xmath151 mev ) there is no clear feature emerging over noise , the lower panel , fig .
[ fig : g2](b ) ( strong surface gravity , @xmath152 mev ) shows several features that can be interpreted as a direct consequence of the analog hawking radiation  @xcite : the two bogoliubov excitations forming the hawking pair are simultaneously emitted from the horizon . at later times ,
the temporal correlation between the emission time of the hawking partners reflects into a long - distance correlation in the density fluctuations due to each member of the pair .
while in the atomic case this correlation extends at long times to indefinitely large distances from the horizon , the finite decay rate of bogoliubov excitations in polariton fluids predicted in ( [ eq : sonic ] ) restricts the correlation signal to a finite region of size @xmath153 around the horizon , @xmath154 being the group velocity of the bogoliubov branch under consideration .
away from this region around the horizon , the hawking phonons have decayed and the corresponding density fluctuations have disappeared .
the dashed lines in fig .  [ fig : g2](b ) indicate the location at which one would expect the maximal density correlations due to the quantum vacuum emission process into the @xmath155 ( white ) , the @xmath156 ( purple ) , the @xmath157 ( yellow ) pair of modes .
the first process is the traditional hawking emission process where the two quanta are emitted into opposite directions from the horizon .
the main properties of the density correlation pattern due to this process is well captured by the non - dispersive quantum field theory on curved space - time as proposed in ref .  .
for @xmath158 and @xmath159 , the correlation is peaked along the half - line @xmath160 here , @xmath161 and @xmath162 are the flow speed and speed of sound in the upstream and downstream regions , so that @xmath163 and @xmath164 are the group velocity of the @xmath146 and @xmath144 modes in the low-@xmath135 , sonic region .
while in the atomic bec case of ref .
their value did not depend on space in the asymptotic regions away from the horizon , the driven - dissipative nature of polaritons is responsible for the spatial dependence of the flow @xmath82 and sound @xmath165 speeds that is visible in fig .
[ fig : bh_scheme](b ) , mostly in the downstream region .
however , as the space dependence remains moderate within the regions of thickness @xmath166 and @xmath167 on either side of the horizon where correlations are significant , no appreciable curvature of the hawking tongues is visible in fig .  [
fig : g2](b ) .    the orientation of the purple and yellow dashed lines in fig .
[ fig : g2](b ) is analogously obtained from the group velocity of the corresponding modes , namely @xmath168 and @xmath169    the excellent agreement of the geometrical location of the numerically observed features with the analytical predictions ( [ eq : ud2][eq : d1d2 ] ) confirms our interpretation . remarkably  , the peak value of the @xmath155 hawking correlation signal is on the order of @xmath170 , which is not too far from the analytical prediction of ref .  .
as the surface gravity is here comparable to the interaction energy @xmath171 , the quantitative discrepancy can be easily traced back to the breakdown of the hydrodynamical approximation underlying the analytical model  @xcite .
thanks to the continuous - wave nature of the proposed experimental setting , the quantitatively low value signal to be observed can be overcome by a sufficiently long integration time .
while agreement with the atomic case is good also for the @xmath156 feature , the @xmath157 one has a different sign .
a possible explanation of this behavior can be traced back to the non - universality of the @xmath172 back - scattering effect that is responsible for the conversion at the horizon of the standard @xmath155 hawking correlation into the @xmath157 one .    before concluding , it is important to note that throughout the whole discussion , we have implicitly assumed that the @xmath173 correlation function is assumed to be measured at the same time : a rough estimate of the required temporal resolution of the detector is given by @xmath174 , where @xmath175 is the speed of sound in the polariton gas ( on the order of @xmath176m / ps from figs .
[ fig : bh_scheme ] ) and @xmath177 is the spatial width of the correlation signal . using the value @xmath178 m taken from fig .
[ fig : g2 ] , one can estimate the needed temporal resolution to be on the order of @xmath179ps , which is within the state - of - the - art of optical technology .
in this article we have presented a comprehensive study of classical and quantum hydrodynamic properties of acoustic black holes in superfluids of exciton - polaritons in semiconductor microcavities . inspired from on - going experimental research ,
we have identified and characterized one - dimensional polariton wire devices as model systems , where the polaritons injected by a single monochromatic laser beam with a suitable spot profile generate a flow configuration showing an acoustic black hole horizon .
inserting a repulsive potential defect provides a dramatic enhancement of the analog surface gravity at the horizon . even if our calculations have been performed for the specific case of a semiconductor microcavity in the strong - coupling regime , all conclusions straightforwardly extend to a generic planar cavity filled by a nonlinear optical medium , as originally suggested in ref .  , or photonic crystal polariton states with engineered mode dispersion  @xcite .    in analogy with previous works on acoustic black holes in atomic condensates and classical hydrodynamics of surface waves ,
two experiments have been proposed and numerically simulated to assess the efficiency of hawking wave - conversion processes at the horizon . when a coherent wavepacket of bogoliubov phonons is incident on the horizon , the hawking effect is visible as the production of an extra emerging wavepacket on the negative - norm branch .
this is the classical evidence for the hawking process taking place at the horizon . on the other hand , the usual hawking
radiation is a purely quantum effect consisting of the conversion of zero - point fluctuations into a stream of observable phonons at the horizon : the correlation between the emission times of a phonon and its hawking partner is visible as a long - distance correlation between the density fluctuations on either side of the horizon .
the density correlations provide the clearest experimental signature of hawking radiation in superfluids : in the present polariton case , in - cavity density correlations directly transfer into the secondary emission from the cavity photon , and can be observed as analogous correlations in the intensity fluctuations of the near - field emission .
as soon as a sufficient signal over noise ratio is achieved , an experiment along the proposed lines will hopefully provide a clear evidence of a fundamental effect of quantum field theory that so far has eluded all experimental observation , overcoming the barrier to direct observation imposed by the negligibly small hawking temperature of astrophysical black - holes .
we are grateful to a. amo , r. balbinot , j. bloch , a. bramati , s. finazzi , s. robertson , d. sanvitto for continuous stimulating discussions , and d. sarchi for a fruitful collaboration at an early stage of the work .
dg warmly acknowledges the friendly hospitality of the ino - cnr bec center in trento , where most of this work has been conceived and carried out , and l.c .
andreani for supporting this research .
r. schtzhold and w. g. unruh , phys .
lett . * 107 * , 149401 ( 2011 ) ; f. belgiorno , s. l. cacciatori , m. clerici , v. gorini , g. ortenzi , l. rizzi , e. rubino , v. g. sala , and d. faccio , phys . rev . lett . * 107 * , 149402 ( 2011 ) .
a. amo , d. sanvitto , f.p .
laussy , d. ballarini , e. del valle , m. d. martin , a. lematre , and j. bloch , d. n. krizhanovskii , m. s. skolnick , c. tejedor , and l. via , nature * 457 * , 291 ( 2009 ) ; a. amo , j. lefrere , s. pigeon , c. adrados , c. ciuti , i. carusotto , r. houdr , e. giacobino , and a. bramati , nat . phys . * 5 * , 805 ( 2009 ) ; a. amo , s. pigeon , d. sanvitto , v. g. sala , r. hivet , i. carusotto , f. pisanello , g. lemenager , r. houdr , e. giacobino , c. ciuti , and a. bramati , science * 332 * , 1167 ( 2011 ) .
o. el daf , a. baas , t. guillet , j .- p .
brantut , r. idrissi kaitouni , j. l. staehli , f. morier - genoud , and b. deveaud , appl .
. lett . * 88 * , 061105 ( 2006 ) ; d. lu , j. ahn , s. freisem , d. gazula , and d. g. deppe , appl .
. lett . * 87 * , 163105 ( 2005 ) .
d. f. walls and g. j. milburn , _ quantum optics _ ( springer verlag , berlin , 2006 ) ; c. cohen - tannoudji , j. dupont - roc , and g. grynberg , _ atom - photon interactions . basic processes and applications . _
( wiley science paperback series , new york , 1998 ) .


fu orionis is the prototype of a small , but quite remarkable class of low - mass young stellar objects ( ysos ) normally referred to as fu ori objects ( fuors ) .
for the first members of this class an outburst in optical light of up to 4 - 6 magnitudes over short time scales , followed by a decrease in luminosity over several years or decades , was observed @xcite .
other objects were included in the class as they shared common specific spectroscopic features , e.g. , double - peaked line profiles , a spectral type varying with wavelength and often co bandhead absorption features . as for one object ( v1057 cyg ) the pre - outburst spectrum is known and as it resembles that of a classical ttauri star @xcite , it is commonly assumed that fuors should be low - mass ysos .
most observational data can be explained by the presence of an accretion disk surrounding the young stars ( * ? ? ?
* however , found that some spectral properties are well explained in the context of rapidly rotating late - type stars with strong stellar winds ) . a dramatic temporal increase in the accretion rate , where the disk outshines the star by several orders of magnitude , can account for the observed outbursts in luminosity .
several scenarios , possibly triggering such an increased accretion rate , have thus far been proposed .
they include ( a ) interactions of binary or multiple systems where tidal forces disturb the circumstellar disk @xcite , ( b ) planet - disk interactions , where thermal instabilities in the disk are caused by the presence of a massive planet @xcite , or ( c ) thermal instabilities in the disk alone @xcite . for a detailed overview concerning the fu ori phenomenon we refer to @xcite .
apart from revealing the mechanism leading to the observed outbursts , it is also important to investigate whether all ttauri stars undergo such epochs of enhanced accretion or whether fuors are a special class of ysos .
most observations of classical ttauri stars show that the derived accretion rates of @xmath1 ( e.g. , * ? ? ?
* ) are not sufficient to build up a low - mass star over time scales of a few myr .
even if part of the matter is supposed to be accreted in the very early phases of a yso , fuor - phases might provide an elegant solution to this problem as they would speed up the accretion process .    as so far only a dozen or so well established fuors are known and as the fuor - phase has significant impact on the young star - disk system it is crucial to combine as much observational information as possible for these objects to eventually derive a coherent theoretical picture . based on near - infrared ( nir ) and/or
mid - infrared ( mir ) interferometry as a technique to study the inner few au of the accretion disks new insights to some of the best studied fuors were provided recently .
@xcite found that accretion disks alone can not reproduce the sed and observed low k - band visibilities for v1057 cyg , v1515 cyg and z cma - se simultaneously .
they concluded that additional uncorrelated flux may arise due to scattering by large dusty envelopes .
@xcite presented the first vlti/ midi observations of v1647 ori whose eruptive behavior suggests that it is either an fuor or an ex lupi ( exor ) type object . in this case
it was possible to fit both , the sed and the observed mir visibility with a simple disk model with moderate disk flaring .
for fu ori itself it was @xcite who analyzed a wealth of nir interferometric data .
they showed that the nir visibilities and the sed could be fitted with two models : one consisting simply of an optically thick and geometrically thin accretion disk and a second one consisting of an accretion disk and an embedded `` hot spot '' . from their error statistics
these authors concluded that the latter model was more likely .    in summary
it becomes clear that up to now no coherent picture can be derived from the interferometric observations of fuors and the group of objects seems to be rather inhomogeneous .
in this paper we present the first mir interferometric measurements of fu orionis .
the data are thus complementary to the nir observations of malbet et al .
( 2005 ) . in section 2
we briefly describe the observations and the data reduction process . in section 3
we discuss the findings derived from the n - band acquisition images .
the 8 - 13@xmath0 m spectrum of fu ori is analyzed in section 4 . in sections 5 and 6
we discuss the results from the interferometric measurements , i.e. the visibilities and the correlated flux spectra , respectively . a simple analytical disk model and its fit to the sed and the visibilities is presented in section 7 .
finally , we summarize our conclusions and mention some future prospects in section 8 .
the observations were carried out between october 31@xmath2 and november 4@xmath3 2004 with the mid - infrared interferometric instrument ( midi ) at eso s very large telescope interferometer ( vlti ) on paranal / chile .
together with the keck interferometer nuller that recently started to produce first scientific results @xcite , midi is currently worldwide the only instrument able to conduct spectrally resolved interferometric observations in the mid - infrared @xcite . for the observations midi was used in high - sens mode using the nacl prism as dispersive element yielding a spectral resolution of 30 .
the maximum projected baseline was 86.25 m and the minimum projected baseline 44.56 m leading to an angular resolution at 10@xmath0 m of 0.@xmath4029 and 0.@xmath4056 , respectively .
for an assumed distance of 450 pc these values correspond to 13.1 au and 25.2 au .
a journal of the observations including projected baselines , position angles and calibrator stars is given in table  [ journal ] .
for completeness we also mention the observations from december 2004 although they are disregarded in the following sections .
these observations had almost exactly the same baseline and position angle as the october observations but showed in general a lower level of total and correlated flux due to rather poor observing conditions in terms of seeing and transparency . the calibrator stars for the data reduction
were chosen by analyzing all calibrator stars observed over the whole night .
we selected those showing a good agreement in their transfer functions , i.e. their instrumental visibilities after the assumed sizes were taken into account .
the data reduction was carried out with the software package mia+ews-1.3 .
this software consists of two independent reduction programs ( mia and ews ) which were both applied for comparison . the software , further information and manuals can be downloaded from the internetnevec / midi / index.html ] . a general description of the basic data reduction steps is also given in @xcite and we refer the reader to this paper . as the results derived with mia and ews agreed quite well ( overall differences in the correlated flux @xmath5 ) , we decided to show only plots resulting from the mia package .
fu ori was found to be a binary system by @xcite .
@xcite confirmed the detected companion and concluded that it was a young star of spectral type k showing clear nir excess . in the same paper it was
also stated that it was very unlikely that the binary component ( fu ori s ) triggered the outburst of fu ori observed in the late 1930s .
knowing about the existence of the fainter companion the integration time of some midi acquisition images was increased in order to derive n - band photometry for both components . in three images
fu ori s was clearly visible ( figure  [ image ] ) and aperture photometry could be applied to the observations .
the results are summarized in table  [ fluxes ] .
interestingly , fu ori s shows relatively a higher n - band excess than fu ori itself .
this can be explained in terms of differences in the geometry of an assumed circumstellar disks ( e.g. , larger flaring angle or smaller disk inclination ) .
we furthermore derive a separation between the components of 0.@xmath4484@xmath60.@xmath401 and a position angle of 162.5@xmath74.1@xmath8 ( measured from north eastwards ) which is in good agreement with the results from @xcite and @xcite .
our errors are standard deviations based on gaussian fitting of 3 independent images . for an assumed distance of 450 pc
the separation corresponds to 217.8@xmath64.5 au .
figure  [ spectrum_short ] shows the mir spectrum of fu ori .
the spectrum obtained with midi is the average over the first three observing nights using both unit telescopes ( uts ) . for comparison _
spitzer irs _ data are plotted additionally .
those data were available from the _ spitzer _ data archive .
for the data reduction we used the _ droopres _ intermediate data product processed through the ssc pipeline s12.0.2 . the smart reduction package developed by the irs instrument team at cornell @xcite
was used to extract the spectrum . within
the error bars the spectra agree quite well .
a broad and weak silicate emission feature is present between 8 - 13@xmath0 m .
however , the averaged midi spectrum appears to be a little flatter than the more roundish _ spitzer _ spectrum .
in contrast to other ysos the observed silicate emission feature of fu ori is rather weak ( see also * ? ? ? * ) . such a broad and flat mir dust spectrum can be explained by grains that already underwent some coagulation process @xcite . to test this hypothesis we derived a possible dust composition by fitting a dust model to the normalized and continuum subtracted _ spitzer _ spectrum from 6 - 13@xmath0 m .
this analysis method is similar to approaches that were successfully applied in previous studies for determining the dust composition of observed dust emission features of other ysos ( e.g. , * ? ? ?
in the model we assume that the observed emission can be computed from the sum of the emission of individual dust species .
the species we used are summarized in table  [ dust_species ] together with references for their optical properties .
we furthermore took into account different grain sizes as for each of the dust species the opacities were calculated for 0.1 , 1.5 , and 6.0@xmath0 m sized particles . to estimate the possible contribution of polycyclic aromatic hydrocarbons ( pahs ) to the emission we used a template spectrum based on profiles that were derived from observations @xcite .
finally , in order to reproduce the observed spectrum we fitted the following emission model to the _ spitzer _ data using a least square minimization @xmath9 here , @xmath10 is the observed flux , @xmath11 denotes the planck function at the temperature of the underlying continuum , @xmath12 the planck function at the temperature of the silicate grains , @xmath13 is the mass absorption coefficient for the silicate species @xmath14 with the grain size @xmath15 , @xmath16 is the template emission spectrum for the pahs , and @xmath17 , @xmath18 , and @xmath19 are the weighting factors for the continuum , the silicate emission and the pah contribution , respectively .
a more detailed description of the dust model and plots showing the different dust opacities is given in @xcite .
the model results are shown in figure  [ dust_spectra ] and table  [ dust_compo ] summarizes the corresponding dust composition . as for the underlying continuum we could fit a temperature of 880 k , we believe that this optically thick emission most likely arises from the inner parts of the disk close to the star
. the outer and thus cooler regions of the disk would then be responsible for the weak silicate emission for which we fitted a temperature of 230 k to the optically thin emission layer .
the results strongly support the idea that most particles are amorphous and much larger in size than those in the ism .
pahs , silica particles and crystalline silicates do not seem to be present at all .
although the derived fractions of different dust species do depend on the model applied to the data ( e.g. , in terms of the precise grain sizes considered and also in terms of particle structures ) the main conclusions , that we do not find evidence for crystalline silicates and that the grains are considerably larger than those found in the ism , remain unaltered . for a more detailed discussion about the influence of different model parameters we refer to @xcite and @xcite .    as the disk is much brighter in the optical and nir than the central star ( section 7 ) ,
stellar radiation is probably not able to produce the observed weak emission feature in an optically thin disk surface layer .
however , it was shown that the silicate feature of fu ori can be reproduced taking into account self - irradiation of the disk @xcite . in this case
the hot inner parts of the accretion disk serve as flux source and illuminate the disk surface layer under a certain angle .
in figure  [ visibilities ] and table  [ visibilities_table ] the visibilities measured for fu ori at three different baselines are shown .
the errors are computed as standard deviations resulting from the use of different calibrators for each night ( see table  [ journal ] ) . in most of the following plots
the position of the atmospheric ozone band is indicated as this part of the spectrum suffered sometimes from imperfect corrections during the data reduction process . here ,
small `` dips '' and `` bumps '' in otherwise flat spectra are remnants from the data reduction .      the calibrated visibility of the shortest baseline ( ut2-ut3 ) remains almost constant at @xmath200.95 over the whole wavelength range and the source is only marginally resolved . as expected , for the longest baseline ( ut2-ut4 ) the lowest visibility
is observed ranging from @xmath200.73 at 8.3@xmath0 m to @xmath200.6 at 13@xmath0 m . for the intermediate baseline ( ut3-ut4 ) the visibility shows a slight increase from @xmath200.73 at 8.3@xmath0 m to @xmath200.85 at 13@xmath0 m .
the fact that the object is clearly resolved with two baselines in the mir supports the argument of @xcite that we are observing an extended circumstellar structure , i.e. a disk , and not only the stellar photosphere @xcite .
for the two baselines ut2-ut3 and ut2-ut4 the observations are consistent with expectations from thermal disk emission as despite the decreasing resolution for longer wavelengths the visibilities indicate larger sizes for the emitting regions ( see also section 5.3 . ) .
the observed increase in visibility for the ut3-ut4 baseline implies that for this baseline , however , the object appears smaller at 13@xmath0 m than at 8@xmath0 m .
this seems difficult to imagine in the context of a circumstellar disk as emitting source and certainly a second observation for this baseline configuration seems eligible . theoretically it is possible that the photometric measurements carried out directly after the interferometric observations are corrupted due to technical problems or different weather conditions .
this in turn might then lead to a change in the visibility function .
however , the calibrators observed at this night did not show any sign of poor photometric measurements and their transfer function was very stable over the whole night .
in addition , also other ysos with circumstellar disks showed an increasing visibility for certain baselines when observed with midi .
thus , since we do not find any evidence for excluding this dataset due to bad quality we decided to keep it in our analyses .      as thus far no midi visibilities for ttauri stars
have been published we are limited to a comparison between fu ori and circumstellar disks around herbig ae / be stars ( haebes ) .
it shows that these stars as well as models applied to them normally show a prominent drop in the visibility between 8@xmath0 m and 10@xmath0 m from where the curve remains almost constant @xcite .
qualitatively , this drop results from the intensity distribution of the passive disks around the haebes : at the short wavelength end the hot inner rim of the disks provides an overproportional contribution to the flux and is at the same time confined to a small spatial region leading to a high visibility .
most of the rest of the mir emission originates from a large area of the hot , illuminated surface layer of the flared circumstellar disk .
the visibilities for fu ori show hardly any wavelength dependence and are very flat from 8 - 13@xmath0 m regardless of baseline .
thus , its seems as if the flux distribution of fu ori is smoother and the visibilities show no sign of a significant contribution from a hot inner rim .
however , as fu ori is surrounded by a heavily active accretion disk where the disk alone produces the majority of the observed flux at almost all wavelengths ( see section 7 ) differences in the intensity distribution and thus in the visibility can be expected .
apart from differences in the shape of the visibility curves most haebes are much better resolved , i.e. show lower visibilities than fu ori @xcite .
this , however , can at least partly be explained by the distance to these objects which is in general less than 200 pc .
fu ori on the other hand has an assumed distance of 450 pc , and if it was closer to the earth we would observe lower visibilities also for this object .      to derive a simple model for the geometry of the emitting regions we assume a simple gaussian brightness distribution for each baseline .
this is a reasonable first approximation for objects showing high visibilities .
the fwhm , and hence the physical size , of this gaussian in arcsec can be computed by @xmath21 where @xmath22 is the measured visibility for a certain spatial frequency @xmath23 ( in arcsec@xmath24 ) derived from the projected baseline @xmath25 and the wavelength @xmath26 .
equation  [ visi_eq ] results from a simple fourier transformation of the assumed brightness distribution .
we computed the fwhm for all three baselines at three different wavelengths ( 9.0 , 11.0 , and 12.5@xmath0 m ) . for this we averaged for each baseline 5 visibility points from table  [ visibilities_table ] centered on the specified wavelengths .
table  [ table_sizes ] gives the resulting sizes of the emitting regions in au for an assumed distance to fu ori of 450 pc .
the results are also visualized in figure  [ figure_sizes ] where the fwhm are shown in their orientation on the sky .
as expected for thermal disk emission the fwhm increases with wavelength for a given baseline .
the only exception is the 12.5@xmath0 m size for the ut3-ut4 baseline which is surprisingly a little smaller than that seen at 11.0@xmath0 m .
in addition to size estimations for each individual baseline our measurements based on different position angles allow to constrain the geometry of the disk .
we fitted the derived fwhm with an ellipse for each of the considered wavelengths in order to derive a simple model for the spatial orientation of an assumed disk - like structure .
the resulting best fit ellipses are overplotted in figure  [ figure_sizes ] and their parameters are summarized in the lower half of table  [ table_sizes ] . as one would expect from a disk - like structure showing a decrease in temperature with radius the semimajor axes of the ellipses increase with wavelength . however , due to the increasing visibility of the ut3-ut4 baseline the semiminor axis and also the total area of the ellipse are slightly smaller at 12.5@xmath0 m than at 11.0@xmath0 m .
an inclination angle for an assumed circular disk can be derived by computing the @xmath27 of the ratio of the semiminor and semimajor axis .
the resulting angles agree very well ( 55.4@xmath72.4@xmath8 ) and are also in good agreement to what was found for fu ori disk models based on nir observations @xcite . by combining the derived inclination angle with rotational velocity measurements in the optical spectral region it is possible to estimate the central mass . based on @xcite
we derive for fu ori @xmath280.36 which is reasonably consistent with typical values for low - mass pre - main sequence stars .
in addition to the inclination the fitted ellipses provide also information on the position angle of the disk . the position angles ( measured from north eastwards )
we find for the semimajor axes at different wavelengths show a larger scatter ( 109.1@xmath711.6@xmath8 ) than the inclinations and differ significantly from the nir findings of @xcite ( see also section 7 ) .
it seems as if the position angle becomes smaller for longer wavelengths , although one has to keep in mind that , again , at least for the ut3-ut4 baseline the increasing visibility is the main reason for the observed rotation . at this point
we leave it to future investigations to derive a 3-d disk model based on a more complex dust distribution possibly able to confirm the apparent changes in the position angle for different wavelength regimes .
the correlated flux is directly linked to the total flux ( i.e. the flux from a single ut telescope ) via the visibility : @xmath29 figure  [ correlflux ] depicts the results for the correlated flux for the three different baselines .
the total spectrum is plotted for comparison .
taking into account the spatial resolution of the three baselines it is possible to estimate from where the correlated 8 - 13@xmath0 m flux originates .
our ut2-ut3 baseline has a spatial resolution of @xmath2025 au at the distance of fu ori , and if we assume an average visibility of 0.95 we hence know that 95% of the 8@xmath0m-13@xmath0 m flux must come from within this 25 au .
similarly , we find that for the ut2-ut4 baseline 65% of the flux is emitted from the inner 13 au . for the ut3-ut4 baseline 75% of the 8 - 13@xmath0 m flux comes from the inner 20 au .
the spatially resolved mir spectra of the midi instrument make it possible to study the radial dependence of the dust composition in the protoplanetary disk @xcite . in section 4
we showed that the spatially unresolved ( total ) spectrum of fu ori can be fitted with mainly large amorphous dust grains .
the shape of the spatially resolved ( correlated ) flux spectra contains information about the dust composition of the inner parts of the disk .    to compare the shape of the correlated spectra to the total spectrum we first subtracted the continuum for which we fitted a straight line between the flux at 8.25 and 12.95@xmath0 m .
then we normalized all continuum subtracted correlated flux spectra and the total _
spectrum by @xmath30    where @xmath31 and @xmath32 denote the continuum subtracted spectra and their mean value , respectively . by applying this formula all correlated spectra conserve their shape and can be compared to the total spectrum ( figure  [ shape_correl ] ) . to identify those parts
where the spectra differ significantly we furthermore computed the deviation from the total _
spectrum for each wavelength in units of the data error : @xmath33    @xmath34 and @xmath35 are the normalized correlated spectra and the normalized total spectrum and @xmath36 denotes the standard deviation of the correlated spectra ( midi errors ) at a given wavelength .
we applied the same normalization factors to the errors as we did to the spectra .
figure  [ shape_correl ] shows that no clear deviations from the total unresolved spectrum are observed for any of the correlated spectra .
this means that with the given spatial resolution and sensitivity no significant chemical difference in the dust composition is observed and the spectra of the inner parts of the disk look very similar to the total unresolved disk spectrum .
this finding is somewhat surprising .
as seen in section 4 the silicate dust grains apparently already underwent noticeable coagulation .
thus , the disk does not consist of purely ism dust anymore and must already have a significant age .
in addition it is known that apart from grain growth also thermal annealing takes place within protoplanetary disks and transform the amorphous silicates into crystalline silicates ( see , * ? ? ?
* ; * ? ? ?
however , neither the total spectrum nor the correlated spectra indicate the presence of crystalline silicates .
this is even more surprising considering the high temperatures in the disk interior due to viscous heating of this actively accreting star ( section 7 ) . annealing processes transforming ism amorphous dust into crystalline dust set in in disk regions where @xmath37 ( e.g. , * ? ? ? * )
in addition to that , high accretion rates directly relate to a faster vertical and horizontal mixing within the protoplanetary disk @xcite .
this effect also eases the detection of crystalline silicates within the accretion disk surrounding fu ori .
but why do nt we find any evidence for their existence ?
a possible explanation might be that they simply can not be detected as the related feature is rather weak and simply dominated by the underlying continuum .
this hypothesis can be checked by analyzing the longer wavelength regime of the _ spitzer _ spectrum as shown in figure  [ spectrum_long ] .
as at longer wavelengths the temperature of the continuum is decreasing the contrast between the continuum and possible emission features due to the presence of crystalline silicate particles increases .
hence , spectral signatures of these particles should be more easily detectable .
figure  [ spectrum_long ] , however , shows that also between 14 and 34@xmath0 m the spectrum of fu ori does not show any clear signs of crystalline silicates .
normally crystalline forsterite particles have emission bands at 16.4 , 23.9 , 27.7 and 33.8@xmath0 m and enstatite grains show features at 18.8 , 21.5 , and 24.5@xmath0 m @xcite .
if one argued that possibly around 27.7 - 27.8@xmath0 m a weak forsterite feature was present one should keep in mind that the features at 23.9 and 33.8@xmath0 m are normally stronger @xcite but are not seen at all .
thus , also the longer wavelength mir spectrum of fu ori shows that crystalline silicates are not present .
this shows that apparently the non - detection is at least not only due to contrast problems between the emission features and the underlying continuum .    what else can prevent us from detecting crystalline silicates ? in case
fu ori underwent its first ( of possible multiple ) outburst(s ) and phase(s ) of enhanced accretion , it is possible that not enough time has yet elapsed for the crystalline silicate particles to reach the disk surface from where they can be detected by means of mir spectroscopy .
this hypothesis is , however , questionable as the large sizes of the amorphous silicates indicate that enough time , at least for dust coagulation , has elapsed .
recently , @xcite found that assuming the crystallization by thermal annealing happens at the very early phases of disk evolution , the level of crystallinity is linked to the properties of the molecular cloud core from which the disk formed .
it was found that rapidly rotating cores produce rather massive disks with a high accretion rate but a low level of crystallinity . as we observe these properties for fu ori the proposed model ,
although yet very simple and with limited predictive power , might provide an interesting way to explain at least partly the lack of crystalline silicates in the fu ori accretion disk . finally , another possibility might be that the vertical and radial mixing within the disk does not work as efficient as expected . in consequence
not enough crystalline dust grains reach the disk surface at radii where also the contrast between the emission feature and the dust continuum allows their detection .      the results from the correlated flux allow a comparison to disks around other ysos that have also been observed with vlti / midi .
not only seem crystalline silicates be commonly present in protoplanetary disks , but for three intermediate - mass haebes stars @xcite found that the degree of crystallization within the disk depends on the distance from the central star .
the inner parts of the disks are more crystalline than the outer disk regions .
this can naturally be explained by assuming that in the inner disk regions where the disk has higher temperatures dust processing and annealing occurs and that over time crystalline particles are transported outwards .
not only for intermediate mass stars but also for the ttauri star ry tau the same radial dependence of the crystallinity was observed @xcite .
thus , crystalline silicates should in general be more easily detectable in the correlated spectra of protoplanetary disks . for reasons mentioned in the previous section it seems thus surprising that the disk around fu orionis does not seem to contain any detectable amount of crystalline dust .    as described in the introduction fuors
are believed to be low - mass or ttauri - like stars in a state of enhanced accretion . from sub - millimeter and millimeter observations
@xcite concluded that fuors are possibly in an evolutionary youger age than classical ttauri stars .
this younger age is also supported by theoretical models where the protostellar disk is continuously replenished by a dusty envelope .
this material infall on the disk can lead to thermal instabilities in the disk and thus provides a possible mechanism for explaining the observed outbursts @xcite .
the puzzling thing is now that apparently in the disk surrounding fu ori the amorphous grains have already grown considerably .
hence , the disk can not be in its native dust composition and might already have a significant age .
furthermore , from the _ spitzer _ feps ( _ formation and evolution of planetary systems _ ) legacy program it is known that apparently the amount of crystalline silicates in ttauri disks is not correlated with the average size of amorphous dust particles @xcite . hence , crystalline dust particles are created in the _ early _ phases of disk evolution where the accretion rate and the disk temperature are thought to be higher .
if thus the fuor - phase is indeed common to most / all ttauri stars then the non - detection of crystalline particles can not be explained in the current picture of protoplanetary disk evolution .
the present paper is more intended to show the observational results rather than to present a sophisticated disk model possibly able to match the observed sed and visibilities with highest accuracy .
thus , instead of applying a complex radiative transfer model we decided to use in a first step a simple analytical disk model . in this model the effective temperature and also the surface density are prescribed by the following broken power law distributions    @xmath38    @xmath39    for the temperature law the temperature at 1 au and the exponent @xmath40 are free parameters . for the surface density
the total disk mass @xmath41 and the exponent @xmath42 are free parameters .
as the integrated surface density must equal the total mass of the disk the radial dependence of the surface density @xmath43 can be computed .
the shape of the total sed ( and finally also the observed visibilities ) require a double power - law to fit the data . at a radius of 3
au the temperature and surface density distribution change . however , we made sure that for both parameters the transition was smooth and continuous .
figure  [ sed ] shows the sed and the resulting fit from the simple disk model .
the model parameters are summarized in table  [ model_parameter ] .
based on the value for the visual extinction @xmath44 we interpolated the extinction for the other wavelengths following the findings of @xcite .
the disk inclination @xmath15 was derived from the ellipses fitted to the fwhm of assumed gaussian brightness distributions ( see section 5.3 . ) .
for the position angle @xmath45 we used two values : one was also derived in section 5.3 . , and
for the second one we fitted the disk model to the observed visibilities with @xmath45 being the only free parameter ( see also section 7.2 . ) .
it is clear that this simple approach is not able to reproduce all features of the sed ( e.g. , the weak silicate emission feature ) . for the moment , however , we restrict ourselves to this simple approach and , as mentioned above , a detailed modeling is left for further investigations .
it is striking that rather high disk temperatures are required to fit the sed . in the inner 0.35 au
no dust can withstand the high temperatures @xmath461500 k even on the disk surface and only the gaseous disk component can survive .
it is interesting to note that the power law index of @xmath47 of the inner disk region equals exactly that what was derived in analytical disk models for geometrically flat , optically thick _ passive _ disks ( e.g. , * ? ? ?
* ) as well as for flat , steady state , optically thick _ accretion _ disks ( e.g. , * ? ? ?
thus , from the temperature distribution alone one can not distinguish between the two different cases ( see also * ? ? ?
. however , in the context of fu ori the passive disk model where the disk is merely re - processing radiation from the central star is unlikely and only an accretion disk model matches most observations @xcite .
a nice example is the recently detected magnetic field in the innermost regions of the disk surrounding fu ori which can only be explained in the picture of a heavily accreting disk @xcite . as the disk temperature
is directly related to the mass accretion rate @xmath48 in classical accretion disk models it is possible to derive the @xmath48 required to produce the observed high temperatures .
@xcite analyzed the effective disk temperature and also the midplane temperature of an accretion disk as a function of the mass accretion rate .
figure 1 of their paper shows that one needs @xmath49 to get close to the effective temperature we find in our model .
such high values for @xmath48 were also found by @xcite and @xcite . due to the high accretion rate the midplane temperature of the disk in the inner @xmath2010 au
is even much higher than the effective disk temperature @xcite .
concerning the temperature law farther out than 3 au the value we adopt here for @xmath50 is in good agreement to what can be found in the literature for isothermal flared disks @xcite . in their model
the index value diverges from the flat solution where @xmath51 at a certain point in the inner disk and approaches @xmath52 for radii much larger than the stellar radius . at intermediate radii
the temperature index might be approximated by values in between these most extreme cases depending on the actual flaring of the disk .    in summary
, the temperature distribution indicates that in the innermost disk regions the disk luminosity is accretion dominated whereas further out it is irradiation dominated .
this is expected for an active accretion disk since the luminosity released by accretion is a much steeper function of distance to the star than that by absorption of radiation from the central star and the hot innermost disk regions ( @xmath53 , whereas @xmath54 , roughly ) .        for this simple disk model
it is straightforward to compute the visibilities as a function of projected baseline and wavelength .
assuming that each part of the disk surface is emitting as a blackbody with an effective temperature @xmath55 the corresponding total flux and visibility for an axisymmetric face - on disk is given by ( see also equations ( 2 ) , ( 3 ) , and following in * ? ? ?
* ) : @xmath56    @xmath57dr\ ] ]    with @xmath58 being the flux , @xmath59 the planck function , @xmath60 the zeroth - order bessel function of the first kind , @xmath61 the projected baseline and @xmath62 the distance to the object .    taking into account the inclination of the disk and the position angle
as given in table  [ model_parameter ] we computed the corresponding visibilities for the three different baselines .
the results for two different position angles are shown in figure  [ model_visi ] . in the upper panel the disk inclination and the position angle are derived from the ellipses fitting in section 5.3 . in the lower panel
the position angle was fitted to the observed visibilities with a least squares fit .
while the ellipse fitting yielded a position angle of 109.1@xmath6311.6@xmath64 a least square analysis leads to a value of 93.4@xmath636.8@xmath64 . from the plots in figure  [ model_visi ]
one can see that changing the position angle mostly affects the visibility of the shortest baseline while the curves for the other baselines change only marginally .
fitting the ut3-ut4 baseline proofs quite a challenge as the observed visibility is increasing over the 8 - 13@xmath0 m wavelength interval while all other modeled visibilities are decreasing . especially for this baseline
better fits are expected from more sophisticated disk models . however , given the simplicity of the current models the fits seems in reasonable agreement with the observations .
it clearly shows that interferometric measurements can put additional constraints on the structure and geometry of circumstellar disks that can not be derived from sed fitting alone .
interestingly , @xcite were able to fit their nir visibilities of fu ori with a similar disk model and also with a disk and an embedded `` hot spot '' .
they concluded that the latter model was more likely based on statistical arguments .
the major differences between their models and ours is that ( 1 ) they fit the sed with single power laws for the temperature neglecting the flux longwards of 20@xmath0 m and ( 2 ) their best fits yield disk position angles of 47@xmath65 and 8@xmath66 .
the inclination angles ( 55@xmath67 and 48@xmath68 , respectively ) , are in good agreement with our value . in figure  [ malbetmodels ]
we compare both of their models to our observations .
we find that both models fit the midi visibilities not as good as the model presented here .
even more we find that our mir interferometric measurements rule out their first model of a simple accretion disk .
this model predicts a higher visibility for the ut3-ut4 baseline than for the ut2-ut3 baseline which is not observed .
the model with an embedded `` hot spot '' in the accretion disk shows the right trend for the visibilities , but the values for the ut3-ut4 and ut2-ut4 baseline are higher than for our models which already give slightly too high values compared to the observations .
the main reason for the higher visibilities resulting from the model applied by @xcite is the single power law approach which provides less mir flux compared to our models .
we also compared our models to the measured nir visibilities from @xcite .
we restricted ourselves to those observations where the object was clearly resolved .
the results are given in table  [ table_nirvisi ] .
except for the observations done with north - west baseline of the palomar testbed interferometer ( pti / nw ) both of our models agree with the nir visibilities within the error bars .
concerning the predicted `` hot spot '' in the accretion disk @xcite our observations do not yield new insights into its nature as they do not have the required spatial resolution to confirm its existence .
also the calibrated phase of our observations do not contain any information on which basis we could speculate about this possible second companion .
we presented the first multi - baseline mir interferometric observations of fu orionis with midi / vlti .
the findings can be summarized as follows :    1 .
fu orionis was clearly resolved in the mir with two vlti baselines and marginally resolved with a third baseline indicating the presence of warm dusty material surrounding fu ori out to several tens of au .
the inclination and the position angle of the accretion disk can be inferred from multi - baseline measurements with midi / vlti .
thus , the instrument provides means for deriving disk parameters that are otherwise only poorly constrained from models and sed fits .
the correlated flux indicates that 95% of the 8 - 13@xmath0 m flux comes from within the inner 25au of the disk for the shortest of our baselines while 65% of this flux arise from the inner 13au for our longest baseline .
the shapes and strengths of the total 8 - 13@xmath0 m spectrum and the ( spatially resolved ) correlated spectra indicate that most dust particles within the accretion disk are amorphous and already significantly larger than typical particles in the ism .
no spectra , neither the total nor the spatially resolved correlated spectra bear significant traces of crystalline silicates . given the high accretion rate and disk temperature of the system this is unexpected and requires further investigations .
the sed and the observed mir visibilities can be fitted reasonably well with a simple analytical disk model prescribing a broken power law for the effective temperature distribution of the protoplanetary disk . within the innermost 3au
the temperature decreases with @xmath69 , farther out the temperature goes with @xmath70 .
the derived power law index for the inner regions is equal to what is derived theoretically for flat , steady state , optically thick accretion disks @xcite . for the outer disk region
the derived value is in good agreement to what can be found for isothermal flared disks @xcite .
7 .   based on points 4 . and 5 . and recent spectroscopic observations of protoplanetary disks around ttauri stars and haebes it seems questionable that fu ori is a very young classical ttauri star as the size of amorphous dust grains in combination with the non - detection of crystalline dust contradicts this assumption .
8 .   from the acquisition image n - band aperture photometry could be derived for the companion fu ori s. it shows that this young k - type star has more relative mir excess than fu ori itself hinting at a circumstellar disk with larger flaring angle or smaller disk inclination .      1 .
most of the published nir visibilities of fu ori @xcite do agree with our newly derived disk models .
2 .   our interferometric observations
rule out one out of two disk models for fu ori presented by @xcite . due to the lack of sufficient spatial resolution the presence of the proposed second very close companion
can neither be confirmed nor disproved with the current set of midi data .
3 .   in contrast to nir interferometric observations of other fuors @xcite our data do not seem to require the introduction of an extended dusty structure providing additional uncorrelated flux .
4 .   for the inner 3au
we require a significantly different disk model ( e.g. , @xmath71 , @xmath69 ) to fit our visibilities and sed in comparison to what was found for the outbursting star v1647 ori ( @xmath72 , @xmath70 ) @xcite .
although the present results provide important new insights into the fu ori system new questions arose as well .
future investigations will include a more thorough and detailed modeling of the observations with numerical radiative transfer disk models .
these models will not only need to consider the high accretion rates but also such effects as self - irradiation of the disk .
it will be interesting to see whether these models are able to confirm the apparent rotation of the disk s position angle with wavelength and whether they can reproduce the increasing visibility curve for the ut3-ut4 baseline .
also , these disk models will have to reproduce the observed nir and mir visibilities simultaneously with higher accuracy than the current simple models are able to do . from the observational side re - measuring
the ut3-ut4 baseline might be eligible and any other new baseline configuration will certainly help to constrain the disk geometry ( especially the position angle ) even better . concerning the dust structure
the lack of crystalline silicates needs to be analyzed in greater detail .
is it still just a contrast effect that we do not see any crystalline features in the spectra or are they simply not there ?
eventually we will have to face the problem that currently all fuors that have been observed with high spatial resolution techniques draw a rather inhomogeneous picture of the group .
apparently , some are still surrounded by nearby dense envelopes while for others the circumstellar disks suffice to explain the observations . however , even those disk - systems ( i.e. , fu ori and v1647 ori ) seem to be different in terms of physical properties as different temperature and surface density profiles are derived .    without any doubt interferometric observations of circumstellar disks in the nir and mir provide an unprecedented means for deriving geometrical disk properties and put new constraints on the physical ( and also chemical ) processes taking place within protoplanetary disks . and as more and more interferometric data of circumstellar disks are published our view on the cradle of planetary systems will be once more refined .    s. p.
quanz kindly acknowledges financial support by the german _ friedrich - ebert - stiftung_. th .
henning and j. bouwman acknowledge financial support by the ec - rtn on `` the formation and evolution of young stellar clusters '' ( hprn - ct-2000 - 00155 ) .
we are most grateful to r. khler and f. lahuis for their support during the data reduction and to r. lachaume and f. malbet for providing mir data derived from their disk models ( section 7.2 . ) .
we thank our anonymous referee for a detailed and thorough review that improved the manuscript .
this research has made use of the simbad database , operated at cds , strasbourg , france , and was partly based on observations made with the _ spitzer space telescope _ , which is operated by the jet propulsion laboratory , california institute of technology , under nasa contract 1407 .
lllll 31.10.2004 & fu ori & 44.56 m ( ut2-ut3 ) & 46.54@xmath8 & - + & hd37160 & & & calibrator + 02.11.2004 & fu ori & 86.25 m ( ut2-ut4 ) & 84.24@xmath8 & - + & hd 31421 & & & calibrator + & hd 37160 & & & calibrator + & hd 50778 & & & calibrator + 04.11.2004 & fu ori & 56.74 m ( ut3-ut4 ) & 106.64@xmath8 & - + & hd 25604 & & & calibrator + & hd 20644 & & & calibrator + & hd 37160 & & & calibrator + & hd 50778 & & & calibrator + 29.12.2004 & fu ori & 44.80 m ( ut2-ut3 ) & 46.61@xmath8 & - + & hd 37160 & & & calibrator + & hd 94510 & & & calibrator +    llr j [ mag ] & 6.30@xmath73 & 10.75@xmath74 + h [ mag ] & 5.64@xmath75 & 9.92@xmath76 + k [ mag ] & 5.25@xmath77 & 9.15@xmath78 + l [ mag ] & 4.18@xmath79 & 8.09@xmath80 + n [ mag ] & 2.75@xmath81 & 5.28@xmath82 +    llll amorphous olivine & mgfesio@xmath83 & homogeneous & ( 1 ) + amorphous pyroxene & mgfesi@xmath84o@xmath85 & homogeneous & ( 1 ) + crystalline forsterite & mg@xmath84sio@xmath83 & inhomogeneous & ( 2 ) + crystalline enstatite & mgsio@xmath86 & inhomogeneous & ( 3 ) + amorphous silica & sio@xmath84 & inhomogeneous & ( 4 ) +    lllll amorphous olivine & 0.17 & @xmath87  0.01 & 0.18 & 0.35 + amorphous pyroxene & @xmath87  0.01 & 0.07 & 0.57 & 0.64 + crystalline silicates & @xmath87  0.01 & @xmath87  0.01 & @xmath87  0.01 & @xmath87  0.01 +    cccc 12.95 & 0.98@xmath60.02 & 0.60@xmath60.02&0.84@xmath60.02 + 12.83 & 0.93@xmath60.02 & 0.61@xmath60.03&0.83@xmath60.02 + 12.70 & 0.94@xmath60.03 & 0.63@xmath60.03&0.84@xmath60.02 + 12.57 & 0.96@xmath60.02 & 0.62@xmath60.02&0.84@xmath60.02 + 12.44 & 0.96@xmath60.02 & 0.62@xmath60.02&0.81@xmath60.02 + 12.31 & 0.94@xmath60.02 & 0.60@xmath60.01&0.81@xmath60.02 + 12.17 & 0.92@xmath60.01 & 0.61@xmath60.02&0.79@xmath60.02 + 12.03 & 0.94@xmath60.02 & 0.61@xmath60.02&0.81@xmath60.02 + 11.89 & 0.93@xmath60.03 & 0.62@xmath60.02&0.79@xmath60.02 + 11.75 & 0.94@xmath60.03 & 0.62@xmath60.01&0.78@xmath60.02 + 11.60 & 0.94@xmath60.03 & 0.63@xmath60.01&0.77@xmath60.02 + 11.45 & 0.91@xmath60.02 & 0.63@xmath60.01&0.79@xmath60.02 + 11.30 & 0.94@xmath60.03 & 0.62@xmath60.01&0.78@xmath60.02 + 11.14 & 0.92@xmath60.03 & 0.62@xmath60.01&0.76@xmath60.02 + 10.98 & 0.93@xmath60.02 & 0.62@xmath60.01&0.77@xmath60.02 + 10.82 & 0.92@xmath60.02 & 0.62@xmath60.01&0.76@xmath60.01 + 10.66 & 0.92@xmath60.02 & 0.62@xmath60.01&0.75@xmath60.01 + 10.49 & 0.91@xmath60.01 & 0.62@xmath60.01&0.75@xmath60.02 + 10.32 & 0.92@xmath60.02 & 0.63@xmath60.01&0.75@xmath60.02 + 10.15 & 0.92@xmath60.01 & 0.63@xmath60.01&0.75@xmath60.02 + 9.97 & 0.94@xmath60.02 & 0.65@xmath60.01&0.73@xmath60.01 + 9.80 & 0.97@xmath60.02 & 0.64@xmath60.01&0.73@xmath60.01 + 9.62 & 0.97@xmath60.03 & 0.64@xmath60.01&0.74@xmath60.01 + 9.43 & 0.93@xmath60.02 & 0.61@xmath60.01&0.76@xmath60.02 + 9.25 & 0.93@xmath60.01 & 0.67@xmath60.01&0.76@xmath60.02 + 9.06 & 0.93@xmath60.01 & 0.69@xmath60.01&0.74@xmath60.02 + 8.86 & 0.94@xmath60.01 & 0.69@xmath60.01&0.74@xmath60.02 + 8.67 & 0.95@xmath60.01 & 0.70@xmath60.01&0.73@xmath60.02 + 8.47 & 0.95@xmath60.01 & 0.70@xmath60.01&0.74@xmath60.02 + 8.27 & 0.94@xmath60.01 & 0.73@xmath60.01&0.73@xmath60.01 +    llll ut2-ut3 ( 44.56 m ) & 2.47@xmath88 & 3.27@xmath89 & 3.43@xmath90 + ut2-ut4 ( 86.25 m ) & 3.25@xmath91 & 4.34@xmath92 & 4.93@xmath93 + ut3-ut4 ( 56.74 m ) & 4.19@xmath94 & 5.00@xmath95 & 4.83@xmath96 + semimajor axis @xmath97 [ au ] & 2.28 & 2.51 & 2.54 + semiminor axis @xmath98 [ au ] & 1.21 & 1.54 & 1.42 + ellipse area [ au@xmath99 & 8.67 & 12.14 & 11.33 + position angle @xmath45 & 122.2@xmath64 & 111.1@xmath64 & 93.9@xmath64 + inclination @xmath15 & 58.0@xmath64 & 52.2@xmath64 & 56.0@xmath64 +    lll inner disk radius & @xmath100 & @xmath101 + outer disk radius & @xmath102 & 100 au + disk inclination  & @xmath15 & @xmath103 + disk position angle  & @xmath45 & @xmath104  /  @xmath105   + extinction & a@xmath106 & @xmath107 + temperature at 1 au & @xmath108 & @xmath109 + power law index for temperature & @xmath110 & @xmath111 + power law index for surface density & @xmath112 & @xmath113 + disk mass & @xmath114 & @xmath115 +      lcllcc pti / ns & h & 103.7 & 61.8 & 0.83@xmath60.04 & 0.80 / 0.82 + pti / nw & h & 85.6 & 14.1 & 0.79@xmath60.05 & 0.92 / 0.94 + pti / ns & k & 102.7 & 63.6 & 0.72@xmath60.08 & 0.78 / 0.80 + pti / nw & k & 84.2 & 14.0 & 0.79@xmath60.05 & 0.91 / 0.93 + pti / sw & k & 82.5 & -55.6 & 0.82@xmath60.08 & 0.85 / 0.87 + vlti / ut1-ut3 & k & 89.6 & 54.2 & 0.87@xmath60.05 & 0.83 / 0.85 +