outflows are ubiquitous among quasars . they manifest themselves as blue - shifted broad and narrow absorption lines in uv and x - rays ( e.g. , boksenberg et al .
1978 ; turnshek et al .
1980 ; weyman et al .
1991 ) . the broad absorption lines ( bals ) with widths of around a few thousands to a few ten thousands
km  s@xmath0  appear in 15 - 25% of all quasars ( e.g.,weymann et al .
1991 ; knigge et al .
2008 ) , and intrinsic narrow absorption lines ( nals ) with widths around hundreds to one thousand km  s@xmath0  in 40 - 50% of quasars ( nester et al .
almost all known bal qsos show broad high ionization absorption lines , such as c  iv  , n  v  ( weymann et al . 1991 ; cf , mrk 231 , veilleux et al . 2013 ) , while only a small fraction ( about 25% ) of bal qsos display also low ionization broad absorption lines ( lobals , weymann et al .
1991 ; zhang et al . 2010 ; c.f .
urrutia et al . 2009 ) , commonly mg  ii  . a yet smaller subset of lobal qsos ( @xmath1% of all quasars ) exhibit broad fe  ii  absorption lines ( hereafter , felobals ; trump et al .
2006 ; c.f .
, dai et al .
2012 ) .
there are two different scenarios for bal and non - bal quasars or for a subclass of bal qsos . in the first scenario , the presence or the absence of bals , or more generally ,
the manifestation of different subclasses of bals , are attributed to the different viewing angles of the same intrinsic population ( weymann et al .
1991 ) . while in the second scenario ,
lobal and non lobal qsos are at different evolutionary stages ( boroson & meyers 1992 ; voit et al .
1993 ; hall et al .
in other words , only a small fraction of quasars possess low - ionization outflows .
early studies showed that the emission lines and broad - band spectral energy distributions ( sed ) are very similar for hibal and non bal quasars ( weymann et al .
1991 ) , suggesting the unification scheme .
hamann , korista & morris ( 1993 ) constrained the global covering factor to less 20% from the fraction of resonant scattering light . however , systematic trends in the sed and emission lines have been revealed with a large sample of bal , lobal in particular , and a non - bal quasar sample from sdss ( e.g. , zhang et al .
2010 ; baskin , laor & hamann 2013 ; farrah et al .
2010 ; wang et al . 2013 ; and reference therein ; c.f .
, lazarova et al . 2012 ) , suggesting that the global covering factor and properties of lobals depend on the intrinsic properties of the system .
feedback through these outflows is now widely believed to be a key physical process in the formation of massive galaxies ( fabian 2012 ) .
however , the mass loss rate and kinetic power , as well as the size , of outflows are still poorly constrained for the majority of bal quasars ( lucy et al .
2014 ) . to estimate these parameters from absorption lines
, one needs to know the global covering factor , column density and the gas density or distance .
the gas column density can be inferred from the column densities of a set of ions through photoionization modeling .
however , the gas density diagnostics for absorption lines are available for only a rare subset of bal qsos with measured absorption lines from low - lying excitation levels of heavy elements ( de kool et al . 2001 ; korista et al . 2008 ; dunn et al . 2012 ) .
previous measurements of density in a dozen of bals , mostly felobals , quasars , indicate that bal region ( balr ) is extended from 100 parsec to several kpc from quasars ( de kool et al .
2001 ; korista et al .
2008 ; borguet et al .
2013 ; chamberlain et al .
the large size implies a very large mass outflow rate ( @xmath2 m@xmath3  yr@xmath0 ) and large kinetic power ( a few to 10 percent of the total luminosity ) if the outflows cover a substantial fraction of the sky . however , absorption lines can only probe the gas on the line of sight , so the global covering factor of absorbing gas is virtually unknown in an individual object . because these quasars are so rare , the average global covering factor must be very small , according to the unification scheme , but if they are at a special evolution stage , on the other hand , the covering factor can be substantial .
interestingly , liu et al . ( 2015 ) detected extended outflow components in two low z agns on the physical sizes consistent with those inferred from ultraviolet absorption lines .    in this paper
, we report four quasars , including one that was noted already by ross et al . ( 2015 ) , with peculiar ultraviolet fe  ii  spectra .
they were discovered during the course of studying a sample of heavily obscured quasars in wise : @xmath4 .
we perform an analysis of the emission lines and the broad band sed , and present evidence that the fe  ii  emission lines in these objects are dominated by continuum fluorescence pumping in the extended outflows , so they can be used to constrain the covering factor and size of balr of lofebal .
the paper is organized as follows . in section 2
, we present an analysis of the optical spectra and broad band sed . in section 3
, we discuss the nature of fe  ii  emission and continuum emission , and a brief conclusion is given in section 4 .
spectral models of fluorescent fe  ii emission are presented in an appendix .
the optical spectra are shown in figure 1 .
one of these quasars , sdss j153542.4 + 090341 has already been noted in ross et al .
( 2015 ) ; and another one ( sdss j163246.6 + 340526 ) appeared on the list of unusual quasars in meusinger et al .
the most striking features of these objects are weak continua , strong ultraviolet fe  ii  spikes and mg  ii  emission lines , and a relatively narrow mg  ii  line .
individual objects are described and analyzed as follows .
* sdss j153542.4 + 090341 : the boss spectrum shows three prominent fe  ii  spikes due to uv1,2,3 and 62,63 , and mg  ii  , al  iii  and c  iv  lines .
he  ii  1640+o  iii1663 and c  iii ] are weak with respective to these lines .
we measure an extreme fe  ii  spike / gap ratio of 8.5 .
a narrow fe  ii  78 is also visible .
in fact , the equivalent widths of c  iii ] and he  ii  are consistent with those of i zw 1 ( see figure 1 , also laor et al .
1997 ) , while those of fe  ii  , mg  ii  , al  iii  and c  iv   are a factor of 10 - 20 times larger .
this indicates that the uv spectrum consists of two components , a weak i zw 1 component and another component dominated by strong emission lines , fe  ii  , c  iv  , al  iii  and mg  ii  .
c  iv  shows two components , a narrow component at the systematic redshift , derived from [ o ii ] , and a blue - shifted component .
alternatively , c  iv  may be split into two components due to a blue shifted narrow absorption line .
the mg  ii  doublet can be well modeled with two gaussians .
the best fit converges to a width of @xmath5 km  s@xmath0  and no shift in velocity ( 2@xmath6500 km  s@xmath0  ) .
[ oii ] appears fairly narrow overall , but with weak broader wings .
[ ne iii]@xmath73869 looks similar with a narrow core and extended wings
. however , both lines are severely contaminated by the subtraction of sky lines . * sdss j100552.6 + 493448 : in addition to strong ultraviolet fe  ii  spikes , this object displays very strong high - ionization forbidden and permitted lines : [ o  iii]@xmath74363 , [ neiii]@xmath83869,3987 , [ nev]@xmath83346,3426 , [ neiv]@xmath72422 , [ ciii]@xmath71906,1909 , as well as low ionization lines @xmath9 , @xmath10 , [ oii]@xmath83726,3729 , cii]@xmath72326 , mg ii , and al iii@xmath81855 , 1863 .
in fact , the spectrum can be reproduced by a supposition of broadened sdss j1535 + 0903 emission lines and the narrow lines mentioned above ( see figure 1 ) . the al  iii  strength
can be well accounted for by the sdss j1535 + 0903 component .
to first order , all but [ o ii ] , have the same width , around fwhm@xmath113150 km  s@xmath0  , using a single gaussian fit . [ o ii ]
is considerably narrower than other lines , probably from an extended star - forming region .
there is also a weak broad wing in [ o ii ] , which may be similar to other lines .
the weak line to the red side of mg  ii  can be either fe  ii  or mg i. [ o  iii]@xmath74363 is stronger than h@xmath12 according to the single gaussian fit with their centers tied to other lines .
unfortunately , the boss spectrum does not cover [ o  iii]@xmath84959,5007 , which should impose important constraints on the gas density and gas temperature in combination with [ o  iii]@xmath74363 .
another constraint on the density comes from the fit to 1900 blending .
we find that in the single gaussian model for each line the fit is significantly worse when the lines are assumed to be [ si iii ] and [ c  iii ] with a @xmath13 with respect to models assuming they are si iii ] and c iii ] .
a further fit with a mixture of c  iii ] and [ c  iii ] yields a [ c  iii]/c  iii]=0.11@xmath60.02 with a @xmath14 lower than pure c  iii ] model .
this gives a gas density @xmath15   @xmath16 ( osterbrock & ferland 2006 ) . the outflow has a considerably lower density than that in q 1321 + 058 ( wang et al .
detailed photoionization modeling of the outflow component will be presented elsewhere .
* sdss j163246.6 + 340526 : the spectrum is dominated by narrow ultraviolet fe  ii  lines . the fe  ii  spikes / gap ratio ( 3.2 ) is less extreme than in j1535 + 0903 .
we adopt the redshift of narrow [ o ii ] as the systemic redshift . fe  ii  and mg  ii  appears fairly narrow .
we fit the mg  ii  multiplet using two gaussian and obtain a best fit with an @xmath17 km  s@xmath0  and a doublet ratio of @xmath18 with no systematic redshift .
mg  ii  is slightly narrower than mg  ii  in i zw 1 ( laor et al . 1997 ) .
interestingly , the observed spectrum in the 2100 - 3100   range can be reproduced approximately by the addition of a scaled i zw 1 spectrum with four - times stronger emission lines plus the fe  ii  spectrum of sdss j1535 + 0903 ( see figure 1c ) . [ o iii ] @xmath75007 is visible in the gap of the sky line forest , and has the same redshift as [ o ii ] .
unfortunately , h@xmath19 is seriously affected by residuals of sky lines . * sdss j212546.9 + 004455 : unlike other objects , the continuum emission in this object is relatively strong .
the narrow mg  ii  line has the same redshift as [ o ii ] , which is considered as the systematic redshift .
a comparison with the fe  ii  template of i zw 1 and the spectrum of sdss j1535 + 0903 suggests the presence of two blue - shifted ultraviolet fe  ii  components . a broadened sdss j1535 + 0903-like component
is shifted by 5700 km  s@xmath0   and a i zw 1-like ultraviolet fe  ii  component is blueshifted by @xmath202850 km  s@xmath0  .
the i zw 1-like component has weaker mg  ii  , c  iii , si iii and c  iv  lines but broader line widths
than i zw 1 .
we measure the fwhm of mg ii using the above two - component fit ( see also figure 1 ) . there appear narrow absorption lines of c  iv  at velocities of 900 and 2300 km  s@xmath0  .
but the signal to noise ratio is low .
table 1 presents the emission - line parameters , including the width of mg  ii   and the equivalent widths of fe  ii  spikes , mg  ii  and al  iii  . to summarize , sdss j1535 + 0903 represents a prototype peculiar fe  ii   spectrum , while the spectra in other objects can be modeled as a combination of this and a i zw 1-like fe  ii  component .
broad band seds of these quasars are shown in figure 2 .
they are characterized by strong and steep mid - infrared emission , and a relatively weak ultraviolet and optical continuum .
the photometry is strongly affected by the ultraviolet fe  ii   and mg  ii  .
we corrected the emission - line flux using the sdss spectrum .
we fit these seds using a model consisting of a clumpy torus plus a galaxy template , accounting for the host galaxy light .
we use all galaxy templates , excluding agns , that are in the swire library ( polletta et al .
the purpose of this fit is to give an estimate how much of the infrared emission comes from the agn heated dust and how this dust is distributed .
the seds of clumpy tori were calculated for a wide range of parameters by hnig & kishimoto ( 2010 ) . in their models ,
the torus consists of optically thick clouds with a small total filling factor .
the number density of clouds decreases outward in a power law form , and has a gaussian distribution in the vertical direction .
the parameters of the models are the number of clouds along the line of sight on the equator ( @xmath21 ) , the radial dust - cloud distribution index ( @xmath22 ) , the subtending angle of the torus ( @xmath23 ) , the optical depth of individual clouds , and the inclination ( refer to hnig & kishimoto 2010 for details ) . the inner radius is fixed to the dust sublimation radius and the outer radius of the torus to 300 dust sublimation radii .
the broad - band sed can be fitted fairly well with such a model .
the best fits are shown in figure [ sedfit ] and the parameters are listed in table [ emissionline ] .
we also compute the distribution of the agn fraction in the rest frame 3 - 15 @xmath24 m using all models weighted by the @xmath25 probability of each model and by the solid angle of the viewing angle , i.e. , @xmath26 .
the results suggest that the infrared emission is dominated by a torus component , and that the galaxy component contributes less than 5% in all cases .
the power - law index for the dust distribution is quite broad from @xmath27 .
thus the steepness in the mid - infrared spectrum is due to weak host galaxy light in the near - infrared , indicating a high eddington ratio of the system , and an edge - on viewing angle .
this result is probably a color selection effect , where we only choose objects with a very red infrared color .
second , the half opening angle of the torus is small ( 45 to 30 degree ) , and the system is viewed at a relatively large inclination angle ( 65 - 75 degree ) , consistent with the large infrared to optical ratio .
finally , torus models predict a uv flux higher than the observed ones if there is no internal extinction .
the reddening to the reflected light can be due to diffuse interstellar dust further out in the host galaxy , which may be partially covering .
note that our model does not include other optical or uv continuum sources , such as direct or electron scattered accretion disk light .
if such components are present , then the contribution from the host galaxy should be even lower and the optical and uv continuum from the dusty torus should be highly extinguished .
the large infrared to optical / uv ratio and the steep near to mid - infrared ratio indicates that the optical to uv , and even the near - infrared , emission from the accretion disk is obscured . assuming the mid - infrared emission , which is dominated by the dusty torus emission as suggested by sed fitting , is nearly isotropic , the continuum emission at rest frame 3000   is a factor of 30 - 300 weaker than in normal quasars ( refer to figure 1 ) .
if the 3000  flux is attributed to the transmitted light , the dust extinction would be @xmath28 for an smc - like extinction curve . applying such dust extinction to the average uv spectrum of quasars
, we will get a slope @xmath29 ( @xmath30 ) between 2000 and 3000 , which is much steeper than the observed one .
adopting an average galactic extinction curve will yield only slightly smaller slope .
this indicates that the observed uv continuum is reflected / scattered light . in fact , the residual ultraviolet flux can be well accounted for by reflected light from the clumpy dust torus itself as in the sed fitting .
alternatively , the continuum can be due to electron scattering , which is known to be important in seyfert 2 galaxies ( e.g. , antonucci & miller 1985 ) .
this suggests that the optical depth to the continuum emission region is much larger than indicated by the @xmath31 derived above . in either case
, the continuum should be polarized , and future polarization observation could verify this .    in the sed fit , we found that the dusty torus has a large subtending angle , i.e. , covering 71 to 87% of the sky .
the large inclination angle suggests that our line of sight to the continuum source is blocked by the torus , consistent with weak optical and uv emission . in this case , the hot dust in the illuminated face of inner dust torus is also blocked , leading to a steep rise from the near to mid - infrared .
so the dusty torus can naturally explain the obscuration and infrared sed .
the sed fitting also suggests that the contribution of stellar light to the near and mid - infrared is small , which implies a high eddington ratio .
once the direct continuum is obscured , the large line equivalent width can be understood if the emission - line region is more extended than the dusty torus . in a few nearby seyfert galaxies , the warm few hundred k dust which is responsible for mid - infrared emission
is resolved into a thick structure of a few parsecs with a bolometric luminosity around @xmath32 erg  s@xmath0 ( tristam et al .
2007 ; lpez - gonzaga et al . 2014 ; tristram et al . 2014 ) .
note that this is not the outer radius of the torus , which might extend to several tens of parsecs .
since the size of the torus is proportional to @xmath33 , with an infrared luminosity of @xmath34 erg  s@xmath0 for our quasars , the true size will be at least several tens of parsecs . in either case
, the emission line region should be larger than this .    before making a detailed discussion the size of the obscurer ,
let us consider the constraint on the size of emission line region first .
the emission line width ( fwhm ) , as fitted with mg  ii  , is in the range 900 - 8000 km  s@xmath0  , which is broader than narrow lines in seyfert galaxies , but is comparable with or narrower than the broad lines . naively , this can be interpreted as lines that come from an intermediate region between the traditional blr and nlr ( hu et al .
2008 ; li et al .
however , note that sdss j1005 + 4934 also displays strong forbidden lines of similar width ( @xmath203000 km  s@xmath0  ) , different from traditional intermediate broad lines .
thus , we consider it more likely that these emission lines come from the inner region of nlr , and that the large width is due to outflows ( e.g. , komossa et al . 2008 ; zhang et al . 2011 ) . in this picture
the traditional blr is completely obscured .    while the dusty torus provides the most naturally explanation for the obscuration , it is worth considering other alternatives .
the obscurer can be a single dusty cloud in the clumpy torus model or the distant interstellar medium .
the obscuration of the blr requires a minimum size cloud to be @xmath35 parsecs ( bentz et al .
2013 ) . using the infrared luminosity as an approximation for the bolometric luminosity and a bolometric correction to @xmath36
, we estimate blr sizes in the range of 0.1 - 0.3 parsecs . assuming a galactic dust to gas ratio ( @xmath37 @xmath38 , bohlin et al .
1978 ; cox et al . 2006 ) , the gas column density of the obscurer should be much greater than @xmath39 @xmath40 as discussed above .
note that in the clumpy torus , the gas density should be of the order of a few @xmath41 @xmath16 ( krolik & begelman 1988 ) , so a blr - sized cloud would have a column density two orders of magnitude larger than the above minimum value , or @xmath42 mag for a galactic dust to gas ratio . on the other hand , assuming a gas density more typical of the interstellar medium @xmath43 @xmath16 , the minimum size of cloud would be several tens of parsecs to meet the minimum column density .    in the case of the single - cloud obscurer
, one would expect to see long - term very large amplitude variability as the cloud moves in and out of the line of sight .
the typical time scale would be the size of the cloud divided by the transverse speed , which is of order the local keplerian velocity . in a typical galaxy ,
the orbit velocity is order of 500 - 1000 km  s@xmath0  , depending on the distance to the center , and a minimum size of a cloud should be blr sized , of order of 0.1 to 1 parsecs .
this would give a variability time scale @xmath44
so it is difficult to observe the cloud transit event fully .
however , the continuum source is much smaller ( a hundred to thousand times smaller ) , so it is possible to observe transient obscuration of the continuum source with small changes in broad lines on time scales of years if the cloud has a sharp edge .
there can also be cases where the cloud blocks the continuum source fully but only partially blocks the broad line region as the obscuring cloud moves in and out .
these sources will show very large equivalent widths of broad lines .
the probability of detecting such a source depends on the relative size of the obscuring cloud relative to the blr size . in the case of comparable sizes of the two
, one would expect that in most sources only the continuum is fully covered , while if the cloud is much larger than the blr , most likely both continuum and blr are obscured .
similarly , we will observe many more such partially obscured sources in the clumpy torus than in the smooth dust torus .
in reality , the cloud may have a smooth edge so the contrast between emission line and continuum should be much moderate .
we note that a single cloud occultation event is likely to give a smaller polarization degree than the dusty torus obscuration due to different symmetries ( marin et al .
2015 ) .    the weakness of the optical to uv continuum and the strong mid - infrared emission may be due to a temporarily shut - off of the central engine .
a comparison of infrared , optical , and uv fluxes requires a drop in the uv flux from the inner disk by a factor of 30 - 300 in these sources , compared to typical agn ( see figure 2 ) .
the size of the mid - infrared emission region is of order of tens to hundreds of dust sublimation radii at the infrared luminosity according to our model of the sed ( hnig et al .
2010 ) . at the infrared luminosity of these sources ,
the light - crossing time of the mid - infrared emission region is of order ten to a hundred years , so the central engine would have to have shutdown in less than this time .
a comparison of sdss photometric data with palomar observatory sky survey ( poss ) on a time scale of 50 years in the observer frame shows a small fraction of quasars with a variability amplitude of 3 mags in @xmath45-band ( de vries et al .
although these amplitudes are still smaller than that is required , since we do not have light curves of agn on such long time scales , it is impossible to tell whether such a large amplitude variability happened on time scale of a hundred years and how frequent they are .
thus this remains a possibility .
there are several tests that could distinguish among these possibilities .
first , scattered light should be highly polarized , while direct light from a much less active accretion disk is likely not .
the observed scattered light in seyfert 2 galaxies is polarized by a few to tens percent , while the polarization of direct emission from quasars is usually less than 1% ( marin et al . 2015 ) .
it should be noted that if a quasar resides in the gas - rich environment , there is inevitable scattered light from surrounding gas . as the central continuum source dims quickly , the fraction of scattered light increases , and so does the polarization degree .
so a more critical test would be to observe the transmitted hard x - rays .
second , direct emission from accretion disk is likely to be variable on relatively short time scales , while the scattered light is not .
thus , the detection of variability on time scales of less than a year would support the intrinsic weak continuum emission model .
we now turn to the properties of the emission region .
previous studies suggested that quasar emission lines consist of distinct kinematic components : broad , intermediate and narrow components ( brotherton et al .
1994 ; marziani et al .
( 2008 ) found that optical fe  ii  lines are formed in the intermediate region , that is redshifted with respective to the system .
the fe  ii  spectrum of our sources does not look like those in normal quasars or from most theoretical simulations aimed at reproducing observed fe  ii  in i zw 1 ( baldwin et al .
2004 ; sigut & pradhan 2003 ; ferland et al .
2009 ; bruhweiler & verner 2008 ) .
we will mainly focus on the fe ii emission in sdss j1535 + 0903 first , because in other sources the spectra can be reproduced with a combination of sdss j1535 + 0903 and i zw 1-like components .
the uv fe  ii  emission is dominated by a few bumps , namely multiplets uv 1 , uv 2,3 , uv 62,63 , with weak emission in between and very weak or no emission at wavelengths longer than 2900 . in sdss j1005 + 4934
, the sdss spectrum extends to 4400  in the quasar rest frame , and optical fe  ii  lines appear relatively weak
.    for thermal excitation , the extreme spike / gap ratio of 8.5 in sdss j1535 + 0903 may indicate very high densities ( @xmath46 @xmath16 ) and very low ionization parameters ( @xmath47 ) , as shown in figure 3 of baldwin et al .
but their models are only for the broad line region , and the ionizing photon flux is too large for our case
. however , in their figure , the curves of iso - spike / gap ratio run diagonally , so the dimensionless ionization parameter may be more crucial than the density .
unfortunately their models do not extend to the low density and very low ionization regime , so it is not clear whether extreme spike / gap ratios can be reproduced there .
but if this is indeed the case , the fe  ii   emission region must be far from the quasar so that it is not obscured by the dusty torus . at very low values of the ionization parameter , some other lines , such as c  ii and ni  ii , might become strong .
detailed photo - ionization modeling would be required .
also the presence of strong c  iv  and al  iii  , and the lack of c  iii  and si iii can not be explained by this model .
alternatively , the extended emission may come from resonant scattering in low - ionization quasar outflows .
previous observations suggest that felobals and hibals , at least in some quasars , are formed on the scale of 100 pc to a few kpc ( see references in  1 ) , so it is likely that outflows are not obscured by the putative dusty torus that blocks the line of sight to the central continuum in type 2 quasars .
the outflows may produce fe  ii  via collisional excitation , as well as via continuum fluorescence . the contribution of fluorescence relative to collisional excitation depends on the optical depth of fe  ii  as well as the density and temperature of the gas . in the case of felobal outflows , the density is usually low ( order of @xmath48 @xmath16 , de kool et al .
2001 ; lucy et al . 2014 ; c.f .
, zhang et al . 2015 ) , and the optical depth of most fe  ii  lines is usually not large ( order of a few to tens for strong lines ) .
thus , continuum fluorescence may be the dominant process in the outflows for fe  ii  lines .
some other emission lines may be also dominated by resonant scattering from such outflows .
consequently , these lines should appear strong in comparison with non - resonant lines .
the presence of strong al  iii  , but very weak or missing similar component of c  iii ] and si iii ] , which are usually much stronger , in sdss j1535 + 0903 , indicates that resonant scattering is the dominant excitation mechanism in this object .
al  iii  can be produced by either resonant scattering or collisional excitation , while c  iii ] and si iii are only collisionally excited .
although the relative strength of al  iii  increases with gas metallicity , this does not explain the fact that al  iii  is much stronger than c  iii ] . in table 2 ,
the al  iii  equivalent width is a factor of 3.7 larger than the c  iii ] equivalent width . as discussed in
[ 3.1 ] , the observed uv continuum is likely scattered light , and , as seen in figure [ specfit ] , the equivalent widths of c  iii ] and si iii is similar to that of i zw 1 .
so c  iii ] emission may be fully accounted for by scattered component and the contribution from the outflow component may be negligible .
if fe  ii  is produced by continuum fluorescence in outflows , a simple idea is that its spectrum may resemble those of the fe  ii  absorption line spectra in felobal qso sdss j080957.39 + 181804.42 .
we note that there are usually very weak fe  ii  absorption lines above 2800 , and no optical fe  ii   absorption lines , consistent with the observed fe  ii  in these spectra .
to test this idea , we extract the absorption line spectrum of a felobal qso .
this spectrum is shown in figure 1 , and is similar , to first order .
it should be noted that the spectrum of fluorescent emission is not exactly the same as the fe  ii  absorption spectrum , because an fe  ii  ion in an excited level , after absorbing a photon , may decay to a different lower level , thus emitting a photon with an energy different from the absorbed one , especially when resonance lines become optically thick .
detailed calculations are presented in  [ model ] .    in appendix
a , we demonstrate that a pure fluorescent model can indeed reproduce the observed fe  ii  band ratios in sdss j1535 + 0903 very well . in order to make fe  ii  the dominant ion of fe , we let the incident ionizing continuum be filtered by an intervening absorbing gas of column density 10@xmath49   @xmath40 .
we adopt a gas density @xmath50 @xmath16 , consistent with those in the felobal qsos , and a turbulent velocity of 1000 km  s@xmath0 to mimic the outflow velocity .
we find that at the optical depth of the resonant fe  ii  line at @xmath51 is 900 , and that the model produces uv feii band ratios that are consistent with the observed ones .
as far as we know , this is the only model that gives such an extreme fe  ii  ratio at a low ionizing continuum flux , which is consistent with its large line formation region .
furthermore , the model gives weak fe  ii  lines at wavelength longer than 2800 , which is also consistent with the observation .
the best matched model predicts an equivalent width of fe ii at 2355  band of 99 , which is about one 6th of the observed value for sdss j1535 + 0903 , for a covering factor of one .
however , as we discussed in
[ 3.1 ] , the continuum incident to the gas is obscured and the observed continuum is attenuated by two orders of magnitude , so only a few percentages covering factor would be sufficiently to explain the observed fe ii strength .
this covering factor is consistent with the estimate in next subsection .
the gas producing the fe ii emission is cool , as shown in the appendix , so collisional excitation is unimportant .
this accounts for the weak optical fe ii emission .
note that this type of model will also predict a number of other low ionization lines , such as ti ii , ni ii and si ii , but their atomic data are very sparse .
so a full treatment of these lines has to wait for a major update in the atomic data .
our models in appendix a do not produce al iii and c iv , because ionizing photons above 13.6 ev are required to produce these ions while in our scenario these photons are intentionally filtered out to prevent destruction of fe@xmath52 ions .
both al  iii  and c  iv  are observed in absorption in a felobal qso , and their average emission spectrum should not be different from the absorption line spectrum in a felobal qso because they are scattered resonantly .
so our purpose is only to demonstrate the effect of fe  ii  radiation transfer after fe  ii  resonant absorption .
note also that the observed fe  ii  ratios from different fe  ii  bands depends on the reddening , because felobal qsos often show a large internal dust extinction .
we do not apply such a reddening correction because we do not know the location of the dust .
if dust lies outside of the balr , then an extinction correction should be applied . on the other hand ,
if dust is within the balr , no correction is needed .
anyway , to first order the absorbed line spectrum of felobal quasars are similar to the quasars observed in the sample .    in sdss j1005 + 4934
, we also detect strong high ionization forbidden or semi - forbidden lines that are excited by collisional process , while the fe  ii  spectrum as well as fe  ii  /al  iii  ratio , are nearly the same as sdss j1535 + 0903 .
we also obtained an electron density of ( 3 - 4)@xmath53 @xmath16 from c  iii  and si iii .
these lines are rather broad with widths of 3000 km  s@xmath0  .
the presence of such lines does not necessarily contradict our interpretation of the strong fe  ii  , al  iii   and mg  ii  lines as resonantly scattering from an outflow . according the the recent model of faucher - gigure et al .
( 2012 ) , the felobal may be produced in the shocked cooling gas during interaction of outflows with the cold gas .
the gas producing strong high ionization forbidden lines may be the warm phase of such shocked gas , and the large width of the line may be explain in this way as well .
blue shifted fe  ii  and mg  ii  lines are seen in the sdss j1632 + 3405 , while emission lines in other quasars do not display systematic blueshift .
the blue shifted lines can be easily explained as from an outflow since one expects that the far side of the outflow is subject to large dust extinction and would not be seen .
however , emission lines from outflows are not necessarily systematically blueshifted if the farside of the outflows is not particularly obscured or the system is seen nearly edge - on . in the continuum fit , we found that sdss j1632 + 3405 has the smallest inclination angle among these objects .
if the fe  ii  in these objects is due to fluorescent excitation , we can constrain the global covering factor of the outflows in these objects .
to first order , we can assume an isotropic angular redistribution function , so the scattered fe  ii  flux will be @xmath54f_c(\lambda)*c\ ] ] where @xmath55 is the un - absorbed continuum flux . assuming a maximum @xmath56 at certain @xmath7 , one can then estimate a minimum covering factor @xmath57_{max}$ ] . to estimate @xmath55 , we assume that the intrinsic sed is similar to the mean quasar sed , and scale the latter to match the wise luminosity at 22@xmath24 m , the peak wavelength of sed . by doing so
, we obtain covering factors of 0.05 , 0.20 , 0.01 , and 0.06 for sdss j2125 + 0044 , sdss j1632 + 3405 , sdss j1005 + 4934 and j1535 + 0903 . note that if mg ii is dominated by resonant scattering in the balr , then the covering factor of mg ii balr would be a factor of 2 - 3 times larger than the above value .
there are several uncertainties in this estimate .
first , due to large number of energy levels , fe ii can not be rigidly considered as resonant scattering , but rather as the fluorescent process .
the absorbed photons may emerge at other wavelength , so the strong resonant lines are usually weaker than in the case of resonant scattering .
second , it is possible that infrared emission of these quasars is enhanced because we selected objects with very red infrared colors . in this case
, the unabsorbed continuum is over - estimated . since sed modeling suggested that infrared emission is dominated by the agn component , while it was shown that , on average , thermal dust emission accounts for 35 percent of the total power of quasars ( e.g. , ma & wang 2013 ) , thus at most , the intrinsic continuum is over - estimated by a factor of 2 - 3 .
third , the scattered photons may be destroyed by dust within or outside the balr , so we may underestimate the covering factor .
it has been known for quite a long time that felobal qsos are reddened . but
the location of the dust is still unknown .
recently , there are arguments that dust may be mixed within the balr because near - infrared emission is correlated with the outflow velocity ( zhang et al .
2014 ) , and the variability of bal is correlated with the uv spectral slope ( he et al .
further , dunn et al . (
2015 ) showed that reddening must be caused by large scale dust because both the continuum and narrow lines are reddened .
if dust is outside of the balr and the scattered light suffers the same extinction as the continuum source , and adopting a median value for the continuum extinction for lobal @xmath58 ( dunn et al .
2015 ) , the correction will raise the covering factor by about 40% . despite these uncertainties
, we note that in 3 of 4 objects , the uncorrected covering factors are already larger than 0.05 , while according to statistics , felobal qsos account for no more than 2% of total quasars ( trump et al .
2006 ; dai et al .
this suggests strongly that felobal qsos must be in a special phase of quasar evolution , rather than representing a very small covering factor .
we presented an analysis of four quasars with peculiar uv fe  ii  spectra .
they had extremely strong uv1,3,4,62,63 and very weak fe  ii  at wavelengths longer than 2800 .
their seds rise very steeply from the optical to the mid - infrared , suggesting they are obscured .
the obscurer is most likely the dusty torus , since that also provides a self - consistent explanation for the infrared emission although other alternatives can not be ruled out .
the characteristics of the fe  ii  spectra can be understood as due to resonant scattering in an extended outflow since the size of the scattering region must be outside the obscuring torus and its fe  ii  spectrum similar to absorption spectrum of felobals .
our fluorescent models can successfully produce the observed uv fe  ii  band ratios and the weakness of fe  ii  lines at wavelengths longer than 2800with a low ionizing photon fluxes .
this is because the strong uv lines are excited by continuum pumping within gas that is too cool for collisional excitation to be significant .
strong al  iii  and weak c  iii ] and si iii emission may be explained as resonant scattering , as is observed in felobal qsos .
the obscuration of the direct continuum also explains the very large equivalent widths of fe  ii  , mg  ii  and al  iii  lines .
these objects are misaligned felobals . by considering the amount of scattered light we set a lower limit on the covering factor of felobals to be 5% to 20% in three objects , supporting
the scenarios that felobal qsos are a special stage of quasar evolution .
we thank the referee for constructive comments that help to improve the presentation of the paper .
we acknowledge the financial support by the strategic priority research program `` the emergence of cosmological structures '' of the chinese academy of sciences ( xdb09000000 ) , nsfc ( nsfc-11233002 , nsfc-11421303 , u1431229 ) and national basic research program of china ( grant no .
2015cb857005 ) . funding for sdss - iii has been provided by the alfred p. sloan foundation , the participating institutions , the national science foundation , and the u.s .
department of energy office of science .
the sdss - iii web site is http://www.sdss3.org/. sdss - iii is managed by the astrophysical research consortium for the participating institutions of the sdss - iii collaboration .
gjf acknowledges support by nsf ( 1108928 , 1109061 , and 1412155 ) , nasa ( 10-atp10 - 0053 , 10-adap10 - 0073 , nnx12ah73 g , and atp13 - 0153 ) , and stsci ( hst - ar- 13245 , go-12560 , hst - go-12309 , go-13310.002-a , hst - ar-13914 , and hst - ar-14286.001 ) .
40 antonucci , r.  r.  j. , & miller , j.  s.  1985 , , 297 , 621 baldwin , j.  a. , ferland , g.  j. , korista , k.  t. , hamann , f. , & lacluyz ,  a. 2004 , , 615 , 610 baskin , a. , laor , a. , & hamann , f.  2013 , , 432 , 1525 bentz , m.  c. , denney , k.  d. , grier , c.  j. , et al .
2013 , , 767 , 149 bohlin , r.  c. , savage , b.  d. , & drake , j.  f.  1978 , , 224 , 132 boksenberg , a. , carswell , r.  f. , smith , m.  g. , & whelan , j.  a.  j.  1978 , , 184 , 773 borguet , b.  c.  j. , arav , n. , edmonds , d. , chamberlain , c. , & benn , c.  2013 , , 762 , 49 boroson , t.  a. , & meyers , k.  a.  1992 , , 397 , 442 brotherton , m.  s. , wills , b.  j. , francis , p.  j. , & steidel , c.  c.  1994 , , 430 , 495 bruhweiler , f. , & verner , e.  2008 , , 675 , 83 cicone , c. , maiolino , r. , sturm , e. , et al .  2014 , , 562 , a21 chamberlain ,
c. , arav , n. , & benn , c.  2015 , , 450 , 1085 cox , n.  l.  j. , cordiner , m.  a. , cami , j. , et al .
2006 , , 447 , 991 dai , x. , shankar , f. , & sivakoff , g.  r.  2012 , , 757 , 180 de kool , m. , arav , n. , becker , r.  h. , et al .
2001 , , 548 , 609 dunn , j.  p. , arav , n. , aoki , k. , et al .  2012 , , 750 , 143 dunn , j.  p. , wasik , b. , holtzclaw , c.  l. , et al .
2015 , , 808 , 94 farrah , d. , urrutia , t. , lacy , m. , et al .
2010 , , 717 , 868 fabian , a.  c.  2012 , , 50 , 455 faucher - gigure , c .- a .
, quataert , e. , & murray , n.  2012 , , 420 , 1347 feruglio , c. , maiolino , r. , piconcelli , e. , et al .
2010 , , 518 , l155 ferland , g.  j. , peterson , b.  m. , horne , k. , welsh , w.  f. , & nahar , s.  n.  1992 , , 387 , 95 ferland , g.  j. , hu , c. , wang , j .-
, et al .  2009 , , 707 , l82 ferland , g.  j. , porter , r.  l. , van hoof , p.  a.  m. , et al .
2013 , , 49 , 137 hall , p.  b. , anderson , s.  f. , strauss , m.  a. , et al .
2002 , , 141 , 267 hamann , f. , korista , k.  t. , & morris , s.  l.  1993 , , 415 , 541 he , z .- c . ,
bian , w .- h . ,
jiang , x .-
, & wang , y .- f .
2015 , , 454,3962 hnig , s.  f. , & kishimoto , m.  2010 , , 523 , a27 hu , c. , wang , j .- m . , ho , l.  c. , et al .
2008 , , 683 , l115 jenkins , e.  b.  2009 , , 700 , 1299 knigge , c. , scaringi , s. , goad , m.  r. , & cottis , c.  e.  2008 , , 386 , 1426 krolik , j.  h. , & begelman , m.  c.  1988 , , 329 , 702 komossa , s. , xu , d. , zhou , h. , storchi - bergmann , t. , & binette , l.  2008 , , 680 , 926 korista , k. , bautista , m.a.,arav , n. et al .
2008 , , 688 , 108 laor , a. , jannuzi , b.  t. , green , r.  f. , & boroson , t.  a.  1997 , , 489 , 656 lazarova , m.  s. , canalizo , g. , lacy , m. , & sajina , a.  2012 , , 755 , 29 li , z. , zhou , h. , hao , l. , et al .
2015 , , 812 , 99 liu , g. , arav , n. , & rupke , d.  s.  n.  2015 , , 221 , 9 lpez - gonzaga , n. , jaffe , w. , burtscher , l. , tristram , k.  r.  w. , & meisenheimer , k.  2014 , , 565 , a71 lucy , a.  b. , leighly , k.  m. , terndrup , d.  m. , dietrich , m. , & gallagher , s.  c.  2014 , , 783 , 58 ma , x. , & wang , t.g .
2013 , , 430 , 3445 marin , f. , goosmann , r.  w. , & gaskell , c.  m.  2015 , , 577 , a66 marziani , p. , sulentic , j.  w. , dultzin - hacyan , d. , calvani , m. , & moles , m.  1996 , , 104 , 37 meusinger , h. , schalldach , p. , scholz , r .- d . , et al .
2012 , , 541 , a77 nestor , d. , hamann , f. , & rodriguez hidalgo , p.  2008
, , 386 , 2055 osterbrock , d.  e. , & ferland , g.  j.  2006 , astrophysics of gaseous nebulae and active galactic nuclei , 2nd .
ed .  by d.e .
osterbrock and g.j .
ferland .
sausalito , ca : university science books , 2006 , 126 polletta , m. , tajer , m. , maraschi , l. , et al .
2007 , , 663 , 81 richards , g. t. , lacy , m. , storrie - lombardi , l. j. , et al . 2006 , , 166 , 470 rllig , m. , abel , n.  p. , bell , t. , et al .
2007 , , 467 , 187 ross , n.  p. , hamann , f. , zakamska , n.  l. , et al .
2015 , , 453 , 3932 sigut , t.  a.  a. , & pradhan , a.  k.  2003 , , 145 , 15 tristram , k.  r.  w. , burtscher , l. , jaffe , w. , et al .
2014 , , 563 , a82 tristram , k.  r.  w. , meisenheimer , k. , jaffe , w. , et al .
2007 , , 474 , 837 o trump , j.  r. , hall , p.  b. , reichard , t.  a. , et al .
2006 , , 165 , 1 turnshek , d.  a. , weymann , r.  j. , liebert , j.  w. , williams , r.  e. , & strittmatter , p.  a.  1980 , , 238 , 488 urrutia , t. , becker , r.  h. , white , r.  l. , et al .  2009 , , 698 , 1095 verner , e.  m. , verner , d.  a. , korista , k.  t. , et al .  1999 , , 120 , 101 voit , g.  m. , weymann , r.  j. , & korista , k.  t.  1993 , , 413 , 95 veilleux , s. , trippe , m. , hamann , f. , et al .
2013 , , 764 , 15 de vries , w. h. , becker , r. h. , white , r. l. , & loomis , c. 2005 , , 129 , 615 wang , t. , zhou , h. , yuan , w. , et al .
2009 , , 702 , 851 wang , h. , xing , f. , zhang , k. , et al .
2013 , , 776 , l15 wang , t.g .
, yang , c.w . ,
wang , h.y . , & ferland , g. j.  2015 , apj , in press ( arxiv:1510.04945 ) weymann , r.  j. , morris , s.  l. , foltz , c.  b. , & hewett , p.  c.  1991 , , 373 , 23 zhang , s. , wang , t .-
g . , wang , h. , et al .
2010 , , 714 , 367 zhang , k. , dong , x .- b . ,
wang , t .-
, & gaskell , c.  m.  2011 , , 737 , 71 zhang , s. , wang , h. , wang , t. , et al .
2014 , , 786 , 42 zhang , s. , zhou , h. , wang , t. , et al .
2015 , , 803 , 58    ccccccccccccccc    100552.6 + 493448 & 1.1217 & 3150 & 173 & 6.4 & 310 & 35 & 103 & + 153542.4 + 090341 & 1.5330 & 1230 & 545 & 8.5 & 782 & 186 & 49 & 330 + 163246.6 + 340526 & 0.8306 & 880 & 280 & 3.2 & 140 & & & + 212546.9 + 004455 & 1.4240 & 8300 & 82 & 2.9 & 103 & 3 & 6 & 2 +    ccccccc
100552.6 + 493448&0.048 & 15 & 60 & -1.5 & 30 + 153542.4 + 090341 & 0.003 & 15 & 45 & -0.5 & 80 + 163246.6 + 340526 & 0.086 & 35 & 60 & -1.0 & 80 + 212546.9 + 004455 & 0.005 & 15 & 60 & -2 & 50 +
here we examine the properties of a cold cloud where iron is mainly singly ionized .
fe ii emission will be produced by continuum fluorescence under these circumstances .
our scenario is that an intervening cloud has extinguished the agn continuum and that this attenuated sed strikes an outer dust - free cloud .
the intervening cloud is assumed to have a column density of n(h ) = 10@xmath49 @xmath40 and an unattenuated ionization parameter of log @xmath59 ( although this has little affect on our predictions ) .
the sed transmitted through this cloud , which then strikes the outer cloud , is shown in fig a1 .
we assume that the outer cloud is dust free .
this is to allow fe to exist in the gas phase .
fe is strongly depleted in dusty gas ( jenkins 2009 ) so little fe ii emission occurs when grains are present .
it is possible that the cloud was at one time too close to the agn for dust to exist .
the cloud would be similar to a classical pdr ( rollig et al .
2007 ) were grains present . in a pdr
grains heat the gas while absorbing the balmer continuum of the incident sed .
since grains are , by hypothesis , absent in our cloud , the gas will be heated mainly by photoelectric absorption by the gas , mainly of second and third row elements with ionization potentials smaller than hydrogen .
a low gas kinetic temperature results .
the radiation striking the gas is bright at wavelengths longward of the lyman limit , so the main species that will be present in the cloud will be h i , he i , c ii , n i , o i , ne i , na ii , mg ii al ii , si ii , s ii , and fe ii .
typical gas temperatures in our resulting models are around 300 k ( figure [ figa3 ] ) , too cool to collisionally excite optical or uv lines .
uv line emission is produced by continuum fluorescent excitation in this case .
note that it is quite likely that the sed striking the cloud is not visible from our line of sight .
this would occur if parts of the torus lie between us and the central object .
we use version c13.03 of cloudy , a spectral simulation code last described by ferland et al .
verner et al . ( 1999 ) describe the fe ii model in this version of cloudy .
fe ii emission is complex ( baldwin et al . 2004 ) and we report it using the sums over wavelength bands described in that paper .
results are mainly sensitive to the column density of the outer cloud and we will show predictions resulting from varying this parameter below . in this cool environment , fe ii
is excited mainly by absorption of the incident sed by uv lines whose lower levels are close to , or within , the ground term .
the level populations within the ground term are affected by the density which then affects which uv lines are efficiently pumped
. the density does not affect the emission from high levels so long as it is below the critical density of the excited levels , @xmath60 @xmath16 . at lower densities the emission spectrum
is mainly determined by the branching ratios as electrons decay to lower levels .
we experimented with a range of densities but settled on @xmath61 @xmath16 .
lines absorb mainly over their doppler core so the amount of pumping will depend on the line width ( ferland 1992 ) .
we assume 1000 km  s@xmath0  in the examples we show although this does not greatly affect the details .
the fe ii grows stronger as the column density increases and fe ii lines absorb more of the incident continuum .
optical depths in the strong uv lines do affect the branching ratios .
as they grow more opaque decays from their upper levels will favor subordinate lines .
this is quite analogous to the transition from case a to case b in h i emission ( osterbrock & ferland 2006 ) .
finally the intensity of the radiation field give in figure [ figa1 ] must be set .
we assume an unabsorbed ionization parameter of @xmath62 .
results are shown in figure [ figa2 ] .
the lower axis gives the total hydrogen column density .
an iron abundances of log fe / h = -4.55 was assumed .
a good fit is obtained at @xmath63 @xmath40 .
the upper axis shows the optical depth in one of the strong uv lines that is responsible for pumping the atom .
the net fe  ii  emission is shown in [ figa4 ] .
fe  ii  uv 1 , 3 , 4 , 63 , 64 are the strongest lines , while those at wavelengths longer than 2800are much weaker .
the spectrum looks similar to the observed fe  ii  .    with these parameters ,
the thickness of the cloud is one parsec , and the distance to the continuum source is estimated about 1 kpc , following the definition of the ionization parameter , and the inferred intrinsic continuum luminosity of j1535 + 0903 .
the total volume of fe ii region is about @xmath41 parsec@xmath64 , and a total gas mass of about @xmath65 m@xmath3 , which is comparable to the reported @xmath66 m@xmath3 molecular gas in outflows of infrared luminous quasars , such as mrk 231 ( feruglio et al .
2010 ; cicone et al .
however , the gas density is poorly constrained so these numbers are only very rough estimate .


it has become a widespread practice to study the physics of systems in various numbers of spatial dimensions , not necessarily @xmath0 .
for instance , graphene with @xmath1 , and string or membrane theory with @xmath2 as high as @xmath3 , are two examples that immediately come to mind for both their experimental and theoretical interest .
moreover , geometric ideas provide a common framework used to pursue such studies .
a pedagogical goal of this paper is to encourage students to think along these lines in the context of a familiar subject  electrostatics .
for example , in three spatial dimensions a uniformly charged straight - line segment gives rise to an electric potential @xmath4 whose equipotential surfaces are prolate ellipsoids of revolution about the segment , with the ends of the segment providing the foci of the ellipsoid .
as an immediate consequence of the geometry for these prolate ellipsoidal equipotentials , the associated electric field  @xmath5  always normal to surfaces of constant @xmath4  has at any observation point a direction that bisects the angle formed by the pair of lines from the observation point to each of the two foci of the ellipsoid .
this beautiful electrostatic example was presented by george green in 1828 @xcite , and it has been discussed in many books since then @xcite-@xcite including at least two texts from this century @xcite .  while the straight - line segment is an idealization , nevertheless it provides insight into the behavior of real thin - wire conductors , especially upon approximating those real wires as very narrow , needle - like  ellipsoids .
but as it turns out , the line segment problem in three dimensions is a very special case , in some sense the most ideal of all possible electrostatic worlds .  in any _ other _ dimension of space , uniformly charged segments
do * not * produce ellipsoidal equipotentials .
conversely , in any other dimension , if the equipotentials are ellipsoidal about a linearly distributed straight - line segment of charge , then that charge distribution can * not * be uniform .
one need look no farther than two dimensional systems to see clearly that there are differences between uniformly charged segments and those with charge distributed so as to produce ellipsoidal equipotentials .
indeed , @xmath1 is the only ellipsoidal equipotential case which is intuitive in the sense that the associated linear charge distribution has maxima at the ends of the segment , as one might naively expect from the repulsive force between like charges placed on a segment of a real , thin conductor at a finite potential .
in contrast , the distribution of charge needed to produce ellipsoidal equipotentials is counter - intuitive in higher dimensions .
to produce such equipotentials for @xmath6 the charge distribution must have a maximum at the center and vanish at the ends of the segment .
we begin our discussion in section 2 in two dimensions , where the two types of charged segments are readily analyzed .
then we compare and contrast uniformly charged line segments with those admitting ellipsoidal equipotentials for any number of spatial dimensions .  as a preliminary , we first discuss briefly in section 3 the potential of a point charge in @xmath2 dimensions .
we then use this information in section 4 to compute the potentials and electric fields for uniformly charged line segments .
we continue in section 5 by considering systems with ellipsoidal equipotentials in @xmath2  dimensions .
we then determine the linear charge distributions that produce such potential configurations , and we find the remarkable result that the linear charge density giving rise to ellipsoidal equipotentials is peaked at the center of the segment , for @xmath6 .
the discussion of general @xmath2 affords the opportunity to illustrate how continuous @xmath2 can be used as a mathematical device to regulate singular behavior .
this too is a widespread practice in theoretical physics .
we use @xmath2 in this way in section 5 of the paper to interpolate continuously between intuitive and counter - intuitive charge distributions for segments with ellipsoidal equipotentials .     finally , in section 6 , we invoke well - known methods to indicate how the various line segment results also give solutions to a class of electrostatic boundary value problems where the charge is moved outward and distributed on one of the equipotential surfaces that surrounds the original segment .
the point - particle electric potential in 2d is well - known to be logarithmic .  for a point charge @xmath7 located at the origin , up to a constant @xmath8 that sets the distance scale , @xmath9 where @xmath10 is the two dimensional analogue of coulomb s constant .
by linear superposition , a finite length , uniformly charged line segment , with constant charge / length @xmath11 distributed on the @xmath12 axis for @xmath13 , produces a potential @xmath14 this result may be established by integrating the contributions of infinitesimal point - like bits of charge that make up the segment , using ( [ 2dpointcharge ] ) and the indefinite integral @xmath15 written out in full,@xmath16 a plot of the potential surface shows the essential features .
eurojphyslinesegmentsv2__1.pdf + figure 1 :  potential surface for a uniformly charged line segment in 2d .
the top of the potential surface is curved and not a straight line , indicating that the charged segment itself is _ not _ an equipotential .
this follows analytically from ( [ uniformline2d ] ) .
although points on the segment are not at the _ same _ potential , they are all at _ finite _ values of the potential , for @xmath1 .
explicitly , for @xmath17 and all @xmath12,@xmath18 we plot @xmath19 versus @xmath12 to show the shape of the potential along the @xmath12-axis , for @xmath20 and @xmath21 .
eurojphyslinesegmentsv2__2.pdf + figure 2 :  uniformly charged line segment potential for @xmath1 , along the @xmath12-axis .    on the other hand , transverse to the @xmath12-axis the potential has a discontinuous slope across the line segment .
that is to say , the electric field normal to the line segment , @xmath22 , is discontinuous due to the presence of the charge density on the segment .
for example , for @xmath23 this transverse profile is given by@xmath24   \ .\ ] ] we also plot @xmath25 versus @xmath26 to show the shape of the potential along the @xmath26-axis , for @xmath20 and @xmath21 .
eurojphyslinesegmentsv2__3.pdf + figure 3 :  uniformly charged line segment potential for @xmath1 , along the @xmath26-axis .    moreover , the equipotentials are _ not _ ellipses surrounding the segment on the @xmath27-plane , as is especially clear for points close to the segment .
the potential contours actually _ intersect _ the segment .  a view of the potential contours from _ below _ the potential surface shows these features graphically .
( see the plot to follow .
but note the view in that plot is actually an orthogonal projection of the contours onto the @xmath27-plane , and not the true perspective of an observer on the potential axis a finite distance below that plane . )
eurojphyslinesegmentsv2__4.pdf + figure 4 :  potential contours for a uniformly charged line segment in 2d .
the components of the electric field produced by the segment are given by@xmath28 where we have used @xmath29 .
note that @xmath30 is independent of the scale @xmath8 , since changing @xmath8 just amounts to adding a constant to the potential .
but of course @xmath30 does depend on @xmath31 as this sets the physical length scale for the system .    as a check , in the units
we have chosen the charge density along the line follows from @xmath32 , integrated over a horizontal rectangle containing an infinitesimal portion of the @xmath12-axis , as the height of the rectangle is taken to zero .  thus , using @xmath33,@xmath34{cl}0 & \text { \ \ if \ \ } x > l/2\\ \pi & \text { \ \ if \ \ } -l/2<x < l/2\\ 0 & \text { \ \ if \ \ } x<-l/2 \end{array } \right .   \ . \label{uniformcharge2dderivation}\ ] ] that is to say , the charge / length between @xmath35 and @xmath36 is the constant @xmath11 , as expected .
elsewhere , the charge density vanishes , as follows from @xmath37 for all points not coincident with the segment .
it is instructive to make a vector field plot of @xmath38 , especially near the charged segment .
again it is evident graphically that the segment itself is not an equipotential , since the electric field lines are not perpendicular to the segment as the @xmath12-axis is approached for @xmath39 , except at the single point @xmath23 .  here
are vector field plots for @xmath20 .
eurojphyslinesegmentsv2__5.pdf    eurojphyslinesegmentsv2__6.pdf + figures 5 & 6 :  vector field plots of @xmath40 for a uniformly charged line segment situated between @xmath41 and @xmath42 .
the field is evaluated at the center of each arrow .
in contrast to the uniformly charged segment in 2d , consider a distribution of charge along the segment such that equipotentials _ are _ ellipsoidal .
the relevant charge distribution turns out to be@xmath43{cl}0 & \text { \ \ if \ \ } x > l/2\\ 1 & \text { \ \ if \ \ }
-l/2<x < l/2\\ 0 & \text { \ \ if \ \ } x<-l/2 \end{array } \right .   \ , \label{lambda(x)2d}\
] ] as we shall confirm in the following .
since @xmath44 , the total charge on the segment is still @xmath45 , the same as for the uniformly charged case .
the corresponding potential is now @xmath46 where @xmath47 is a sum of two distances , from the observation point @xmath48 to each of the two ends of the segment .
that is,@xmath49 because the positional dependence of the potential is given entirely by @xmath47 , the equipotentials are ellipses on the @xmath27-plane , with the ends of the segment serving as the foci of each equipotential ellipse .
again , we plot the potential surface .
eurojphyslinesegmentsv2__7.pdf + figure 7 :  potential surface for a non - uniformly charged line segment in 2d , with @xmath50 .
the top of the potential surface is now a straight line , indicating that the charged segment itself _ is _ an equipotential , and the points on the segment are at a _
potential for @xmath1 , namely , @xmath51 .
this follows analytically from ( [ ellipsoidalline2d ] ) .
explicitly , for @xmath17 and all @xmath12,@xmath52 we plot @xmath53 versus @xmath12 to show the shape of the potential along the @xmath12-axis , for @xmath20 and @xmath21 .
eurojphyslinesegmentsv2__8.pdf + figure 8 :  non - uniformly charged line segment potential for @xmath1 , along the @xmath12-axis .
transverse to the @xmath12-axis , the potential again has a discontinuous slope across the line segment due to the presence of the charge density on the segment .
for example , for @xmath23 this transverse profile is given by@xmath54 we plot @xmath55 versus @xmath26 to show the shape of the potential along the @xmath26-axis , for @xmath20 and @xmath21 .
eurojphyslinesegmentsv2__9.pdf + figure 9 :  non - uniformly charged line segment potential for @xmath1 , along the @xmath26-axis .    a view of the potential contours from below the potential surface shows the equipotential ellipses .
( but again note the view in that plot is actually an orthogonal projection of the contours onto the @xmath27-plane , and not the true perspective of an observer on the potential axis a finite distance below that plane . )
eurojphyslinesegmentsv2__10.pdf + figure 10 :  ellipsoidal potential contours for a non - uniformly charged line segment in 2d .
the components of the electric field are now@xmath56 whose final form in rectangular coordinates follows from @xmath57 in terms of @xmath58 and @xmath59,@xmath60 writing everything out in terms of @xmath12 and @xmath26 gives@xmath61 the components of @xmath62  are not elegant in these coordinates @xcite but their properties are fully encoded and amenable to machine computation.@xmath63 although the geometric features of the electric field may not be transparent from these expressions , nevertheless the fact that the equipotentials are ellipsoidal allows one to immediately visualize the direction of @xmath62  as normal to those surfaces of constant @xmath4 .    to confirm the charge density along the line , we again use @xmath64 , and integrate over a horizontal rectangle containing an infinitesimal portion of the @xmath12-axis , as the height of the rectangle is taken to zero .
the crucial features here are @xmath65 as @xmath66 for @xmath67 ( away from the segment ) , but @xmath68 as @xmath66 for @xmath39 ( on the segment ) .  more precisely , as points on the segment are approached transversely , @xmath69 thus for @xmath39 , @xmath70 combining this with the more obvious @xmath71 for @xmath39 gives@xmath72{cl}0 & \text { \ \ if \ \ } x > l/2\\ 1 & \text { \ \ if \ \ } -l/2<x < l/2\\ 0 & \text { \ \ if \ \ } x<-l/2 \end{array } \right .   \ . \label{nonuniformcharge2dderivation}\ ] ] that is to say , the charge / length between @xmath35 and @xmath36 is@xmath73 as anticipated above in ( [ lambda(x)2d ] ) .  elsewhere , the charge density vanishes ,
again as follows from @xmath74 for all points not coincident with the segment .
once more it is instructive to make a vector field plot of @xmath75 , especially near the charged segment .
the fact that the segment itself is an equipotential is evident graphically since the electric field lines are perpendicular to the segment as the @xmath12-axis is approached for @xmath39 .  here
are vector field plots for @xmath20 .
eurojphyslinesegmentsv2__11.pdf    eurojphyslinesegmentsv2__12.pdf + figures 11 & 12 : vector field plots of @xmath40 for a non - uniformly charged line segment situated between @xmath41 and @xmath42 .
the field is evaluated at the center of each arrow .
as previously stressed , the geometry of the ellipsoidal equipotentials ensures that the electric field at any observation point always has a direction that bisects the angle formed by the two lines from the end points of the segment to the observation point @xcite .
this result is also manifest in an integral expression for @xmath62 that follows from linear superposition of the field contributions from infinitesimal @xmath76 charges along the segment , upon choosing an appropriate integration variable .  that is,@xmath77 the unit vector @xmath78 points from the infinitesimal charge on the segment to the observation point , with angle @xmath79 measured from the @xmath12-axis in the usual counterclockwise sense on the @xmath27-plane .
thus @xmath80 .
correspondingly ,  @xmath81 are the angles from the ends of the line segment to the observation point , as given by @xmath82 .
since the integration over @xmath83 weights @xmath78 by an _ even _ function of @xmath83 , and the range of integration is _ symmetric _ about @xmath84 , it follows that the resulting direction of @xmath30 will be proportional to @xmath85 which points in direction @xmath86 .
that is to say , @xmath87 the change of variables needed to obtain ( [ efieldpunchline ] ) will be discussed more fully below , for non - uniformly charged segments giving rise to ellipsoidal equipotentials in any number of dimensions .
a direct graphical comparison of the uniformly charged segment and the non - uniformly charged segment in 2d is obtained by superimposing their equipotentials in a true orthogonal projection of the @xmath4 surface contours onto the @xmath27-plane .
eurojphyslinesegmentsv2__13.pdf + figure 13 :  2d equipotential contours for a uniformly charged segment , with constant @xmath11 for @xmath88 , for @xmath89 ,  @xmath90 ,  @xmath91 ,  @xmath92 , @xmath93 , and @xmath94 , as inner to outer black curves , and ellipsoidal contours for a coincident non - uniformly charged segment , for @xmath95 , @xmath92 , and @xmath94 , as inner to outer red curves .    for large distances from the line segments , whether uniformly charged or otherwise , the potential approaches that of a point charge as given in ( [ 2dpointcharge ] ) .
so asymptotically both sets of contours become circles .  regarding this ,
recall the total charge on either segment under consideration is the same , namely , @xmath96 .
thus for large distances from the segment the equipotentials in figure 13 will coalesce , although it is perhaps surprising how rapidly this occurs .  in figure 13
the equipotential contours for both the uniform and non - uniform charge distributions , for the same value of @xmath4 , are very nearly coincident if @xmath97  the exact locations of the two sets of contours never differ by more than a few percent if @xmath97  well before the contours reach their asymptotic circular form .
for a point charge @xmath7 located at the origin of coordinates , in @xmath98 dimensions , the scalar potential at an observation point @xmath99  is hypothesized to be @xcite @xmath100 where @xmath10 is the @xmath2-dimensional analogue of coulomb s constant .
integrating the radial gradient of @xmath101 over the surface of a sphere fixes the normalization , by gauss law.@xmath102 where the total solid angle ( i.e. the area of the _ unit radius _ sphere , @xmath103 , embedded in @xmath2 dimensions ) is given by @xmath104 for example , @xmath105 , @xmath106 , @xmath107 , @xmath108 , etc .    to put it differently , in terms of a dirac delta in @xmath2-dimensions,@xmath109 and in fact
, this gives the correct result even for @xmath1 , by taking a limit :     @xmath110 , hence @xmath111 .
so for @xmath1 the point particle potential is logarithmic , as previously noted in section 2 .
perhaps the simplest convention would be to set @xmath112 which would require , for @xmath98,@xmath113 this is singular at @xmath1 where the potential is a logarithm , not a power .
in that case the corresponding choice would be @xmath114 , so @xmath115 and again @xmath116 .
the line of charge and observation point are shown in red in the following figure .
eurojphyslinesegmentsv2__14.pdf + figure 14 :  segment coordinates .    from the figure ( @xmath117 and @xmath118 in the formulas correspond , respectively , to the bottom point b  and the top point t  in the figure )    @xmath119    the relations in the third line above follow from the law of sines.@xmath120note that the last expression must be used with care since
@xmath121  is multi - valued .  in particular , if @xmath122 , then as the segment is approached it is always true that @xmath123 , no matter how the segment is approached .
more generally , when @xmath124 decreases to cross the surface of the sphere for which the segment is a diameter , the value of @xmath121  increases through @xmath125 .
this follows from @xmath126 so that the argument of the @xmath121 in ( [ arcsin ] ) is just @xmath127 for @xmath128 , and @xmath129 at this radius .
reducing @xmath124 below @xmath36 increases @xmath130 above @xmath125 , i.e. @xmath121 has moved onto another branch of the function .
this change of branch can be explicitly taken into account through the use of heaviside step functions , @xmath131 , to write    @xmath132    where the @xmath121 in this expression is the principal branch of the function .    assuming linear superposition for the potential , a uniformly charged line segment as shown above in figure 14 , of length @xmath31 , centered on the origin , with top ( t ) and bottom ( b ) ends at @xmath133 ( i.e. @xmath134 ) and @xmath135 ( i.e. @xmath136 ) , will produce in @xmath98 dimensions a potential @xmath137 given by @xcite @xmath138 here @xmath11 is the constant charge / length on the line segment , and @xmath139 are the polar angles for vectors from the top and bottom endpoints of the line segment to the observation point @xmath140 , as in figure 14 .
the potential has no dependence on the additional @xmath141 angles needed to specify the location of a point in @xmath2 dimensions using spherical polar coordinates .
that is to say , in @xmath2 dimensions the equipotentials are always higher dimensional surfaces of revolution about the line segment .
the corresponding electric field is given by@xmath142 and it depends manifestly on @xmath124 and @xmath143 .  the only dependence of the electric field on the additional @xmath141 angles in @xmath2 dimensional spherical polar coordinates is carried by the unit vectors @xmath144  and @xmath145 .    to see that ( [ phid ] ) is correct , we need only sum the potential contributions for infinitesimal pieces of the line segment with charge @xmath146 located on the vertical axis of figure 14 at position @xmath147 , and at a distance @xmath148 from the observation point @xmath149 as given by @xmath150 but writing @xmath151 and @xmath152 , for fixed @xmath153 and @xmath154 , we have @xmath155 , and @xmath156 on the other hand , @xmath157 .  hence the result ( [ phid ] ) .
for instance , in the special case @xmath0 the potential of the uniformly charged segment is @xmath158 note that the potential is infinite for all points on the charged segment itself , for which points @xmath159 and @xmath160 .
be that as it may , after a bit of algebra @xmath161 can be reduced to @xmath162 where @xmath163 is the sum of the distances from the end points of the segment to the observation point .
thus equipotentials in this case are given by constant @xmath47 , and as is common knowledge , this defines an ellipsoid of revolution about the line segment with focal points @xmath164 and @xmath165 .
a view  from below the potential surface again clearly shows the equipotentials are ellipses .
eurojphyslinesegmentsv2__15.pdf + figure 15 : @xmath166 for @xmath20 , where @xmath167 and @xmath168 .
eurojphyslinesegmentsv2__16.pdf + figure 16 :  contours of constant @xmath169 for @xmath20 , plotted versus @xmath167 ( vertical axis ) and @xmath168 ( horizontal axis ) .    for other dimensions , however , the geometrical shapes of the equipotentials are not so easily discerned .  for @xmath6 , it remains to evaluate the angular integral in ( [ phid ] ) .
define the indefinite integral @xmath170 , and compute@xmath171{ll}i\left (   4\right )   = \vartheta\ , & i\left (   5\right )   = -\cos\vartheta\ , \\ & \\ i\left (   6\right )   = \frac{1}{2}\vartheta-\frac{1}{2}\cos\vartheta\sin \vartheta\ , & i\left (   7\right )   = \frac{1}{3}\cos^{3}\vartheta-\cos \vartheta\ , \\ & \\ i\left (   8\right )   = \frac{3}{8}\vartheta+\left (   \frac{1}{4}\cos^{3}\vartheta-\frac{5}{8}\cos\vartheta\right )   \sin\vartheta\ , \ \ \ \ \   & i\left (   9\right )   = -\frac{1}{5}\cos^{5}\vartheta+\frac{2}{3}\cos^{3}\vartheta-\cos\vartheta\ , \end{array}\ ] ] etc .  for odd @xmath172
the integral is always a polynomial of order @xmath173 in @xmath174 , while for even @xmath175 the integral always has a term linear in @xmath176  plus , for @xmath177 , a term with @xmath178 multiplying a polynomial of order @xmath179 in @xmath174 .  for any @xmath2
we therefore obtain the angular integral in ( [ phid ] )  in terms of @xmath180 , @xmath181 , @xmath182 , @xmath183 , and @xmath184 .
these quantities may then be expressed in terms of @xmath143 and @xmath124 upon using the relations in ( [ trig]-[arcsin ] ) .
for example , in @xmath185 the potential of the uniformly charged segment is@xmath186 unlike @xmath187 , the potential @xmath188 is _ not _ a function solely of the variable @xmath47 as defined in ( [ s ] ) .
consequently , the equipotentials are _ not _ ellipsoidal for this four dimensional example .
indeed , a plot now shows that the equipotentials are not ellipsoidal , especially for points close to the uniform line of charge , although the shape of the potential surface is similar to that for @xmath0 .
eurojphyslinesegmentsv2__17.pdf + figure 17 : @xmath189 for @xmath20 , where @xmath167 and @xmath168 .
note that @xmath190 is infinite for all points on the segment .  also note , to obtain the correct shape of the potential surface it is important to interpret @xmath121 in the formula for @xmath190 as a multi - valued function , as given explicitly by ( [ arcsinmulti ] ) .
a view from below the potential surface shows that the potential contours are _ not _ ellipses in this case , although the difference is somewhat subtle .
this next graph should be compared closely to the ellipsoidal case for @xmath185 as presented below in section 5 .
eurojphyslinesegmentsv2__18.pdf + figure 18 :  contours of constant @xmath191 for @xmath20 .
the electric field is given by ( [ linee ] ) and in this case has @xmath124 and @xmath143 components@xmath192 all the complication in these expressions lies in the derivatives @xmath193  and @xmath194 whose geometrical significance is not yet easy to visualize .  on the other hand , the last terms involving @xmath195 in ( [ uniformer ] ) and ( [ uniformetheta ] ) have simple geometrical interpretations since they define a vector that is perpendicular to the axis of the segment and points away from that axis for positive @xmath11  ( i.e. in the direction of @xmath196 if we were to use cylindrical coordinates ) .    in any case , it is worthwhile to compute @xmath197  and @xmath198 since these derivatives appear for all even @xmath175 and always have the same form .
after some algebra we find @xmath199 note the difference in sign for @xmath200 .
once again @xmath196 , so the bulk of the contributions from these derivatives again gives a vector that is perpendicular to the axis of the segment .
but there remains a contribution that is entirely in the @xmath144 direction if @xmath201 is written as @xmath202 , or else entirely in the @xmath145 direction if @xmath203 is written as @xmath204 .
choosing the first of these options gives@xmath205 to compare to the three dimensional ellipsoidal case , it suffices to express @xmath144 and @xmath145 in terms of @xmath206 and @xmath207 .
although the latter two unit vectors are not orthogonal , they are independent except at @xmath133 and @xmath135 , and at those particular angles either one of @xmath206 and @xmath207  will suffice to give the direction of the electric field , since @xmath208 .  we find@xmath209@xmath210
consider the same geometry as in figure 14 , but rather than assuming uniform charge density on the line segment , suppose the potential produced by the line segment depends positionally only on @xmath47 in any number of dimensions .  including some convenient numerical and dimensionful factors , suppose@xmath211 in this case , equipotentials are ellipsoidal by assumption , since the positional dependence is only on @xmath47 .
the only issue is to determine @xmath212 .
if the only charge present is on the line segment , then for other points the potential function must be harmonic , @xmath213 .
so , using @xmath214 and @xmath215 ( see eqns(5)-(7 ) in @xcite  for more details ) , we compute@xmath216 the potential will then be harmonic for points not lying on the segment if and only if@xmath217 for @xmath218 .  for general @xmath2 _ the relevant solution _ of this second order , ordinary differential equation ,
is given by a gauss hypergeometric function represented by the standard series around @xmath219 , namely,@xmath220 the various factors in the last line are  useful to determine the exact expression for the charge distribution along the segment , as presented below ( see ( [ chargedistribution ] ) ) .
note the large @xmath47 behavior of @xmath221 implies@xmath222 since @xmath223 .
the normalization of @xmath221 and the @xmath224 factors in ( [ ellipsoidalanyd ] ) were chosen so that the total charge on the segment is the same as in the uniformly charged case , namely , @xmath96 .
this total charge appears in the limiting form of the potential for @xmath225 , i.e. @xmath226 .
alternatively , closed - form expressions for @xmath212 in terms of more elementary functions are sometimes more easily obtained by evaluating the corresponding hypergeometric series around @xmath227 , as given by @xcite@xmath228 and then analytically continuing to @xmath218 .  here
@xmath229 and @xmath230 are constants chosen so that the result ( [ alternate3pointpotential ] ) is real - valued for @xmath218 , and so that any constant term ( the trivial harmonic ) is eliminated from @xmath221  as @xmath231 .    for integer
@xmath2 the hypergeometric functions that appear in the solutions ( [ 3pointpotential ] ) and ( [ alternate3pointpotential ] ) always reduce to elementary functions .
for example , in various dimensions the relevant harmonic functions are:@xmath232 by construction , the equipotentials are always ellipsoidal for any @xmath2 , except for the trivial one dimensional case where equipotentials consist of just pairs of points on the line .
in contrast to the uniformly charged line segment for @xmath233 , if the potential depends only on @xmath47 then the charged line itself  where @xmath234 for all points on the segment  is _ always _ an equipotential , albeit with _ infinite _
@xmath4 when @xmath235 , while for @xmath236 and @xmath1 the potential along the charged line is _
finite_.    for example , consider graphically the case for @xmath185 .
qualitatively the potential surface is similar to that for the uniformly charged segment in three dimensions .
potential contours are ellipsoids that surround @xmath237 and _ never _ intersect the charge segment .  also , the top of the potential surface is just an exact copy of the segment itself , albeit at an infinite value of @xmath4 .  that is to say ,
the charged segment is itself an equipotential , but in fact the potential is infinite for points on the segment .  again , a view from below the potential surface shows more clearly the ellipsoidal potential contours .
eurojphyslinesegmentsv2__19.pdf + figure 19 : @xmath238 for @xmath20 , versus @xmath167 and @xmath168 .
eurojphyslinesegmentsv2__20.pdf + figure 20 :  contours of constant @xmath239 for @xmath20 .
admittedly , in 4d it takes a discerning eye to see differences in the shape of the ellipsoidal equipotential surface compared to that for the uniformly charged segment , if graphs of the two cases are viewed separately .
but if viewed side - by - side the difference _ is _ evident in contour plots for the potentials .  perhaps even more clearly , the difference is highlighted by plotting both cases in the same graph .
eurojphyslinesegmentsv2__23.pdf + figure 21 :  comparison of  @xmath4 contours in 4d for a uniformly charged segment ( black ) and a non - uniformly charged segment with ellipsoidal equipotentials ( red ) for @xmath20 .    upon doing
so , it is also evident that the non - uniform charge distribution must be greater near the center of the segment rather than at the ends of the segment , a counter - intuitive feature .
the ellipsoidal equipotentials are shorter in the direction of the line segment , and wider transverse to the segment , than those of the uniformly charged line for the same value of @xmath240 .
this is exactly the opposite of what happens in 2d , where the ellipsoidal equipotentials were longer and narrower than those of the uniformly charged segment for the same value of @xmath4 .     the non - uniform charge distribution on the segment that produces ellipsoidal equipotentials
can now be determined using the integral form of gauss law in @xmath2 dimensions , ( [ integralgausslawind ] ) , in complete parallel to the calculation leading to ( [ nonuniformcharge2d ] ) in 2d .
here we only give the results of that calculation .    in @xmath2 dimensions the charge density for a line segment that produces ellipsoidal equipotentials is@xmath241 where we have used the ratio of total solid angles  in @xmath2 and @xmath242 dimensions,@xmath243 for example , @xmath244 to give @xmath245 in agreement with ( [ lambda(x)2d ] ) , while @xmath246 to give the uniform distribution , @xmath247 .
these charge densities are all normalized so that the total charge on the segment is the same for any @xmath2 , namely,@xmath248 this follows from@xmath249    a geometrical construction to obtain the distribution ( [ chargedistribution ] ) is to consider a uniformly charged @xmath103 , i.e. a hypersphere embedded in @xmath2 dimensions , with radius @xmath250 and with _ constant _ hypersurface charge density @xmath251 .
an orthogonal projection of the charge , @xmath252 , from a hyper - cylindrical ribbon  of revolution about a diameter of the hypersphere onto the underlying surrounded segment @xmath253 of the diameter , gives precisely @xmath254 for @xmath255 .
this follows directly from @xmath256 , where the hypersurface area  @xmath257 of the hyper - ribbon that surrounds @xmath253 is @xmath258 with @xmath259 , @xmath260 , and @xmath261 .
this construction is easily visualized for @xmath1 and @xmath262 .
we plot the densities for various integer @xmath2 , from @xmath1 up to @xmath263 , shown respectively as the lower to upper curves ( for @xmath264 ) in the following graph .
eurojphyslinesegmentsv2__24.pdf + figure 22 :  linear charge density @xmath265 leading to ellipsoidal equipotentials , in various dimensions .    in this graph
the 2d case is the only one which is intuitive in the sense that the infinitesimal pieces of charge repel each other so we would expect excess charge / length to be pushed towards the ends of a finite length , equipotential segment of real conductor with small but nonzero transverse size , all parts of which are at the same finite @xmath4 .  for the idealized line segment in 3d
this is not so  the charge density is uniform  but then the segment itself is at an infinite @xmath4 so physical intuition based on real conductors is perhaps difficult to apply in this situation @xcite .  in higher dimensions , the result is even _ more _ counter - intuitive for idealized line segments .
the charge distribution is peaked at the center of the segment .
of course , in making these statements , we are assuming the physical properties of the idealized segment itself can be completely understood mathematically by taking the limit of the surrounding equipotential ellipsoids of revolution as their girth goes to zero .    to gain more insight about the behavior of the charge distribution and the resulting potential
, we may think of @xmath2 in ( [ chargedistribution ] ) as a continuous variable  not just an integer  to be used as a regulating parameter by which the @xmath0 case can be approached as a limit .
this point of view shows that all @xmath266 behave intuitively , with @xmath267 peaked at the ends of the segment .
the extreme case is @xmath236 where all charge is located _ only _ at the ends of the segment .
moreover , for all @xmath266 the potential along the segment , as given by ( [ 3pointpotential ] ) , is a _
finite _ constant , so there is nothing pathological about the potential that might obviate physical intuition based on real conductors .
on the other hand , all @xmath6 behave counter - intuitively , with @xmath268 peaked at the center of the segment and vanishing at the segment ends .
and for all @xmath6 the potential along the segment , as given by ( [ 3pointpotential ] ) , is _ infinite _ , so real - world physical intuition is not guaranteed to be reliable for these idealized situations .
the uniformly charged 3d case , also with infinite @xmath4 along the segment , acts as a separatrix between intuitive and counter - intuitive charge distributions .
a plot of @xmath265 as a surface over the @xmath269 plane helps to visualize these features .
eurojphyslinesegmentsv2__25.pdf + figure 23 :  the @xmath265 surface as a continuous function of both @xmath147 and @xmath2 for a straight - line segment from @xmath270 to @xmath271 .
constant density for a uniformly charged @xmath0 segment is represented by a thick red line .
we close this section with a discussion of the direct calculation of the potential and the electric field as linear superpositions of the @xmath272 s and @xmath273 s due to infinitesimal bits of charge along the segment , @xmath274 , for a line segment in @xmath2 dimensions .  in terms of the obvious rectilinear coordinates ,
as shown in figure 14 of section 4 , the results are@xmath275@xmath276 where @xmath154 is the transverse distance from the segment and @xmath153 is the @xmath147-coordinate of the observation point , and where @xmath277 is a unit vector pointing from the location @xmath278  of the bit of charge on the segment to the observation point @xmath99 .
now , using the charge distribution ( [ chargedistribution ] ) that produces ellipsoidal equipotentials , we eventually obtain@xmath279@xmath280 where we have changed integration variables from @xmath147 to an angle @xmath176 according to@xmath281 note that @xmath176 here is the polar angle _ measured from the location of the bit of charge _ , and in general this is not the @xmath143 shown in figure 14 of section 4 .
nonetheless , @xmath282 and @xmath283 _ are _ the angles shown in figure 14 .
these last relations , in particular ( [ crucialid ] ) and ( [ electricfieldthetaintegral ] ) , are useful to establish _ the direction _ _ of the electric field _ _ without actually performing the integration_.  to verify this statement , let @xmath284 in ( [ electricfieldthetaintegral ] )  to find@xmath285 the integration here weights @xmath78 by an _ even _ function of @xmath83 , with the range of integration symmetric about @xmath84 , so it follows that the resulting direction of @xmath286 will be proportional to @xmath287 .  that is to say , at any observation point around the segment the direction of @xmath288 bisects the angle formed by the two lines from the end points of the segment to the observation point , thereby confirming that the equipotentials are ellipsoids of revolution about the segment .
note that this direct calculation of @xmath288 using the integral expression ( [ elineind ] ) generalizes prior calculations  for @xmath0 where the charge density is uniform , as given in @xcite , to cases where the charge distribution on the segment is _ not _ uniform for @xmath233 .
we leave it as an exercise for the reader to evaluate the integral in ( [ elineind ] ) to obtain explicit expressions , and to show that these are in agreement with @xmath289 .  even before attempting to evaluate the integral in ( [ elineind ] ) , however , it should be evident that a determination of @xmath288 by taking the gradient of the potential is easier to carry through .
having found the equipotentials for various charged line segments , we have in hand a set of solutions for a variety of electrostatic boundary value problems .
if all the charge from the line segment is moved outward along electric field lines and placed on a single surface selected from among the equipotentials , in such a way as to preserve @xmath4 on that surface , then the surrounding equipotentials are unchanged as well .
this is widely known , and was in fact a feature emphasized by green @xcite .
the charge distribution on the selected surface , as required to carry out this feat , is of course given by the electric field normal to that surface as provided by the original charged line segment solution .
in particular , ellipsoidal equipotentials in any @xmath2 yield hypersurface charge densities that are expressible as elementary functions .
the results are straightforward generalizations of the well - known 3d ellipsoidal case .
in fact , it is not even necessary for the hyperellipsoids to be surfaces of revolution .  for a general ellipsoid , as defined in @xmath2 dimensions by @xmath290 the charge distribution on that hyperellipsoid such that it is an equipotential is given by@xmath291 this can be established by a straightforward adaptation of the argument that applies to the 3d case , as given for example in smythe @xcite , chapter 5 ,  5.00-5.02 .
if this charged hyperellipsoidal surface is squashed  to obtain an equipotential hyperdisk ( see smythe  5.03 for the 3d case ) , say by letting @xmath292 while maintaining the constraint ( [ ellipsoidd ] ) , then the hypersurface charge density on the resulting hyperdisk ( counting charge on both sides of the disk ) becomes@xmath293 this charge distribution may now be projected along any of the remaining unsquashed principal axes of the original ellipsoid by elementary integrations .
the final result is exactly the linear charge distribution given by ( [ chargedistribution ] ) , namely , @xmath294 where the projection has been made onto the @xmath10th axis .
alternatively , the surface charge distribution ( [ ellipsoidsurfacecharge ] ) of the unsquashed ellipsoid may be projected directly without being squashed , again by elementary integrations , to obtain the same result for @xmath295 .
in the context of a familiar subject  electrostatics  we have carried out several elementary calculations in various numbers of spatial dimensions , @xmath2 , to encourage students to think more critically about the role played by @xmath2 .
although for @xmath0 a uniformly charged straight - line segment gives rise to an electric potential @xmath4 whose equipotential surfaces are prolate ellipsoids of revolution about the segment , with the ends of the segment providing the foci of the ellipsoid , this is a very special case .
if @xmath233 , uniformly charged segments do not produce ellipsoidal equipotentials .  if @xmath233 , a non - uniform distribution of charge is required to produce equipotentials that are ellipsoidal about a straight - line segment of charge .
we have illustrated these @xmath233 features in detail for @xmath1 and @xmath185 , and we have provided a framework as well as explicit formulas to carry through the same level of detail for any @xmath2 .  we have shown how @xmath1 is the only ellipsoidal equipotential case which is intuitive in the sense that the associated linear charge distribution has maxima at the ends of the segment , as one might naively expect .
in contrast , we have shown that the distribution of charge needed to produce ellipsoidal equipotentials is counter - intuitive for @xmath6 , becoming all the more so as @xmath2 is increased , in the sense that the requisite charge distribution has an absolute maximum at the center of the segment and vanishes at its ends .     for @xmath0 and @xmath185 ,
the potential surface plots were found to have similar features that can be distinguished only by careful inspection .
differences in @xmath4 between uniformly charged line segments and those giving rise to ellipsoidal equipotentials were shown to be somewhat subtle .
distinctions between the two cases were more easily drawn by examining the corresponding electric fields close to the segment , or , relatedly , by directly comparing the charge densities .
we also took the opportunity to illustrate , albeit briefly , how continuous @xmath296 can be viewed as a mathematical tool to regulate singular behavior and to interpolate between intuitive and counter - intuitive charge distributions .  in the course of our discussion
, we employed some basic geometrical ideas at an elementary level , such as in the determination of the direction of the electric field associated with ellipsoidal equipotentials .    finally , we indicated how the line segment results may be generalized to solve a class of electrostatic boundary value problems by distributing the charge on hypersurfaces , rather than straight lines , while maintaining the line segment equipotentials outside the charged hypersurface .                            for more on the history of the 3d line segment problem , including citations of some recent pedagogical papers on the subject , see t s van kortryk , on the fields due to line segments  arxiv:1410.6832 [ physics.hist-ph ] .    throughout our discussion
we use rectangular cartesian , spherical polar , and cylindrical coordinates , in view of student familiarity with these coordinate systems , even though prolate ellipsoids of revolution are more naturally described in terms of prolate spheroidal coordinates .     it is tempting to introduce the variable @xmath297 ,
the coordinate locally orthogonal to @xmath47 in a prolate spheroidal coordinate system .
this would make the result for @xmath62  just as compact as that for @xmath190 when the equipotentials are ellipsoidal
.  however , the results are sufficiently manageable in terms of rectangular cartesian coordinates , at least for graphical purposes , that there is no real numerical advantage achieved by using prolate spheroidal coordinates .
beginning students can appreciate why the normal to a prolate ellipsoidal surface bisects the angle formed by the two vectors from the foci to any point on the ellipsoid by thinking about the law of reflection applied to a light ray that strikes an ellipsoidal mirror while en route from one focus to the other .
alternatively , a few banked shots on any good ellipsoidal billiard table should clarify the situation .
the form in ( [ alternate3pointpotential ] ) is also more useful to obtain the correct result for @xmath1 since the form in ( [ 3pointpotential ] ) requires taking a limit @xmath298 to obtain a nontrivial result .
just how real conductors of finite length and girth reduce to uniformly charged ideal line segments in 3d is a venerable problem with a dignified history , beginning with j c maxwell , on the electrical capacity of a long narrow cylinder , and of a disk of sensible thickness  proc .
london math .
ix ( 1878 ) 94101 , reproduced as paper xcii , pp 672 - 680 , in _ the scientific papers of james clerk maxwell _ , vol .
ii , w d niven , editor , cambridge university press ( 1890 ) .  for a guide to the literature since then
, see j d jackson , charge density on a thin straight wire : the first visit  am .
70 ( 2002 ) 409 - 410 .