a classical nova eruption plays an important role in the evolution of cataclysmic variable stars ( cvs ; see the books by * ? ? ? * ; * ? ? ?
it represents the primary mechanism for mass - loss in these systems , preventing ( most of ) them to end up as supernovae ia .
irradiation of the donor star by the eruption - heated white dwarf is furthermore suspected to cause an increase in the mass - transfer rate , which would have a strong effect on the evolutionary time - scale @xcite . still , the connection between the nova eruption and the physical parameters of cvs is not very well known , mainly due to the fact that a large fraction of the post - nova systems remain unidentified @xcite .
we are carrying out a systematic observing campaign in order to recover and study classical novae that erupted at least 30 yr ago ( earlier than 1986 ) .
this implies that sufficient time has passed since the nova eruption for the optical signal to become dominated by the underlying cv instead of representing mainly the ejected material .
the targets presented here were selected from @xcite and refer to fields with hitherto ambiguous identification of the post - nova ( el aql , v606 aql , v908 oph , v1149 sgr and v1583 sgr ) .
they were investigated by means of @xmath0 photometry to identify potential candidates for the post - nova , based on their colours .
the selected candidates were subsequently observed spectroscopically to confirm the post - nova via their spectroscopic characteristics ( mainly the presence of characteristic emission lines ) .
in addition , we also obtained spectrograms of candidates for the post - novae v1301 aql , v1151 sgr and v3964 sgr that were included in @xcite .
other collections of post - nova spectra can be found in @xcite , @xcite and parts of our own series tappertetal12 - 1,tappertetal14 - 1,tappertetal15 - 2 .
@llllllll object & ra ( 2000.0 ) & dec .
( 2000.0 ) & rms ( arcsec ) & date & filter / grating & @xmath1 ( s ) & mag + el aql & 18:56:01.87 & @xmath203:19:18.8 & 0.21 & 2012 - 07 - 18 & @xmath3/@xmath4/@xmath5/@xmath6 & 1576/420/160/120 & 20.9v + & & & & 2015 - 04 - 27 & b600 ( 1.5 arcsec ) & 4760 & 19.9r + v606 aql & 19:20:24.29 & @xmath200:08:07.8 & 0.30 & 2012 - 06 - 24 & @xmath3/@xmath4/@xmath5/@xmath6 & 1576/420/160/120 & 20.4v + & & & & 2015 - 06 - 15 & b600 ( 1.5 arcsec ) & 3600 & 20.1r + v908 oph & 17:28:04.58 & @xmath227:43:04.4 & 0.33 & 2014 - 04 - 03 & @xmath3/@xmath4/@xmath5/@xmath6 & 1655/495/200/170 & 20.5v + & & & & 2015 - 04 - 28 & b600 ( 1.5 arcsec ) & 3600 & 18.5r + v1149 sgr & 18:18:29.85 & @xmath228:17:24.9 & 0.21 & 2012 - 08 - 15 & @xmath3/@xmath4/@xmath5/@xmath6 & 1576/420/160/120 & 18.4v + & & & & 2015 - 06 - 06 & b600 ( 1.5 arcsec ) & 1350 & 18.1r + v1583 sgr & 18:15:26.53 & @xmath223:23:17.5 & 0.27 & 2012 - 07 - 24 & @xmath3/@xmath4/@xmath5/@xmath6 & 1576/420/160/120 & 20.5v + & & & & 2015 - 04 - 21 & b600 ( 1.5 arcsec ) & 3600 & 19.5r + v3964 sgr & 17:49:42.42 & @xmath217:23:34.8 & 0.17 & 2015 - 04 - 26 & b600 ( 1.5 arcsec ) & 3600 & 19.0r +    the @xmath0 photometric data for the objects with ambiguous identification were taken in two service mode runs in 2012 ( el aql , v606 aql , v1149 sgr and v1583 sgr ) and in 2014 ( v908 oph ) at the antu 8 m telescope that is part of the very large telescope of the european southern observatory ( eso ) , at cerro paranal , chile , using the focal reducer / low dispersion spectrograph ( fors2 ; * ? ? ?
* ) system with the high - throughput broad - band filters u_high , b_high , v_high and r_special .
series of four and three frames per filter were obtained in the 2012 and the 2014 observations , respectively . in service mode
, fors2 usually operates with a mosaic of two 2k@xmath74k mit ccds . for the data obtained on 2012 august 18 ,
exceptionally , the blue - sensitive e2v ccds were used .
the field of view of fors2 is 6.8@xmath76.8 arcmin@xmath8 .    after combining the two mosaic ccd frames into one with the fsmosaic routine from the fors instrument mask simulator package
, the further reduction was performed using the ccdred package of iraf .
the bias - subtracted and flat - fielded data were then corrected for the telescope offsets between them and averaged to a single frame per filter .
a 3@xmath9 clipping algorithm was applied to exclude bad pixels and cosmic rays .
subsequent photometry was performed using iraf s daophot package and the stand - alone daomatch and daomaster routines @xcite .
the aperture radius for the photometry was chosen to be slightly smaller than the full width at half - maximum ( fwhm ) of the point spread function ( psf ) .
the conversion to calibrated magnitudes was obtained using observations of standard fields @xcite .
the @xmath3 passband is not a part of eso s standard calibration plan , and in most nights , the standard field observations in this filter covered only a very limited range of airmasses that did not allow us to derive a value for the extinction from them .
thus , a standard value for the extinction was used that was taken from the la silla observatory website .
while this observatory is located roughly 480 km south of cerro paranal , comparison of the @xmath10 passbands and also the @xmath3 data for the few nights , where a sufficiently large airmass range was covered , showed identical values within the uncertainties .
similarly , for some nights , also the @xmath10 standard data proved insufficient .
for those filters , median values for a 60 d range are available at the observatory s health check monitor .
the spectroscopic data were obtained using the gemini multi - object spectrograph ( gmos ; * ? ? ?
* ) mounted at the gemini - south 8 m telescope situated on cerro pachn , chile .
three spectra were taken per object , in most cases subsequently .
since 2014 , gmos is equipped with a row of three 2048@xmath74176 hamamatsu ccds .
grating b600 was employed with a 1.5 arcsec slit and a two and four pixels binning in spatial and dispersion direction , respectively .
this combination yielded a typical wavelength range of 4100  7100  and a resolving power of @xmath11800 .
the grating was operated at a central wavelength of 5600  to put the two 0.915 mm wide gaps between the chips at regions of comparatively low scientific interest .
the typically affected wavelength ranges are 5030  5080  and , in principle , 6100  6050  .
however , all observations suffered from a bright large - scale distortion in spatial direction that severely affected the signal at a much broader wavelength range of roughly 5900  6200  and could not be properly corrected for .
the data were reduced using the gemini  gmos extension for iraf as contained in the ureka software package .
after subtraction of an averaged bias frame and subsequent division by a flat - field that was normalized by fitting a high - order cubic spline , a two - dimensional wavelength correction was applied that was obtained from spectra taken with a cuar comparison lamp . after this transformation , the individual spectra were combined employing a 2@xmath9 clipping algorithm to correct for cosmic rays .
the sky was subtracted by fitting a chebychev function of typically third order in spatial direction .
finally , the spectra were extracted and corrected for the instrumental response function by comparison with a spectrum of the standard star ltt 3864 taken on 2015 february 25 .
since all calibrations are based on this single observation of a standard , we emphasize that , while the spectral energy distribution ( sed ) of our objects should be sufficiently well reproduced , the zero - point is not well determined , and thus the data can not be regarded as flux calibrated . above steps were performed both for the combined and the individual spectra to examine them for variability .
the details of the observations of the confirmed post - novae are summarized in table [ obslog_tab ] .
the coordinates in columns two and three were determined by calculating an astrometric correction on the photometric @xmath6 fors2 images or the @xmath12 gmos acquisition frames using starlink s gaia tool ( version 4.4.6 ) with the us naval observatory ccd astrograph catalogue version 4 @xcite .
the associated root mean square ( rms ) is given in column four .
the fifth column states the date of the observations referring to the start of the night , while the sixth column gives details on the instrumental setup , i.e.  filter or grating and slit , and column seven contains the total exposure time per object .
finally , the brightness of the object at the time of the observations is presented in column eight , with the corresponding passband indicated by a letter . for the spectroscopic data ,
the magnitudes were determined via differential photometry on the acquisition frames with respect to five comparison stars in the vicinity of the target .
these values were calibrated with respect to previous photometric data , where available , or with the guide star catalogue , version 2.2 @xcite .
@lllll object & @xmath5 & @xmath13 & @xmath14 & @xmath15 + el aql & 20.88(03 ) & @xmath20.10(15 ) & 0.62(08 ) & 0.70(08 ) + v606 aql & 20.41(03 ) & @xmath20.28(15 ) & 0.40(07 ) & 0.44(07 ) + v908 oph & 20.54(11 ) & @xmath20.10(15 ) & 1.02(15 ) & 1.23(17 ) + v1149 sgr & 18.36(05 ) & @xmath21.44(15 ) & 0.30(12 ) & 0.30(07 ) + v1583 sgr & 20.54(04 ) & @xmath20.16(18 ) & 0.96(09 ) & 0.83(06 ) +     data . in each field , the post - nova is indicated by a circle . ]
@llllllllllll object & & & bowen & & civ + & 4340 & 4861 & 6563 & 4922 & 5016 & 6678 & 7065 & 4645 & 4686 & 5412 & 5812 + el aql & 26(1 ) & 35(1 ) & 47(1 ) & 4(1 ) & 3(1 ) & 8(1 ) & 5(1 ) & & 3(1 ) & 4(1 ) + v606 aql & 18(1 ) & 21(1 ) & 26(1 ) & 4(1 ) & 2(1 ) & 4(1 ) & 2(1 ) & &  &
 + v908 oph & 1(1)@xmath8 & 1(1)@xmath8 & 3(1 ) &  &  &  &  &  &  &  &
+ v1149 sgr & 3(1 ) & 2(1 ) & 6(1 ) &  &  & 1(1 ) &  & 5(1 ) & 3(1 ) & 1(1 ) & 3(1 ) + v1583 sgr & 6(1 ) & 9(1 ) & 26(5 ) & 3(1 ) &  & 4(1 ) & 2(1 ) & 5(1 ) & 2(1 ) &  &  + v3964 sgr & 3(1 ) & 5(1 ) & 12(1 ) & 1(1 ) &  & 2(1 ) &  & 3(1 ) & 2(1 ) &  &  +     + @xmath16 bowen/ blend .
+ @xmath8 emission cores in absorption troughs .
colour  colour diagrams for the fields of reported novae with ambiguous identification for the post - nova are presented in fig .
[ ubvr_fig ] .
we have excluded objects with photometric errors @xmath170.1 mag or @xmath180.005 mag to minimize contamination by very faint sources and artefacts , respectively .
candidates for subsequent spectroscopy were selected based on their position in the colour  colour diagram with respect to the other stars in the field .
an additional selection criterion was the distance to the reported coordinates of the nova , taking into account their detection history .
for example , the search area for novae that were identified on photographic plates was limited to the central 1 arcmin@xmath8 .
subsequent long - slit spectroscopic data were used to confirm the post - nova .
the spectra of the confirmed novae are presented in fig .
[ allspec_fig ] , their photometric data are summarized in table [ ubvr_tab ] .
we briefly describe the methods used to analyse the spectroscopic properties that apply to all targets .
the emission line strength is characterized via the equivalent width @xmath19 , which is calculated with respect to three manually selected wavelength ranges .
two of them are used to define the continuum and the noise on either side of the line , while the third one contains the line itself .
we then integrate over the range of the latter divided by a linear fit to the continuum .
the uncertainty of the resulting @xmath19 is estimated via a monte carlo simulation by repeating the measurement a thousand times to a variation of the data with respect to the addition of a random value within the noise interval .
the resulting equivalent widths are presented in table [ eqw_tab ] .
we make an attempt at recovering the intrinsic sed by correcting for the interstellar extinction in the line of sight .
we use the interface on nasa s infrared science archive ( irsa ) web pages to obtain values @xmath20 from @xcite .
these serve as input for the dereddening routine implemented in iraf and based on the relations from @xcite .
a standard value for the ratio of the total to the selective extinction @xmath21 is employed . in order to characterize the corrected continuum slope ,
we fit the wavelength range between 5000 and 7000  with a power law @xmath22 . for the intrinsic sed
, the exponent @xmath23 then depends on the mass - transfer rate and the system inclination .
larger values for @xmath23 correspond to higher mass - transfer rates ( hotter and thus bluer discs ) and lower inclinations ( a larger fraction of the hotter inner disc is visible ) . however , there are large uncertainties involved in the reddening correction , because these are large - scale extinction maps with a resolution of 2@xmath72 deg and local fluctuations may be present . furthermore , the exact depth position of the nova with respect to the absorbing material is also unknown .
these uncertainties propagate to the derived @xmath23 , which should thus be taken only as an approximate indicator of the intrinsic sed .
line profiles for the four individual spectra of el aql in sequence from bottom to top .
the vertical line indicates the rest wavelength of h@xmath23 . ]
the nova was discovered by @xcite as a star of ninth magnitude on photographic plates taken on 1927 july 30 and 31 .
examination of pre - discovery harvard plates by @xcite revealed the object having reached its peak brightness before that date .
the brightest recorded value is 6.4 mag on 1927 june 15 , which is thus taken as maximum brightness , although observational gaps before and after that date leave room for an even brighter peak value .
the decline of the nova is comparatively well covered , especially in the first 200 d after eruption .
description of the spectral evolution during the decline is found in a number of publications .
the most detailed data can be found in @xcite , while @xcite present the only drawing of the spectrum .    in spite of the good photometric coverage
, there is no complete representation of the long - term light curve in the literature .
the plot in @xcite is limited to the first 100 d , and @xcite does not include the harvard data . in fig.[elaqllc_fig ] , we collect the observations from @xcite , @xcite , @xcite , @xcite , @xcite , @xcite and @xcite .
the dates and brightness values from @xcite were extracted from his plot with a precision of 0.25 d and 0.1 mag , respectively .
all other data were available in tabulated form .
for data where only the start of night was given , ut = 0 h was assumed for the observation .
we omitted from the plot the first four data points from @xcite , where he describes the target as ` blurred ' , because they differ drastically from those of other authors .
the upper limits in the plot are from @xcite and @xcite . for reasons of visualization ,
we omit another upper limit ( @xmath1716 mag ) from 1927 may 28 @xcite .
after the initial decline , the nova experienced a series of minor rebrightenings and then settled on a prolonged plateau phase .
there also appear to be further fluctuations around 400 and 500 d after eruption .
these are based on the @xcite data only , but since in the overlapping regions his magnitudes usually agree within 0.2 mag to those of other authors , there is no reason to doubt the reality of these rebrightenings . to our knowledge , there are no further data on the nova after the seasonal gap between 540 and 720 d after eruption . assuming that this is not due to the lack of observations , the final decline must have occurred sometime during that period .
@xcite observed the field again on 1940 june 5 and reported a brightness of 19 mag for the nova , but , as we argue below , this is a likely misidentification .
the decline light curve is classified in @xcite as ` ba ' that includes novae with ` standstills or other minor irregular fluctuations during decline ' .
although el aql is not included in the @xcite catalogue , its light curve likely corresponds to the ` p ' ( plateau ) class , and in particular , shows some similarity to hs sge , especially with respect to the rebrightenings . in el aql ,
the major peak around day 63 makes it difficult to define the end of the initial decline and the start of the plateau phase . if we select the end of this event as the corresponding point , we find the following approximate values for the light - curve parameters as defined by @xcite : the time of the start of the plateau phase , and thus the duration of the initial decline , @xmath2487 d , the duration of the plateau phase @xmath25450630 d , the brightness difference between the start of the plateau phase and the peak magnitude @xmath263.9 mag , and the slope of the decline during the plateau phase to @xmath27 .
these are typical values for class p novae and justify the inclusion of el aql into this category .
@xcite find that this class contains a large fraction ( 6090% ) of the recurrent novae and suggest that some of the other systems in this class are yet unrecognized recurrent novae .
the remaining members of the p class novae appear to be magnetic systems . a candidate for the post - nova
was selected from our @xmath0 photometric data on the basis of its colours ( fig .
[ ubvr_fig ] , table [ ubvr_tab ] ) and its proximity ( 3 arcsec ) to the coordinates reported in @xcite .
the spectrum of this object presents emission lines typical for a cv and we conclude that it is indeed the post - nova .
its brightness is @xmath28 mag and thus considerably fainter than the 19th magnitude star observed by @xcite in 1940 .
@xcite gives a value of @xmath29 mag for el aql for observations made in 1989 .
it is , in principle , possible that both observations above caught the post - nova in a bright state ( or our data were taken in a low state ) .
however , we consider this unlikely , because cvs in bright states are also expected to have bluer colours than in low states , and the object observed by @xcite is significantly redder in the visual range ( @xmath30 mag compared to our 0.62 mag , see table [ ubvr_tab ] ) .
we thus conclude that those observations did not concern the post - nova , but a different object .
similarly , @xcite gives infrared data @xmath31 mag , @xmath32 mag , @xmath33 mag for el aql .
this led @xcite to suspect the nova to contain a giant secondary star .
however , the spectral appearance of the post - nova ( fig .
[ allspec_fig ] ) and its position in the colour  colour diagram does not support this interpretation , and it thus again appears more likely that @xcite observed a different object .
we furthermore note that the resulting eruption amplitude ( defined as the difference between the brightness of the post - nova and the maximum brightness ) of 14.5 mag basically excludes el aql as a recurrent nova candidate .
el aql presents a rich emission line spectrum ( fig .
[ allspec_fig ] , table [ eqw_tab ] ) .
the balmer lines are exceptionally strong for an old nova , which might be indicative of a low mass - transfer rate @xcite .
the strength of heii @xmath344686 and the presence of heii @xmath345412 points to a magnetic system , but since the former does not quite match the strength of h@xmath35 , el aql is more likely to be an intermediate polar rather than a discless system .
the continuum has a slightly red slope , but as can be seen from the comparison of the colour  colour diagrams ( fig .
[ ubvr_fig ] ) , the field of el aql is significantly affected by interstellar reddening . from the irsa website
, we find @xmath20 = 1.11(08 ) mag , and when corrected for that value , an inverse power law with an exponent @xmath36 provides an excellent fit to the continuum slope .    the spectroscopic data on el aql consist of four individual spectra . in fig .
[ elaqllprof_fig ] , we examine the variation of the h@xmath23 line profile over the course ( @xmath111 h ) of our observations .
we find that the line as a whole shows a significant doppler shift .
estimating the shift in radial velocity units by fitting a single gaussian function to the line profiles yields a variation of @xmath37@xmath11100 km s@xmath38 over the observed time range .
furthermore , also the shape of the line profile appears to be variable , pointing to the presence of an additional emission source in the system , such as from an irradiated secondary star or a bright spot ( e.g. * ? ? ?
* ) .    from the gathered evidence ,
we conclude that el aql is a likely intermediate polar seen at a moderately high inclination . in spite of its faintness , its orbital period might therefore be accessible with time series photometry .
line profiles for the three individual spectra of v606 aql in sequence from bottom to top .
the vertical line indicates the rest wavelength of h@xmath23 . ]
this is a little observed nova that was discovered by fleming as an object of seventh magnitude on a photographic plate from 1899 april 21 , as reported by @xcite .
the brightness of the nova on this plate was later more precisely determined to 6.75 mag by @xcite .
there is a @xmath11200 d gap of observations prior to the discovery , and it is thus possible that the real maximum has been missed .
@xcite estimates the peak brightness to 5.5 mag .
the observed decline light curve as shown in @xcite covers @xmath11160 d , and suggests a classification as type ` p ' @xcite . in 1991 , @xcite obtained spectroscopy of the suspected post - nova , but found a g7k4 type spectrum for a @xmath5 = 16.3 mag object .
the south - western component of an unresolved visual binary that was marked in @xcite as the possible post - nova showed promising colours ( fig .
[ ubvr_fig ] , table [ ubvr_tab ] ) .
subsequent spectroscopy of the @xmath5 = 20.4 mag object confirmed the nova ( fig .
[ allspec_fig ] ) .
the spectrum presents comparatively strong balmer emission lines , although not as strong as in el aql .
the high excitation lines are even considerably weaker , and we do not detect heii @xmath345412 . as discussed in section [ elaql_sec ] ,
most members of the ` p'-class appear to be recurrent or magnetic novae .
even if we use the actually observed brightness maximum to calculate the eruption amplitude , the result of @xmath39 mag appears much too large for a recurrent nova @xcite . on the other hand ,
the spectroscopic properties do not present particularly strong evidence for a magnetic nature either .
apparently , v606 aql is an unusual member of the p class .
the spectrum presents a slightly blue continuum slope , and the colour  colour diagram shows that the field is only moderately affected by interstellar extinction .
after dereddening with a value @xmath20 = 0.35(03 ) mag , it became apparent that the continuum was distorted by a number of large ` bumps ' . after examining the three individual spectra
, we concluded that the bumpy structure is not intrinsic but an artefact .
however , the presence of the bumps rendered large parts of the continuum unusable for fitting .
we thus chose a number of small regions of continuum between the blue emission lines as well as a couple of regions around 6000  to determine an inverse power - law exponent of @xmath40 2.02(02 ) .
when examining the h@xmath23 line profile of the individual spectra , the changes are not as pronounced as in el aql , and the dopper shift of the whole profile is not as large ( fig .
[ v606aqllprof_fig ] ) . still , the shape of the line profile is clearly variable , and the width of the emission lines is basically identical to those in el aql .
we thus tentatively conclude that also v606 aql is likely seen at a moderately high inclination , which makes it a potential target for follow - up time series observations .
very few data are available for v908 oph .
its discovery was reported by @xcite who detected oi emission in a near - infrared spectrum taken on 1954 , july 2 .
he estimated the brightness in this spectral range to 9 mag , but suspected that the nova was observed about a week after its maximum brightness .
the only other observations consist of two objective prism spectra taken on 1954 july 27 by @xcite . apart from strong h@xmath23 emission , they could not distinguish more detailed spectral features due to the object having declined to a blue magnitude ` probably less than 14 ' .
if the object was indeed observed several days after having reached its maximum brightness , and accounting for the colour differences between the two observations and for the strong interstellar reddening , the time @xmath41 to have declined by 3 mag below maximum brightness is likely to be in the order of 3040 d , which would classify the object as a moderately fast nova .
the coordinates of the nova had only been recorded with a precision in the order of arcmin , but fortunately , a candidate for the post - nova can be easily identified in the colour  colour diagram ( fig .  [ ubvr_fig ] ) .
its spectrum presents a few weak hydrogen emission lines , with those bluewards of h@xmath23 being embedded in absorption troughs , and a red continuum slope ( fig .
[ allspec_fig ] ) . from the @xmath0 photometry ,
it is clear that the field is strongly affected by interstellar reddening , and the respective catalogue yields an accordingly high @xmath20 = 1.26(02 ) mag .
once corrected for this value , the sed corresponds to a steep blue power law with a negative exponent @xmath42 .
thus , the spectral characteristics are consistent with a high mass - transfer system . in fig .
[ v908ophlprof_fig ] , we show the soectral regions around h@xmath35 and h@xmath23 . in more detail .
we note that high excitation emission lines like the bowen blend and heii are absent or , at least , very weak .
this is in stark contrast to other likely high mass - transfer post - novae like gr sgr . towards the red side of h@xmath23 emission
there appears to be another peak that could be identified with the cii @xmath436578,6583 doublet .
this line has been previously detected in a few old novae with sufficiently narrow h@xmath23 emission lines @xcite , and the latter also applies to v908 oph .
however , the signal - to - noise ratio ( s / n ) in that part of the spectrum is low , and there are a number of similar large features close - by , which we were unable to identify and thus suspect to be artefacts . in consequence ,
the detection of cii in this object remains ambiguous .
we finally remark that we did not detect any significant positional variation of the h@xmath23 line in the three individual spectra .
together with the narrowness of that line , this suggests that the system is seen at low inclination , and that it will be difficult to determine its orbital period .
the nova was discovered by @xcite on an objective prism plate at a photographic brightness of about 9 mag .
later research by @xcite revealed the maximum brightness to 7.4 mag .
the light curve is described as a ` very steep increase in brightness [ ... ] followed by a rapid initial decrease ' , reaching a plateau by 1946 may 28 .
the nova was visible from 1945 may 16 to 1948 may 30 , and afterwards fainted below 14.7 mag .
large gaps in the observations prevent a more detailed analysis .    from our @xmath0 photometry ( fig .
[ ubvr_fig ] ) , we selected two potential candidates for the post - nova based on their colours and coordinates .
spectroscopy showed one of those to present emission lines and a blue continuum , and we thus identify this object as the post - nova .
it is located @xmath1111 arcsec south - west of the previously reported position .
we note that these coordinates do not coincide with the objects labelled as v1149 sgr in @xcite and in @xcite .
the spectrum shows a blue continuum , that becomes somewhat steeper ( @xmath44 ) when corrected for the catalogued interstellar reddening of @xmath20 = 0.39(03 ) mag , and weak emission lines ( fig .
[ allspec_fig ] ) .
in contrast to v908 oph ( section [ v908oph_sec ] ) , v1149 sgr counts with a strong bowen blend and heii @xmath344686 emission ( fig .
[ v1149sgrlprof_fig ] ) . since no other heii lines can be identified , the system is not a strong candidate for a magnetic cv .
the balmer lines are comparatively narrow , and a bump on the red wing of h@xmath23 could tentatively identified with cii ( see also section [ v908oph_sec ] ) . the three individual spectra do not show any significant radial velocity variations , and thus , like v908 oph , this object is likely to be seen at a rather low system inclination .       line profiles for the three individual spectra of v1583 sgr in sequence from bottom to top .
the vertical line indicates the rest wavelength of h@xmath23 . ]
an examination of archival photographic plates by @xcite led to the discovery of this nova almost 30 yr after its eruption .
still , in contrast to other novae , the maximum brightness of 8.9 mag on 1928 june 24 appears to be well established .
the authors report that two days before that night , the nova was fainter than 15.6 mag and that it fell again below the plate detection limit 75 d after maximum brightness .
the system is classified as a moderately fast nova with a decline time @xmath45 d @xcite .    from its colours ( fig .
[ ubvr_fig ] , table [ ubvr_tab ] ) and its position only 3.5 arcsec off the reported coordinates , we identified a candidate for the nova .
subsequent spectroscopy show a slightly reddish continuum ( fig .
[ allspec_fig ] ) , which is likely the consequence of a comparatively strong interstellar reddening ( @xmath20 = 1.09(04 ) mag ) .
like for most of the other systems , the corrected sed presents a blue slope , with @xmath46 .
superposed are moderately strong emission lines ( table [ eqw_tab ] ) of the balmer and hei series .
similar to v908 oph ( section [ v908oph_sec ] ) , the bowen blend and heii @xmath344686 lines are weak and hardly detectable .
the emission lines are also comparatively narrow , but still the line profile shows some clear asymmetric structure and marked differences between the individual exposures ( fig .  [ v1583sgrlprof_fig ] ) .
we especially notice that the line strength in the third spectrum of that series is significantly diminished , which indicates that the emission source is being obscured .
a possible explanation for this apparent discrepancy between narrow lines and an obscuration effect is that the system is seen at high inclination , but that the doppler - broadened component of the line profile produced by the disc is weak and that the line profile is dominated by a narrow component that originates on the secondary star . in any case
, the observed obscuration makes v1583 sgr an attractive target for time series photometry .
line profiles for the three individual spectra of v3964 sgr in sequence from bottom to top .
the vertical line indicates the rest wavelength of h@xmath23 . ]    in , we had observed a candidate for this nova that turned out to be a k1v star .
however , examination of the 2d spectra showed the possibility of h@xmath23 emission from a different star in the slit , and a likely source of this emission was identified . here
, we present further spectroscopic observations of that target .
the spectrum in fig .
[ allspec_fig ] clearly reveals it as the post - nova system , presenting the typical hydrogen and helium emission lines , as well as a weak bowen blend .
the continuum has a red slope , but like in the other objects , this is due to interstellar reddening .
still , even after the correction for @xmath20 = 0.56(01 ) mag , the slope appears significantly flatter than in the other systems , being represented by an exponent @xmath23 = 1.43(03 ) .
the emission lines are not particularly strong ( table [ eqw_tab ] ) , but broad , and even show , at times , a hint of a double - peaked disc profile ( fig .
[ v3964sgrlprof_fig ] ) .
furthermore , like in v1583 sgr ( section [ v1583sgr_sec ] ) , the emission line in the third spectrum is significantly diminished compared to the other two spectra . on the other hand ,
a possible doppler shift of the line is less clear .
still , the width of the line , the potential obscuration effect and the flat slope of the corrected continuum suggest that v3964 sgr is seen at a comparatively high inclination .
thus , it should be feasible to determine the orbital period with time series photometry .      together with v3964 sgr ( section [ v3964sgr_sec ] )
, we had also reported on candidates for the novae v1301 aql and v1151 sgr that turned out to be misidentifications .
for both targets , we had indicated alternative candidates in the immediate vicinity of the original targets . in the current observing run
, we incorporated some of these new candidates , namely the objects marked as ` w1 ' and ` w2 ' for v1301 aql , and the one marked as ` 1 ' in the case of v1151 sgr .
all objects show rather unremarkable spectra with reddish continua and balmer lines in absorption . `
1 ' turned out to be a close visual binary , but neither component presented any emission lines .
these two post - novae thus remain unidentified .
likely , additional photometric data will be needed to select further candidates for spectroscopy .
@lllllllllllllllll object & @xmath47@xmath16 & @xmath48@xmath16 & @xmath49 & @xmath50 & @xmath41 & @xmath20 & @xmath23 & @xmath51 & & @xmath52 + & ( mag ) & ( mag ) & ( mag ) & ( yr ) & ( d ) & ( mag ) & & (  ) & (  ) & ( km s@xmath38 ) & ( arcsec ) + el aql & 6.4@xmath53 & 20.9@xmath5 & 14.5 & 88 & 25 & 1.11 & 3.12 & 47 & 27 & 1230 & 3 + v606 aql & 6.8@xmath53 & 20.4@xmath5 & 13.6 & 116 & 65 & 0.35 & 2.02 & 26 & 28 & 1280 & 2 + v908 oph & @xmath189@xmath54 & 20.5@xmath5 & @xmath1711.5 & 61 & @xmath5530 & 1.26 & 3.38 & 3 & 12 & 550 & 234 + v1149 sgr & 7.4@xmath53 & 18.4@xmath5 & 11.0 & 70 & @xmath18270 & 0.39 & 2.41 & 6 & 21 & 960 & 11 + v1583 sgr & 8.9@xmath53 & 20.5@xmath5 & 11.6 & 87 & 37 & 1.09 & 3.27 & 26 & 17 & 780 & 4 + v3964 sgr & @xmath188.5@xmath53 & 19.0@xmath6 & @xmath1711.5 & 40 & 32 & 0.56 & 1.43 & 12 & 29 & 1320 & 3 +     + @xmath16 @xmath53 : photographic , @xmath5 : @xmath5 band , @xmath6 : @xmath6 band , @xmath54 : near - infrared .
+    in table [ prop_tab ] we have collected several distinguished properties of the here presented post - novae .
these are the maximum brightness @xmath47 , the brightness of the post - nova @xmath48 at the time of our observations and the resulting eruption amplitude defined as @xmath56 . in the next two columns , we give the time @xmath50 that has passed between the nova eruption and our spectroscopic data and the time @xmath41 by which the nova had declined to 3 mag below maximum brightness .
the parameters in the subsequent two columns are associated with the determination of the continuum slope , representing the catalogued value for the interstellar reddening @xmath20 and the negative exponent @xmath23 of the fitted power law ( section [ results_sec ] ) .
the next three columns give selected properties of the h@xmath23 emission line : the equivalent width @xmath51 and the fwhm of a gaussian fit to the line , the latter both in wavelength and in velocity units .
the final column states the deviation of the position of an object from the originally recorded coordinates of the nova eruption as given by @xcite .    since the absolute maximum brightness of novae shows a scatter of typically only 23 mag @xcite , but the absolute brightness of cvs in quiescence is distributed over a range of @xmath119 mag ( e.g. , * ? ? ?
* ) , differences in @xmath49 will mostly reflect the difference in the brightness of the post - nova .
the latter principally depends on the brightness of the accretion disc ( or , more generally , on the emission from the accretion process ) , and thus on the mass - transfer rate @xmath57 and on the system inclination @xmath58 @xcite . because @xmath47 can be assumed to be largely independent of these two quantities , @xmath49 can be taken as an indicator of the absolute brightness of the post - nova and of the involved parameters .
apart from the spread in the absolute maximum brightness , additional uncertainties are introduced by the possibility that the recorded @xmath47 does not correspond to the real maximum brightness , but that the nova was observed at some later point , and by the different passbands used for @xmath47 and for @xmath48 . while in general , the colour differences in the optical range for cvs will amount to less than 1 mag ( e.g. , * ? ? ?
* ) , these differences will be enhanced in regions with strong interstellar reddening .
a thorough analysis of @xmath49 will thus have to correct for these effects .
however , here we only aim at a rough classification of our post - novae within the cv subtypes , and we will do so by also considering other parameters , thus reducing the weight of above uncertainties .
further quantities of table [ prop_tab ] that are indicative of @xmath57 are @xmath23 and @xmath51 .
the significance of the former and the involved uncertainties have already been presented in section [ results_sec ] .
the indicative power of the latter stems from low-@xmath57 , optically thin discs producing stronger emission lines than high-@xmath57 , optically thick ones @xcite .
however , in cvs , typically the accretion disc is not the sole contributor to the emission line profile @xcite . because post - novae contain an eruption heated white dwarf , irradiation of the secondary star is likely to produce an additional emission component .
furthermore , there is also the possibility that the ejected material still contributes to the line profile especially in ` younger ' , i.e.  more recent post - novae . finally , because the continuum emission depends on @xmath58 , so does @xmath51 @xcite , and thus in our interpretation , we will also use the width of the line fwhm , which is indicative of @xmath58 due to the larger projected velocity values @xmath59 sampled in the line profile .
el aql and v606 aql are the novae with the largest eruption amplitudes in the present sample . indeed , a comparison with fig.7 in shows that there are only @xmath1110 novae in our sample of pre-1986 eruptions with similar large @xmath49
the idea of this parameter being an indicator for the disc brightness in these two systems fits well with the comparatively strong emission lines . on the other hand , the lines are also broad , and thus additionally , the inclination will have contributed to the high value of @xmath49 . while , for v606
aql , also @xmath23 is in agreement with a comparatively low-@xmath57 disc , this is not the case for el aql .
however , we note that the steepest continuum slopes are found for the systems with the largest interstellar reddening .
this appears suspicious and suggests a tendency to overcorrect for this effect .
v908 oph and v1149 sgr have the smallest @xmath49 and are also spectroscopically consistent with an optically thick disc and thus high @xmath57 in that they present only weak emission lines . in v908 oph ,
the bluer hydrogen lines are even embedded in absorption troughs .
still , also here the line width suggests that the inclination is such that it amplifies that effect . in conclusion , because of the strong dependence of @xmath49 and @xmath19 on @xmath58 , those four post - novae could have very similar intrinsic properties with respect to @xmath57 in spite of their different spectroscopic appearance .
the picture is somewhat different for the remaining two systems .
v1583 sgr presents narrow , but comparatively strong emission lines with a variable line profile and medium @xmath49 .
we interpret this as the combination of a broad but weak disc line profile and a narrow component that likely originates in the secondary star . in this way
, all parameters would agree with v1583 sgr having a medium to high @xmath49 and a medium to high inclination .
assuming that the narrow component is produced by irradiation , we would also expect that either the white dwarf in this system is still very hot or that the orbital period is comparatively short .
because of the ` age ' of the post - nova ( @xmath50 = 87 yr ) , the latter appears more likely . finally , v3964 sgr has broad and weak lines and shows a possible obscuration effect in one spectrum .
this points to a high-@xmath57 system seen at high inclination .
the recorded medium @xmath49 represents only a lower limit that , if the speculations by @xcite are correct , could be up to 3 mag larger .
thus , it would also be in agreement with a high @xmath58 .
last , not least , the continuum slope appears rather flat for the suspected high @xmath57 , but again might reflect a high system inclination .    in , we have argued that there is a high probability that the cvs detected by us do , in fact , represent the post - novae corresponding to the respective nova eruption and not different cvs that coincidentally have similar coordinates .
our first point was that the low space density of cvs ( * ? ? ? * and references therein ) makes it unlikely that more than one cv could be present in a field of a few square arcmin and a limiting magnitude of @xmath6024 mag , although a proper quantitative treatment would have to take into account a detailed model of the galactic structure .
as a second point of support to our argumentation , we estimated our detection probability of a specific post - nova system .
this is based on the limiting magnitude @xmath61 of our photometry , calculated as the faintest object that is detected in @xmath0 and has s / n = 5 in @xmath5 , corresponding to an uncertainty in the @xmath5 brightness of 0.198 mag .
this gives us the maximum amplitude @xmath62 that a nova with a brightness at maximum @xmath47 can have to still be detected in that specific photometry .
finally , the fraction of novae with @xmath63 translates into the detection probability .
this probability is not fixed , but depends on our current knowledge of the @xmath49 distribution . as of 2016
may 15 , we count 92 pre-1986 post - novae with known @xmath49 , with the shape of the distribution suggesting that it is biased to bright , small-@xmath49 post - novae .
our estimate furthermore ignores the possibility that a nova could have very unusual colours and thus is not recognized as such .
table [ prop_tab ] shows that four of our objects were detected within 5 arcsec of the reported eruption and thus in its immediate vicinity .
we estimate the possibility of a different cv being present within this radius to have a very low probability , making it unnecessary to additionally calculate the detection probability for them .
v1149 sgr is slightly farther off , but only v908 oph shows a very large discrepancy , due to the low precision of the coordinates of the nova eruption . for the photometries of v1149 sgr and of v908 oph
we find @xmath61 = 23.0 and 20.5 mag , and thus @xmath64 = 15.6 and 11.5 mag , respectively .
the much lower brightness and @xmath64 limit for v908 oph is not due a lower quality of the photometry but a result of the strong interstellar reddening causing less than 1 per cent of the objects detected in @xmath10 being also detected in @xmath3 .
the corresponding detection probabilities then calculate to 99 and 55 per cent . additionally taking into account the small positional deviation to the reported coordinates
, we can thus be confident that v1149 sgr is indeed the post - nova . for v908 oph
, the case is less clear , and both parameters leave room for speculation . here ,
basically the only argument supporting the identification of this system with the post - nova of nova oph 1954 is the absence of other candidates in the field .
1 .   we have used @xmath0 photometry and optical spectroscopy to recover the post - novae el aql , v606 aql , v908 oph , v1149 sgr , v1583 sgr and v3964 sgr .
we furthermore report on spectroscopy of candidates for the post - novae v1301 aql and v1151 sgr , 2 .
historic observations of el aql were collected to build a complete eruption light curve that allows it to be classified as a class p nova @xcite .
its eruption amplitude and colours indicate that the previously suspected possibility that el aql is a recurrent nova , is based on source confusion .
furthermore , its spectroscopic features suggest that the object belongs to the group of intermediate polars .
3 .   comparison of several properties of the novae indicates that the differences in the spectroscopic appearances of el aql , v606 aql , v908 oph and v1149 sgr regarding the aspect of the mass - transfer rate are likely to be an effect of the system inclination and thus do not allow for a conclusive classification .
the shape of the h@xmath23 emission line profile in v1583 sgr indicates that it is dominated by a non - disc emission component , likely originating in the irradiated donor star .
inspection of the emission line profiles of the individual spectra proposes el aql , v606 aql , v1583 sgr and v3964 sgr as good candidates for time series observations in order to determine the orbital period .
we have calculated the probability that the detected objects are , in fact , the post - novae corresponding to the reported eruptions .
we have found that the only potentially ambiguous case is that of v908 oph , where the low precision of the reported nova coordinates and the strong interstellar reddening strongly reducing the number of objects detected in @xmath3 leave some room for doubt .
still , in the absence of other candidates , the here presented object should be regarded as the post - nova .
this paper is based on observations with gemini - south , program ids gs-2015a - q-54 and gs-2015a - q-75 , and eso telescopes , proposal numbers 089.d-0505(a ) and 092.d-0225(a ) .
we are indebted to the gemini and eso astronomers who performed the service observations .
eso furthermore provided office space and infrastructure for the two week stay of ae in chile , during which large parts of the work were discussed and prepared .
many thanks for the kind hospitality !    this research was supported by fondecyt regular grant 1120338 ( ct and nv ) .
ae acknowledges support by the spanish plan nacional de astrononoma y astrofsica under grant aya2015 - 66211-c2 - 1 .
ct , nv and ifm acknowledge support by the centro de astrofsica de valparaso .
we gratefully acknowledge ample use of the simbad data base , operated at cds , strasbourg , france , and of nasa s astrophysics data system bibliographic services .
iraf is distributed by the national optical astronomy observatories .
the guide star catalogue - ii is a joint project of the space telescope science institute and the osservatorio astronomico di torino .
space telescope science institute is operated by the association of universities for research in astronomy , for the national aeronautics and space administration under contract nas5 - 26555 .
the participation of the osservatorio astronomico di torino is supported by the italian council for research in astronomy .
additional support is provided by european southern observatory , space telescope european coordinating facility , the international gemini project and the european space agency astrophysics division .
gemini observatory is operated by the association of universities for research in astronomy , inc .
, under a cooperative agreement with the nsf on behalf of the gemini partnership : the national science foundation ( united states ) , the national research council ( canada ) , conicyt ( chile ) , the australian research council ( australia ) , ministrio da cincia , tecnologia e inovao ( brazil ) and ministerio de ciencia , tecnologa e innovacin productiva ( argentina ) .
58 natexlab#1#1    i. et  al .
, 1998 , the messenger , 94 , 1    n. , strashny g. , 1928 , astron .
, 234 , 213    m. , 1929 , astronomische nachrichten , 235 , 427    v.  m. , 1954 , , 1461 , 1    f. , 1931 , astronomische nachrichten , 244 , 207    a.  j. , 1927 , harvard college observatory bulletin , 851 , 9    j.  a. , clayton g.  c. , mathis j.  s. , 1989 , , 345 , 245    m. , livio m. , 1995 , , 452 , 704    j. , hoffleit d. , 1955 , , 60 , 259    r.  a. , webbink r.  f. , shara m.  m. , ritter h. , kolb u. , duerbeck h.  w. , 2005 , journal of astronomical data , 11 , 2    h.  w. , 1981 , , 93 , 165    h.  w. , 1987 , , 45 , 1    t.  e. , 1992 , , 259 , 17p    m. , 1927 , harvard college observatory bulletin , 851 , 10    c. , 2001 , cataclysmic variable stars .
springer    i.  m. , jrgensen i. , allington - smith j.  r. , davies r.  l. , metcalfe n. , murowinski r.  g. , crampton d. , 2004 , , 116 , 425    m.  m. , cenko s.  b. , kulkarni s.  r. , ofek e.  o. , quimby r. , rau a. , 2011 , , 735 , 94    a. , prialnik d. , shara m.  m. , 1988 , , 325 , 828    a.  u. , 1983 , , 88 , 439    a.  u. , 1992 , , 104 , 340    b.  m. et  al . , 2008 , , 136 , 735    h.  s. , 1933 , annals of harvard college observatory , 84 , 121    i. , stenholm b. , 1976 , , 2997 , 2    i. , stenholm b. , 1977 , information bulletin on variable stars , 1351 , 1    m.  w. , 1949 , , 54 , 191    p. et  al . , 2015 , , 219 , 26    a. , schaefer b.  e. , 2014 , , 788 , 164    j. , 1984 , , 54 , 443    j. et  al . , 2013 , , 434 , 1902    e.  c. , 1900 , , 12    m. , 2014 , in the x - ray universe 2014 , p. 164
d. , shara m.  m. , 1986 , , 311 , 172    f.  a. , naylor t. , mukai k. , 1996 , , 281 , 192    r.  k. et  al .
, 2013 , , 554 , a123
g. , nikonoff w. , 1928 , astronomische nachrichten , 233 , 221    e.  f. , finkbeiner d.  p. , 2011
, , 737 , 103    c.  k. , yoss k.  m. , 1954 , harv .  coll .
card , 1260    m.  m. , livio m. , moffat a.  f.  j. , orio m. , 1986 , , 311 , 163    p.  b. , 1992 , in astronomical society of the pacific conference series , vol .  25 ,
astronomical data analysis software and systems i , worrall d.  m. , biemesderfer c. , barnes j. , eds .
, p. 297
p.  b. , 2000 , , 112 , 925    r.  j. , schaefer b.  e. , henden a.  a. , 2010 , , 140 , 34    p. , 1994
, , 108 , 639    c. , ederoclite a. , mennickent r.  e. , schmidtobreick l. , vogt n. , 2012 , , 423 , 2476    c. , hanuschik r. , 2001 , lnp vol .  573 : astrotomography , indirect imaging methods in observational astronomy , 573 , 119    c. , schmidtobreick l. , vogt n. , ederoclite a. , 2013 , , 436 , 2412    c. , vogt n. , della valle m. , schmidtobreick l. , ederoclite a. , 2014 , , 442 , 565    c. , vogt n. , schmidtobreick l. , ederoclite a. , 2015 , , 450 , 943 ( paper v )    t. , swierczynski e. , mikolajewski m. , ilkiewicz k. , 2015 , , 576 , a119    r. , 1981 , , 31 , 127    j. , 1927 , , 4 , 106    b. , 1986 , , 222 , 11    b. , 1987 , , 227 , 23    b. , 1995 , cataclysmic variable stars .
cambridge astrophysics series , cambridge , new york : cambridge university press    j. , 1965 , aavso abstracts , 8    m. , 1927 , astronomische nachrichten , 230 , 421    a.  b. , 1940 , publications of lick observatory , 14 , 217    n. , finch c.  t. , girard t.  m. , henden a. , bartlett j.  l. , monet d.  g. , zacharias m.  i. , 2013 , , 145 , 44    w. , 1928 , astronomische nachrichten , 231 , 383
[ lastpage ]


classical spin liquids @xcite , such as the coulomb phase @xcite in spin ice @xcite and related systems , are examples of phases whose behavior is not captured by the standard landau picture of broken symmetries @xcite .
their two defining characteristics are fractionalization , the emergence of excitations not constructed from finite combinations of the elementary degrees of freedom , and topological order , the presence of structure that can only be discerned by observing the system globally @xcite .    in the case of the classical spin ices
@xcite , a family of magnetic pyrochlore oxides , the consequences of these properties have been of sustained interest , from both theoretical and experimental perspectives .
the spins in these materials carry large magnetic moments , and their fractional excitations take the form of magnetic monopoles @xcite , acting as sources for the physical magnetic field .
furthermore , transitions between phases where these excitations are confined and deconfined have particularly interesting properties that can not be captured by the landau  ginzburg theory of phase transitions @xcite .
spin ice also has the unusual feature that the magnetization is a topologically constrained quantity , and so topological order can be probed directly @xcite . in particular , within the low - energy configuration space relevant at low temperatures , any local dynamics conserves the uniform magnetization @xcite .
this fact is known to have interesting consequences for critical properties at certain confinement transitions , such as the kasteleyn transition in an applied magnetic field @xcite .    in this work ,
we investigate the interplay between these two aspects of the coulomb spin liquid by studying confinement phase transitions in spin ice across the full range of magnetization . as noted by henley @xcite , one can distinguish three categories of phase transition from the coulomb phase : the first , where the magnetization is zero across the transiton , was studied systematically in ref .
@xcite . a second , where the ordered state has saturated polarization ,
was considered in refs .
here we consider the phase structure and transitions more generally , including the third case , where the magnetization is nonzero but unsaturated , and may or may not change across the transition .    to do so , we introduce a mean - field theory ( mft ) designed to capture the distinction between confined and deconfined phases .
standard mean - field approaches , based on the landau picture of ordered phases , are not well suited to describing confinement transitions .
instead , we take inspiration from the well - known mft for the bose  hubbard model @xcite , making use of a mapping from spin ice to a quantum model of bosons @xcite .
we apply the method to a simplified model of classical spin ice , in order to illustrate the general approach and its physical interpretation .
we also briefly discuss extensions to the method that could be used to treat more physically realistic perturbations .
in addition , we present arguments for the generality of our results , beyond mean - field theory and our particular choice of hamiltonian .
we argue that the phase diagram generically consists of a set of confined phases at low temperature in which the magnetization is , to a very good approximation , fixed .
this plateau structure is strongly reminiscent of the lobes present in the phase diagram of the bose
hubbard model @xcite , a connection that helps clarify the general phase structure .
while the focus here is on classical spin ice in a @xmath0 $ ] field , the distinction between transitions in different flux sectors is more general .
for example , spin ice in a field along the @xmath1 $ ] direction also exhibits magnetization plateaux @xcite , while it has been suggested that the zero - field ground state of quantum spin ice may have nonzero magnetization @xcite .    in section  [ secmodel ]
, we introduce the model to be used and briefly outline the general properties of the coulomb phase and confinement transitions .
the mean - field method is then described and applied to the model in section  [ secmeanfieldtheory ] .
we present general arguments for the phase structure and the nature of the phase transitions in section  [ secphasestructure ] , before concluding in section  [ secconclusions ] .
our presentation will be based on the case of classical spin ice in a magnetic field applied along the @xmath0 $ ] crystallographic direction .
we first introduce the nearest - neighbor model of spin ice ( nnsi ) that describes the physics of the coulomb spin liquid , before discussing perturbations , such as an applied field , and the resulting ordering transitions .
the spin ices @xcite are a family of frustrated magnetic materials with moments on the sites @xmath2 of a pyrochlore lattice , a network of corner - sharing tetrahedra .
prominent examples are the `` classical '' spin ices , such as @xmath3  and @xmath4 , which are well described in terms of classical spins @xmath5 .
each spin is subject to strong easy - axis anisotropy along the local @xmath6 direction @xmath7 which joins the centres of the two tetrahedra sharing site @xmath2 , as shown in fig .
[ figtetrpair ] .
( shown as green arrows ) lie at the sites @xmath2 of the lattice and are constrained to point along the local @xmath6 axes ( dashed lines ) that join the centers of adjacent tetrahedra.,scaledwidth=35.0% ]    the anisotropy results from a large ( of order ) crystal - field splitting between states with maximum projection along @xmath7 and all other states in the single - moment hilbert space .
this large gap ensures that tunneling processes between the low - energy states are strongly suppressed , rendering the moments effectively classical @xcite .
for the low temperatures of interest , it is therefore sufficient to treat the spins as effectively constrained to @xmath8 , where @xmath9 is a classical ising variable .
each tetrahedron @xmath10 in the pyrochlore lattice can be labeled according to its orientation @xmath11 , and the lattice structure is such that all neighbors of @xmath10 have orientation @xmath12 .
the fixed unit vector @xmath7 for each pyrochlore site @xmath2 , directed along the line joining the centers of its two tetrahedra , can be chosen to point towards the tetrahedron with @xmath13 , and so @xmath14 indicates whether @xmath5 points out of ( @xmath15 ) or into ( @xmath16 ) tetrahedron @xmath10 . with this convention , @xmath17 and hence @xmath18 , for all nearest neighbors @xmath2 and @xmath19 .
the minimal hamiltonian that captures the spin - liquid physics of the classical spin ices @xcite contains only interactions between nearest - neighbor pairs @xmath20 , @xmath21 the coupling between moments @xmath5 is ferromagnetic ( @xmath22 ) , incorporating the net effect of the dipolar interactions between nearest neighbors and ( antiferromagnetic ) exchange @xcite . the effective interactions between the ising variables @xmath23 are therefore antiferromagnetic , and hence frustrated .
further - neighbor interactions , both dipolar and exchange , are significant in spin ice , but can be treated as a relatively small perturbation to @xmath24 @xcite ; we will return to these in section  [ secperturbations ] .    because two sites are nearest neighbors if and only if they share a tetrahedron , the interactions can be rewritten as @xmath25 where @xmath26 includes all tetrahedra ( of both orientations ) and @xmath27 is referred to as the ( fictional ) `` charge '' on tetrahedron @xmath10 in terms of a sum over its sites @xmath28 .
the energy is therefore minimized by configurations with @xmath29 , which are those where two spins point into and two out of each tetrahedron .
configurations that satisfy this condition are said to obey the _ ice rule _
@xcite ; all configurations that do so at every tetrahedron are degenerate minimal - energy states of @xmath24 .    tetrahedron configurations with three spins in and one out , or vice versa , have energy @xmath30 above the ice - rule configurations ; they have charge @xmath31 and are hence referred to as `` monopoles '' .
( note that @xmath32 , and so all configurations are globally charge - neutral . ) for @xmath3 , the effective nearest - neighbor coupling is @xmath33 @xcite , corresponding to a monopole cost of @xmath34 .
the corresponding boltzmann weight @xmath35 is called the monopole fugacity . while the significance of the term `` monopole '' is here
only that they carry effective `` charge '' @xmath36 , it has been argued that monopole excitations in the classical spin ice materials in fact carry true magnetic charge @xcite .    for temperature @xmath37 significantly below @xmath38 ,
the system is effectively constrained by the ice rule , and the ensemble of ice - rule configurations therefore determines its behavior @xcite .
the number of such configurations grows exponentially with the total number of spins @xmath39 , and so the entropy density remains nonzero even well below @xmath38 , the temperature scale at which an unfrustrated system would order . in spite of the absence of order , the ensemble exhibits correlations that decay only algebraically with distance .
the resulting correlated paramagnet is referred to as a coulomb phase @xcite , and is an example of a classical spin liquid @xcite .
the long - wavelength properties of the coulomb phase can be captured by coarse - graining the vector spins @xmath40 to give a continuum field @xmath41 referred to as the effective `` magnetic field '' @xcite .
the ice - rule constraint applied to @xmath40 implies that @xmath42 is divergenceless , while monopoles act as sources or sinks , according to an effective gauss law .
the resulting predictions for spin
spin correlation functions are in good agreement with neutron - scattering experiments in classical spin ice @xcite
.    the ice - rule configurations also display interesting topological properties @xcite . for our purposes ,
the most important is that any pair of ice - rule configurations is connected by flipping all spins along one or more loops . for a system with periodic boundary conditions ( which we will assume ) ,
such a loop can only change the total magnetization @xmath43 if it has nontrivial winding number . for other systems that host coulomb phases @xcite , such as classical dimer models , an analogous quantity , called the `` flux '' ,
can be defined , but does not necessarily correspond to a quantity that is so directly accessible experimentally .
a second important property of the coulomb phase is that monopoles , defects in the ice rule , are deconfined @xcite . to make this statement precise ,
imagine imposing a pair of monopoles with opposite charges @xmath44 at tetrahedra @xmath10 and @xmath45 in a background that otherwise obeys the ice rule .
let @xmath46 denote the set of spin configurations @xmath47 compatible with this charge configuration .
their number @xmath48 is a function of the separation @xmath49 that decreases with @xmath50 but reaches a finite limit as @xmath51 .
the result is an effective entropic interaction @xmath52 between the charges that is bounded , and hence allows them to be separated to infinity at finite free - energy cost .
in fact , since monopoles are charges in the continuum field @xmath42 , coarse - graining predicts that the interaction obeys the coulomb law , @xmath53 , for large @xmath54 @xcite .
more generally , consider the partition function @xmath55 where @xmath56 is a perturbation that splits the degeneracy of the ice - rule configurations .
this will tend to suppress fluctuations , and can eventually , as the temperature is reduced compared to the scale of the perturbation , drive the system out of the coulomb phase . a transition into a phase where @xmath57 increases without limit as @xmath51 , and so @xmath58 ,
is referred to as a confinement transition @xcite .
( in terms of the effective continuum field @xmath42 , gauss law implies that an imposed pair of charges must be joined by a fixed quantity of flux .
linear confinement , where @xmath59 for large @xmath60 , occurs when the flux forms a narrow tube connecting the monopoles , with finite tension . )
while confinement may occur simultaneously with appearance of conventional symmetry - breaking order @xcite , there is no local order parameter that provides a signature of confinement .
the confinement criterion , that @xmath61 approaches zero as @xmath51 in a confined phase and conversely has a finite limit in a deconfined phase , instead resembles the characterization of a conventional symmetry - broken phase in terms of long - range order .
this connection can be pursued further by considering the partition function with an unconstrained charge distribution and a position - dependent fugacity @xmath62 , @xmath63 = \sum_{\{\sigma_i\ } } { \mathrm{e}}^{-v / t } \prod_t \zeta_t^{\lvert q_t \rvert}{\,{\text{,}}}\ ] ] in terms of which @xcite @xmath64 \right\rvert_{\zeta_\cdot = 0}{\,{\text { , } } } \label{eqmonopoledistributionfunction}\ ] ] where @xmath65\right\rvert_{\zeta_\cdot = 0}$ ] is the partition function restricted to ice - rule configurations .
the relationship between the monopole distribution function @xmath66 and the fugacity @xmath67 expressed by eq .
( [ eqmonopoledistributionfunction ] ) is exactly analogous to the fluctuation
dissipation theorem relating , for example , a spin
spin correlation function and an applied field .
this connection motivates a mean - field theory , analogous to standard weiss theory @xcite , that describes the confinement transition in terms of an effective self - consistent fugacity .
it should be noted that the confinement distinction is only sharp in the limit @xmath68 , because any nonzero density of monopoles screens the effective interaction @xmath69 at large separation .
there are nonetheless implications for the critical properties away from this limit @xcite .
if the phase transition also involves a different type of order ( such as spontaneous symmetry breaking ) , then it may survive for @xmath70 , but with a different universality class @xcite .      within the nearest - neighbor model @xmath24 , all configurations that obey the ice rules are degenerate and the system exhibits a deconfining coulomb phase .
to drive a confinement transition , one must add perturbations that split the degeneracy of the ice - rule configurations .
here we consider a hamiltonian of the form @xmath71 where @xmath72 is the coupling to an applied magnetic field and @xmath73 contains additional interactions described below .
coupling to a magnetic field @xmath74 is described by a zeeman term @xmath75 where a factor dependent on the size of the magnetic moments has been included in the definition of @xmath74 .
the case that we consider here has @xmath76 where @xmath77 is a unit vector along @xmath0 $ ] , and so @xmath78 . the coupling can therefore be written as @xmath79 where @xmath80 . using the choice described in section  [ secnnsi ] , where @xmath7 points towards the tetrahedron with @xmath13 , @xmath81 alternates in sign on successive @xmath82 planes of spins .
the component of the magnetization along the field can likewise be expressed as @xmath83    the perturbation @xmath72 alone can drive a confinement transition , referred to as a kasteleyn transition @xcite . for @xmath84 ,
all spins align with the field ( to the extent allowed by the easy - axis anisotropy ) , an arrangement that satisfies the ice rules .
any excitation within the ice - rule states involves flipping a loop of spins , but starting from the fully polarized state , the only available loops span the system in the @xmath0 $ ] direction .
such excitations therefore have energy cost proportional to the linear system size @xmath85 and are suppressed in the thermodynamic limit .
the result is a strictly saturated magnetization @xmath86 even for nonzero @xmath87 ( but in the limit @xmath88 ) . above a temperature @xmath89 ,
however , the energy cost is outweighed by the gain in entropy , also @xmath90 , associated with the multitude of available paths for a spanning loop , and a kasteleyn transition occurs from the saturated paramagnet to the coulomb phase .
the transition is possible only because the ice rule restricts to loop - like excitations ; away from the limit @xmath91 , the transition is replaced by a crossover @xcite .    to study the interplay between the topological constraints on the magnetization , inherent in the ice - rule configurations , and conventional ordering transitions , we also consider an additional interaction @xmath73 that induces spontaneous symmetry breaking
. a natural choice would be the further - neighbor coupling present in the spin - ice compounds , due to a combination of further - neighbor exchange @xcite and long - range dipolar interactions .
the latter have been shown using monte carlo @xcite to lead to a phase transition into the ordered configuration shown in fig .
[ figmdg ] , which we will call the melko  den hertog  gingras ( mdg ) state @xcite .     and ordering wavevector along @xmath0 $ ] ( vertical ) .
the spins in each @xmath82 plane are aligned ( to the extent allowed by the local easy axes ) , and between successive planes the magnetization rotates in the right - handed sense by @xmath92 .
the perturbation @xmath73 reduces the ferromagnetic coupling of the red bonds by @xmath93 , and hence favors the configuration shown , or the one with all spins inverted.,scaledwidth=85.0% ]    the actual perturbation @xmath73 that we will use here is one that is simpler to treat within our mean - field theory but that is expected to lead to an ordered phase of the same type .
it is given by @xmath94 where @xmath95 is the set of nearest - neighbor bonds highlighted in fig .
[ figmdg ] , which form a set of left - handed screw chains .
the effect of @xmath96 is to reduce the strength of the interactions on these bonds and to favor the mdg configuration shown , or the one with all spins inverted .
[ figtetrgrid ] shows the energy for each of the six configurations of a single tetrahedron that obey the ice rule , including both the applied field and the perturbation @xmath73 .    ,
defined in eq .
( [ eqhu ] ) , as well as an applied field @xmath97 along @xmath0 $ ] ( upwards ) .
the red bonds , as in fig .
[ figmdg ] , are those where the ferromagnetic coupling is reduced by @xmath93 .
the configurations favored by this perturbation are therefore those in the right column where the spins joined by the red bond are antialigned . ]    unlike dipolar interactions , the perturbation @xmath73 reduces the symmetry of the lattice by picking an axis and a chirality .
the 12-fold degeneracy of the mdg state is therefore reduced by a factor of @xmath98 , but the symmetry under inversion of all spins remains , and so an ordering transition into this state spontaneously breaks an ising symmetry .
this transition ( at @xmath99 ) is in fact a type considered in ref .
@xcite , where it was argued to be most likely of first order .
as we will show , this simplified model has a phase diagram that illustrates general features emerging from the interplay between deconfinement , ordering , and topological constraints . in section  [ secconclusions ] , we will briefly discuss the prospects for including more realistic perturbations in the mean - field theory .
some aspects of the phase structure can be elucidated using a mapping between classical spin ice and an effective quantum model of hard - core bosons on a lattice @xcite .
this proceeds by identifying the @xmath0 $ ] direction , along which the external field is applied , with the imaginary - time axis in a path - integral representation of the quantum problem .
the thermodynamic limit of the original classical system therefore corresponds to the zero - temperature limit in the quantum problem .
this mapping has previously been studied in detail for the hamiltonian @xmath100 @xcite .
the fully polarized state ( along the @xmath101 direction , say ) is identified with the vacuum ; flipped spins relative to this configuration form system - spanning loops , which are viewed as the world lines of bosons .
the result is an effective model of hard - core bosons on a lattice with density @xmath102 and hence chemical potential @xmath103 . in general
, one expects all other local terms compatible with symmetry , such as hopping between nearby sites and further - neighbor interactions , with parameters related in nontrivial ways to those of the spin model @xcite .    as @xmath104
is increased , the system passes from the vacuum through a quantum phase transition @xcite to a superfluid , which corresponds to the coulomb phase of spin ice .
( off - diagonal long - range order in the superfluid can be associated with deconfinement in the coulomb phase @xcite . ) for sufficiently large positive @xmath104 , the system is driven into a fully occupied mott insulator , with @xmath105 , corresponding to a saturated paramagnet with @xmath106 in an inverted applied field .
critical properties at the phase transitions in spin ice are , _ mutatis mutandis _ , equivalent to those in the quantum model @xcite .
the presence of interactions between bosons on different lattice sites also makes possible mott insulators with fractional filling .
these can be identified with confining phases of spin ice ( because they lack superfluid order ) in which the magnetization is fixed to a value less than @xmath107 .
the mdg phase , for example , is an ordered phase with @xmath108 , and corresponds to a mott insulator at half filling .    in general ,
more exotic phases such as supersolids are also possible in the effective quantum model .
a supersolid would correspond to a coulomb phase with broken spatial symmetry ; such phases have been proposed in spin ice @xcite , but are not captured by the mean - field theory used here .
the phase structure of @xmath109 in eq .
( [ eqhamiltonian ] ) is expected to contain : the coulomb phase , characterized by deconfinement of monopoles ; the saturated paramagnet , in which fluctuations are suppressed by the ice - rule constraints ; and the mdg phase , with conventional order that breaks spin - inversion and spatial symmetries .
a mean - field theory that describes all three must therefore capture the confinement  deconfinement distinction , unlike mean - field approaches based on finite clusters @xcite , and also respect the spatial structure , unlike bethe - lattice calculations @xcite .
here we apply an approach that is analogous to the standard mean - field theory for the bose
hubbard model @xcite , but is applied directly to the classical partition function , written in terms of a transfer matrix .
( for other applications of the variational method to the transfer matrix , see refs .
we start by partitioning the lattice into chains that span the system , which are the minimal units that can obey the ice rule .
an effective model is then defined for each chain , with the coupling between chains treated self - consistently .    in standard mean - field theory for a symmetry - breaking transition @xcite , coupling between lattice sites is replaced by an effective field related to the order parameter . here , in the absence of a local order parameter
, we capture the confinement  deconfinement transition by introducing an effective self - consistent fugacity @xmath110 for monopoles .
since , according to eq .
( [ eqmonopoledistributionfunction ] ) , the monopole distribution function can be interpreted as the correlation function corresponding to a monopole fugacity , a nonzero @xmath110 implies deconfinement .
this decoupling is closely analogous to that used for the bose
hubbard model , where superfluid order is captured by introducing a self - consistent source term for bosons . following the mapping between this quantum problem and the classical statistical mechanics of spin ice ,
we apply the method at the level of the transfer matrix .
it should be emphasized that this differs from the standard mean - field approach for classical statistical systems , where one chooses a trial distribution to minimize the free energy .
( the latter method can not be used , because a trial configuration with nonzero fugacity does not enforce the ice rule and therefore has infinite energy . )
while we treat deconfinement using an effective mean - field fugacity @xmath110 , we set the true monopole fugacity @xmath35 to zero .
nonzero @xmath67 could straightforwardly be included in the approximation and is analogous to standard weiss mean - field theory in the presence of an applied magnetic field @xcite .
we choose to decompose the pyrochlore lattice into right - handed screw chains with axes parallel to the @xmath0 $ ] direction , as illustrated in fig .
[ figmdgchains ] .    $ ] crystallographic direction . in the mdg state ,
the chains are ordered in a checkerboard pattern , with spins along each chain aligned and neighboring chains antialigned.,scaledwidth=85.0% ]    for our purposes , this choice has two advantages : first , the chains span the system in the imaginary - time direction under the mapping to a bosonic model , and so the connection to the bose
hubbard mean - field theory is transparent .
second , the ordered states that we aim to describe have natural interpretations in terms of chain order : in the mdg state , shown in fig .  [ figmdgchains ] , the chains order in a checkerboard pattern , while in the fully polarized paramagnet , all chains are aligned .
other orderings , such as the partially polarized state observed by lin and kao @xcite , can be described in terms of ordered configurations of the same chains , while the general approach can likely be adapted for ordered states that are not captured by this decomposition of the lattice .
( for example , mcclarty et al .  @xcite consider ordered states described in terms of linear chains of spins . )
the division of the lattice into screw chains has the property that each chain @xmath111 has @xmath112 neighbors , which it meets in turn at tetrahedra in successive @xmath82 layers .
the set of tetrahedra at a given layer , at vertical position @xmath113 , amounts to a partition of the chains into pairs , which we denote @xmath114 .
we also define @xmath115 as the chain that meets chain @xmath111 in a tetrahedron at layer @xmath113 .
the pyrochlore lattice is symmetric under translations by @xmath116 layers in the @xmath0 $ ] direction @xcite , and so @xmath117 and @xmath118 .    in the quantum mapping ,
the vertical direction is treated as imaginary time .
applying the mapping at a microscopic level would therefore lead to a time - evolution operator that couples different pairs of lattice sites at successive time steps . this `` stroboscopic '' @xcite implementation of the hopping in the effective quantum hamiltonian has been shown to be irrelevant for the critical properties @xcite , but must be handled with some care when implementing the mean - field theory .      the partition function @xmath119 can be expressed exactly in terms of a set of operators acting on degrees of freedom associated with the chains . to each chain @xmath111
is associated a pair of basis vectors , @xmath120 and @xmath121 , representing the orientation , up or down ( i.e. , @xmath122 ) , of the spin at a given vertical position @xmath113 along the chain . denoting by @xmath123 the space spanned by these vectors ,
the full hilbert space for a layer of the system is @xmath124 .
( note that there is a single two - dimensional hilbert space for each chain , and so the basis vectors are labeled only by their chain @xmath111 , and not by their position @xmath113 along the chain . )
the partition function can then be written as @xmath125 where @xmath85 is the system size ( number of layers ) in the @xmath113 direction and @xmath126 is the transfer operator for the layer of tetrahedra at vertical position @xmath113 .
( we assume periodic boundary conditions in the @xmath113 direction . )
the tetrahedron transfer operator @xmath127 , acting on @xmath128 , accounts for the possible configurations of the chains @xmath111 and @xmath129 on entering and leaving this tetrahedron . using the symmetry of the pyrochlore lattice under translations by @xmath116 layers in the @xmath113 direction , the product in eq .
( [ eqpartitionfunction ] ) can be rewritten as @xmath130^{l_z/\lambda_z}{\,{\text{,}}}\quad\text{where}\quad{\mathcal{t}}^{(\lambda_z)}=\prod_{z=\lambda_z}^{1}{\mathcal{t}}_z { \,{\text{.}}}\ ] ] in the thermodynamic limit , the partition function is therefore given by @xmath131^{l_z/\lambda_z}$ ] , where @xmath132 is the spectral radius ( largest absolute eigenvalue ) .
it should be noted that , even if @xmath133 is hermitian , in general @xmath134\neq 0 $ ] for @xmath135 , and so @xmath136 is not hermitian .    the approximation we propose is to replace the system with an effective model of decoupled chains , whose parameters are determined self - consistently .
the partition function for chain @xmath111 is taken as @xmath137 where the effective transfer operator for chain @xmath111 is @xmath138 an operator on @xmath123 .
the normalized vector @xmath139 , which is interpreted as an ansatz for the configuration of chain @xmath111 immediately above layer @xmath113 , is determined self - consistently as the thermal state of the one - dimensional system defined by eq .
( [ eqpartitionfunctionmf ] ) . assuming a self - consistent solution exists with at least periodicity under translation in @xmath113 by @xmath140 , it is given by the normalized ( right ) eigenvector with largest absolute eigenvalue of @xmath141    for the types of phases that we aim to describe
, it will in fact be sufficient to consider @xmath142 uniform along each chain .
this leads to the simplification that @xmath143 which is therefore hermitian .
for such a solution to exist , @xmath144 must be a simultaneous eigenvector of @xmath145 for all @xmath113 .
for the case with an applied field @xmath74 and a chain - favoring interaction @xmath93 , the tetrahedron transfer operator can be written explicitly as @xmath146 & { \mathrm{e}}^{\frac{2}{\sqrt{3}}h / t } \ket{\uparrow}_c \ket{\uparrow}_{c'}\bra{\uparrow}_{c}\bra{\uparrow}_{c ' } + { \mathrm{e}}^{-\frac{2}{\sqrt{3}}h / t}\ket{\downarrow}_c \ket{\downarrow}_{c'}\bra{\downarrow}_c \bra{\downarrow}_{c'}\\ & { } + { \mathrm{e}}^{2u / t}\ket{\uparrow}_c\ket{\downarrow}_{c'}\bra{\uparrow}_c \bra{\downarrow}_{c ' } + { \mathrm{e}}^{2u / t}\ket{\downarrow}_c\ket{\uparrow}_{c'}\bra{\downarrow}_c \bra{\uparrow}_{c'}\\ & { } + \ket{\uparrow}_c\ket{\downarrow}_{c'}\bra{\downarrow}_c \bra{\uparrow}_{c ' } + \ket{\downarrow}_c\ket{\uparrow}_{c'}\bra{\uparrow}_c \bra{\downarrow}_{c'}{\,{\text{. } } } \end{aligned}\ ] ] each term corresponds to one of the six ice - rules configurations of the tetrahedron , illustrated in fig .
[ figtetrgrid ] , and is associated with a boltzmann weight implementing the terms @xmath72 and @xmath73 in the hamiltonian .
( terms breaking the ice rule could also be included ; we are here treating such configurations as having zero boltzmann weight . )
the first four terms correspond to cases where the two chains maintain their configurations after passing through the tetrahedron  i.e . , the two spins on each chain are aligned  while in the last two terms the chains exchange orientations .
taking the inner product in @xmath147 gives @xmath148 & \left[\frac{1}{2}(1+m{^{\text{mf}}}_{c'}){\mathrm{e}}^{\frac{2}{\sqrt{3}}h / t } + \frac{1}{2}(1-m{^{\text{mf}}}_{c'}){\mathrm{e}}^{2u / t}\right]\ket{\uparrow}_c\bra{\uparrow}_c \\ & { } + \left [ \frac{1}{2}(1-m{^{\text{mf}}}_{c'}){\mathrm{e}}^{-\frac{2}{\sqrt{3}}h / t } + \frac{1}{2}(1+m{^{\text{mf}}}_{c'}){\mathrm{e}}^{2u / t}\right]\ket{\downarrow}_c\bra{\downarrow}_c\\ & { } + \zeta{^{\text{mf}}}_{c'}\ket{\uparrow}_c\bra{\downarrow}_c + \left({\zeta{^{\text{mf}}}_{c'}}\right)^*\ket{\downarrow}_c\bra{\uparrow}_c{\,{\text { , } } } \end{aligned } \label{eqteff}\ ] ] where the expectation values @xmath149 are mean fields describing the magnetization and density of monopoles on chain @xmath150 .
the effective model for chain @xmath111 , described by the transfer operator @xmath151 , involves an applied field that is a function of @xmath152 and an effective monopole fugacity @xmath153 .
this is the crux of the approximation : the exchange of orientation between chains , described by the last two terms of eq .
( [ eqtetrahedrontransferoperator ] ) , has been replaced by terms where a single chain flips orientation , thus creating or destroying a monopole .
the criterion for deconfinement , that a distant pair of monopoles can be inserted with finite free - energy cost , is trivially satisfied whenever @xmath154 .
( in the analogous mean - field theory for the bose  hubbard model @xcite , the expectation value of the bosonic annihilation operator is used as a mean field . )    in general , the solution of the self - consistent approximation requires finding simultaneous eigenvectors of @xmath151 , given in eq .
( [ eqteff ] ) , for all the values of @xmath113 . for the types of ordered phases expected in our model of spin ice , however , some further simplifications are possible . in particular , to describe the mdg and saturated phases , it is sufficient to divide the chains into two `` chain sublattices '' , which we label @xmath155 and @xmath156 , with one chain of each type meeting in every tetrahedron , as illustrated in fig .
[ figmdgchains ] .
the operator @xmath151 is therefore independent of @xmath113 , and it is sufficient to find @xmath157 and @xmath158 such that @xmath144 is the eigenvector of @xmath159 with largest eigenvalue @xcite , for both @xmath160 , @xmath161 and vice versa . because @xmath162 is hermitian , this is equivalent to maximizing @xmath163 with respect to normalized @xmath157 and @xmath158 .    as an aside
, we note that this amounts to maximizing the matrix element of the layer transfer operator @xmath133 in the subspace of factorizable vectors @xmath164 . in these terms
, it is clear how this is related to the mean - field theory for the bose
hubbard model , where one minimizes the matrix element of the hamiltonian @xmath109 .
( recall the standard connection between classical and quantum statistical mechanics , where @xmath165 is the evolution operator in imaginary time @xmath166 . ) in both cases , the operator is hermitian and so the matrix element is bounded by its extremal eigenvalue , but only for the bh model is the hamiltonian constant in time and the matrix element bounded by the true ground - state energy . in the present case , the product of extremal matrix elements @xmath167 provides an approximation to the exact partition function , but this approximation is not variational .    to maximize @xmath168 , we express both vectors as @xmath169 where @xmath170 and @xmath171 , in terms of which @xmath172 the distinct maxima of this expression correspond to phases of the mean - field theory , shown in fig .
[ figphasediagram ] .    , the melko  den hertog
gingras ( mdg ) phase @xcite with @xmath173 , and a saturated paramagnet with @xmath174 . at higher - temperatures is a coulomb spin liquid whose magnetization is a smooth function of the applied field @xmath97 .
( dotted lines are contours of constant magnetization . ) ]    the phase boundaries can be found by expanding in small deviations from the ordered phases and identifying the points where they cease to be stable maxima .
( we have verified numerically that the ordered states maximize the matrix element globally if and only if they are local maxima . )    at low temperature , there are three confined phases with different , fixed , values of the magnetization @xmath175 for @xmath176 , the maximum of eq .
( [ eqtv1v2 ] ) has @xmath177 ( @xmath178 ) or @xmath179 ( @xmath180 ) for all @xmath111 , representing a fully polarized paramagnet with @xmath181 . for @xmath182
, it has @xmath177 for @xmath160 and @xmath179 for @xmath183 , or vice versa .
the two chain sublattices therefore have opposite orientations , @xmath184 , and the total magnetization vanishes .
comparison with fig .
[ figmdg ] confirms that this corresponds to the mdg phase .    at higher temperature
, the maximum has @xmath185 for all @xmath111 and @xmath186 , corresponding to a paramagnet with continuously varying magnetization .
this last phase is deconfined , having @xmath187 , and is therefore identified as the coulomb phase .
[ figphasediagram ] also shows the magnetization @xmath188 evaluated at the maximum of @xmath168 .
the resulting phase diagram consists of a set of low - temperature confined phases with fixed magnetization , surrounded by the coulomb phase in which the magnetization varies smoothly with the parameters .
the mdg phase consists of a lobe , whose width , i.e. , the range of fields for which it is stable , decreases with increasing temperature .    as a benchmark for the method
, we note that our model reduces for @xmath189 to that considered by jaubert et al .
@xcite ( in the limit of zero monopole fugacity ) . in this case
, one can restrict to confined phases with all chains equivalent ( i.e. , @xmath190 ) , and it is straightforward to show that the magnetization is given by @xmath191 where @xmath192 .
our method therefore agrees , in this case , with the results of the bethe - lattice calculation @xcite .
we have shown that applying mean - field theory to the model defined by @xmath109 in eq .
( [ eqhamiltonian ] ) produces a phase diagram containing the coulomb phase @xcite as well as a fully - polarized paramagnet @xcite and the mdg phase @xcite . in this section
, we argue that the general features in the phase diagram of fig .
[ figphasediagram ] , and in particular the lobe structure , can be understood by considering the interplay of confinement and flux in the coulomb phase .
we first consider the limit of vanishing monopole density , approximately valid at temperatures well below @xmath38 , before discussing the qualitative effect of including monopoles .
the crucial observation underlying the phase structure is that confinement implies suppression of fluctuations in the magnetization @xmath193 , and hence vanishing uniform magnetic susceptibility @xmath194 .
by contrast , phases where monopoles are deconfined , such as the coulomb phase , are magnetizable ( i.e. , have nonzero susceptibility @xmath194 ) .
this relationship between confinement and vanishing susceptibility ( or , more generally , flux variance ) occurs in any local model with a confining phase @xcite .
the lobe structure of the phase diagram follows immediately from this observation . starting within a confining phase and changing the applied field
will eventually tip the balance in favor of the coulomb phase , whose free energy is quadratic in the field .
the higher entropy in the coulomb phase implies that the range of applied fields for which a given ordered state is favored narrows as the temperature increases .
the confined phases therefore comprise lobes with fixed magnetization , each centered around favorable values of the applied field .
( the fully polarized phases are exceptions , extending to arbitrarily large @xmath195 as a result of their saturated magnetization . )
this qualitative structure is confirmed in the mean - field phase diagram , fig .
[ figphasediagram ] . with more realistic interactions ,
we expect more lobes with other commensurate values of magnetization to appear in the gaps between the lobes shown , as indeed seen in the mc results of ref .
@xcite .
the general phase structure can also be understood through the mapping to a quantum model of bosons , described in section  [ secbosonmapping ] and ref .
@xcite . for an extended hard - core bose
 hubbard model in the grand canonical ensemble , the phase diagram contains a set of mott - insulator lobes at small hopping , consisting of different possible ordered states @xcite . returning to spin ice , each of these corresponds to a confined phase with order in the spin degrees of freedom and fixed magnetization .
the susceptibility strictly vanishes in the ordered phases only in the thermodynamic limit at zero monopole fugacity @xmath35 . with finite monopole cost , confinement of monopoles
is no longer precisely defined and the susceptibility becomes nonzero , though remains exponentially small at low temperature .
the plateaux are therefore rounded for nonzero @xmath67 .
those phases that have conventional order , such as the mdg phase , remain separated from the paramagnet by a symmetry - breaking transition .
certain properties of the phase transitions can also be determined on general grounds , following similar arguments to those for the bose
hubbard model @xcite , but with density replaced by magnetization .
we first note that phase transitions can be classified by whether the magnetization changes across the transition .
the magnetization necessarily changes in the case where the ordered phase has saturated magnetization , but it can also change across a transition into a magnetically ordered phase .    at the tip of the lobe ,
the magnetization is identical in the two phases ( ordered within and coulomb outside ) , and so this is the only point at which the magnetization stays the same across the transition .
the argument is that , inside the ordered phase but sufficiently close to the tip , an arbitrarily small change in the applied field pushes the system into the coulomb phase , regardless of its sign .
if the two phases had different magnetizations at this point , their linear coefficients in the free energy would be different and this could not be the case .
( this argument is essentially identical to the corresponding one for the bose
hubbard phase diagram @xcite . )
for any continuous transitions with fixed magnetization , critical properties can be calculated using the approach of ref .
@xcite , but with a modified projective symmetry group due to the nonzero magnetization . while mdg found a strongly first - order transition @xcite and the transition into the partially polarized phase
was seen to be of first order in ref .
@xcite , it remains possible that transitions in the presence of different perturbations could be continuous .
( the possibility of unconventional transitions between a superfluid and fractional mott insulators has been discussed in the extended bose  hubbard model @xcite . )    when the magnetization changes across the transition , i.e. , everywhere except at the tip , the transition is in the same universality class as the kasteleyn transition @xcite .
it should be noted , however , that only the transitions into the saturated phases , with @xmath181 , are true kasteleyn transitions , which are distinguished by the complete absence of fluctuations on the ordered side .
we have presented a mean - field theory designed to study confinement transitions , based on the analogy between the confinement criterion and long - range order in conventional ordering transitions .
the method has been applied to a simplified model of classical spin ice and shown to produce a phase diagram containing , besides the coulomb spin liquid @xcite , a set of magnetization plateaux .
these include the saturated paramagnet @xcite as well as the mdg phase expected to occur in dipolar spin ice at low temperature @xcite .
we have argued that the general properties of the phase diagram follow from the interplay between confinement and magnetization that characterizes the ice - rule states and pointed out the analogy with the phase diagram of the extended bose  hubbard model @xcite .
the simplified model that was studied here reproduces the mdg phase and is particularly amenable to the mean - field approach .
it would be desirable to extend the approach to allow for more realistic further - neighbor interactions between spins , particularly the dipolar interactions that are known to be significant in the classical spin - ice materials .
further - neighbor interactions can also stabilize the partially magnetized phase observed using mc in ref .
( this phase can be described in terms of the chains shown in fig .
[ figmdgchains ] , but distinguishes four chain sublattices , rather than two . )    a transfer operator written in the form of eq .
( [ eqtetrahedrontransferoperator ] ) has the limitation , however , that it can represent interactions only within a single tetrahedron .
( one can not even define the layer transfer operator for general interactions . )
a potential route to include further - neighbor interactions is through an additional mean - field decoupling , replacing a pairwise interaction @xmath196 with terms such as @xmath197 , representing a mean field acting on site @xmath2 .
this extension , and others such as including nonzero monopole fugacity @xmath67 , are left to future work .
k. matsuhira , z. hiroi , t. tayama , s. takagi , and t. sakakibara , _ a new macroscopically degenerate ground state in the spin ice compound @xmath3  under a magnetic field _ , j. phys . : condens .
matter * 14 * , l559 ( 2002 ) .
p. henelius , t. lin , m. enjalran , z. hao , j. g. rau , j. altosaar , f. flicker , t. yavorskii , and m. j. p. gingras , _ refrustration and competing orders in the prototypical @xmath3  spin ice material _ , , 024402 ( 2016 ) .
this nomenclature follows lin and kao @xcite .
mcclarty et al .
@xcite call the same state a `` cubic antiferromagnet '' , while henelius et al .
@xcite call it the `` single - chain state '' . in ref .
@xcite , it was called a `` spin spiral '' .
mdg themselves @xcite refer to it by its ordering wavevector , @xmath199 , where @xmath200 is the side length of the cubic unit cell .
one can imagine creating a monopole
antimonopole pair , winding one of the pair around the periodic boundaries , and then recombining them @xcite .
the result of this charge transport is to increase the flux through the system ; since the separation required diverges in the thermodynamic limit , it is possible with finite energy cost only in a deconfined phase .