compact radio sources in nearby starburst galaxies often appear to be supernova remnants ( snr ) .
snr in our galaxy and the large magellanic cloud ( lmc ) are sources of bright near - infrared [ feii ] line emission with luminosities ranging from @xmath5w to @xmath6w ( for the [ feii]@xmath0 m line , e.g. , oliva , moorwood , & danziger 1989 ; burton & spyromilio 1993 ) .
this behavior suggests a connection between supernova activity and the [ feii ] emission in starbursts .
this is supported by near - infrared spectroscopy ( moorwood & oliva 1988 ; greenhouse et al .
1991 ; vanzi & rieke 1997 ) , by high - resolution [ feii ] imaging studies ( forbes et al .
1993 ; greenhouse et al . 1997 ; alonso - herrero et al .
2000 ) , and by the tight correlation between the radio 6 cm and the [ feii ] emissions ( forbes & ward 1993 ) .
further determination of the relation of the infrared [ feii ] emission with the supernova activity in nearby starforming galaxies would be helpful both for the understanding of these galaxies and for using the near - infrared [ feii ] luminosity as estimator of the supernova rate in more distant galaxies where individual radio supernovae can not be identified .
there are indications of some complexities in this relation .
for example , the bright [ feii ] emission may not be spatially coincident with the radio snr .
lumsden & puxley ( 1995 ) conducted a spectroscopic survey for [ feii ] emission from radio snrs in m33 and found that in some cases the [ feii ] emission does not peak at the position of the radio snr , as well as that some radio snr show no [ feii ] emission . recently ,
morel , doyon , & st - louis ( 2002 ) have carried out an [ feii ] imaging survey of optically selected snr in m33 and found a similar result .
both works have concluded that the [ feii ] emission from snr may be produced over a time scale of @xmath7 10@xmath8 years .
greenhouse et al .
( 1997 ) found that none of the brightest sources of [ feii ] emission in m82 were coincident with the position of the radio snrs .
m82 and ngc  253 are the nearest starburst galaxies ( distances @xmath9mpc and @xmath10mpc respectively ) .
both galaxies show a large number of radio compact sources ( e.g. , @xmath11 in m82 kronberg et al .
1985 and huang et al . 1994 ; @xmath12 in ngc  253 turner & ho 1985 and ulvestad & antonucci 1991 , 1997 and references therein ) .
their proximity makes these two galaxies ideal to study the relation between supernova activity and the [ feii ] emission in starburst galaxies .
greenhouse et al .
( 1997 ) presented a fabry - perot [ feii]@xmath0 m map for m82 with a spatial resolution of @xmath13 .
their main conclusions were that the compact [ feii ] sources contribute only @xmath14 of the total [ feii ] emission in m82 and that they trace a population of snr older than the radio snr .
their work is somewhat inconclusive because of the limited angular resolution and also because strong extinction could hide [ feii ] emission from some of the radio snrs .
this paper presents _
hst_/nicmos images of m82 and ngc  253 .
the high spatial resolution of the nicmos [ feii]@xmath0 m images ( @xmath15 ; this corresponds to 5pc and 4pc for m82 and ngc  253 , respectively ) allows a detailed study of the location of the bright sources of [ feii ] and of the rle of supernova remnants in this line emission .
_ hst_/nicmos observations were obtained on april 7 , 8 , 12 and 22 1998 ( for m82 ) and on august 13 and 20 1998 ( for ngc  253 ) using the nic2 and nic3 cameras .
table  1 lists details of the observations .
the observational strategy consisted of taking images in a spiral dither with a 5.5 pixel spacing with four positions for a given field .
the plate scales for the nic2 and nic3 cameras are 0.076pixel@xmath16 and 0.20pixel@xmath16 , respectively .
the fwhm of the point sources in the nic3 f166n continuum images is @xmath17 .
lcclccc nic2 & f160w & @xmath18 & m82 & 192 & 4 & @xmath19 + & & & ngc  253 & 192 & 2 & @xmath20 + nic2 & f222 m & @xmath21 & m82 & 448 & 4 & @xmath19 + & & & ngc  253 & 512 & 2 & @xmath20 + nic2 & f187n & pa@xmath1 & m82 & 384 & 4 & @xmath19 + & & & ngc  253 & 384 & 2 & @xmath20 + nic2 & f190n & continuum & m82 & 384 & 4 & @xmath19 + & & & ngc  253 & 384 & 2 & @xmath20 + nic3 & f164n & [ feii]@xmath0 m & m82 & 4096 & 2 & @xmath22 + & & & ngc  253 & 2560 & 1 & @xmath23 + nic3 & f166n & continuum & m82 & 3840 & 2 & @xmath22 + & & & ngc  253 & 2560 & 1 & @xmath23 + nic3 & f215n & continuum & m82 & 2304 & 2 & @xmath22 + & & & ngc  253 & 2304 & 1 & @xmath23 +    the reduction of the nicmos images used routines from the package nicred ( mcleod 1997 ) .
a master dark image was produced by combining between 10 and 20 darks for a given sample sequence after the subtraction of the first readout .
the darks used were part of other proposals observed close in time , and the flatfield images were on - orbit data .
the data reduction was performed using the following steps : subtraction of the first readout , dark current subtraction on a readout basis , correction for linearity and cosmic ray rejection ( using fullfit ) , and flatfielding .
the individual dithered galaxy images were registered to a common position using fractional pixel offsets and cubic spline interpolation , and combined to produce the final image of each field .
prior to the flux calibration of the images , the background needs to be subtracted ( only for wavelengths longwards of @xmath24 m ) . due to the large projected size of both galaxies ,
the field of view of the nic2 images is not large enough to enable measurements of the background on empty corners of the images , so we used background measurements taken during the servicing mission observatory verification ( smov ) program for the filter nic2 f222 m .
the f215n images cover a larger field of view , and an estimate of the thermal background was obtained from the corners of the images .
the flux calibration of the broad - band , on - line and off - line images was performed using the conversion factors ( from adu / s to mjy ) based on measurements of the standard star p330-e during smov ( marcia rieke 1999 private communication ) .    ideally one would like to measure the continuum on both sides of the emission line to estimate the continuum at the emission line .
however , nicmos only provides continuum bands to the red of the emission line . in cases of objects such as m82 and ngc  253 , where the extinction to the nuclear regions is very high ( even at infrared wavelengths ) , a straight subtraction of the longer wavelength off - line image may result in an over correction of the continuum at the wavelength of the emission line . to estimate the continuum for the nic3 f164n images , we fitted the continuum on a pixel - by - pixel basis between @xmath25 m and @xmath26 m by linear regression of the flux as a function of the wavelength for the nic3 f166n and nic3 f215n line - free images .
the extrapolated continuum was then subtracted from the nic3 f164n image to produce the final continuum - subtracted [ feii]@xmath0 m line image .
this procedure does not correct for extinction , but it rather ensures the subtraction of the correct _ obscured _ continuum at each position of the galaxy . for the nic2 f187n image ( pa@xmath1 ) ,
since we did not have other narrow - band images taken with camera 2 , we used the nic2 f190n continuum image .
fortunately , the continuum need not be subtracted as accurately for pa@xmath1 as for the [ feii ] emission line because the pa@xmath1 emission is much brighter .
the fully - reduced images of each individual field were combined to produce the mosaics .
the images presented in this paper were taken during different orbits .
this resulted in slightly different orientations ( from @xmath27 to @xmath28 for m82 and from @xmath29 to @xmath30 for ngc  253 ) . before combining the images into mosaics , they were rotated by linear interpolation to north up , east to the left .
this provided the same orientation for all the images , and more importantly allows comparisons with previously published data .    to construct extinction maps to correct the [ feii ] emission , and [ feii]/pa@xmath1 maps ,
the nic2 f187n , f190n , f160w and f222 m images were demagnified to match the plate scale of the nic3 images and filtered with a gaussian with @xmath31 to match the point spread function ( psf ) of the nic3 images .
figure  1a and 1b are mosaics ( after rotation ) of the continuum at @xmath25 m and the continuum - subtracted [ feii]@xmath0 m line emission , respectively .
both images are displayed on a logarithmic scale to stress the contrast between the diffuse extended continuum and the compact sources .
the [ feii]@xmath0 m line map has been smoothed with a gaussian filter with @xmath32 .
the line emission image shows some residuals which are due to both the under - sampling of the nic3 camera psf and the variation in the shape of the psf with wavelength .
in addition , the right hand side of the line image shows a slightly lower s / n because the f164n and f166n images at that position had different orientations .
the [ feii]@xmath0 m line emission exhibits two distinct components : a number of relatively compact regions distributed along the disk of the galaxy and an extended diffuse component in the form of filaments extending above and below the plane of the galaxy .
the filaments seen in the [ feii]@xmath0 m emission line show a remarkable similarity with the chimneys at the base of the superwind detected in radio ( wills et al .
1999 ) , which are interpreted by these authors as expelled material from the central region of m82 .    in figures  2a and 2b we present contour maps of the @xmath25 m continuum emission and the [ feii]@xmath0 m line emission , respectively . comparison of the nicmos [ feii]@xmath0 m line map with that of greenhouse et al .
( 1997 ) reveals that some of the sources seen in the lower resolution map of greenhouse et al .
( 1997 ) break up into individual sources at higher spatial resolution .
in particular , their source fe1 contains some 4 or 5 individual sources .
source fe6 ( the most south - west source in greenhouse et al .
1997 map ) appears as a point source at our resolution in both the nic3 f164n and nic3 f166n images , and the signal to noise on it is compromised by the differing image orientations ( see above ) .
therefore it is not surprising we do not detect line emission from it .
the nicmos images detect a number of compact sources in addition to those reported in greenhouse et al .
( 1997 ) .
in figures  2a and 2b we display the positions of 44 radio snrs with peak flux densities above 1mjy ( list provided by dr .
karen wills , see more details in section  5.1 ) .
the distribution of the radio snrs and their relation with the [ feii]@xmath0 m emission will be discussed in section  5 .    in figures  3a and 3b
we display the mosaics of the @xmath33 m continuum ( nic2 f190n images ) and the pa@xmath1 line emission ( continuum subtracted nic2 f187n images ) respectively .
the pa@xmath1 image is dominated by a ring of ionized gas with an approximate diameter of 400pc .
the resemblance with the mid - infrared [ neii ] line emission map at @xmath34 m ( achtermann & lacy 1995 ) is remarkable .
the south side of this ring of ionized material appears to be incomplete ( see achtermann & lacy 1995 ) , and it also shows a `` hole '' of emission at the center of the galaxy , which approximately coincides with the peak of the @xmath35 m surface brightness .
ccccccc     + & ( ) & ( w m@xmath36 ) & ( w m@xmath36 ) & ( mag ) + pa@xmath1/h@xmath1 & 5.8 & @xmath37 & @xmath38 & @xmath39 & this work , oconnell & mangano ( 1978 ) + pa@xmath1/pa@xmath40 & 6 & @xmath37 & @xmath41 & @xmath42 & this work , mcleod et al . ( 1993 ) + pa@xmath1/pa@xmath43 & 3 & @xmath44 & @xmath45 & @xmath46 & this work , mcleod et al .
( 1993 ) + pa@xmath1/pa@xmath40 & 3 & @xmath44 & @xmath47 & @xmath42 & this work , mcleod et al .
( 1993 ) + br@xmath43/pa@xmath1 & 3.8 & @xmath48 & @xmath49 & @xmath50 & lester et al .
( 1991 ) , this work + br@xmath1/pa@xmath1 ( e1 ) & 6 ? & @xmath51 & @xmath52 & @xmath53 & achtermann & lacy ( 1995 ) , this work + br@xmath1/pa@xmath1 ( w1 ) & 6 ? & @xmath54 & @xmath55 & @xmath56 & achtermann & lacy ( 1995 ) , this work + br@xmath1/pa@xmath1 ( w2 ) & 6 ? & @xmath57 & @xmath58 & @xmath59 & achtermann & lacy ( 1995 ) , this work + @xmath60feii]1.64/[feii]1.25 & 3 & @xmath61 & @xmath62 & @xmath46 & this work , mcleod et al .
( 1993 ) +   + & ( ) & ( mjy ) & ( mjy ) & ( mag ) + f222m / f160w & 3 & 237 & 166 & @xmath56 & this work +    ccccccc + & ( ) & ( w m@xmath36 ) & ( w m@xmath36 ) & ( mag ) + pa@xmath1/br@xmath43 ( center ) & @xmath63 & @xmath64 & @xmath65 & @xmath66 & this work , harrison et al .
( 1998 ) + pa@xmath1/br@xmath43 ( center ) & @xmath67 & @xmath68 & @xmath69 & @xmath70 & this work , puxley & brand ( 1995 ) + pa@xmath1/br@xmath43 ( 3  ne ) & @xmath63 & @xmath71 & @xmath72 & @xmath73 & this work , harrison et al .
( 1998 ) + pa@xmath1/br@xmath43 ( 3.9  ne ) & @xmath67 & @xmath74 & @xmath75 & @xmath76 & this work , puxley & brand ( 1995 ) +   + & ( ) & ( mjy ) & ( mjy ) & ( mag ) + f222m / f160w ( center ) & @xmath67 & 213 & 142 & @xmath56 & this work +    it has been suggested that higher extinction in the south half of m82 may be causing the uneven distribution of ionized gas in the ring .
however , the ring appears to be complete in [ feii]@xmath0 m line emission .
if extinction were the cause , we should expect to see a similar morphology in the [ feii ] emission , which is more affected by extinction than the pa@xmath1 emission .
it is therefore tempting to suggest that the south part of the ring of ionized gas is actually _ broken _ , perhaps due to sn explosions .
this hypothesis would also explain the excess of [ feii ] emission in this region .
the presence of the hole of gas emission at the center of m82 along with the anti - correlation between the morphology of the hydrogen recombination lines and the co index has been invoked as evidence for an outwards propagating starburst ( e.g. , satyapal et al .
the diffuse component of pa@xmath1 emission is also in the form of `` fingers '' or `` filaments '' above and below the plane of the galaxy .
it is now known from ground - based h@xmath1 imaging that this emission along the minor axis of the galaxy extends over 11kpc ( devine & bally 1999 ) .
figure  1b and figure  3b show how there is a very close spatial correspondence between the morphologies of the extended [ feii]@xmath0 m and pa@xmath1 emissions .
figures  4a and 4b are the @xmath25 m continuum and the [ feii]@xmath0 m line emission maps of ngc  253 .
as seen in m82 , the [ feii]@xmath0 m emission is in the form of a number of compact regions along the disk of the galaxy superimposed on a diffuse emission component .
figures  5a and 5b are contour plots of these two emissions .
we also show the positions of the 6 cm radio sources with @xmath77mjy from ulvestad & antonucci ( 1991 ) , excluding those sources that are classified as hii regions in ulvestad & antonucci ( 1997 ) .
ground - based [ feii]@xmath0 m and br@xmath43 maps at lower spatial resolution ( @xmath78 ) were previously observed by forbes et al .
they found that the central [ feii ] emission of ngc  253 is dominated by three peaks of emission .
our higher resolution nicmos image shows a larger number of [ feii ] compact sources .
for instance , source b ( in the forbes et al .
1993 notation ) breaks up into three individual sources .
in addition to the nucleus and sources a and b detected by forbes et al .
( 1993 ) , there are other sources along the disk of the galaxy , as well as diffuse emission perpendicular to the disk , more prominent in the south half of the galaxy .
we believe that forbes et al . (
1993 ) did not detect the diffuse component above and below the disk of the galaxy because of the lower s / n of their [ feii ] map .
we also mark in these figures the position of the peak ( the filled square ) of the near - infrared emission ( also coincident with the peak of the pa@xmath1 emission ) which was identified as the nucleus of the galaxy by forbes et al .
other authors ( e.g. , keto et al .
1999 ) claim that a bright non - thermal radio source located ne of the peak of the infrared continuum emission ( source 2 in turner & ho 1985 ) may be the true nucleus of ngc  253 .
it is shown in these figures as a filled star symbol .
the nicmos pa@xmath1 line map ( figure  6b ) shows very bright emission stemming from the peak of infrared continuum emission , together with more diffuse emission in a plume - like morphology along the plane of the galaxy , primarily to the ne of the nucleus . as for m82
, there is an extended component of the pa@xmath1 emission perpendicular to the disk , better seen in our nic3 images ( not shown here ) which are more sensitive to diffuse emission .
this component is however not as strong as that in m82 .
there is a remarkable resemblance of the central pa@xmath1 emission with both the [ neii]@xmath34 m emission line and the nearby continuum images presented in bker , krabbe , & storey ( 1998 ) and keto et al .
both m82 and ngc  253 are known to be highly obscured ( see for instance the comparison between the optical _ hst_/wfpc2 f555w and the _ hst_/nicmos f222 m images of m82 in alonso - herrero et al .
2001 ) . to calibrate the supernova rate in terms of the [ feii]@xmath0 m emission and to obtain the total number of ionizing photons ,
it is necessary to correct the observed quantities for extinction .
we can estimate the extinction to the gas in m82 using observations of the ratios of hydrogen recombination and near - infrared [ feii ] lines .
we used published values of the line fluxes for the smallest apertures and matched those apertures on the nicmos line images when measuring the pa@xmath1 and [ feii]@xmath0 m fluxes .
the intrinsic ratios of hydrogen recombination lines are from hummer & storey ( 1987 ) for an electron temperature of @xmath79k and electron density of @xmath80@xmath81 .
the ratio between [ feii]@xmath82 m and [ feii]@xmath0 m does not depend upon the physical conditions of the gas .
we used a polynomial fit to the near - infrared extinction law given in he et al .
( 1995 ) .    since we do not have nicmos images of two hydrogen recombination lines , to measure extinction at the full nicmos angular resolution we made use of the broad - band filters nic2 f160w and nic2 f222 m to construct a map of extinction to the stars in m82 .
we assumed a color for the unreddened stellar population of @xmath83 ( which applies almost independently of the age of the population ) , and , for computational convenience , a foreground dust screen model following the extinction law of he et al .
the contribution of dust emission at the @xmath21-band was assumed to be negligible .
figure  7a is the extinction map ( @xmath84 ) from the @xmath85 color map of m82 , with the positions of the radio sources superimposed .
the values of the extinction in the @xmath21-band derived from this map are between @xmath86mag and @xmath87mag .
the m82 extinction map qualitatively looks very similar to those of satyapal et al .
( 1995 ) , but quantitatively the values of the extinction are different .
the estimates of the extinction to the gas for m82 are presented in the first part of table  2 along with the line fluxes .
the errors in table  2 are the uncertainties associated with the aperture size , aperture centering and flux calibration . for comparison
, we also estimated the extinction through a small aperture ( 3 ) using the f222m / f160w continuum flux ratio , which gives an estimate of the visual extinction to the stars .
as can be seen from this table the agreement between the extinctions to the gas and to the stars is excellent at infrared wavelengths .
satyapal et al . (
1995 ) reached a similar conclusion
. a much lower level of extinction is indicated by comparing the h@xmath1 and pa@xmath1 line fluxes .
this behavior is known to arise from the different optical depths probed by optical and near - infrared lines .
this behavior is confirmed by the much higher extinction of @xmath88mag inferred by puxley et al .
( 1989 ) from the radio h53@xmath1 to mid - infrared br@xmath1 line ratio .
we also computed the extinctions to the e1 , w1 and w2 regions in achtermann & lacy ( 1995 ) notation ( see table  2 for the positions ) .
these were found to be bright in [ neii]@xmath34 m and br@xmath1 .
we used their br@xmath1 measurements along with the mid - infrared extinction curve in rieke & lebofsky ( 1985 ) . since achtermann & lacy ( 1995 )
did not quote aperture sizes , we assumed that their fluxes were measured over the extent of these regions ( approximately @xmath89  in diameter ) .
the values of the extinction for these three regions are slightly larger than the values we would measure from the @xmath85 extinction map , perhaps indicating minor differences in the optical depth .
a further test of the extinction to the gas was based on the total number of ionizing photons from the extinction corrected pa@xmath1 line emission map .
the number of ionizing photons obtained from the extinction - free h53@xmath1 radio recombination line is @xmath90s@xmath16 ( puxley et al . 1989 ) , and from the @xmath91 mm continuum is @xmath92s@xmath16 ( carlstrom & kronberg 1991 ) .
we measured through a large aperture the pa@xmath1 flux on the observed and the extinction - corrected images .
we obtained @xmath93s@xmath16 ( observed ) and @xmath94s@xmath16 ( corrected for extinction ) , for electron temperature @xmath95k .
the total number of ionizing photons corrected for extinction is in excellent agreement with a number of estimates including those of puxley et al .
( 1989 ) , carlstrom & kronberg ( 1991 ) and frster schreiber et al .
( 2001 ) .      as for m82
, the extinction to the stars was computed from the @xmath85 color map of ngc  253 ( figure  7b ) .
it shows values between @xmath96mag and @xmath97mag and is very similar to that of sams et al .
( 1994 ) derived from a @xmath98 color map .
we find a similar value of the extinction ( @xmath99mag ) for the point with the highest value of the extinction which was identified by sams et al .
( 1994 ) as the true nucleus of ngc  253 .    to calculate the extinction to the gas in ngc  253 we used pa@xmath1 and br@xmath43 fluxes ( the latter from puxley & brand 1995 and harrison et al . 1998 ) .
both works centered their apertures on the approximate coordinates of the secondary near - infrared peak and therefore we will assume that the largest br@xmath43 fluxes in these works correspond to the primary near - infrared peak . if their ( 0,0 ) position were to be identified with the peak of the continuum emission , their fluxes would not be consistent with forbes et al .
( 1993 ) br@xmath43 line emission map .
this map and our pa@xmath1 line emission map ( figure  6b ) show that the brightest hydrogen line emission arises from the peak of the continuum emission . in table
3 we give estimates for the extinction for the primary near - infrared peak ( center ) and secondary near - infrared peak ( positions 3  and 3.9  to the ne of the primary peak for puxley & brand 1995 and harrison et al .
1998 works respectively ) .
the errors take into account 10% and 20% combined flux calibration / centering uncertainties for the nucleus and ne source , respectively .
the values of the gas extinction to these two regions as derived from the pa@xmath1/br@xmath43 line ratio are consistent ( to within the errors ) with the extinction to the stars from the @xmath85 color map . despite this apparent good agreement
we caution the reader of the aforementioned centering uncertainties .
we also measured the observed and corrected for extinction pa@xmath1 fluxes which yielded a number of ionizing photons of @xmath100 ( observed ) and @xmath101 ( corrected for extinction ) , for the field of view shown in figure  6b .
the latter value is consistent with engelbracht et al .
( 1998 ) and keto et al .
( 1999 ) estimates of the number of ionzing photons from their large aperture br@xmath43 and [ neii ] measurements , but smaller than puxley et al .
( 1997 ) estimate ( @xmath102 ) from millimeter - wavelength hydrogen recombination lines .
to determine whether there is a spatial correspondence between the compact radio sources and the [ feii ] sources in both m82 and ngc  253 , very precise coordinates are needed .
unfortunately , due to problems with the _ hst _ telemetry during the observations of m82 , and the uncertainties when producing and rotating the mosaics , the coordinates of the nic3 images may not be completely reliable .
in addition , it is expected that there are minor discrepancies between the radio and the optical / infrared coordinate reference systems .
our approach was to determine the coordinates of the [ feii ] image by trying to find the highest possible number of coincidences between the positions of the radio snrs and the [ feii ] sources .
we use daofind within the daophot package in iraf to find compact sources . the term `` compact '' is not used in the sense of point sources .
in fact , most of these sources seem to be slightly resolved , although due to the under - sampling of the nic3 psf , some uncertainties remain .
the sources were selected to be easily distinguishable from the extended diffuse emission detected above and below the plane of the galaxy .
in m82 we identified 65 `` compact '' sources in the [ feii]@xmath0 m map ( not corrected for extinction ) . for the radio sources in m82 we used the positions of snrs identified by searching for compact features within radio maps at 408mhz , 1.4ghz , 5ghz and 8.4ghz , with peak flux densities above 1mjy ( dr .
karen wills , private communication 2000 ) .
we excluded from this list radio sources that are classified as hii regions based on their radio spectral index ( see mcdonald et al .
2002 for details ) , resulting in 44 radio snrs .
we began the comparison of positions using the coordinates in the header of the nicmos nic3 f164n images and made small offsets of the [ feii ] image , cross - correlating the positions of the [ feii ] sources with the radio snrs .
based on the resolution of the nic3 images and the average angular size of the radio snrs ( @xmath103 , huang et al .
1994 ) , we counted a coincidence when the angular separation between a radio source and an [ feii ] source was @xmath104 . with a small shift , we found 21 coincidences between radio snrs and [ feii ] sources ,
whereas if we used the original coordinates we could only find 7 coincidences .
monte carlo simulations showed that the probability of finding 21 coincidences by pure chance is @xmath105 ( see whitmore & zhang 2002 for details on the method ) , so we believe that now the infrared and the radio images are properly aligned .
for ngc  253 we found discrepancies between the coordinates derived for the near - infrared peak of emission from the nic2 and nic3 images ( which were taken during two different orbits ) .
although no problems were reported for these two sets of observations , some error must be affecting the coordinates given in the headers of the images . again as with m82 , we cross - correlated the positions of the [ feii]@xmath0 m ( 36 sources ) and 35 radio sources taken from ulvestad & antonucci 1991 excluding those sources that are classified as hii regions in ulvestad & antonucci ( 1997) using the coordinates of the nic3 images .
we find a total of 9 coincidences for a separation between the radio sources and [ feii ] sources @xmath106 .
the probability of finding 9 coincidences in by pure chance in ngc  253 is approximately @xmath107 .
one important conclusion , despite the uncertainties of the image alignment , is that we can not assign an [ feii ] source to each radio snr or vice versa .
we find approximately @xmath3% of radio snrs in m82 and ngc  253 have an [ feii ] .
this point is illustrated in the contour images of the [ feii]@xmath0 m line emission and continuum emission at @xmath25 m in figures  2a and 2b for m82 and figures  5a and 5b for ngc  253 . in all figures
we show the positions of the radio snrs as asterisks .
we also tried to determine if there is a relation between the presence of a coincidence and the age of the radio snr .
such a correlation would be expected if the [ feii ] emission were brighter in older supernova remnants .
for m82 we used the angular sizes of the radio snr as an indication of their age ( see huang et al .
in general , we did not find that older radio snrs tend to have a coincidence with an [ feii ] source .
we performed aperture photometry on the [ feii ] sources identified in m82 image registered to match the field of view and position of the pa@xmath1 image , the number of detected sources ( 71 ) is slightly different from that in section  5.1 ] and ngc  253 using a fixed @xmath108 diameter aperture .
this aperture corresponds to a physical size of 12pc and 10pc for m82 and ngc  253 , respectively .
figures  9a and 9b show the distributions of the measured [ feii]@xmath0 m luminosities ( not corrected for extinction ) .
some of the sources may be larger than the aperture used for the photometry ; on the other hand , patchy extinction will make some of the sources appear smaller . from inspection of typical source sizes , however , the luminosities would be accurate on average .
recent studies of galactic snr show that the [ feii ] emission is filamentary and have sizes of the order of 11pc ( e.g , ic443 rho et al .
2001 ; see also morel et al .
2002 for snr in m33 ) .    for the assumed distances to m82 and ngc  253 the observed [ feii]@xmath0 m luminosities of these sources range from @xmath109w to @xmath110w .
as a consistency check for the photometry , we measured the [ feii]@xmath0 m fluxes for the fe1 , fe2 and fe3 sources in m82 identified by greenhouse et al .
( 1997 ) using the aperture diameters listed in their table  1 .
the agreement with greenhouse et al .
( 1997 ) fluxes for these sources is excellent , to within @xmath111 .
however , when we compared our photometry for the sources reported in forbes et al .
( 1993 ) for ngc  253 we found some discrepancies .    the corrected [ feii]@xmath0 m luminosities for compact sources in m82 and ngc  253 range from @xmath112w to @xmath113w . as a comparison ,
galactic snrs yield [ feii]@xmath0 m luminosities of @xmath114w to @xmath115w ( oliva et al .
1989 ; keller et al . 1995 ) , snrs in the lmc yield luminosities from @xmath116w to @xmath112w ( oliva et al . 1989 ) and snrs in m33 ( assuming a distance of 720kpc ) @xmath117 to @xmath112w ( lumsden & puxley 1995 and morel et al .
2002 ) .
the individual sources account for some 22% and 27% of the total [ feii]@xmath0 m emission ( extinction corrected ) measured in m82 and ngc  253 , respectively .
given the estimates that the [ feii ] emitting phase of a snr has a duration of @xmath118yr ( lumsden & puxley 1995 ; morel et al .
2002 ) and the supernova rates in m82 and ngc  253 are @xmath119yr@xmath16 ( huang et al . 1994 and ulvestad & antonucci 1997 ) , the detected [ feii ] sources account for @xmath120yr worth of supernovae .
presumably much of the remaining unresolved [ feii ] arises from additional supernovae that are packed so densely that they are confused in our images .
also , snrs in these two galaxies may expand and merge into a general hot ism within a few thousand years losing their identity as individual sources .
photometry on the [ feii ] emission arising from the disk of m82 shows that as much as 70% of the total emission ( disk plus component perpendicular to the disk ) is associated with snrs .
this together with the good agreement between the radio and [ feii ] morphologies of the components perpendicular to the disk gives further indication for a snr origin of the [ feii ] emission in starburst galaxies .
lumsden & puxley ( 1995 ) pointed out that one possibility is that the [ feii ] lines reach maximum luminosity when the remnant as a whole undergoes the transition from adiabatic expansion to radiative expansion . according to shull & draine ( 1986 )
, this happens when radiative cooling occurs faster than dynamical times , around @xmath121yr after the supernova explosion .
therefore if the brightest [ feii ] emission occurs at approximately @xmath121yr , it is not surprising that we do not find a correspondence for all the [ feii ] sources and all the radio snrs .
muxlow et al .
( 1994 ) found that half of the radio sources in m82 are unresolved indicative of a young age , but the other resolved half are probably intermediate age .
this together with the fact that we found that @xmath122 of the radio snr have an [ feii ] counterpart implies that the [ feii ] sources may be @xmath123yr .
one possibility for the lack of one - to - one coincidences between the radio snrs and the [ feii ] sources would be that most of the [ feii ] emission in these two galaxies is produced in hii regions ( photoionization ) . to check this
, we measured the pa@xmath1 fluxes for the sources detected in the [ feii ] extinction corrected map of m82 , and assumed an [ feii]@xmath0m / pa@xmath1 ratio for hii regions as that observed in orion .
we found that only @xmath124 of the total [ feii]@xmath0 m arises from hii regions .
this further supports the hypothesis that most of the [ feii]@xmath0 m emission along the plane of m82 is produced by snrs .
the [ feii]/br@xmath43 ratios in hii regions and snr differ by a factor of a hundred or more , and therefore can be used to indicate the excitation mechanism of the [ feii ] line within our images . in alonso - herrero et al .
( 1997 ) we showed that the behavior of the [ feii]@xmath0m / br@xmath43 line ratio in starbursts and active galaxies can be understood as a progression from pure photoionization to pure shock excitation going from pure hii regions through starbursts , transition objects , seyferts to supernova remnants .    figures  8a and 8b are the observed [ feii]@xmath0m / pa@xmath1 maps for m82 and ngc  253 , respectively .
the values of the [ feii]@xmath0m / pa@xmath1 ratio range from 0.02 in black ( mostly along the disk of the galaxy ) to 0.8 in white ( some of the bright [ feii ] sources ) .
these ratios translate into [ feii]@xmath0m / br@xmath43 line ratios of between 0.24 and 10 ( see also the histograms in figure  10 ) .
since these maps are not corrected for extinction , these numbers are actually lower limits . when the positions of radio sources are plotted in the [ feii]@xmath0m / pa@xmath1 maps we find again that there is a number of coincidences between the radio snrs and the regions of enhanced [ feii]@xmath0 m emission .
this indicates that the relative contribution of the snrs is nearly as large in the diffuse background as it is in the radio snr .
interestingly , we find that most of the coincidences with the radio snrs occur in regions of enhanced [ feii]@xmath0m / pa@xmath1 line ratio , or conversely we do not find a large number of coincidences in regions with high pa@xmath1 flux ( that is , hii regions ) .
the supernova rate in m82 has been estimated from number counts of the radio supernova remnants and their ages to be @xmath125yr@xmath16 ( huang et al .
the total [ feii]@xmath0 m luminosity in m82 is @xmath126([feii]@xmath127w .
our value for the [ feii ] luminosity of m82 is approximately 1.5 and 2 times higher than the previous estimates by greenhouse et al .
( 1997 ) and vanzi & rieke ( 1997 ) , respectively .
this difference is most likely due to better extinction estimates made possible by our nicmos images .
if the supernova rate in m82 is @xmath128yr@xmath16 , then we can derive the following calibration for the supernova rate in terms of the [ feii]@xmath0 m luminosity :    @xmath129}}{10^{34 } { \rm w}}\,\,\ , ( { \rm yr}^{-1})\ ] ]    the total [ feii]@xmath0 m luminosity ( corrected for extinction ) measured in ngc  253 is @xmath130 .
an upper limit on the supernova rate from radio supernova counts is @xmath131 ( ulvestad & antonucci 1997 ) .
a rate of @xmath132 has been calculated assuming that the non - thermal radio flux is entirely due to synchrotron emission ( van buren & greenhouse 1994 ) ; this value is consistent with the similar number of compact radio sources in this galaxy compared with m82 .
the calibration of supernova rate in terms of [ feii]@xmath0 m luminosity for this galaxy is then    @xmath133}}{10^{34 } { \rm w}}\,\,\ , ( { \rm yr}^{-1})\ ] ]    since m82 is the more thoroughly studied galaxy in the radio , we give higher weight to the calibration for it in averaging these results to estimate    @xmath134}}{10^{34 } { \rm w}}\,\,\ , ( { \rm yr}^{-1})\ ] ]    in the above calculations , we have included the entire [ feii ] luminosity , since that is the parameter usually measured in a distant galaxy .
however , in these nearby examples , we can compare the [ feii ] emission produced in hii regions with that from supernovae .
we assume that the ratio for a typical hii region is @xmath135([feii]@xmath0m@xmath136 ( e.g. , orion ) .
we obtained that the [ feii ] luminosity produced in hii regions accounts for less than 8% the total [ feii ] luminosity in m82 .
we have presented _ hst_/nicmos observations of the starburst galaxies m82 and ngc  253 .
the morphology of the [ feii]@xmath0 m line emission in both galaxies displays two components : ( 1 ) a disk component with a number of superimposed compact sources and ( 2 ) a diffuse component in the form of filaments above and below the disk of the galaxy . in m82
we find a remarkable similarity between the extended component , the fingers or chimneys above and below the disk , of the radio and [ feii ] emissions .    in m82 ,
we have detected in both pa@xmath1 and [ feii ] a nuclear ring of star formation .
this ring , with an approximate diameter of 400pc , has also been observed in other near- and mid - infrared hydrogen recombination lines and in the [ neii]@xmath34 m line .
the peak of the @xmath35 m continuum brightness of m82 lies very close to ( within 1 ) the dynamical center , and corresponds to the center of the ring of star formation . in ngc  253 ,
bright pa@xmath1 emission stems from the near - infrared continuum peak and extends in a plume - like morphology along the disk of the galaxy .
there is also a diffuse pa@xmath1 component above and below the plane of galaxy , however , much fainter than that in m82 .    prior to calibrating the surpernova rate in terms of the [ feii ] emission ,
we have derived the extinction to the central regions of ngc  253 and m82 . to assess the extinction correction ,
we have compared the number of ionizing photons infered from the extinction - corrected pa@xmath1 maps with other estimates based on longer wavelength indicators , and found a good agreement .    since the [ feii ]
emission in starburst galaxies is thought to be associated with supernova activity , we have cross - correlated the positions ( with a spatial resolution of a few tenths of arcseconds ) of the radio snrs and the identified compact [ feii ] sources in m82 and ngc  253 .
we have found that approximately @xmath3% of the radio snrs have an [ feii ] counterpart .
however , we can not associate an [ feii ] source to each radio snr or vice versa .
this is well understood in terms of the different ages traced by the radio snrs ( a few hundred years ) and the [ feii ] stage of snrs ( @xmath4yr ) .
the compact [ feii ] emitting sources in m82 and ngc  253 only contribute to some @xmath137 of the total [ feii]@xmath0 m emission . however , much of the remaining emission in the disk of the two galaxies is most likely produced by snrs that expanded and merged into a general ism a few @xmath138yr ago .
photometry on the [ feii]@xmath0 m emission arising from the disk of m82 shows that as much as @xmath139 of the total [ feii]@xmath0 m emission ( corrected for extinction ) is associated with snrs . the remaining diffuse [ feii ] emission perpendicular to
the disk of m82 seems to be related to the superwind ( via expulsion of material from the central starburst ) , as argued from the radio emission morphology ( wills et al .
1999 ) .
all these arguments support the supernova origin for the observed near - infrared [ feii ] line emission in starburst galaxies .
thus we have recalibrated the supernova rate in terms of the [ feii ] luminosity relation for starburst galaxies using our newly measured ( corrected for extinction ) [ feii]@xmath0 m luminosities of m82 and ngc  253 .
we would like to thank drs .
karen wills , t. muxlow and a. pedlar for kindly providing us with the radio images for comparison with the nicmos images , and the list of radio snr .
we are also grateful to dr .
ernesto oliva for enlightening discussions , and an anonymous referee for constructive comments that helped improve the paper .
during the course of this work aa - h was supported by the national aeronautics and space administration on grant nag 5 - 3042 through the university of arizona .
the work was also partially supported by the national science foundation under grant ast-95 - 29190 .
satyapal , s. , watson , d. m. , pipher , j. l. , forrest , w. j. , coppenbarger , d. , raines , s. n. , libonate , s. , pich , f. , greenhouse , m. a. , smith , h. a. , thompson , k. l. , fisher , j. , woodward , c. e. , & hodge , t. 1995 , apj , 448 , 611


the description of quantum mechanical systems by means of the schrdinger equation has been developed more than 80 years ago . at that time
only a few resonance states in nuclei and atoms were known which are well separated from one another .
the energies of these states are well described by means of the schrdinger equation with a hermitian hamiltonian . in order to describe the finite lifetimes of these states ,
the r - matrix theory has been developed which is however too complicated for practical calculations .
the finite lifetimes of the individual states are calculated usually perturbatively , see e.g. @xcite .
later , systems at high level density were in the center of interest .
under these conditions , a statistical description of the system turned out to be very efficient @xcite .
the resonances are well isolated from one another , and the lifetimes of the states are very long and do not play any role .    in the course of time , experimental studies have been performed for different quantum systems with a much improved accuracy . also the theoretical calculations are carried out today not only for states being well separated in energy from one another but also for resonance states in the regime of overlapping . here
the single individual states can no longer be identified what results in problems of their interpretation .
contradictions between experimental results and conventional hermitian quantum physics appeared in different small quantum systems .
for example , heiblum et al .
@xcite found experimentally a crossover from the mesoscopic to a universal phase for electron transmission in quantum dots about than 10 years ago .
these results could not be explained in conventional quantum physics in spite of much effort @xcite . recently
koehler et al .
@xcite observed that neutron resonance data exclude random matrix theory .
deviations between experimental data and random matrix theory in nuclear physics studies were observed also earlier , e.g. @xcite .
these and other experimental results show that the schrdinger equation originally introduced with a hermitian hamilton operator for the description of well isolated resonances , has to be extended . above all , the lifetimes of the resonance states have to be calculated also in the regime of overlapping resonances and the justification of statistical assumptions has to be proven for small quantum systems .
the lifetimes of the resonance states can be calculated when the system is explicitly considered to be open and the calculations are performed quantum mechanically for both the system and the environment of scattering wavefunctions into which the system is embedded . using the feshbach projection operator technique @xcite , first the _ energy - independent _ many - body problem of the system ( with the hermitian hamiltonian @xmath4 ) is solved in the standard manner .
these solutions provide the energies @xmath5 and wavefunctions @xmath6 of the discrete states with inclusion of the so - called _ internal _ interaction . in a second step ,
the _ energy - dependent _ scattering wavefunctions @xmath7 of the environment are calculated and , further , the ( energy - dependent ) coupling matrix elements @xmath8 between the discrete states of the system and the environment are evaluated ( see @xcite , section 2.1 ) .
the corresponding schrdinger equation contains an energy dependent ( nonlinear ) source term describing the coupling between system and environment .
the hamiltonian @xmath9 of this equation is non - hermitian and provides the lifetimes of the states , but without any additional interaction of the states via the environment .
the schrdinger equation with @xmath9 and source term can be rewritten into a schrdinger equation without source term but with a non - hermitian hamiltonian @xmath10 that contains the interaction of the states via the environment ( the so - called _ external _ interaction ) in the nondiagonal matrix elements @xcite .
the hamiltonian @xmath11 reads @xmath12 where @xmath13 and @xmath14 stand for the coupling between system and environment and @xmath15 is the green function in the subspace of scattering states .
the external interaction of the states via the continuum is complex , generally .
the principal value integral is @xmath16 and the residuum reads @xmath17 the interaction of the states of the system via the environment is involved in the eigenvalues @xmath18 and eigenfunctions @xmath19 of the hamiltonian @xmath10 .
it is therefore relatively easy , in this formalism , to study the influence of the environment onto the states of the system .
that means , the non - hermitian quantum physics is  in contrast to the widely spread meaning  _ not _ a further approximation introduced in the theory .
it is an expression of the fact that the system considered is really open and its properties are influenced by the coupling to the environment .
this influence is unimportant at low level density where it can be described by perturbation theory .
it becomes , however , decisive in the regime of overlapping resonances .
meanwhile , many calculations for open quantum systems are performed with high accuracy . without using any perturbation theory or statistical assumptions , the schrdinger equation with the non - hermitian hamilton operator @xmath10
is solved and the eigenvalues @xmath20 and eigenfunctions @xmath19 are obtained .
that means , not only the energies @xmath21 of the states of the system are evaluated but also their lifetimes @xmath22 .
the control of the eigenvalues by a parameter allows to draw conclusions on the dynamics of open quantum systems .
the limits of the applicability of the standard quantum theory with hermitian hamilton operator can be seen best at high level density , where the individual resonance states overlap and interact with one another via the continuum of scattering wavefunctions . here , the dynamics of the system is determined by singular points , the so - called _
exceptional points_. they cause level repulsion in energy as well as a bifurcation of the widths ( inverse proportional to the lifetimes ) of the states .
results are obtained theoretically as well as experimentally , which are counterintuitive at first glance : due to width bifurcation , coherent short - lived states are formed together with states that are almost decoupled from the continuum of scattering states @xcite .
an example is the above - mentioned crossover from the mesoscopic to a universal phase in the transmission in quantum dots @xcite and its qualitative explanation on the basis of a schrdinger equation with non - hermitian hamilton operator @xcite .
this phenomenon is observed in different systems and is called usually _ dynamical phase transition _ @xcite .
the short - lived states are analog to the coherent superradiant states considered by dicke @xcite in 1954 as shown by means of a toy model @xcite .
moreover , it could be shown in an experiment on non - locally coupled pairs of quantum point contacts that discrete states undergo a robust interaction that is achieved by coupling them to each other through the continuum @xcite .
most of these results are of high value for fundamental questions of quantum mechanics as well as for applications .
it is the aim of the present paper to show some generic ( mostly numerical ) results for open quantum systems in order to receive a deeper understanding for the dynamical phase transitions occurring in the regime of overlapping resonances .
a toy model is used in order to receive conclusions on the role played by exceptional points and avoided level crossings at high level density in open quantum systems .
of special interest are the effects arising from the imaginary part ( [ form12 ] ) of the coupling term via the environment .
the calculations are performed with respectively two , four and ten resonance states coupled mostly to one open decay channel ( corresponding to the common continuum of scattering wavefunctions ) under the condition that they avoid crossing or even cross in a certain parameter range .
the coupling of @xmath23 resonance states to @xmath24 channels corresponds to the general situation of open quantum systems .
the formalism used in the calculations is the same as that discussed in our earlier paper @xcite .
the results show very clearly that the imaginary part of the coupling term between system and environment causes width bifurcation and , finally , the formation of different time scales in the regime of overlapping resonances . due to width bifurcation ,
the system splits into two parts that exist at different times .
this process is irreversible .
the wavefunctions of the short - lived _ aligned _ states are mixed coherently in relation to the open decay channel while those of the long - lived _ trapped _ states are mixed incoherently such that the states are almost decoupled from the open decay channel .
the results confirm the characteristic features of the dynamical phase transitions appearing in open quantum systems at high level density .
they show that the number of states of the system is reduced during the dynamical phase transition due to ejecting the aligned short - lived state and , furthermore , that the wavefunctions of the states before and beyond the dynamical phase transition are completely different from one another .
while the states of the original system ( without interaction via the continuum ) have individual spectroscopic features , those of the system consisting of the long - lived states beyond the dynamical phase transition show chaotic features .
these results prove once more that the experimental results @xcite on the crossover from the mesoscopic to a universal phase for electron transmission in quantum dots can be explained by means of the non - hermiticity of the hamiltonian .
additionally it should be mentioned here that the aligned state corresponds to the dicke superradiant state and the trapped states to the subradiant states which are considered in optics , e.g. @xcite .    in section [ form ] , we sketch the formalism used in the present calculations . for @xmath25 states with equal decay widths ,
the analytical and numerical solutions of the problem show clearly the effect of width bifurcation .
numerical results for eigenvalues and eigenfunctions obtained for @xmath25 , 4 and 10 states are provided in sections [ 2 ] to [ v10 ] .
the results are discussed in section [ disc ] while conclusions on the relation between exceptional points and width bifurcation are drawn in the last section [ concl ] .
we consider an @xmath26 matrix @xmath27 the diagonal elements of which are the @xmath23 complex eigenvalues @xmath28 of a non - hermitian operator .
the @xmath29 are the so - called selfenergies of the states arising from their coupling to the environment of scattering wavefunctions into which the system is embedded . in atomic physics , these values are known as lamb shift .
our calculations are performed with coupling matrix elements @xmath29 the values of which do not depend on the parameter considered . in such a case , the @xmath29 can considered to be included into the diagonal matrix elements , which read @xmath30 .
the @xmath31 and @xmath32 denote the energies and widths , respectively , of the @xmath23 states ( including their selfenergies ) without account of the interaction of the different states via the environment .
the internal interaction of the two states @xmath33 and @xmath34 ( appearing in the closed system ) _ as well as _ their external interaction ( via the environment ) are contained in the @xmath35 .
the internal interaction can be caused only by some part of re@xmath36 while the external interaction contains complex @xmath35 , see equations ( [ form11 ] ) and ( [ form12 ] ) .
most interesting part of the external interaction is therefore im@xmath36 .
it becomes important at high level density where the corresponding resonance states overlap .    when the number @xmath23 of states is equal to the number @xmath37 of common channels ( i.e. equal to the number @xmath37 of different common continua of scattering wavefunctions ) all the coupling matrix elements @xmath35 are different from zero and the matrix ( [ form1 ] ) is full . in the case with only one open decay channel @xmath38 ,
all @xmath35 different from @xmath39 and @xmath40 are zero .
an example of @xmath41 states coupled to only the fourth channel is the following matrix @xmath42 with @xmath43 for @xmath44 and @xmath45 .
the eigenvalues of @xmath11 will be denoted by @xmath46 where @xmath21 and @xmath47 stand for the energy and width , respectively , of the eigenstate @xmath33 .
the eigenfunctions of the non - hermitian @xmath11 are biorthogonal ( see sections 2.2 and 2.3 of @xcite ) , @xmath48 it follows @xmath49 and @xmath50 at an exceptional point , @xmath51 and @xmath52 .
the @xmath53 and @xmath19 contain global features that are caused by many - body forces induced by the coupling @xmath35 of the states @xmath33 and @xmath34 via the environment , see equations ( [ form11 ] ) , ( [ form12 ] ) and the corresponding discussion in @xcite .    in the case
@xmath25 , the two eigenvalues of @xmath10 are @xmath54 according to this expression , two interacting discrete states ( with @xmath55 ) avoid always crossing since @xmath0 and @xmath56 are real in this case . resonance states with nonvanishing widths @xmath47 repel each other in energy according to the value of re@xmath57 while the widths bifurcate according to the value of im@xmath58 .
the two states cross when @xmath59 .
this crossing point is an exceptional point according to the definition of kato @xcite .    at the exceptional point @xmath59 ,
the eigenfunctions of ( [ form1 ] ) of the two crossing states are linearly dependent from one another , @xmath60 according to analytical as well as numerical and experimental studies , see appendix of @xcite and section 2.5 of @xcite .
that means , the wavefunction @xmath61 of the state @xmath62 jumps , at the exceptional point , via the wavefunction  @xmath63 of a chiral state to  @xmath64 @xcite . from ( [ eif5 ] ) follows :    _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ \(i ) when two levels are distant from one another , their eigenfunctions are ( almost ) orthogonal , @xmath65 .
+ ( ii ) when two levels cross at the exceptional point , their eigenfunctions are linearly dependent according to ( [ eif5 ] ) and @xmath66 .
+ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _    these two relations show that the phases of the two eigenfunctions relative to one another change when the crossing point is approached .
this can be expressed quantitatively by defining the _ phase rigidity _
@xmath67 of the eigenfunctions @xmath68 , @xmath69 it holds @xmath70 .
the non - rigidity @xmath67 of the phases of the eigenfunctions of @xmath11 follows also from the fact that @xmath71 is a complex number ( in difference to the norm @xmath72 which is a real number ) such that the normalization condition ( [ int3 ] ) can be fulfilled only by the additional postulation im@xmath73 ( what corresponds to a rotation ) .    when @xmath74 , an analytical expression for the eigenfunctions as a function of a certain control parameter can , generally , not be obtained .
an exception is the special case @xmath75 for which @xmath76 . in this case
, the condition @xmath59 can not be fulfilled if @xmath77 is real due to @xmath78 the exceptional point can be found only by analytical continuation into the continuum @xcite and the two states avoid crossing .
this is analogue to the avoided level crossings of discrete states .
the condition @xmath59 can be fulfilled however if @xmath79 is imaginary , @xmath80 and two exceptional points appear .
it holds further @xmath81 independent of the parameter dependence @xmath82 . in the first case ,
the eigenvalues @xmath83 differ from the original values @xmath84 by a contribution to the energies and in the second case by a contribution to the widths .
the width bifurcation starts at one of the exceptional points and becomes maximum in the middle between the two exceptional points .
this happens at the crossing point @xmath85 where @xmath86 .
some years ago , the case @xmath25 with @xmath87 , fixed real @xmath88 , and different fixed values of @xmath32 , including @xmath89 , is studied as a function of the parameter @xmath90 in the neighborhood of avoided and true crossings of the two levels @xcite . the results for the @xmath25 case @xcite show further that the wavefunctions of the two states @xmath61 and @xmath91 are mixed in a finite range of the parameter @xmath90 around the critical value @xmath92 at which the two states avoid crossing .
this holds true not only for resonance states but also for discrete states .
the eigenfunctions @xmath19 of @xmath10 can be represented in the set of basic wavefunctions @xmath93 of the unperturbed matrix ( corresponding to the case with vanishing coupling matrix elements @xmath94 ) , @xmath95 also the @xmath96 are normalized according to the biorthogonality relations of the wavefunctions @xmath97 .    in our calculations ,
the mixing coefficients @xmath96 of the wavefunctions of the two states due to their avoided crossing are not calculated .
we simulate the fact that the two wavefunctions are mixed in a finite parameter range around the critical value of their avoid crossing @xcite by assuming a gaussian distribution @xmath98 for the coupling coefficients .
the results reproduce very well those discussed in @xcite for 2 levels and real coupling @xmath0 .    of special interest
is the situation at high level density where the ranges of avoided crossings , defined by ( [ int7 ] ) , of different levels overlap .
some generic results obtained with 2 , 4 and 10 resonance states will be presented in the following sections .
we start our calculations with the two - level case which is studied mostly in literature .
we choose the matrix ( [ form1 ] ) with @xmath87 , different fixed values of @xmath32 and fixed @xmath0 .
the functional dependence of the energies over the parameter @xmath90 is similar as in @xcite .
however , the @xmath0 are real in @xcite while they are mostly complex in the present paper as well as in @xcite . according to ( [ int6 ] ) to ( [ int6d ] )
, the exceptional points appear at different values of the parameter @xmath90 when the ratio re@xmath1 to im@xmath1 is varied , see also figs
. 1 to 3 in @xcite .    in fig .
[ fig1 ] we show the numerical results for the avoided level crossing phenomenon as a function of the parameter @xmath90 beyond , at and below the exceptional point . in all cases , the value of @xmath0 is fixed to @xmath100 while the @xmath101 are different in the different subfigures ( with @xmath102 in all cases ) .
the exceptional point appears at @xmath103 , see figs .
[ fig1].c , d. when @xmath104 ( figs .
[ fig1].a , b ) , the eigenvalue trajectories cross in energy while the widths @xmath105 and @xmath106 are always different from one another .
the situation is another one when @xmath107 . here
, the eigenvalue trajectories avoid crossing in energy while the trajectories @xmath105 and @xmath106 cross at certain values of @xmath90 ( figs .
[ fig1].e to h ) .    the comparison of fig .
[ fig1 ] with complex @xmath108 and fig .
2 in @xcite with real @xmath109 shows the influence of im@xmath1 .
while the @xmath110 vary symmetrically around the critical value @xmath111 when @xmath0 is real , this is not the case when @xmath0 is complex . in the last case ,
the difference @xmath112 blows up in the critical region .
furthermore , the @xmath47 approach their asymptotic values in a relatively small parameter range of @xmath90 when @xmath0 is real in contrast to the case with complex @xmath0 . these results will be discussed further in the following sections .     and widths @xmath113 ( full lines ) of @xmath25 states beyond  ( a , b ) , at  ( c , d ) and below  ( e to h ) the exceptional point .
the parameters of the subfigures are @xmath114  ( a , b ) ;  @xmath115  ( c , d ) ;  @xmath116  ( e.f ) ;  @xmath117
( g , h ) . further parameters :  @xmath118 .
the dashed lines show @xmath82 . ,
width=529,height=529 ]     of @xmath25 states beyond ( a ) , at ( b ) and below ( c , d ) the exceptional point .
the parameters are the same as in fig .
[ fig1 ] . ,
width=264,height=529 ]    in fig .
[ fig2 ] , we show the mixing coefficients @xmath119 of the eigenfunctions as a function of the parameter @xmath90 corresponding to the eigenvalue figures shown in fig . [ fig1 ] . the mixing coefficient is complex , @xmath120 .
this figure can be compared with fig .
5 in @xcite calculated exactly with real @xmath109 .
in both cases , the difference between the two curves @xmath119 is always 1 when @xmath104 .
this result corresponds to the fact that the two states will not be exchanged when @xmath104 .
the energy trajectories of the eigenstates cross freely and the widths of the two states are always different from one another , see figs .
[ fig1].a , b.    at the exceptional point ( where @xmath121 ) , @xmath122 in both cases . when @xmath107 , the two curves @xmath119 coincide at a certain value of @xmath90 .
this happens for real ( see @xcite ) as well as for complex ( fig .
[ fig2 ] ) @xmath0 , but is symmetrically only for real @xmath0 . at these points ,
the two states are exchanged .    in both cases @xmath107 and @xmath104 , the asymptotic values 1 and 0 are reached in a smaller parameter range when @xmath0 is real than in the case with complex @xmath0 . within this parameter
range the two eigenfunctions are mixed . the range of mixing shrinks to one point when the two states really cross ( at the exceptional point ) .
it is especially large when two discrete states ( with @xmath123 ) avoid crossing , as can be seen from fig . 5 in @xcite ( where the calculations are performed exactly with real @xmath0 for all @xmath96 ) .    the parameter range in which two states are mixed due to the existence of an exceptional point in the neighborhood plays an important role in realistic systems
it will be discussed in the following sections with more than one avoided crossing for both real and complex @xmath0 .
and widths @xmath113 ( full lines ) of @xmath25 states ( a , b , g , h ) and @xmath41 states ( c to f and i , j ) coupled to @xmath38 channel .
the parameters of the subfigures are @xmath125 ( a to d ) ;  @xmath126 ( e , f ) ;  @xmath127 ( g , h ) ;  @xmath128 ( i , j ) . further parameters : @xmath129 .
the dashed lines show @xmath82 . , width=491,height=642 ]     but @xmath130 ( a , b , e , f ) ; @xmath109 ( c , d , g , h ) ;  @xmath41 ;  @xmath131 ( a to d ) ;  @xmath132 ( e to h ) . , width=529,height=529 ]     but @xmath41 and @xmath133 ;  @xmath130  ( a , b ) ;  @xmath134  ( c , d ) ;  @xmath109  ( e , f ) ;  @xmath135 . + ,
width=529,height=396 ]    let us first compare the calculations with @xmath25 and @xmath41 states and imaginary coupling @xmath134 between discrete states and environment of scattering states ( fig .
[ fig3 ] ) . in the calculations ,
the eigenvalue trajectories @xmath136 and @xmath110 are traced as a function of the parameter @xmath90 which is varied in such a manner that one state crosses in energy one and three , respectively , other states .
the results show two exceptional points in both cases appearing at the two different parameter values @xmath137 and @xmath138 , see figs .
[ fig3].a , b for two states and figs .
[ fig3].c to f for four states . in between these two values , the widths bifurcate : the differences @xmath139 blow up .
the two - level case is described analytically by equations ( [ int6 ] ) to ( [ int6d ] ) .
the width bifurcation is completely reproduced in the numerical results .
the parameter range @xmath140 is larger in the case with four states than in the other one due to the larger region in which avoided level crossings take place .
furthermore , the difference @xmath141 is larger in the case of 4 states than in the case of 2 states . in the case that all states have equal widths @xmath142 , the widths @xmath113 of only two states bifurcate also in the case with four states .
the widths of the two other states remain unchanged ( figs .
[ fig3].c , d ) .
this result holds also when the different @xmath142 differ slightly from one another ( figs .
[ fig3].e , f ) . in this latter case ,
the two states avoid crossing at the two critical values @xmath137 and @xmath138 and the width bifurcation is almost the same as in the case with equal widths @xmath32 .
the situation changes completely when the widths @xmath32 differ stronger from one another . when @xmath143.g to j ) .
the eigenvalue trajectories cross freely in energy and the widths do ( almost ) not bifurcate .
the dynamics of the system is therefore completely different from that determined by the results shown in figs .
[ fig3].a to f.    in fig
. [ fig4 ] , we show the eigenvalues @xmath21 and @xmath113 for the case with four levels and one channel with complex and real coupling coefficients ( @xmath120 and @xmath109 , respectively ) .
width bifurcation can be seen when @xmath0 is complex and the widths @xmath32 are equal ( figs .
[ fig4].a , b ) , while the @xmath144 are independent of @xmath90 when @xmath0 is real ( figs .
[ fig4].c , d ) .
the results with different @xmath32 differ from those obtained with imaginary coupling @xmath145 ( figs .
[ fig3].g to j ) : the widths of two states cross at parameter values at which their energies avoid crossing ( figs . [
fig4].e to h ) .    in fig .
[ fig5 ] , we show the results of calculations with equal number of states and channels , @xmath41 and @xmath133 .
although usually @xmath24 in realistic systems , the results allow us to receive a deeper understanding for the spectroscopic redistribution processes taking place in the critical region .
the calculations are performed with the same parametric dependence of the energies as in the foregoing calculations with @xmath41 and @xmath38 .
further assumptions : @xmath146 and @xmath147 , and 0.05 , respectively .
the results show width bifurcation when @xmath0 is complex or imaginary ( figs .
[ fig5].a to d ) which is however smaller than in the corresponding cases with @xmath38 ( figs .
[ fig4].a , b and figs .
[ fig3].c , d , respectively ) . instead , the widths @xmath47 of the states are changed asymptotically due to the coupling to the different channels . when @xmath0 is real ( figs .
[ fig5].e , f ) , the widths @xmath113 are independent of @xmath90 and are equal to @xmath148  ( @xmath149 1 to 4 ) .
of @xmath25 and @xmath41 states the eigenvalues of which are shown in fig .
[ fig3 ] ;  @xmath38 channel ;  @xmath134 . , width=264,height=680 ]     of @xmath41 states the eigenvalues of which are shown in fig .
[ fig4 ] ; @xmath38 channel ;  @xmath130  ( a , c ) ;  @xmath109  ( b , d ) . , width=302,height=585 ]     of @xmath41 states the eigenvalues of which are shown in fig .
[ fig5 ] ; @xmath133 channels ;  @xmath130  ( a ) ;  @xmath134  ( b ) ;  @xmath109  ( c ) . + , width=302,height=453 ]    the coefficients @xmath96 defined in ( [ int20 ] ) determine the mixing of the eigenfunctions @xmath19 in relation to the basic wavefunctions @xmath150 .
the mixing is nonvanishing in the neighborhood of the critical points of avoided level crossings as has been shown in exact numerical calculations for two crossing states and real coupling coefficients @xmath94 @xcite . in the calculations of the present paper ,
the mixing of the wavefunctions is simulated by assuming a gaussian distribution for the @xmath94 , equation ( [ int7 ] ) .
it reproduces very well the exact results obtained in @xcite .    in the following ,
we consider the mixing coefficients @xmath96 for all wavefunctions @xmath19 the eigenvalue trajectories of which are studied in figs .
[ fig3 ] to [ fig5 ] as a function of the parameter @xmath90 . in fig .
[ fig6 ] , the coefficients @xmath119 are shown for two and four states with imaginary coupling @xmath145 ( see the corresponding eigenvalue trajectories in fig .
[ fig3 ] ) . when the widths @xmath32 of all states are equal ( figs .
[ fig6].a , b ) , @xmath151 at the two exceptional points @xmath137 and @xmath138 . in the parameter range @xmath152 ,
the wavefunctions are mixed while @xmath153 and @xmath154 for @xmath155 beyond the two values @xmath156 .
the states are completely ( 1:1 ) mixed in the range @xmath152 in the two - level - case ( fig . [ fig6].a ) .
the picture is more complicated in the four - level - case ( fig . [
fig6].b ) . here
, all states are involved in the redistribution taking place in the critical parameter range @xmath152 .
the exchange of the two states the widths of which remain unchanged according to fig .
[ fig3].d can be seen from the energy eigenvalue trajectories , fig .
[ fig3].c , as well from the mixing coefficients , fig .
[ fig6].b .
the figures are similar when the widths @xmath32 of the states differ slightly from one another . instead of exceptional points ,
there are avoided level crossings at @xmath137 and @xmath138 , see fig .
[ fig6].c for the four - level - case .
the mixing coefficients show a dependence on @xmath90 which is similar to that obtained for equal widths ( fig . [
fig6].b ) .
when the widths @xmath32 of the states @xmath33 differ stronger from one another than in fig .
[ fig6].c , the eigenfunctions remain almost unmixed for all parameter values @xmath90 ( figs .
[ fig6].d , e ) .
this result corresponds to the almost constant width trajectories @xmath47 as a function of @xmath90 ( figs .
[ fig3].h , j ) .
an exchange of states does not take place ( figs .
[ fig3].g , i ) .
these results show the great influence of the exceptional points onto the dynamics of the system considered .    the calculations in figs .
[ fig3 ] and [ fig6 ] are performed with imaginary coupling @xmath145 . for comparison ,
we show in fig .
[ fig7 ] the mixing coefficients @xmath119 when @xmath0 is complex and real , respectively , for the four - level - case ( corresponding to the eigenvalue trajectories in fig .
[ fig4 ] ) . when all @xmath32 are equal to one another and @xmath0 is complex , the mixing coefficients @xmath119 point to the participation of all states in the redistribution process in the critical parameter range ( fig .
[ fig7].a ) .
this corresponds fully to the corresponding eigenvalue trajectories that show level exchange as well as width bifurcation ( figs .
[ fig4].a , b ) . when the coupling is real , however , the wavefunctions of only two states are mixed while the wavefunctions of the other two states remain almost unchanged in the critical region ( fig .
[ fig7].b ) .
the mixing of the eigenfunctions is completely different from that discussed above when the widths @xmath32 differ more strongly from one another .
a mixing of the wavefunctions appears only when the energy trajectories @xmath136 avoid crossing and the width trajectories @xmath110 cross , compare figs .
[ fig4].e to h with the corresponding figs .
[ fig7].c , d. a width bifurcation does not take place .    the mixing coefficients @xmath119 for the case with @xmath41 states coupled to @xmath133 channels are shown in fig .
[ fig8 ] ( parameters the same as in fig .
[ fig5 ] ) .
the differences between the cases with coupling to 4 channels to those with coupling to only 1 channel can be seen by comparing fig .
[ fig8].a with fig .
[ fig7].a ,  fig .
[ fig8].b with [ fig6].b and fig .
[ fig8].c with fig .
[ fig7].b .
the mixing of the wavefunctions by coupling the system to four channels is more complicated than that by coupling to only one channel .
this statement holds true also when @xmath0 is imaginary .
furthermore , an additional mixing is caused by the avoided level crossing at small @xmath90 in the four - channel - case when @xmath0 is complex .
we continue our studies with the matrix ( [ int2 ] ) by choosing @xmath158 states coupled to @xmath62 channel .
the eigenvalue trajectories @xmath136 and @xmath110 are shown in figs .
[ fig9 ] and [ fig10 ] . in these cases , one of the states crosses in energy the remaining nine states one by one .
the nine states have the same parametric energy dependence of @xmath90 and are shifted equidistantly relative to one another .     and
widths @xmath113 ( full lines ) of @xmath158 states coupled to @xmath38 channel .
the parameters of the subfigures are @xmath125 ( a , b ) ;  @xmath159 ( c to f ) ;  @xmath160 ( a to d ) ;  0.07 ( e , f ) .
further parameters :  @xmath161 .
the dashed lines show @xmath82 . , width=529,height=415 ]     but  @xmath134 ( a , b ) ;  @xmath162 ( c , d ) ;  @xmath163 ( e , f ) ;  @xmath164 ( g , h ) . , width=529,height=529 ]     but @xmath165 ;  @xmath166 ( a , b ) ;  @xmath167 ( c , d ) and  @xmath168 ( e , f ) . ,
width=529,height=415 ]    in fig .
[ fig9].a , b we show the results obtained with @xmath169 and @xmath170 for all states .
the results are similar to those for the 4-level case ( figs .
[ fig3].c , d ) .
two exceptional points appear as well as width bifurcation : the width of one of the states is much larger while that of another one is much smaller than the widths of all the other states in the whole critical parameter range between the two exceptional points .    in figs .
[ fig9].c to f , the widths @xmath32 differ from one another .
we compare the results obtained with an imaginary and with a real coupling constant , @xmath171 and @xmath172 , respectively .
when @xmath0 is real ( figs .
[ fig9].e , f ) , the energy trajectories @xmath136 are very regular : avoided level crossings appear one by one with all nine states . an exceptional point can be found only by analytical continuation into the continuum ( see @xcite ) .
the width trajectories @xmath110 are less regular : they cross at several parameter values @xmath90 . according to the results obtained for a smaller number of states in sections [ 2 ] to [ f4 ]
, the widths do not bifurcate when @xmath0 is real .
when @xmath0 is imaginary , the energy trajectories @xmath136 are also regular ( fig . [
fig9].c ) .
[ fig9].d shows nicely how widths continue to bifurcate such that the broad state remains broad and the narrow state remains narrow for all @xmath90 in the critical region .
mostly , the widths of two states are different from one another when the energy trajectories cross .
the states can , nevertheless , be exchanged as can be seen from fig .
[ fig9].c .    in fig .
[ fig10 ] , we show results obtained with the imaginary coupling constant @xmath145 ( being smaller than in figs .
[ fig9].c , d ) as well as those obtained with the larger values @xmath173 ( figs .
[ fig10].a , b and figs .
[ fig10].c to h , respectively ) .
we see width bifurcation similar as in fig .
[ fig3].d for @xmath41 states and in fig .
[ fig9].b for 10 states with equal widths @xmath32 . in any case
, width bifurcation appears in the whole critical parameter range between the two exceptional points .
the width bifurcation occurring in the case of 10 levels is stronger than that with a smaller number of levels ( compare section [ v4 ] ) .
the 9 states crossed by the state 10 are exchanged in the critical parameter region as can be seen from the corresponding figs .
[ fig10].a , c , e , g.    the results shown in figs .
[ fig9 ] and [ fig10 ] follow from calculations with a toy model being symmetrical in relation to the crossing states . in the critical parameter range ,
the states are exchanged in all cases but 8 of them do not contribute actively to the width bifurcation when @xmath0 is imaginary .
this result follows from the high symmetry of the exceptional points in relation to the crossing states , i.e. from the linear dependence of the energy of the crossing state on the parameter @xmath90 ,  @xmath174 .
the 8 states play the role of observers similarly as discussed in a toy model with 4 states @xcite and in a realistic model with 3 states @xcite .    in order to exclude the high symmetry of the exceptional points in relation to nearby states we show in fig .
[ fig11 ] the results of another version of the toy model . instead of @xmath174 in figs .
[ fig9 ] and [ fig10 ] , we use @xmath165 in fig .
[ fig11 ] .
the width @xmath47 of one of the states becomes large in the critical parameter range when @xmath0 is imaginary or complex ( figs .
[ fig11].b and d , respectively ) , while the widths of all the other states are much smaller .
this result reflects the situation observed in many realistic cases ( see e.g. the review @xcite and also @xcite ) .
a separation of a short - lived state from the other ones does not appear when @xmath0 is real ( fig . [
fig11].f ) .
in sections [ 2 ] to [ v10 ] , we showed results obtained numerically for the eigenvalues @xmath175 and eigenfunctions @xmath19 by using the matrix ( [ form1 ] ) ( or ( [ int2 ] ) ) with 2 , 4 and 10 levels , respectively . only the energies @xmath31 of the states are varied as a function of a certain parameter @xmath90 .
the widths @xmath32 are assumed to be constant , the interaction @xmath94 between the levels @xmath33 and @xmath176 is either zero or simulated by the gaussian distribution ( [ int7 ] ) where @xmath0 is the same for all states ( and also independent of @xmath90 ) , and the self - energy terms are considered to be included into the diagonal matrix elements .
most interesting is the one - channel case for different @xmath23 , e.g. the matrix ( [ int2 ] ) for @xmath41 .
this toy model describes well the generic features of open quantum systems .
the @xmath94 are complex .
they stand for the interaction of the states via the environment ( continuum of scattering wavefunctions ) : re(@xmath94 ) arises from the principal value integral ( [ form11 ] ) and im(@xmath94 ) is the residuum ( [ form12 ] ) .
both parts play an important role for the dynamics of the open quantum system at high level density where the resonance states overlap .
the influence of re(@xmath94 ) onto eigenvalues @xmath18 and eigenfunctions @xmath19 is studied numerically exact for @xmath25 levels in an earlier paper @xcite .
it causes the avoided crossing of the states which is well known for discrete as well as for narrow resonance states .
equations ( [ int6a ] ) to ( [ int6d ] ) show analytically the main difference between re(@xmath94 ) and im(@xmath94 ) in the 2-level case when the widths of the two states are equal , @xmath177 . according to ( [ int6a ] )
, an exceptional point can be found only by analytical continuation into the continuum when @xmath94 is real .
thus , the two states avoid crossing as it is very well known for discrete as well as for narrow resonance states @xcite . in contrast , im(@xmath94 ) causes two exceptional points according to equation ( [ int6b ] ) .
most interesting is the width bifurcation arising in between the two exceptional points according to ( [ int6d ] ) .    the analytical results ( [ int6a ] ) to ( [ int6d ] ) are well reproduced in our numerical calculations for the eigenvalues @xmath18 , figs .
[ fig3].a , b. also the eigenfunctions @xmath19 show the two exceptional points : @xmath122 in approaching them , fig .
[ fig6].a . in between the two exceptional points , the wavefunctions of the two states
are strongly mixed .
the mixing is 1:1 in the middle between the two singular points . beyond the critical region between the two exceptional points , the @xmath119 approach 1 when @xmath178 and 0 when @xmath155 . here
, the two states may be exchanged , at most .
the figures 3.a , b and 6.a represent the basic process of width bifurcation according to the analytical expressions ( [ int6c ] ) and ( [ int6d ] ) .    in sections
[ v4 ] , [ f4 ] and [ v10 ] of the present paper , numerical results for 4 and 10 states are shown under similar conditions , i.e. with equal or nearly equal @xmath32 for all states : figs . [ fig3].c to f and figs .
[ fig6].b and c for the eigenvalues and eigenfunctions , respectively , of 4 states , and figs .
[ fig9].a , b for the eigenvalues of 10 states .
the development of width bifurcation as a function of increasing imaginary coupling vector @xmath0 can be seen in fig .
[ fig10 ] for @xmath158 states when the single widths @xmath32 are different from one another . the uniform width bifurcation in fig .
[ fig9].b appears at a smaller value of im@xmath1 than in fig .
[ fig10].h due to the different values of the single @xmath32 in the last case . generally , the width bifurcation is stronger when the number of states is larger .    the influence of re(@xmath94 ) onto the eigenvalue trajectories and eigenfunctions is also shown for @xmath41 and 10 ( figs . [
fig4].c , d for the eigenvalues of 4 states , fig .
[ fig7].b for the eigenfunctions of 4 states and figs .
[ fig9].e , f for the eigenvalues of 10 states ) . under the influence of re(@xmath94 ) , the states avoid crossing in energy and the widths do not bifurcate . when the calculation is performed with complex coupling vector ( re@xmath179 _ and _ im@xmath179 ) , the widths bifurcate with some shift of the position of the maximum relative to that of the minimum ( figs .
[ fig4].a , b ) and a larger parameter range of mixed wavefunctions ( fig . [ fig7].a ) .    according to the results presented in sections [ v4 ] to [ v10 ]
the eigenvalue trajectories are strongly influenced by the values of the external ( fixed ) parameters when the resonances overlap and exceptional points determine the dynamics of the system .
however , the eigenvalue trajectories are almost independent of one another when the degree of resonance overlapping is small .
examples are figs .
[ fig3].g to j and [ fig4].e to h for the eigenvalues of , respectively , two and four resonance states and the corresponding figs .
[ fig6].d , e and [ fig7].c , d for the eigenfunctions of these states .
these results show the strong influence of external parameters onto the dynamics of the system at high level density which is known from the study of realistic cases .
an example is the enhanced transmission through microwave cavities when the formation of whispering gallery modes is supported by the manner the leads are attached to the cavity @xcite .
width bifurcation is directly related to the _ alignment _ of one of the @xmath23 resonance states of the open quantum system to a decay channel ( @xmath38 ) with the consequence that it becomes short - lived while other states become _ trapped _ ( long - lived ) , i.e. more or less decoupled from the environment .
mathematically , the alignment of a resonance state is possible since the eigenfunctions of a non - hermitian operator are biorthogonal and their phases ( relative to those of the eigenfunctions of the other states ) are not rigid in approaching an exceptional point , see equation ( [ eif11 ] ) .
width bifurcation appears therefore in our calculations when @xmath38 and @xmath41 or 10 .
this corresponds to realistic situations in which usually @xmath24 .
width bifurcation does , however , not appear when @xmath41 states are coupled to @xmath133 channels ( fig .
[ fig5 ] ) although the wavefunctions are strongly mixed in the critical parameter region also in this case ( fig .
[ fig8 ] ) .
an artifact of our model calculations is the assumption of equidistant energies @xmath31 of all but one state which are crossed by one state the energy of which depends linearly on the parameter @xmath90 .
thus , the critical points are symmetrically in relation to the two neighboring states , and @xmath180 states do not contribute actively to the width bifurcation taking place in the whole system . in fig .
[ fig11 ] , we show results obtained when the symmetry is somewhat disturbed ( as it is usually the case in realistic systems ) . in such a case , the width of one state separates from those of _ all _ the other states .
that means , all nine states remain trapped ( decoupled ) in the whole critical parameter range .
this result corresponds to results obtained theoretically as well as experimentally in different realistic systems ( see section 4 in the review @xcite ) .
it should be added here , that asymptotically ( beyond the critical parameter region ) the states can be exchanged at most . this is the case , indeed , in all our calculations .
all redistribution processes caused by the exceptional points take place only in the critical parameter region .
this holds true for real as well as for imaginary coupling coefficients @xmath94 and can be seen also in the mixing coefficients @xmath96 which approach 1 or 0 when @xmath178 and @xmath155 , respectively , beyond the critical parameter region .
the two cases with real and imaginary @xmath0 differ , however , by the length of the critical region .
the results of the present paper show that exceptional points cause width bifurcation because the coupling of the states via the environment ( continuum of scattering wavefunctions ) is complex .
an example are the bound states in the continuum appearing at a
_ finite _ value of the coupling strength between system and environment when the interplay between internal and external interaction is taken into account , see section 4.4 of the review @xcite . of particular interest
is the width bifurcation appearing at high level density . here
, many exceptional points are near to one another .
as shown in the present paper , width bifurcation causes , under this condition , a splitting of the system into two parts one of which exists in the short - time scale while the other one appears in the long - time scale .
the two parts can not be observed together .
the long - lived states occur as fluctuations in the short - time scale while the short - lived states perform a smooth background in the long - time scale .
this phenomenon is known in literature from theoretical as well as experimental studies on different small quantum systems and systems equivalent to them ( see examples in the review @xcite and in @xcite ) . in many - body open quantum systems ,
it is called mostly _ dynamical phase transition_. it is known also in optics where it is called _ superradiant phase transition _ according to dicke @xcite . in pt - symmetric systems ,
the phase transition is called mostly _ pt - symmetry breaking _ @xcite .
common to all these studies is that the dynamical phase transition is very robust when the necessary conditions are fulfilled , i.e. when the level density is high and the dynamics of the system is determined by exceptional points . although the existence of exceptional points is decisive ( as shown in the present paper ) , the dynamical phase transition does _ not _ appear at the parameter value which gives the position of the exceptional point itself .
instead , the phase transition occurs in the neighborhood of one ( or several ) exceptional points .
such a result is known from different experimental studies on concrete realistic systems .
it is discussed in e.g. @xcite .
the results of the present paper show that the critical parameter range is determined by the distance between ( at least ) two exceptional points .
the dynamical phase transition appearing in small quantum systems at high level density , causes a new understanding of time as discussed in @xcite .
time which is characteristic of the system , is inverse proportional to the widths @xmath47 and related therefore directly to the non - hermitian part of the hamilton operator .
the widths do not increase limitless as shown in the present paper .
instead , the system is dynamically stabilized : the lifetimes of the states of the system are increased and the system as a whole is stabilized by ejecting the short - lived state from the system at the dynamical phase transition .
thus , the time characteristic of the system is bounded from below in a similar manner as the energy . beyond the dynamical phase transition ,
all states have lost their original spectroscopic features and the number of states is reduced .
this fact is very well known in nuclear physics ( although very seldom interpreted in this manner ) .
the properties of compound nuclei are described well by the _ unified theory of nuclear reactions _ developed by feshbach @xcite , which contains both the short - lived direct reaction part and the long - lived compound nuclear reaction part .
however , feshbach introduces statistical assumptions in order to describe the long - lived compound nuclear states .
this is in contrast to the calculations in the present paper where the long - lived states are described _ without _ any statistical assumptions .
instead the long - lived states are shown by us to arise from width bifurcation causing a dynamical stabilization of the system at high level density , a phenomenon that is caused by exceptional points .
many of these results seem to be counterintuitive .
they allow us , however , to explain some unexpected experimental results .
for example , the crossover from the mesoscopic to a universal phase for transmission in quantum dots observed experimentally @xcite can be explained qualitatively by the formation of the short - lived resonance state at the dynamical phase transition .
the results of the present paper support this interpretation @xcite of the experimental results .
another example is the experimental observation @xcite of non - statistical effects in nuclear reactions on middle - heavy nuclei .
the data show directly the formation of different time scales in the system at high level density .
other examples are discussed in @xcite .
the long - lived states beyond the dynamical phase transition are strongly mixed and show chaotic features .
the random matrix theory is surely applicable only after averaging over different decay channels ( as usually done ) , i.e. it is not applicable to the description of small systems coupled to one channel ( and , respectively , to two channels in the case of transmission through the system ) . a theoretical analysis of the spectra ( without any statistical assumptions ) by restricting to the results obtained from one decay channel according to the recent experimental data on compound nuclei @xcite
is not performed up to now . on the basis of the results of the present paper ,
it can be stated today only the following .
the mixing of the long - lived states ( beyond the dynamical phase transition ) is caused by complex many - body forces via the continuum of scattering wavefunctions ( simulated by the @xmath94 in our calculations ) and not by two - body forces .
this may provide an explanation of the fact why compound nucleus spectra ( after averaging over different channels and beyond decay thresholds ) may be described by a gaussian orthogonal ensemble ( containing many - body forces ) , but not by a two - body random ensemble .
this point should be studied in future in more detail .    in order to prove the role of exceptional points in an open many - level quantum system
it is highly interesting to study experimentally time symmetry breaking caused by the influence of a nearby state onto an exceptional point .
theoretical studies with a symmetrical influence ( as in most calculations of the present paper ) can surely not describe the properties of realistic systems with broken time symmetry .
symmetry breaking influences not only width bifurcation ( as shown in fig .
[ fig11 ] ) but will prove , above all , the irreversible processes caused by exceptional points in open quantum systems .
these processes are decisive , among others , for the dynamical stabilization of quantum systems and the formation of quantum chaos .
finally we remark that the sensitive dependence of the dynamics of an open quantum system on the value of external parameters ( as shown in the present paper ) can be used in order to construct small systems with desired properties .
this point is important for applications .
a. yacoby , m. heiblum , d. mahalu , and h. shtrikman , phys .
lett . * 20 * , 4047 ( 1995 ) ; r. schuster , e. buks , m. heiblum , d. mahalu , v. umansky , and h. shtrikman , nature * 385 * , 417 ( 1997 ) ; m. avinun - kalish , m. heiblum , o. zarchin , d. mahalu , and v. umansky , nature * 436 * , 529 ( 2005 )    m. goldstein and r. berkovits , new j. phys .
* 9 * , 118 ( 2007 ) ; d.i .
golosov , and y. gefen , new j. phys .
* 9 * , 120 ( 2007 ) ; y. oreg , new j. phys .
* 9 * , 122 ( 2007 ) ; p.g .
silvestrov and y. imry , new j. phys .
* 9 * , 125 ( 2007 )                i. rotter , contribution to the special issue _ quantum physics with non - hermitian operators : theory and experiment _ , fortschritte der physik - progress of physics * 61 * , 178 ( 2013 ) ; doi : 10.1002/prop.201200054      g.a .
lvarez , e.p .
danieli , p.r .
levstein , and h.m .
pastawski , j. chem . phys . * 124 * , 194507 ( 2006 ) ; g.a .
lvarez , e.p .
danieli , p.r .
levstein , and h.m .
pastawski , phys .
rev . a * 75 * , 062116 ( 2007 ) ; h.m .
pastawski , physica b * 398 * , 278 ( 2007 )        y. yoon , l. mourokh , t. morimoto , n. aoki , y. ochiai , j.l .
reno , and j.p .
bird , phys .
lett . * 99 * , 136805 ( 2007 ) ; y. yoon , m.g .
kang , t. morimoto , l. mourokh , n. aoki , j.l .
reno , j.p .
bird , and y. ochiai , phys .
b * 79 * , 121304(r ) ( 2009 ) ; y. yoon , m.g .
kang , t. morimoto , m. kida , n. aoki , j.l .
reno , y. ochiai , l. mourokh , j. fransson , and j.p .
bird , phys .
x * 2 * , 021003 ( 2012 )    h. eleuch and i. rotter , contribution to the special issue _ quantum physics with non - hermitian operators : theory and experiment _ , fortschritte der physik - progress of physics * 61 * , 194 ( 2013 ) ; doi : 10.1002/prop.201200062          in studies by other researchers , the factor @xmath33 in ( [ eif5 ] ) does not appear .
this difference is discussed in detail and compared with experimental data in the appendix of @xcite and in section 2.5 of @xcite .
a. guo , g.j .
salamo , d. duchesne , r. morandotti , m. volatier - ravat , v. aimez , g.a .
siviloglou , and d.n .
christodoulides , phys .
. lett . * 103 * , 093902 ( 2009 ) ; c.e .
rter , k.g .
makris , r. el - ganainy , d.n .
christodoulides , m. segev , and d. kip , nature physics * 6 * , 192 ( 2010 ) ; t. kottos , nature physics * 6 * , 166 ( 2010 )