ngc 1232 is a nearly face  on , gas  rich , sc galaxy .
face  on galaxies provide excellent targets for studies of radial properties , such as the chemical enrichment of the interstellar medium ( e.g. , mccall , rybski , & shields 1985 ; vila  costas & edmunds 1992 ; zaritsky , kennicutt , & huchra 1994 ; ferguson , gallagher , & wyse 1998 ) or the velocity dispersion of the neutral gas ( e.g. , van der kruit & shostak 1982 , 1984 ; shostak & van der kruit 1984 ; dickey , hanson , & helou 1990 ; boulanger & viallefond 1992 ; kamphuis & briggs 1992 ; rownd , dickey , & helou 1994 ) .
ngc 1232 first came to our attention as part of a study of radial trends in elemental abundances ( van zee , salzer , & haynes 1998a ; van zee 1998b ) .
the outermost region observed by van zee ( 1998a , b ) has an abnormally low oxygen abundance and a higher n / o ratio than expected at that radius .
such features may arise from asymmetries in the disk , or from perturbations induced by tidal interactions .
we thus decided to investigate the neutral gas distribution and kinematics of ngc 1232 to determine if there were kinematic peculiarities associated with this particular region , or more generally in the outer gas disk .    with a relatively low inclination angle , approximately 30 , ngc 1232
is also a prime candidate for studies of the neutral gas velocity dispersion as a function of radius .
since the hi distribution typically extends a factor of 1.3 to 2 times the size of the optical disk in sc galaxies ( e.g. , broeils & rhee 1997 ) , one might naively expect a decrease in the observed velocity dispersion beyond the radius of active star formation [ star formation activity may not be truncated in spiral galaxies , however ( ferguson 1998b ) ] .
previous studies of the neutral gas velocity dispersion in face  on spiral galaxies have been inconclusive .
while some observations indicate that the velocity dispersion may decrease with radius ( e.g. , ngc 6946 , boulanger & viallefond 1992 ) , others indicate that the velocity dispersion is constant in the outer gas disk ( e.g. , ngc 1058 , dickey 1990 ) .
such studies are hampered by two main problems : ( 1 ) rotational broadening and ( 2 ) insufficient spectral resolution .
the observations of ngc 1232 presented in this paper have adequate spectral resolution , but suffer from significant rotational broadening , in part due to a severe warp in the gas disk .
nonetheless , these observations do provide additional insight on radial trends of the neutral gas velocity dispersion in spiral galaxies .
the physical parameters of ngc 1232 are summarized in table [ tab : global ] . throughout this paper ,
we adopt a distance of 21.5 mpc ( van zee 1998b ) , based on an assumed h@xmath1 of 75  mpc@xmath2 and a virgocentric infall model .
optically , ngc 1232 is a rather nondescript sc galaxy with well defined spiral arms and optical colors ( de vaucouleurs 1991 ) that are typical for late  type spiral galaxies .
its total hi mass and m@xmath3/l@xmath4 ratio are also quite typical of the class . on the other hand , it s far infrared luminosity ( soifer 1989 ) is significantly lower than typical for sc galaxies ( see , e.g. , roberts & haynes 1994 for trends as a function of hubble type ) .    this paper is organized as follows .
the hi synthesis observations and data reduction are presented in section [ sec : obs ] . the neutral gas distribution and kinematics of ngc 1232
are discussed in section [ sec : ana ] .
a brief discussion and summary are presented in sections [ sec : disc ] and [ sec : conc ] .
finally , the results of optical spectroscopy of the spatially nearby system ngc 1232a are summarized in an appendix .
hi synthesis imaging observations of ngc 1232 were conducted by the nrao summer students with the dnc configuration of the very large array ( vla ) on 1992 july 4 .
the observations spanned 5 hours , with a total on  source integration time of 227 minutes .
the correlator was used in 2ad mode with the right and left circular polarizations tuned to 1684 .
the total band width was 1.56 mhz .
the on  line hanning smoothing option was selected , producing final spectral data cubes of 127 channels , each 2.6  wide .
standard tasks in aips ( napier 1983 ) were employed for calibration and preliminary data reduction .
calibration of this data set was complicated by the fact that the west arm of the array was shadowed during the observations of the primary flux density calibrator , 3c 286 . however , the observations of 3c 286 were sufficient to set the flux density scale ; phase and bandpass calibration was derived from observations of a nearby continuum source , b0237233 . with an observed flux density of 6.21 jy at 1412.5 mhz , the observations of b0237233 had adequate signal  to  noise for calibration of both the phase and the bandpass .
after calibration , the line data was transformed to the _ x  y _  plane .
a robust weighting technique was employed by the aips task imagr to optimize the beam shape and noise levels of the final data cubes ( briggs 1995 ) .
the `` robustness parameter '' in imagr controls the weighting of the _ u  v _  data , permitting a fine  tuning between sensitivity and resolution .
as currently implemented , a robustness of 5 corresponds to natural weighting of the _ u  v _  data ( maximizing sensitivity ) while a robustness of 5 corresponds to uniform weighting ( lower sensitivity , but better spatial resolution ) .
the relevant imagr parameters for three data cubes are listed in table [ tab : maps ] . throughout this paper
we will refer to the robustness of 1.0 cube as the `` low resolution cube , '' to the robustness of 0.5 cube as the `` intermediate weight cube , '' and to the robustness of 0.5 cube as the `` high resolution cube . ''
all subsequent analysis of the data cubes was performed within the gipsy package ( van der hulst 1992 ) .
line emission from ngc 1232 filled most of the bandpass ( channels 7 to 121 ) . since the few line  free channels were at the extreme edge of the bandpass , they had poor noise characteristics and thus were insufficient to create an accurate model of the continuum emission in either the _ u  v _  or _ x  y _  planes .
continuum subtraction was a necessary step in the data reduction process , however , since ngc 1232 has 58.5 mjy of continuum emission at l  band ( condon 1987 ) .
continuum subtraction was conducted in the _ x  y _  plane , with a continuum image created separately for each data cube . to create the continuum image , the mean of the inner 100 channels ( 1553.8 @xmath5 1814.3 )
was computed .
this image was severely contaminated by the line emission of the galaxy ; the region of galaxy emission was blotted from the image and replaced with the value 0 .
next , the mean was computed for channels with velocities between 1538.2 and 1561.6  ( line emission from only the eastern half of the galaxy ) and for channels with velocities between 1811.7 and 1835.2  ( line emission from only the western half of the galaxy ) .
the continuum sources coincident with the galaxy emission were recovered from these two maps .
the final continuum image for each resolution was inspected to verify that all of the continuum sources were recovered ; an example of one of these continuum images ( from the intermediate weight data cube ) is shown in figure [ fig : cont ] .
the continuum image for each resolution was subtracted from every channel of the data cube to produce the final continuum  free data cubes . selected channels from the continuum
subtracted intermediate weight data cube are presented in figure [ fig : chans ] . the continuum subtraction appears to have been successful ( see , for instance , the baseline of the integrated flux density profile , figure [ fig : flux ] ) .
note , however , that even if the continuum subtraction had not been as successful , the subsequent analysis of the gas kinematics and dynamics would still be valid . in particular , the velocity field and velocity dispersion measurements are fairly robust , even in the presence of continuum emission ; only the total gas distribution would have been significantly effected by poor continuum subtraction .    to determine if the total hi flux density was recovered in the hi synthesis observations , the flux density was measured in each channel of the continuum subtracted data cubes .
the flux density profile , corrected for primary beam attenuation , is shown in figure [ fig : flux ] .
qualitatively , the flux density profile is similar to those obtained from single dish observations ( fisher & tully 1981 ; reif 1982 ; staveley
smith & davies 1988 ) .
the total flux density recovered from the low resolution data cube was 122.4 @xmath0 12.2 jy , in good agreement with previous single dish measurements ( e.g. , 115.5 @xmath0 7.2 jy , green bank 43 m [ fisher & tully 1981 ] ; 127.9 @xmath0 13.8 jy , jodrell bank [ staveley  smith & davies 1988 ] ) .
thus , despite the calibration problems and missing short spacings , the flux density scale appears to be correct , and the total hi flux density has been recovered .    the slight asymmetry between the two horns of the integrated hi profile has been ascribed to emission from the spatially nearby galaxy ngc 1232a ( eso 547g016 ) ( e.g. , reif 1982 ; becker 1988 ) .
however , there is no apparent gas density enhancement or perturbation in the velocity field near this object ( see figure [ fig : mom ] ) .
furthermore , a heliocentric optical velocity of @xmath6  was measured from several emission lines with the palomar 5 m telescope ; further details of the emission lines observed in this galaxy are presented in the appendix .
this newly derived redshift is in basic agreement with the velocity listed in ned ; thus , this object does not appear to be physically associated with ngc 1232 .
both moment analysis and gaussian fitting were used to characterize the gas distribution , velocity field , and velocity dispersion of ngc 1232 .
moment maps of each data cube were computed in the following manner .
first , the low resolution cube was smoothed to a resolution of twice the beam ; second , the smoothed cube was clipped at the 2@xmath7 level ; the resultant clipped cube was then interactively blanked to remove spurious noise spikes .
signal was identified based on spatial continuity between channels .
a conditional transfer was used to blank the corresponding locations in the low , intermediate , and high resolution data cubes ( corrected for primary beam attenuation ) based on the blanked , smoothed cube .
moment maps of the blanked data cubes were created with the gipsy task moments .
in addition , gaussian fits to the observed line profiles at each position in the unclipped data cubes were obtained with the gipsy task gaufit . in regions of high signal  to  noise , gaussian fitting is preferred over moment analysis since it is less susceptible to beam smearing .
the resultant maps were qualitatively similar to the moment maps in both gas distribution and velocity field .
however , significant differences were found in the velocity dispersion maps , probably due to beam smearing .
thus , the moment maps were used for investigation of the gas distribution and kinematics , while interactive gaussian fitting was employed for study of the velocity dispersion as a function of position .
the zeroth and first moments are presented in figure [ fig : mom ] ; the optical images used in this figure are described more fully in van zee ( 1998b ) .
there are several notable features in the gas distribution of ngc 1232 .
first , the hi distribution extends well beyond the optical radius . at a column density limit of 10@xmath8 atoms @xmath9 , the hi  to  optical diameter ratio is 1.6 , quite typical of spiral galaxies ( e.g. , cayatte 1994 ; broeils & rhee 1997 ) .
second , as has been seen in many other massive galaxies ( e.g. , brinks & bajaja 1986 ; broeils & van woerden 1994 ) , the neutral gas column density is lower in the center of the galaxy .
the azimuthally averaged radial profile of the gas density is shown in figure [ fig : rot]d and tabulated in table [ tab : rot ] .
as seen in figure [ fig : rot]d , the eastern and western sides of the galaxy are not perfectly symmetric .
the eastern half of the galaxy has a higher column density , corresponding to the prominent spiral arms on that side . on the western side
, the gas distribution extends to larger radii , perhaps due to a warp of the gas disk in the outer regions .
the observed asymmetry between the horns of the global hi profile is the result of the convolution of this asymmetric gas distribution with the rotation of the disk .
such asymmetries are not expected to be long lived , but also are not uncommon in spiral galaxies ( e.g. , haynes 1998 ) .
closed contours in the velocity field ( figure [ fig : mom]d ) suggest that the outer gas disk is warped , causing an apparent decline in the rotation curve . a model for the rotation curve of ngc 1232
was derived using tilted  ring analysis of the intermediate resolution velocity field .
the process of fitting a rotation curve is iterative ; we began by fitting both sides of the velocity field in order to determine the kinematic center and systemic velocity ( table [ tab : global ] ) .
next , the center and systemic velocity were held fixed and each side was fit separately to determine the position and inclination angles .
the position angle was well constrained by the models .
however , the fit to the inclination angle was relatively unconstrained . based on the optical axial ratios ,
the inclination angle was initially assumed to be approximately 30 .
since the rotcur fits resulted in wild variations of the inclination angle , the inclination was set to 30  in the inner regions of the galaxy and then permitted to slowly decrease in the outer regions .
finally , the rotation curve for each side was determined by fixing the inclination and position angles to the average of both sides with the additional constraint that these values must change smoothly from ring to ring .
the result of rotcur s fit to the approaching and receding sides are shown in figure [ fig : rot ] ; table [ tab : rot ] lists the average values for both sides . as seen in the upper panel of figure [ fig : rot ] ,
a warped disk model is able to recover a flat rotation curve in the outer gas disk .
however , since the inclination angle was poorly constrained by the tilted  ring models , we can not rule out the possibility that the inclination angle is constant at all radii ( i.e. , the gas disk may not be warped ) .
thus , a rotation curve was also derived using the same parameters as the warped disk model , with the exception that the inclination angle was held fixed at 30  throughout the disk ; not unexpectedly , this model results in a falling rotation curve in the outer disk .
figure [ fig : pv ] shows the position  velocity diagram for a slice along the major axis of ngc 1232 . the rotation curve derived from the warped disk model is superposed .
as frequently found in spiral galaxies , the warp appears to begin at the edge of the optical disk ( see e.g. , ngc 4013 , sparke & casertano 1988 ; ngc 4138 , jore , broeils , & haynes 1996 ) .
both the warped disk and falling rotation curve models adequately reproduce the major features of the velocity field .
however , several kinks in the iso  velocity contours are found in the outer regions of the gas disk which can not be well fit by either of these models .
such kinks are usually associated with spiral streaming motions , but these features occur in regions well beyond the prominent spiral arms .
note , however , that such features would be difficult to detect in the inner regions due to the presence of a steeply rising rotation curve .
thus , their apparent absence in the inner region of the galaxy does not eliminate the possibility that they arise from the continuation of spiral density waves to the edge of the disk .
due the the steep rotational gradient , higher spatial resolution observations will be needed to investigate fully the pattern speeds in the outer gas disk .
the derived rotation curves were used to determine the dark matter content of ngc 1232 .
rotation curve decomposition requires knowledge of the gaseous and stellar mass distributions ( assumed to be disk  like ) and an assumed dark matter distribution .
the contribution of the gaseous and stellar components to the gravitational potential were calculated using the gipsy routine rotmod .
the input surface densities were estimated from a series of ellipses with position angles of 270  and inclinations of 30 .
the stellar distribution was taken as the observed azimuthal average to the extent of the optical emission ( figure [ fig : rot]e ) , followed by an extrapolated exponential disk . in fitting the rotation curve ,
the hi distribution was multiplied by 4/3 to account for the he content .
the gas content was held fixed during the fitting process .
note that for the observed neutral gas distribution of ngc 1232 , the net acceleration is outwards for the inner part of the gas disk ; following the usual convention , this outward acceleration is denoted by a negative velocity .
the stellar disk was fit with one free parameter , the mass  to  light ratio ( m / l@xmath10 ) .
note that since the available r  band image was obtained under non
photometric conditions , the observed surface brightness distribution was scaled to match the total blue luminosity ( de vaucouleurs 1991 ) . while b  band and r
band scale lengths are typically slightly different for spiral galaxies ( see , e.g. , de jong 1996 ) , this minor difference should be insignificant for the present mass models .
a simple spherical halo model was used for the dark matter contribution ; this model has only two free parameters , the core radius , r@xmath11 , and the central density , @xmath12 : @xmath13 where @xmath14 is the dark matter density at radius @xmath15 . in practice , the maximum velocity of the halo is fit rather than the central density : @xmath16,\ ] ] where v@xmath17 is the contribution of the halo to the gravitational potential at a radius @xmath15 and v@xmath18 is the asymptotic velocity , equal to @xmath19 .    to determine the parameters of the dark matter halo , the derived rotation curves were fit using a maximum  disk method .
in this method , the disk mass  to  light ratio is fit using the inner part of the rotation curve without overshooting the derived circular velocities .
initially , the maximum mass  to  light ratio ( in the b  band ) is found by assuming no dark matter is present .
for the warped disk model , the rotation curve derived from only the observable mass rapidly underestimated the observed rotation curve at large radii .
thus , the presence of a dark matter halo was assumed and the dark matter halo core radius , maximum velocity , and disk mass  to  light ratio were fit in an iterative manner .
in contrast , the falling rotation curve model was well fit by the luminous matter alone and did not require a dark matter halo ( i.e. , the falling rotation curve is nearly keplarian in nature ) .
the derived rotation curve decompositions are shown in figure [ fig : rotcur ] . in both cases , the derived disk mass  to  light ratios ( 3.0 for a warped disk , 3.7 for a falling rotation curve )
are reasonable for the observed colors ( see , e.g. , broeils 1992 ) .
one significant difference between the two models is the derived dynamical mass . for a spherical mass distribution ,
the dynamical mass within radius @xmath15 is : @xmath20 where @xmath21 is in km s@xmath2 and @xmath15 is in kpc .
however , if no dark matter halo is required , the system is disk  like and the dynamical mass is significantly less than that calculated assuming a spherical halo . for the falling rotation curve model , the dynamical mass is simply the sum of the gaseous and stellar mass , @xmath22 . for the warped disk model
, the dark matter halo dominates at large radii , so the dynamical mass is given by the above equation , @xmath23  at the last measured point of the rotation curve .
this latter value gives a dark  to  luminous ratio of 1.9 , quite typical of sc galaxies ( broeils 1992 ) .
as mentioned in section [ sec : obs ] , interactive gaussian fitting was used to determine the velocity dispersion throughout the galaxy .
both the second moment map and the velocity dispersion map from gaufit showed an apparent trend in the velocity dispersion as a function of radius , with the inner regions of the galaxy having extremely high velocity dispersions . despite the relatively low inclination of this galaxy , however , the velocity dispersion measurements are significantly contaminated by rotational broadening throughout the galaxy .
for instance , for beams between 40 and 120 arcsec from the center , the average rotational gradient is 0.3  arcsec@xmath2 ; for beams between 210 and 300 arcsec , the average rotational gradient is 0.2
arcsec@xmath2 . with a 40  beam
, this corresponds to a minimum of 8  of rotational broadening per beam . only in the outer regions of the galaxy
is the velocity gradient sufficiently shallow to permit detailed analysis of the line shape and width .
a line profile from gas in the outer disk is shown in figure [ fig : prof ] .
this profile was created by averaging line profiles from regions on the eastern and western sides of the galaxy with radii between 200  and 260 .
the profile is well fit by a gaussian with a velocity dispersion of 10.0  .    despite the large rotational gradient ,
it is possible to investigate radial trends in the velocity dispersion by deconvolving the rotational and turbulent motions .
the contribution from rotational broadening can be calculated either from a rotation curve model ( e.g. , table [ tab : rot ] ) , or from the observed rotational gradient , as measured in the position
velocity diagram .
the latter method is more direct , but is limited to regions along the major axis .
however , the major axis is also the region where the rotational broadening is minimized ; therefore , in the following analysis we limited our study to two independent slices @xmath024  from the major axis of the high resolution data cube .
the velocity dispersion was measured interactively for each beam averaged region and the rotational gradient for each region was obtained directly from the position
velocity diagram .
figure [ fig : disp ] shows the results of the deconvolution of thermal and rotational broadening ; each slice has been folded so that all radii are positive .
the observed line widths were corrected for rotational broadening by assuming that the rotational and thermal components add in quadrature .
as illustrated in figure [ fig : disp ] , the velocity dispersion is approximately constant as a function of radius , with an average value of 9.9 @xmath0 1.8 .
this is in remarkably good agreement with the velocity dispersion measured from the average of profiles in the outer regions of the galaxy .
similar results were also obtained from more detailed rotational models .
an estimate of the rotational broadening throughout the disk was obtained by creating a model system with the gipsy task galmod ; the system was assumed to have a negligible intrinsic velocity dispersion and the observed gas distribution and the derived ( no warp ) rotation curve were given as the input parameters .
the derived velocity dispersion of this model system is thus indicative of the rotational broadening at any given position in the galaxy .
the intrinsic velocity dispersion was recovered by subtracting , in quadrature , the model dispersion from the observed velocity dispersion of the high resolution data cube ( pixel  by  pixel ) .
similar to the above method , no radial trend was seen in the final velocity dispersion map .
the average velocity dispersion was measured to be 10.3 @xmath0 1.7  throughout the gas disk using this technique .
the observed velocity dispersion is much higher than one might expect for thermal motions . for thermal
broadening at 21 cm , the velocity dispersion is related to the thermal temperature by : @xmath24 thus , a velocity dispersion of 9.9  corresponds to a doppler temperature of @xmath2511900 k , which is significantly higher than the expected temperature of the warm neutral medium ( c.f .
kulkarni & heiles 1988 ) .
however , similarly high velocity dispersions have been observed in several other galaxies ( e.g. , ngc 3938 , van der kruit & shostak 1982 ; ngc 1058 , dickey 1990 ; ngc 6946 , boulanger & viallefond 1992 ; m101 , kamphuis 1993 ) . since neither the previous nor the present observations have adequate spatial resolution to resolve individual hi clouds , the observed velocity dispersion is the ensemble average of several clouds per beam .
in fact , the mean motions of the gas clouds , rather than the temperature of the gas within each cloud , could be the dominant motion resulting in the broad velocity dispersion .    for a galaxy with a high star formation rate ,
turbulent motion is obviously a significant contributor to the velocity dispersion in the inner regions of the galaxy , where supernovae and outflows provide large amounts of kinetic energy .
one difficulty , however , is to explain how such high levels of turbulence are sustained throughout the gas disk , particularly in regions well beyond the optical edge or in systems with modest star formation rates ( such as ngc 1232 ) .
recently , sellwood & balbus ( 1999 ) proposed that mhd instabilities in the differentially rotating gas disk could support a moderate amount of turbulence at all radii via dynamical heating of the disk .
their model requires a modest magnetic field ( @xmath25 3 @xmath26 g ) to reproduce a velocity dispersion of 6  in the outer gas disks of spiral galaxies .
while it is still unclear whether this mechanism is the dominant source of turbulence at large radii , it is an intriguing possibility .
several studies report a radial decrease in the velocity dispersion , with the outer regions of a galaxy having narrower and more perfectly gaussian line shapes ( e.g. , ngc 6946 , boulanger & viallefond 1992 ; m101 , kamphuis 1993 ) .
as discussed previously , ngc 1232 also has an observable radial decrease in the velocity dispersion , all of which can be accounted for by rotational broadening in the presence of a steeply rising rotation curve . while most studies acknowledge the complications of rotational broadening , it can be quite difficult to account fully for its effects on the data . for instance , when convolved with a decreasing gas density in the center of the galaxy , beam smearing can result in highly asymmetric , non  gaussian spectra .
a velocity dispersion measured either by moment analysis or gaussian fitting will not accurately reflect the width of the line in such spectra .
thus , it is extremely difficult to study the profile shape and radial trends in the velocity dispersion unless the observations are obtained with high spectral _ and _ spatial resolution .
we have investigated the gas distribution and kinematics of the nearly face  on spiral galaxy ngc 1232 . our results and conclusions are summarized below .
\(1 ) as is typical of spiral galaxies , the neutral gas distribution extends a factor of 1.6 beyond the optical radius .
the gas distribution is slightly asymmetric , with a higher gas column density associated with the prominent spiral arms on the eastern side .
\(2 ) the observed rotation curve has a nearly keplarian decline beyond the edge of the optical disk .
small variations in the inclination angle ( which is poorly constrained due to the low inclination of this galaxy ) can recover a flat rotation curve for the system .
assuming that the system is indeed warped , the dark  to  luminous ratio is 1.9 .
\(3 ) the presence of a steeply rising rotation curve complicates measurements of the velocity dispersion as a function of radius .
nonetheless , after correcting for rotational broadening , the velocity dispersion appears to be constant throughout the galaxy , with an average value of 9.9 @xmath0 1.8 .    as mentioned in the introduction ,
one of the reasons for looking at the gas distribution of ngc 1232 was to see if there were unusual gas kinematics associated with the outermost hii regions .
there are no obviously disturbed kinematic features in the present data set , but such features may have been washed out by the poor spatial resolution of the observations . nonetheless , to first order , the spectra are all well fit by a single gaussian ( with the exception of the central region where the rotational broadening is severe ) and are consistent with a slightly warped , differentially rotating , gas disk extending well beyond the optical radius .
thus , it is unlikely that large scale kinematic peculiarities are the cause of the abnormal abundances measured in the outer hii region .
analysis of this data set emphasizes the need for good quality high resolution ( both spectral _ and _ spatial ) observations of face  on galaxies .
the few canonical systems ( e.g. , ngc 1058 and ngc 3938 ) have been studied extensively with both wsrt and the vla , but more systems are needed to address possible dependences of the velocity dispersion on galaxy morphology and star formation rates . while the requirement for low inclination angles significantly restricts the potential target list , corrections for rotational broadening are less important with high spatial resolution observations .
in fact , these studies will become more feasible within the next few years ; with a new correlator , the upgraded vla will permit high velocity resolution over a wide bandpass , thereby improving the quality of similar observations of moderately inclined spiral galaxies .
we thank the nrao summer students of 1992 who elected to observe ngc 1232 as part of their summer student project : chris de pree , robbie dohm  palmer , eric olson , ben oppenheimer , carlos rabaca , james sweeney , and eric wagoner .
we thank michael rupen , katrina koski , andreea petric , and dave westpfahl for many informative conversations about gas in face  on galaxies .
martha haynes kindly provided access to the palomar 5 m telescope for the redshift determination of n1232a .
we thank the anonymous referee for thoughtful comments on rotational broadening .
robbie dohm  palmer and katrina koski are gratefully acknowledged for providing helpful comments on an early draft of the paper .
jb acknowledges funding through the research experience for undergraduates program run by the nrao and funded by the national science foundation .
as noted in section [ sec : obs ] , optical spectra of ngc 1232a were obtained with the palomar 5 m telescope during the night of 1999 january 23 .
the spectra were of sufficient quality that further analysis ( beyond the redshift determination ) was warranted . here
, we describe the observations in more detail and present the results of oxygen and nitrogen abundance calculations for the two hii regions observed .    during the observations , a long slit ( 2 ) with a 2  aperture
was centered on ngc 1232a for one 1200 sec exposure .
the slit was positioned at an angle of 0  ( the parallactic angle at the time of the observation ) and passed through both the center of the galaxy and an hii region 19  to the south ( figure [ fig : n1232a ] ) .
a 5500   dichroic was used to split the light to the two sides ( blue and red ) , providing complete spectral coverage from 36007600  .
the blue spectra were acquired with the 600 lines / mm diffraction grating ( blazed at 4000  ) .
the red spectra were acquired with the 316 lines / mm diffraction grating ( blazed at 7500  ) .
thinned 1024@xmath271024 tek ccds , with read noises of 8.6 e@xmath28 ( blue ) and 7.5 e@xmath28 ( red ) , were used on the two sides of the spectrograph .
both ccds had a gain of 2 .
the effective spectral resolution of the blue camera was 5.0   ( 1.72  / pix ) ; the effective spectral resolution of the red camera was 7.9   ( 2.47  / pix ) .
the spatial scale of the long slit was 0.62 /pix on the blue and 0.48 /pix on the red side .    while the night was non
photometric , relative flux calibration was obtained by observations of standard stars from the list of oke ( 1990 ) .
wavelength calibration was obtained by observations of arc lamps taken before and after the galaxy observations .
a hollow cathode ( fe and ar ) lamp was used to calibrate the blue spectra ; a combination of he , ne , and ar lamps were used to calibrate the red spectra .    the spectra were reduced and analyzed with the iraf package .
the spectral reduction included bias subtraction , scattered light corrections , and flat fielding with both twilight and dome flats .
the 2dimensional images were rectified based on the arc lamp observations and the trace of stars at different positions along the slit .
the sky background was removed from the 2-dimensional images by fitting a low order polynomial along each row of the spectra .
one dimensional spectra of the two hii regions were extracted from the rectified images using a 3.6  extraction region for the southern hii region and a 3.0  region for the central hii region ( figure [ fig : n1232a ] ) .
reddening corrected line strengths for both hii regions are listed in table [ tab : lines ] .
derivation of the oxygen abundance for both hii regions followed the methods described in van zee ( 1998b ) .
the temperature sensitive line [ oiii ] @xmath294363 was detected in the southern hii region and thus the electron temperature could be determined directly ; for this hii region , the oxygen abundance was calculated using the emissivity coefficients from a version of the fivel program ( de robertis 1987 ) assuming a low density ( 100 @xmath30 ) and a t@xmath31 of 11050 @xmath32 . for the central hii region ,
the oxygen abundance was estimated based on the strong line ratios ( r@xmath33 ( [ oii ] + [ oiii])/h@xmath34 ) using the model grid of mcgaugh ( 1991 ) .
an electron temperature of 6500 @xmath0 500 k reproduced this estimated oxygen abundance , and was thus adopted to derive the n / o ratio for the hii region .
the derived oxygen and nitrogen abundances for both hii regions are listed in table [ tab : lines ] .
not unexpectedly , the southern hii region is relatively low abundance while the central hii region is significantly higher abundance . with measurements of only two hii regions in hand ,
the derivation of an abundance gradient in this system is premature ; nonetheless , to first order , the abundance gradient is similar to those found in other spiral galaxies ( e.g. , vila  costas & edmunds 1992 ; zaritsky 1994 ; van zee 1998b ) .
becker , r. , mebold , u. , reif , k. , & van woerden , h. 1988 , , 203 , 21 boulanger , f. , & viallefond , f. 1992 , , 266 , 37 briggs , d. 1995 , phd thesis , new mexico tech brinks , e. , & bajaja , e. 1986 , , 169 , 14 broeils , a. h. 1992 , phd thesis , university of groningen broeils , a. h. , & rhee , m.h .
1997 , , 324 , 877 broeils , a. h. , & van woerden , h. 1994 , , 107 , 129 cayatte , v. , kotanyi , c. , balkowski , c. , & van gorkom , j. h. 1994 , , 107 , 1003 condon , j. j. 1987 , , 65 , 485 de jong , r. s. 1996 , , 313 , 377 de robertis , m. m. , dufour , r. j. , & hunt , r. w. 1987 , , 81 , 195 de vaucouleurs , g. , de vaucouleurs , a. , corwin , h. g. , buta , r. , paturel , g. , & fouqu , p. 1991
, third reference catalogue of bright galaxies ( springer , new york ) ( rc3 ) dickey , j. m. , hanson , m. m. , & helou , g. 1990 , , 352 , 522 fisher , j. r. & tully , r. b. 1981 , , 47 , 139 ferguson , a. m. n. , gallagher , j. s. , & wyse , r. f. g. 1998a , , 116 , 673 ferguson , a. m. n. , wyse , r. f. g. , gallagher , j. s. , & hunter , d. a. 1998b , , 506 , l19 haynes , m. p. , hogg , d. e. , maddalena , r. j. , roberts , m. s. , & van zee , l. 1998 , , 115 , 62 jore , k. p. , broeils , a. h. , & haynes , m. p. 1996
, , 112 , 438 kamphuis , j. 1993 , phd thesis , university of groningen kamphuis , j. , & briggs , f. 1992 , , 253 , 335 kulkarni , s. r. , & heiles , c. 1988 , in galactic and extragalactic radio astronomy , g.l .
verschuur and k.i .
kellermann , eds . , ( springer - verlag : new york ) mccall , m. l. , rybski , p. m. , & shields , g. a. 1985 , , 57 , 1 mcgaugh , s. s. 1991 , , 380 , 140 napier , p. j. , thompson , r. t. , & ekers , r. d. 1983 , proc .
ieee , 71 , 1295 oke , j. b. 1990 , , 99 , 1621 reif , k. , mebold , u. , goss , w. m. , van woerden , h. , & siegman , b. 1982 , , 50 , 451 roberts , m. s. , & haynes , m. p. 1994
, , 32 , 115 rownd , b. k. , dickey , j. m. , & helou , g. 1994 , , 108 , 1638 sellwood , j. a. , & balbus , s. a. 1999 , , 511 , 660 shostak , g. s. , & van der kruit , p. c. 1984 , , 132 , 20 soifer , b. t. , boehmer , l. , neugebauer , g. & sanders , d. b. 1989 , , 98 , 766 sparke , l. s. , & casertano , s. 1988 , , 234 , 873 staveley  smith , l. , & davies , r. d. 1988 , , 231 , 833 van  der  hulst , j. m. , terlouw , j. p. , begemen , k. , zwitser , w. , & roelfsema , p. r. 1992 , in astronomical data analysis software and systems , ed .
d. worall , c. biemesderfer , and j. barnes ( asp conference series ) , 25 , 131 van der kruit , p. c. , & shostak , g. s. 1982 , , 105 , 351 van der kruit , p. c. , & shostak , g. s. 1984 , , 134 , 258 van zee , l. , salzer , j. j. , & haynes , m. p. 1998a
, , 497 , l1 van zee , l. , salzer , j. j. , haynes , m. p. , odonoghue , a. a. , & balonek , t. j. 1998b , , 116 , 2805 vila
costas , m. b. , & edmunds , m. g. 1992 , , 259 , 121 zaritsky , d. , kennicutt , r. c. , & huchra , j. p. 1994
, , 420 , 87    ll distance ( mpc ) & 21.5kinematic center ( b1950 ) & 03 07 30.2 , 20 46 08 hi properties : integrated h i flux density ( jy km s@xmath2 ) & 122.4 @xmath0 12.2 heliocentric velocity ( km s@xmath2 ) & 1678.9 @xmath0 1.0 profile width ( 50% ) ( km s@xmath2 ) & 234 @xmath0 2 .
profile width ( 20% ) ( km s@xmath2 ) & 256 @xmath0 3 .
velocity dispersion ( km s@xmath2 ) & 9.9 @xmath0 1.8 size at 10@xmath8 atoms @xmath9 ( arcmin @xmath27 arcmin ) & 11.6 @xmath27 9.7 peak h i surface density ( atoms @xmath9 ) & 1.4 @xmath27 10@xmath35 optical properties : apparent blue magnitude & 10.52 @xmath0 0.14 ( u  b)@xmath1 & 0.02 @xmath0 0.10 ( b  v)@xmath1 & 0.60 @xmath0 0.04 r@xmath36 ( kpc ) & 23.2 r@xmath37 ( kpc ) & 6.3 derived properties : hydrogen mass , @xmath38 ( 10@xmath39 ) & 1.33 optical luminosity , @xmath40 ( 10@xmath41 ) & 4.46 fir luminosity , @xmath42 ( 10@xmath41 ) & 0.875 dynamical mass , @xmath43 ( 10@xmath39 ) & 43.6 dynamical mass  no warp , @xmath43 ( 10@xmath39 ) & 18.3 @xmath44 ( @xmath45 ) & 0.3 @xmath46 & 0.03 @xmath47 & 1.9      rcccc 10 . & 1.27 & 30 & 270 & 70.7 @xmath0 20.7 40 . &
2.20 & 30 & 270 & 184.9 @xmath0 1.4 70 . & 4.23 & 30 & 270 & 209.9 @xmath0 1.6 100 . & 5.85 & 30 & 270 & 216.6 @xmath0 0.8 130 . & 6.34 & 30 & 270 & 221.8 @xmath0 0.7 160 . & 5.91 & 30 & 270 & 220.1 @xmath0 0.6 190 . & 5.14 & 29 & 270 & 221.0 @xmath0 0.6 220 . & 3.70 & 28 & 269 & 220.9 @xmath0 0.7 250 . & 2.18 & 26 & 268 & 219.4 @xmath0 1.3 280 . & 1.19 & 23 & 268 & 218.8 @xmath0 1.7 310 . &
0.72 & 21 & 268 & 215.8 @xmath0 1.6 340 . &
0.52 & 20 & 268 & 219.1 @xmath0 1.9 370 . & 0.41 & 20 & 269 & 217.9 @xmath0 2.3 400 . & 0.32 & 20 & 269 & 217.8 @xmath0 6.4    lccc & 3728 & 2.32 @xmath0 0.17 & 2.30 @xmath0 0.14 & 3869 & & 0.19 @xmath0 0.02 hi & 4340 & 0.47 @xmath0 0.05 & 0.47 @xmath0 0.03 & 4363 & & 0.037 @xmath0 0.016 hi & 4861 & 1.00 @xmath0 0.06 & 1.00 @xmath0 0.05 & 4959 & 0.61 @xmath0 0.05 & 1.34 @xmath0 0.06 & 5007 & 1.81 @xmath0 0.11 & 4.33 @xmath0 0.19 hei & 5876 & 0.09 @xmath0 0.01 & 0.10 @xmath0 0.01 & 6300 & 0.09 @xmath0 0.12 & 0.064 @xmath0 0.005 & 6364 & 0.03 @xmath0 0.01 & 0.020 @xmath0 0.003 & 6584 & 0.11 @xmath0 0.01 & 0.045 @xmath0 0.004 hi & 6563 & 2.86 @xmath0 0.22 & 2.84 @xmath0 0.18 & 6584 & 0.38 @xmath0 0.03 & 0.14 @xmath0 0.01 hei & 6678 & & 0.038 @xmath0 0.004 & 6716 & 0.53
@xmath0 0.04 & 0.26 @xmath0 0.02 & 6731 & 0.27 @xmath0 0.02 & 0.15 @xmath0 0.01 & 7136 & & 0.079 @xmath0 0.006 c@xmath48 & & 0.02 @xmath0 0.08 & 0.24 @xmath0 0.06 f(h@xmath34)@xmath49 & & 2.05 & 4.98 ew(h@xmath34 ) [  ] & & 10.8 & 95.3 log(r@xmath50 ) & & 0.676 @xmath0 0.019 & 0.902 @xmath0 0.013 t(o@xmath51 ) & & 6500 @xmath0 500 & 11050 @xmath32 12 + log(o / h ) & & 8.79 @xmath0 0.10 & 8.21 @xmath0 0.10 log(n / o ) & & 1.24 @xmath0 0.14 & 1.42 @xmath0 0.15


wireless powered communication network has arisen as a new system with stable and self - sustainable power supplies in shaping future - generation wireless communications @xcite .
the enabling technology , known as simultaneous wireless information and power transfer ( swipt ) , has particularly drawn an upsurge of interests owing to the far - field electromagnetic power carried by radio - frequency ( rf ) signals that affluently exist in wireless communications . with the transmit power , waveforms , and dimensions of resources , etc .
, being all fully controllable , swipt promises to prolong the lifetime of wireless devices while delivering the essential communication functionality , as will be important for low - power applications such as rf identification ( rfid ) and wireless sensor networks ( wsns ) ( see @xcite and the references therein ) .    on the other hand , privacy and authentication
have increasingly become major concerns for wireless communications and physical ( phy)-layer security has emerged as a new layer of defence to realize perfect secrecy transmission in addition to the costly upper - layer techniques such as cryptography . in this regard , relay - assisted secure transmission was proposed @xcite and phy - layer security enhancements by means of cooperative communications have since attracted much attention @xcite .
in particular , cooperative schemes can be mainly classified into three categories : _ decode - and - forward ( df ) _ , _ amplify - and - forward ( af ) _ , and _ cooperative jamming ( cj ) _
@xcite with cj being most relevant to phy - layer security .
specifically , coordinated cj refers to the scheme of generating a common jamming signal across all single - antenna relay helpers against eavesdropping @xcite , while uncoordinated cj considers that each relay helper emits independent artificial noise ( an ) to confound the eavesdroppers @xcite .
it is expected that in the scenarios where the direct link is broken between the transmitter ( tx ) and the legitimate receiver ( rx ) , some of the relays have to take on their conventional role of forwarding the information while others will perform cj @xcite .
a recent paradigm that generalizes all the above - mentioned cooperation strategies is _ cooperative beamforming ( cb ) mixed with cj _
@xcite , where the available power at each relay is split into two parts : one for forwarding the confidential message and the other for cj .
however , mixed cb - cj approaches may be prohibitive in applications with low power devices because idle relays with limited battery supplies would likely prefer saving power for their own traffic to assisting others communication . in light of this
, swipt provides the incentive for potential helpers to perform dedicated cb mixed with cj at no expense of its own power , but opportunistically earn harvested energy .
motivated by this , our work considers secrecy transmission from a tx to a legitimate rx with the aid of a set of single - antenna wireless energy harvesting ( weh)-enabled af - operated relays in the presence of multiple single - antenna eavesdroppers . as a matter of fact
, cooperative schemes that involve weh - enabled relays was recently investigated in @xcite .
we consider the use of the _ dynamic power splitting ( dps ) _ receiver architecture , initially proposed for swipt in @xcite , which divides the received power with an adjustable ratio for energy harvesting ( eh ) and information receiving ( ir ) .
weh - enabled relays using dps receivers have also been considered in @xcite and @xcite , without ( w / o ) and with secrecy consideration , respectively .
note that there is also interest in addressing the threat that weh receivers may attempt to intercept the confidential messages in swipt - enabled networks @xcite .
nevertheless , we will focus on exploiting the benefits of weh - enabled relays when they are trustful .
in particular , motivated by the strong interest in swipt and the vast degree - of - freedom ( dof ) achievable by cooperative relays , this paper aims to maximize the secrecy rate with the aid of weh - enabled af - operated relays , subject to the eh power constraints of individual relays by jointly optimizing the cb of the relays and the cj covariance matrix . in this paper
, we assume that there is no direct link between the source and destination nodes , and perfect global channel state information ( csi ) is available for the case of centralized optimization .
it is worth pointing out that although our setting may look similar to @xcite , their optimal cb - cj design is not applicable to ours due to the multiplicative nature in beamforming weights incurred by the power splitting ( ps ) ratios that intrinsically poses more intractability to our optimization problem .
further , our work also differs from @xcite where an efficient algorithm was proposed to maximize the secrecy rate for the optimization of the ps ratios and af relay beamforming .
the difference is twofold .
first , an was not considered in the second transmission phase in @xcite and in addition , their algorithm only converged to a local optimum , as opposed to our work that gives the global optimal solutions for cb .
the rest of the paper is organized as follows .
section  [ sec : system model ] describes two types of weh - enabled rx architecture for the af relays and defines the secrecy rate region of the relay wiretap channel .
section  [ sec : problem formulation ] then formulates the secrecy rate maximization problems that jointly optimize the an ( or cj ) and the af - relay cb for the weh - enabled relays operating with the two types of rx .
the problems are respectively solved by centralized schemes in section  [ sec : centralized ] and distributed approaches in section  [ sec : distributed ] .
section  [ sec : numerical results ] provides simulation results to evaluate the performance of the proposed schemes .
finally , section  [ sec : conclusion ] concludes the paper .    _
notations_we use the uppercase boldface letters for matrices and lowercase boldface letters for vectors .
the superscripts @xmath0 , @xmath1 , @xmath2 and @xmath3 represent , respectively , the transpose , conjugate , conjugate transpose operations on vectors or matrices , and the optimum .
in addition , @xmath4 stands for the trace of a square matrix .
moreover , @xmath5_{i , j}$ ] denotes the @xmath6th entry of a matrix , while @xmath7 and @xmath8 represent the euclidean norm and the entry - wise absolute value square of a vector , respectively .
also , @xmath9 denotes a diagonal matrix with its diagonal specified by the given vector and @xmath5_{i=1}^n$ ] represents an @xmath10 vector with each element indexed by @xmath11 .
furthermore , @xmath12 and @xmath13 stand for product and hadamard product , respectively .
@xmath14 denotes the field of complex ( real ) matrices with dimension @xmath15 and @xmath16}$ ] indicates the expectation operation .
finally , @xmath17 is short for @xmath18 .
in this paper , we consider secrecy transmission in a swipt - enabled wsn , where a tx ( alice ) wants to establish confidential communication with the legitimate rx ( bob ) with no direct link but with the aid of @xmath19 weh - enabled sensors operating as af relays , denoted by @xmath20 , in the presence of multiple eavesdroppers ( eves ) , denoted by @xmath21 , all equipped with single antenna .
we assume that there is no direct link from the tx to any of the eves , due to , for instance , severe path loss or shadowing , as illustrated in fig .
[ fig : system model ] .
we consider a two - hop relaying protocol based on two equal time slots and the duration of one transmit - slot is normalized to be one unit so that the terms `` energy '' and `` power '' are interchangeable with respect to ( w.r.t . )
one transmit - slot .    at the receiver of each af relay , we introduce two types of weh - enabled receiver architecture , namely , _ static power splitting ( sps ) _ ( fig .  [ fig : subfig : sps rx ] ) and _ dps _ ( fig .
[ fig : subfig : dps rx ] ) , both of which allow the relay to harvest energy and receive information from the same received signal .
specifically , the receiver first splits a portion of @xmath22 , of the received power for eh and the rest @xmath23 for ir , @xmath24 .
the @xmath22 portion of harvested power is further divided into two streams with @xmath25 used for generating the an to confound eves and @xmath26 used for amplifying the received signal , where @xmath27 is the @xmath11th element of the received signal @xmath28 , and @xmath29 denotes the eh efficiency .
note that dps with adjustable @xmath22 s is presently the most general receiver operation because practical circuits can not directly decode the information from the stream used for eh @xcite and sps is just a special case of dps with @xmath30 , @xmath24 , fixed for the whole transmission duration .
however , sps , advocated for its ease of implementation , is introduced separately in the sequel for its simplified relay beamforming design . in the first transmit - slot , the received signal at each individual relay can be expressed as @xmath31 where the transmit signal @xmath32 is a circularly symmetric complex gaussian ( cscg ) random variable with zero mean and unit variance , denoted by @xmath33 , @xmath34 denotes the complex channel from the tx to the @xmath11th relay , @xmath35 is the transmit power at the tx , and @xmath36 is the additive white gaussian noise ( awgn ) introduced by the receiving antenna of the @xmath11th relay , denoted by @xmath37 . as such , the linearly amplified baseband equivalent signal at the output of the @xmath11th relay
is given by @xmath38 where @xmath39 denotes the complex af coefficient , and @xmath40 denotes the noise due to signal conversion from the rf band to baseband , denoted by @xmath41 .
since @xmath42 is constrained by the portion of the harvested power for forwarding , i.e. , @xmath26 , @xmath39 is accordingly given by @xmath43 where @xmath44 denotes the phase of the af coefficient for the @xmath11th relay .
next , we introduce the cj scheme .
denote the cj signal generated from @xmath19 relays by @xmath45^t$ ] and define its covariance matrix as @xmath46 $ ] .
then the coordinated cj transmission can be uniquely determined by the truncated eigenvalue decomposition ( evd ) of @xmath47 given by @xmath48 , where @xmath49)$ ] is a diagonal matrix with @xmath50 s denoting all the positive eigenvalues of @xmath47 and @xmath51 is the precoding matrix satisfying @xmath52 .
note that @xmath53 denotes the rank of @xmath54 which will be designed later . as a result
, the cj signal can be expressed as @xmath55 where @xmath56 s are drawn from the columns of @xmath57 , and @xmath58 s are independent and identically distributed ( i.i.d . ) complex gaussian variables denoted by @xmath59 , which is known as the optimal distribution for an @xcite . on the other hand , @xmath60 , @xmath24 ,
denotes the power constraint for jamming at the @xmath11th relay , which implies that @xmath61 where @xmath62 is a diagonal matrix with its diagonal @xmath63 ( a unit vector with the @xmath11th entry equal to @xmath64 and the rest equal to @xmath65 ) .
note that the cj scheme proposed above is of the most general form .
for the special case when @xmath66 , i.e. , @xmath67 , each relay transmits a common jamming signal @xmath68 with their respective weight drawn from @xmath69 @xcite .
this case is desirable in practice since it has the lowest complexity for implementation . in summary ,
the transmitted signal at the @xmath11th relay is given by @xmath70    according to together with , , and , the transmit signal from all relays can be expressed in vector form as @xmath71 where @xmath72 and @xmath73 are , respectively , diagonal matrices with their diagonals composed of @xmath74 @xmath75 and @xmath76 .
in addition , @xmath77_{i=1}^n$ ] , @xmath78_{i=1}^n$ ] , and @xmath79_{i=1}^n$ ] .
@xmath80    in the second transmit - slot , the received signal at the desired receiver , i.e. , bob , is given by @xmath81 where @xmath82_{i=1}^n$ ] comprises complex channels from the @xmath11th relay to the rx and @xmath83 is the corresponding receiving awgn . by substituting into
, @xmath84 can be expressed as @xmath85 the received signal at the @xmath86th eve , @xmath87 , is given by @xmath88 where @xmath89_{i=1}^n$ ] denotes the complex channels from the relays to the @xmath86th eve and @xmath90 is the awgn at the @xmath86th eavesdropper .
the mutual information for the rx ( bob ) is given by @xmath91 , and that for the @xmath86th eve is @xmath92 , @xmath93 , where @xmath94 and @xmath95 , which denote their respective signal - to - interference - plus - noise ratios ( sinrs ) , are given at the top of next page .
next , we define the secrecy rate region that consists of all the achievable secrecy rate for the relay wiretap channel given transmit power @xmath35 , denoted by @xmath96 , which is given by @xcite @xmath97
in this section , we consider the secrecy rate maximization problem by jointly optimizing the an beams , relay beam and their power allocations for weh - enabled af relays operating with sps , i.e. , @xmath30 , @xmath24 , is fixed .    by replacing @xmath39 with , @xmath98 in at the next page becomes @xmath99 and @xmath100_i\triangleq h_{sr_i}h_{r_id}\sqrt{\tfrac{\eta\bar\alpha_i(1-\bar\alpha_i)|h_{sr_i}|^2p_s}{(1-\bar\alpha_i)(|h_{sr_i}|^2p_s + \sigma_{n_a}^2)+\sigma_{n_c}^2}}.\ ] ] in addition , @xmath101 and @xmath102 in can also be combined as @xmath103_{i , i},\ ] ] where @xmath104_{i , i}=\frac{\eta\bar\alpha_ip_s|h_{sr_i}|^2|h_{r_id}|^2((1-\bar\alpha_i)\sigma_{n_a}^2+\sigma_{n_c}^2 ) } { ( 1-\bar\alpha_i)(|h_{sr_i}|^2p_s+\sigma_{n_a}^2)+\sigma_{n_c}^2}.\end{aligned}\ ] ] as a result , @xmath105 can be rewritten as @xmath106 where @xmath107_{i=1}^n$ ] .
similarly by letting @xmath108_i\triangleq h_{sr_i}h_{r_ie , k}\sqrt{\frac{\eta\bar\alpha_i(1-\bar\alpha_i)|h_{sr_i}|^2p_s}{(1-\bar\alpha_i)(|h_{sr_i}|^2p_s + \sigma_{n_a}^2)+\sigma_{n_c}^2}}\ ] ] and @xmath109_{i , i}\triangleq\frac{\eta\bar\alpha_ip_s|h_{sr_i}|^2|h_{r_ie , k}|^2((1-\bar\alpha_i)\sigma_{n_a}^2+\sigma_{n_c}^2 ) } { ( 1-\bar\alpha_i)(|h_{sr_i}|^2p_s+\sigma_{n_a}^2)+\sigma_{n_c}^2},\ ] ] we have @xmath110    by some simple manipulation , is reformulated as a per - relay jamming power constraint given by @xmath111    now , the secrecy rate maximization problem w.r.t .
@xmath112 s , @xmath44 s and @xmath47 for sps - based relays can be formulated as @xmath113    @xmath114      here ,
we consider the secrecy rate maximization problem for weh - enabled af relays with adjustable ps ratios @xmath115 by jointly optimizing the an beams , relay beam , weh ps ratios @xmath115 , and an ps ratios @xmath116 .
first , consider the following variable transformation : @xmath117 using this , @xmath118 can then be expressed as @xmath119 , where @xmath120_{i=1}^n$ ] and @xmath121_{i=1}^n$ ] .
moreover , @xmath122 and @xmath123 can be simplified as @xmath124 and @xmath125 , respectively , where @xmath126_{i=1}^n$ ] with @xmath127 , @xmath24 , and @xmath128_{i=1}^n$ ] .
similarly , we have @xmath129 , @xmath130 , and @xmath131 , where @xmath132_{i=1}^n$ ] , @xmath133 .
then we apply the above transformation to @xmath134 ( c.f .  ) and @xmath135 ( c.f .  ) , @xmath93 , to get and ( see next page ) .
now , we can recast the constraints w.r.t .
@xmath47 , @xmath22 s , and @xmath112 s to those w.r.t .
the transformed variables @xmath136 s and @xmath137 s . in accordance with , the optimization variables , @xmath22 s and @xmath112 s , can be alternatively given by @xmath138 where @xmath139 .
replacing @xmath22 s and @xmath112 s with , is reformulated as @xmath140 on the other hand , since @xmath141 and @xmath142 , @xmath24 , after some simple manipulation , it follows from that @xmath143    as such , the secrecy rate maximization problem for dps - based relays becomes @xmath144
in this section , we resort to centralized approaches to solve problem @xmath145 and @xmath146 , respectively , assuming that there is a central optimizer that is able to collect all csis including @xmath147 , @xmath148 and @xmath149 , perform the optimization , and broadcast to relays their individual optimized parameters .      to start with , we recast @xmath145 into a two - stage problem by introducing a slack variable @xmath150 .
first of all , we solve the epigraph reformulation of @xmath145 with a fixed @xmath151 $ ] as @xmath152    defining @xmath153 as the optimum value of @xmath154 and denoting @xmath155 , the objective function of @xmath145 is given by @xmath156 where @xmath157 in the objective function has been omitted and we claim a zero secrecy rate if admits a negative value . as a result , @xmath145 can be equivalently given by @xmath158 note that this single - variable optimization problem allows for simple one - dimension search over @xmath159 $ ] , assuming that @xmath160 is attainable given any @xmath150 in this region .
as the physical meaning of @xmath161 in @xmath154 can be interpreted as the maximum permitted sinr for the best eavesdropper s channel , feasibility for a non - zero secrecy rate implies that @xmath162 where cauchy - schwarz inequality has been applied in @xmath163 and @xmath164 follows from @xmath165 , @xmath166 .
the above epigraph reformulation of non - convex problems like @xmath145 has been widely employed in the literature @xcite , and @xmath167 admits the same optimal value as @xmath145 while @xmath154 with the optimal @xmath150 provides the corresponding optimal solution to @xmath145 .
we summarize the steps for solving @xmath145 here : given any @xmath159 $ ] , solve @xmath154 to obtain @xmath160 ; solve @xmath167 via a one - dimensional search over @xmath150 . before developing solutions to @xmath154
, we have the lemma below .
@xmath160 is a concave function of @xmath150 .
[ lemma : concavity of h_1(tau ) ]    see appendix  [ appendix : proof of concavity of h_1(tau ) ]
.    using lemma [ lemma : concavity of h_1(tau ) ] , it is easy to verify that @xmath168 is also a concave function of @xmath150 according to the composition rule @xcite , which allows for a more effective search for the optimum @xmath150 , e.g. , bi - section method , than the exhaustive search used in @xcite .
moreover , although @xmath160 is not differentiable w.r.t .
@xmath150 , the bi - section method can still be implemented , the algorithm involving which is similarly applied in solving @xmath146 and thus will be given later in section  [ sec : proposed solutions for dps ] .    in the sequel , we focus on solving @xmath154 . by introducing @xmath169 and ignoring the rank - one constraint on @xmath170 , @xmath171 can be alternatively solved by @xmath172 note that the objective function has been multiplied by @xmath150 compared with that of @xmath154 for ready computation of @xmath160 .
although @xmath173 is made easier to solve than @xmath154 by rank relaxation , it is still a quasi - convex problem considering the linear fractional form of the objective function and constraints , for which charnes - cooper transformation @xcite will be applied for equivalent convex reformulation .
specifically , by substituting @xmath174 and @xmath175 into @xmath173 , it follows that @xmath176 problem @xmath177 can now be optimally and efficiently solved using interior - point based methods by some off - the - shelf convex optimization toolboxes , e.g. , cvx @xcite .
we have the following results :    1 .
the optimal solution to @xmath177 satisfies @xmath178 ; 2 .
@xmath179 , where @xmath180 is given by @xmath181 in which @xmath182 is given in appendix  [ appendix : proof of rank - one solution with cj for sps ] ; 3 .
@xmath183 .
[ prop : rank - one solution with cj for sps ]    see appendix  [ appendix : proof of rank - one solution with cj for sps ] .
proposition  [ prop : rank - one solution with cj for sps ] implies that the rank - one relaxation of @xmath173 from @xmath154 is tight for an arbitrary given @xmath150 .
the @xmath184 s and @xmath185 s can thus be retrieved from the magnitude and angle of @xmath186 , respectively , by applying evd to @xmath187 .
similar to section  [ sec : optimal solutions for sps ] , in this section , we aim at solving the two - stage reformulation of @xmath146 by introducing a slack variable @xmath188 $ ] .
first , for a given @xmath150 , we solve @xmath189 next , denoting @xmath190 by @xmath191 , where @xmath192 is the optimum value for problem @xmath193 , we solve the following problem that attains the same optimum value as @xmath146 : @xmath194 where @xmath195 is similarly derived as @xmath196 so that we directly arrive at @xmath197 , denoted by @xmath195 .
we claim that @xmath198 can be solved by bi - section for @xmath150 over the interval @xmath199 $ ] assuming that @xmath191 is valid for any given @xmath150 ( otherwise a zero secrecy rate , i.e. , @xmath200 , is returned . ) , since @xmath191 has the following property .
@xmath191 is a concave function of @xmath150 .
[ lemma : concavity of h_2(tau ) ]    the proof is similar to that for lemma  [ lemma : concavity of h_1(tau ) ] , and thus is omitted .
it is also seen that how to attain @xmath191 forms the main thrust for solving @xmath146 . however , the constraints in , and are not convex w.r.t .
@xmath136 and/or @xmath137 , @xmath24 , due to their high orders and multiplicative structure .
@xmath193 thus turns out to be very hard to solve in general . to cope with these non - convex constraints , we introduce the following lemma .    the restricted hyperbolic constraints which have the form @xmath201 , where @xmath202 , @xmath203 , are equivalent to rotated second - order cone ( soc ) constraints as follows .
@xmath204[lemma : hyperbolic to soc ]    for convenience , denoting @xmath205 , @xmath206 , @xmath207 by @xmath208 , @xmath209 , and @xmath210 , respectively , @xmath24 , can be rewritten as @xmath211 according to and , it is easily verified that @xmath212 . hence , is eligible for lemma  [ lemma : hyperbolic to soc ] , which is reformulated into the soc constraint : @xmath213 similarly , can be simplified as @xmath214 , and after some manipulation , it is recast into a constraint of the restricted hyperbolic form as @xmath215 is thus , in line with lemma  [ lemma : hyperbolic to soc ] , equivalent to an soc constraint given by @xmath216 at last , is a linear constraint w.r.t .
@xmath208 and @xmath209 given by @xmath217    @xmath218    note that  have so far been equivalently transformed into the soc constraints , the linear constraints , as well as , the latter two of which are jointly convex w.r.t .
@xmath208 and @xmath209 , @xmath24 . however , is still not convex w.r.t .
@xmath210 , @xmath24 , yet . to circumvent this , in the sequel we propose to solve problem @xmath146 by alternating optimization .
the upshot of the algorithm is that first we fix @xmath47 by @xmath219 and thus @xmath210 by @xmath220 , @xmath24 , and solve problem @xmath221 ( @xmath193,@xmath198 ) with fixed @xmath47 as @xmath221 ( @xmath222,@xmath223 ) in the sequel . ] to find the optimal @xmath224 , @xmath225 and @xmath226 via @xmath222 and @xmath223 ; then with @xmath227 , @xmath24 , we devise the optimal solution derived in section  [ sec : optimal solutions for sps ] to obtain the optimal cj covariance , viz @xmath228 , and thus @xmath229 , @xmath24 ; finally , by updating @xmath230 and @xmath231 , @xmath24 , problems @xmath221 and @xmath145 are iteratively solved until they converge .
the remaining challenges lie in solving problem @xmath222 now that , and are all made convex w.r.t .
their variables @xmath208 , @xmath209 , @xmath24 .
similar to that for @xmath154 , we introduce @xmath232 and @xmath233 and exempt problem @xmath222 from @xmath234 and @xmath235 as follows : @xmath236 @xmath191 ( c.f .  ) and are given at the top of next page . recalling the procedure to deal with @xmath173
, we now apply charnes - cooper transformation to convert @xmath237 into a convex problem , denoted by @xmath238 , by replacing @xmath239 and @xmath240 with @xmath241 and @xmath242 , respectively .
the solution for @xmath238 is tight and characterized by the following proposition .
we have the following results :    1 .
the optimal solution to @xmath238 satisfies @xmath243 such that @xmath244 ; 2 .
@xmath245 is given by @xmath246 where @xmath247 , @xmath248 , and @xmath249 is given in appendix  [ appendix : proof of rank - one solution with cj for dps ] ; 3 .
@xmath250 , of which the diagonal entries compose a vector denoted by @xmath251 , can be reconstructed by @xmath252 with rank one .
[ prop : rank - one solution with cj for dps ]    see appendix  [ appendix : proof of rank - one solution with cj for dps ] .
the @xmath253 s and @xmath112 s are thus attained according to via evd of @xmath254 and @xmath255 .
the proposed algorithm for solving @xmath146 is presented in table  [ table : algorithm i ] .    ' '' ''    @xmath228 ; @xmath256 that denotes the optimum value for @xmath145 given @xmath257 , @xmath24 @xmath258 , @xmath259 @xmath260 @xmath261 , @xmath262 , @xmath24 , and solve @xmath221 : @xmath263 , @xmath264 , @xmath265 , @xmath266 , @xmath267 @xmath268 , @xmath269 solve @xmath222 and @xmath191 @xmath270 @xmath271 , where @xmath272 and @xmath273 denotes an arbitrary small value .
@xmath274 @xmath275 @xmath276 , @xmath277 , and obtain @xmath278 according to @xmath279 , @xmath24 , and solve @xmath145 via @xmath154 and @xmath167 @xmath280 , @xmath281 , and obtain @xmath282 and @xmath283 according to @xmath284 , @xmath24 update @xmath285 according to @xmath278 , @xmath282 , @xmath283 , and @xmath228    ' '' ''
in this section , we investigate heuristic algorithms to solve problems @xmath145 and @xmath146 in a completely distributed fashion . note that different from the paradigm of _ distributed optimization _ that allows for certain amount of information exchange based on which iterative algorithms are developed to gradually improve the system performance , we herein assume that each individual relay can only make decision based on its local csis , namely , @xmath34 , @xmath286 , @xmath287 , @xmath24 , and there is no extra means of information acquisition for ease of implementation .
the purpose for such an algorithm is twofold : on one hand , we aim to answer the question that in the least favourable situation , namely , no coordination over the relays , how to improve the achievable secrecy rate of the system ? on the other hand , it provides a lower - bound for the centralized schemes proposed in section  [ sec : centralized ] , which sheds light upon the trade - off achievable between secrecy performance and complexity .
besides , we emphasize the jamming scheme that is different from the cj in the centralized schemes . unlike the cj signal coordinately transmitted by all relays , in the distributed implementation , each relay is only able to generate its an locally , i.e. , @xmath288^t$ ] , in which @xmath289 s are i.i.d .
an beams , denoted by @xmath290 .
this type of cj is known to be uncoordinated with the covariance matrix given by @xmath291)$ ] . in this section
, we assume that each relay consumes all of its remaining power from af for an , i.e. , @xmath292 , @xmath166 ( c.f .  ) .
hence , the an design solely depends on @xmath22 s and/or @xmath112 s .
first , we propose a heuristic scheme for the @xmath11th af relay to decide on @xmath112 , @xmath24 , which is given by @xmath293 where @xmath294 is a constant controlling the relay s level of jamming .
for example , a larger @xmath295 indicates that each relay prefers to splitting a larger portion of power for jamming and vice versa .
the intuition behind is that if the @xmath11th relay observes that @xmath296 , which means that a nonnegative secrecy rate is achievable even if there is only itself in the system , it will shut down the an ; otherwise , it will split up to @xmath295 portion of the harvested power for jamming . in an extreme case of @xmath297 , probably when an eve is located within the very proximity of this relay , it allocates the maximum permissible portion of power , i.e. , @xmath295 , for an .
next , since an individual relay can not evaluate the secrecy performance of the whole system , @xmath44 s are simply chosen to be the optimum for the multi - af relaying without security considerations , i.e. , @xmath298 , @xmath24 .      following the same designs for @xmath112 s and @xmath44 s in section  [ sec : distributed alg . for sps
] , the remaining task for weh - enabled relays operating with dps is to set proper values for @xmath22 s .
we choose @xmath22 s that maximize the `` hypothetical sinr '' .
this `` hypothetical '' sinr may not be the actual sinr for the destination , but just a criterion calculated based on the `` hypothetical '' received signal seen by the @xmath11th relay , given by @xmath299 the corresponding sinr is thus expressed as @xmath300 where @xmath301 .
consequently , the maximization of w.r.t .
@xmath22 , @xmath24 , is formulated as @xmath302    the optimal @xmath22 , @xmath24 , to @xmath303 is @xmath304[prop : optimal alpha to distr ]    it is easy to verify that problem @xmath303 is convex and the minimum solution of its objective function derived from the first - order derivative happens to fall within the feasible region of @xmath22 , which is seen in .    with @xmath112 s , @xmath44 s and @xmath22 s set
, each af relay is then able to decide its relay weight and an transmission .
in this section we compare our proposed schemes operating with sps or dps with some benchmarks . in the centralized case ,
the optimal solution for sps in section  [ sec : optimal solutions for sps ] is denoted by _ cj - sps _ , while algorithm  [ table : algorithm i ] in section  [ sec : proposed solutions for dps ] is denoted by _ cj - dps_. the distributed schemes in section  [ sec : distributed alg . for sps ] and section  [ sec : distributed alg . for dps ] are referred to as _ distributed - sps _ and _ distributed - dps _ , respectively . to demonstrate the effectiveness of our an - aided secure multi - af relay beamforming algorithms
, we also provide three benchmark schemes : _ nocj - sps _ , _ nocj - dps _ and _ random ps_. for _ nocj - sps _ , we solve problem @xmath145 by replacing @xmath47 with @xmath305 . similarly , for _ nocj - dps _ , we initialize @xmath306 and quit the loop in algorithm  [ table : algorithm i ] after the very first time of solving problem @xmath221 . _ random ps _ , on the other hand , picks up i.i.d .
@xmath22 and @xmath112 uniformly generated over @xmath307 $ ] , respectively , and co - phases @xmath308 , @xmath24 .
consider that @xmath19 weh - enabled af relays and @xmath309 eavesdroppers are located within a circular area of radius @xmath310 .
specifically , we assume that their respective radius and radian are drawn from uniform distributions over the interval @xmath311 $ ] and @xmath312 , respectively .
we also assume that the channel models consist of both large - scale path loss and small - scale multi - path fading .
the unified path loss model is given by @xmath313 where @xmath314 , @xmath315 denotes the relevant distance , @xmath316 m is a reference distance , and @xmath317 is the path loss exponent set to be @xmath318 .
@xmath34 , @xmath286 , and @xmath319 , @xmath166 , @xmath133 , are generated from independent rayleigh fading with zero mean and variance specified by .
the simulation parameters are set as follows unless otherwise specified : the radius defining the range is @xmath320 m ; the transmit power at the source is @xmath321dbm ; the noise variances are set as @xmath322dbm , @xmath323dbm , @xmath324 , and @xmath325 , @xmath93 ; the eh efficiency is set to @xmath326 .
also , numerical results were averages over @xmath327 independent channel realizations .      here
, we evaluate the performance of the proposed centralized designs in section  [ sec : centralized ] .
the efficiency of the alternating optimization that iteratively attains numerical solution to @xmath146 is studied in fig .  [
fig : alternating sec vs iter ] , which shows the increment of the achievable secrecy rate after each round of the iteration .
the most rapid increase is observed after the first iteration , which illustrates that the optimization of the ps ratios , @xmath22 s , accounts for the main factor for the secrecy rate performance gains over a sps scheme that sets @xmath328 .
it is seen that the alternating algorithm converges within the relative tolerance set to be @xmath329 , after an average of @xmath330@xmath331 iterations for several channel realizations , which appears reasonable in terms of complexity .
dbm , @xmath332 , and @xmath333.,width=288 ]    fig .
[ fig : centralized compare vs n ] shows the achievable secrecy rate for the legitimate rx with the increase in the number of af relays by different schemes . as we can see , the secure multi - af relaying schemes assisted by the transmission of an outperforms that without an for both sps and dps . in addition , with the increase in @xmath19 , the role of cj gradually reduces for both schemes of sps and dps
this is because as @xmath19 gets larger , the optimal designs tend to suppress the interception at the most capable eavesdropper more effectively with @xmath19 dof , enforcing the numerator of @xmath334 to a relatively low level , which can also be observed from @xmath335 in fig .
[ fig : subfig : eve s rate vs n ] , and therefore the optimal amount of power allocated to an beams inclines to be little ; otherwise the jamming yielded will be detrimental to the reception at the legitimate rx . besides , given the same number of af relays , the secrecy performance gains brought by the proposed schemes with cj are more prominent in the presence of more eavesdroppers , since it is hard to reduce all the eavesdroppers channel capacity without resorting to cj properly .
db for @xmath332 and @xmath336 , respectively.,width=288 ]    fig .
[ fig : centralized compare sec vs eves ] shows the achievable secrecy rate for the legitimate rx versus the number of eavesdroppers by different schemes .
first , similar to the results shown in fig .
[ fig : centralized compare vs n ] , the proposed an - aided multi - af relaying designs operating with dps - enabled relays , viz _ cj - dps _ , perform best among all the schemes .
secondly , as @xmath309 goes up , the an - aided schemes , _ cj - dps _ and _ cj - sps _ , allow the secrecy rate to drop slowly , in other words , more robust against multiple eavesdroppers , while the secrecy rate of their nocj counterparts almost goes down linearly with @xmath309 .
furthermore , with @xmath309 increasing , for example , more than @xmath337 , the increase in the number of relays , from @xmath332 to @xmath338 , can not replace the role of cj as shown in fig .
[ fig : centralized compare vs n ] and [ fig : centralized compare sec vs eves ] , since in the presence of many eavesdroppers , more relays may also result in improved eavesdroppers decoding ability w / o the assistance of cj .
it is also noteworthy that with @xmath339 , there is little use of cj by the centralized schemes , which was also observed in @xcite .     for @xmath332 and @xmath336 , respectively.,width=288 ]
[ fig : centralized compare sec vs ps ] provides simulation results of different schemes by varying the source transmit power .
it is seen that with more power available at the source , the advantage of cj is more pronounced , since given other variables fixed , larger @xmath35 indicates larger feasible regions for @xmath145 and @xmath146 . furthermore , as similarly seen in fig .
[ fig : centralized compare vs n ] , with a mild number of eavesdroppers ( @xmath333 ) , subject to the same @xmath35 , a large number of cooperative relays enables more dof in designing the optimal @xmath22 s and @xmath44 s , which alleviates the dependence on an beams to combat eves .      here
, we study the performance of the distributed schemes , namely , _
distributed - sps _ and _ distributed - dps _ in section  [ sec : distributed ] .
as mentioned earlier , these heuristics are provided as benchmarks to demonstrate what can be done under the extreme `` no - coordination '' circumstance , in comparison with _
random ps_. note that any other distributed schemes with certain level of cooperation among relays are supposed to increase the secrecy performance up to the proposed centralized algorithms , namely , _ cj - dps _ and _ cj - sps _ , at the expense of extra computational complexity and system overhead .
db and @xmath333.,width=288 ]    fig .
[ fig : distributed compare sec vs n ] provides the results for the achievable secrecy rate of various schemes versus the number of relays . _ distributed - sps _ and _ distributed - dps _ ,
are observed to be outperformed by their centralized counterparts though , they are considerably superior to _
random ps_. it is also seen that the performance gap between the centralized and distributed approaches is enlarged as @xmath19 increases , which is expected , since larger @xmath19 yields more dof for cooperation that is exclusively beneficial for the centralized schemes .
furthermore , compared with the centralized schemes , the distributed ones are more vulnerable to the increase in the eavesdroppers number .
db and @xmath19=8.,width=288 ]    in fig .
[ fig : distributed compare sec vs eves ] , we investigate the relationship between the secrecy rate performance and the number of eavesdroppers by different methods . as can be observed , compared with the centralized schemes , the secrecy rates achieved by _ distributed - sps _ and _ distributed - dps _ both reduce more drastically with the increase in @xmath309 due to the lack of effective cooperation . also , the advantage of dps over sps for the distributed schemes is compromised since @xmath22 s are not jointly designed with other parameters . at last ,
a similar observation has been made as that for fig .
[ fig : distributed compare sec vs n ] , that is , larger @xmath19 yields more visible performance gap between the centralized and distributed approaches .     and @xmath333.,width=288 ]    in fig .
[ fig : distributed compare sec vs ps ] , we examine the effect of increasing the transmit power at the source on the secrecy performance of different schemes under the same settings as those in fig .
[ fig : centralized compare sec vs ps ] . among all the presented designs ,
_ cj - dps _ still achieves the best secrecy rate as observed in other examples . also , the fact that larger @xmath19 benefits more from cooperative designs is corroborated again due to the same reason as that for fig .
[ fig : distributed compare sec vs n ] and [ fig : distributed compare sec vs eves ] .
furthermore , the secrecy rate of _ distributed - sps _ or _ distributed - dps _ is quickly saturated when @xmath340db while that for their centralized counterparts still rises at fast speed .
this paper studied secure communications using multiple single - antenna weh - enabled af relays assisted by an via cj for a swipt network with multiple single - antenna eavesdroppers . using ps at the relays ,
the achievable secrecy rates for the relay wiretap channel were maximized by jointly optimizing the cb and the cj covariance matrix along with the ps ratios for relays operating with , respectively , sps and dps . for dps , reformulating the constraints into restricted hyperbolic forms enabled us to develop convex optimization - based solutions .
further , we proposed an information - exchange - free distributed algorithm that outperforms the random decision .
as the optimal solution to @xmath177 is proved to be optimal to problem @xmath154 , the optimum value for @xmath177 is @xmath341 .
hence , we identify the property of @xmath160 by investigating @xmath177 , with its lagrangian given by @xmath342 where @xmath343 denotes a tuple consisting of all the primal and dual variables .
specifically , @xmath344 , @xmath345 and @xmath346 are lagrangian multipliers associated with @xmath347 , @xmath348 and the first constraint of @xmath177 , respectively ; @xmath349 are the dual variables associated with the sinr constraint for the @xmath86th eve , respectively ; @xmath350_{i=1}^n)$ ] with each diagonal entry @xmath351 denoting the dual variable associated with the per - relay power constraint ; @xmath352 is the lagrangian multiplier associated with @xmath353 . in addition , @xmath354_{i=1}^n)$ ] .
the karush - kuhn - tucker ( kkt ) conditions for are partially listed as follows :      the associated dual problem is accordingly given by @xmath356 where @xmath344 , @xmath345 , and @xmath352 are given by
, respectively .
since it is easily verified that @xmath177 satisfies the slater condition , the strong duality holds @xcite .
this implies that the dual optimum value in @xmath357 is @xmath160 , which turns out to be a point - wise minimum of a family of affine functions and thus concave over @xmath159 $ ] @xcite .
pre- and post - multiplying @xmath359 with the left hand side ( lhs ) and right hand side ( rhs ) of , respectively , and substituting @xmath360 into , @xmath361 can be rewritten as @xmath362 introducing the notation of @xmath5_{\rm offd}$ ] to represent a square matrix with its diagonal entries removed , it follows from   that @xmath363_{\rm offd}=\boldsymbol 0.\ ] ] by subtracting from , @xmath344 can be rewritten as @xmath364_{\rm d}\notag\\ & + \left[\boldsymbol w_0^{\frac{1}{2}}\boldsymbol y_2\boldsymbol w_0^{\frac{1}{2}}\right]_{\rm d}+\sum_{k=1}^k\theta_k p_s\tilde{\boldsymbol h}_{se , k}^\dag\tilde{\boldsymbol h}_{se , k}^t,\label{eq : y1 intermediate2}\end{aligned}\ ] ] where @xmath5_{\rm d}$ ] denotes a square matrix with only the diagonal remained . observing that @xmath365_{\rm d}^{-1}\left[\boldsymbol w_0^{\frac{1}{2}}\boldsymbol h_{re , k}^\dag\boldsymbol h_{re , k}^t\boldsymbol w_0^{\frac{1}{2}}\right]_{\rm d}\nonumber\\ & = { \rm diag}\left([|h_{r_ie , k}|^2\left.\middle/\right.|h_{r_id}|^2]_{i=1}^n\right)\notag\\ & \equiv \boldsymbol r_{ed
, k},\label{eq : r_{ed , k}}\end{aligned}\ ] ] @xmath361 can be finally recast as @xmath366 where @xmath367_{\rm d}-\left(\left[\boldsymbol w_0^{\frac{1}{2}}\boldsymbol h_{rd}^\dag\boldsymbol h_{rd}^t\boldsymbol w_0^{\frac{1}{2}}\right]_{\rm d}-\boldsymbol d_{\hat{sd}}\right ) \left(\lambda\boldsymbol
i-\sum_{k=1}^k\theta_k\left(\frac{1}{\tau}-1\right)\boldsymbol r_{ed , k}\right).\end{aligned}\ ] ]    in the following , we show that @xmath368 is a positive definite matrix . note that since @xmath369 is a diagonal matrix , its definiteness is only determined by the signs of its diagonal entries , for which we commence with the discussion in three difference cases below .
1 .   * case i * : @xmath370 such that @xmath371_{i , i}<0 $ ] . since @xmath372_{\rm d}-\boldsymbol d_{\hat{sd}}\right]_{i , i}=$ ] + @xmath373
, it follows from that @xmath374_{i , i}>0 $ ] in this case .
* case ii * : @xmath370 such that @xmath371_{i , i}>0 $ ] .
we have @xmath375_{\rm d}\right]_{i , i}-[\boldsymbol w_0]_{i , i}\vert h_{r_id}\vert^2(\lambda^\ast-\sum_{k=1}^k\theta_k\left(\frac{1}{\tau}-1\right)[\boldsymbol r_{ed , k}]_{i , i})\ge 0 $ ] in accordance with , which implies that @xmath374_{i , i}=\left[\left[\boldsymbol w_0^{\frac{1}{2}}\boldsymbol y_2\boldsymbol w_0^{\frac{1}{2}}\right]_{\rm d}\right]_{i , i}-[\boldsymbol w_0]_{i , i}\vert h_{r_id}\vert^2(\lambda-\sum_{k=1}^k\theta_k\left(\frac{1}{\tau}-1\right)[\boldsymbol r_{ed , k}]_{i , i})+[\boldsymbol d_{\hat{sd}}]_{i , i}(\lambda-\sum_{k=1}^k\theta_k\left(\frac{1}{\tau}-1\right)[\boldsymbol r_{ed , k}]_{i , i})>0 $ ] ( c.f .  ) .
* case iii * : @xmath370 such that @xmath371_{i , i}=0 $ ] . in this case
, it follows that @xmath374_{i , i}= \left[\boldsymbol w_0^{\frac{1}{2}}\boldsymbol y_2\boldsymbol w_0^{\frac{1}{2}}\right]_{i , i}$ ] + @xmath376 .
it is noteworthy that the number of @xmath11 s such that @xmath371_{i , i}=0 $ ] can not exceed one .
this can be proved by contradiction as follows . if @xmath377 , @xmath378 , such that @xmath371_{i_1,i_1}=0 $ ] and @xmath371_{i_2,i_2}=0 $ ] , it implies that @xmath379_{i_1,i_1}=\sum_{k=1}^k\theta_k[\boldsymbol r_{ed , k}]_{i_2,i_2}$ ] , which contradicts to the fact that for any two independent continuously distributed random variables , the chance that they are equal is zero .    in summary , @xmath374_{i ,
i}\ge 0 $ ] , @xmath24 . if @xmath374_{i , i } > 0 $ ] , then it is obvious that @xmath380 . next , we show that it still holds true in the case that @xmath381 , such that @xmath374_{i^\prime , i^\prime}=0 $ ] , @xmath382 , by definition .
define the null - space of @xmath369 by @xmath383 and multiply @xmath384 and @xmath385 , @xmath386 , on the lhs and rhs of @xmath368 , respectively .
if @xmath387 , it is straightforward to obtain @xmath388 ; otherwise , it follows that @xmath389_{i^\prime}\vert^2>0 $ ] , as @xmath390_{i^\prime}\neq 0 $ ] in probability .
this completes the proof . as a result , @xmath344 in
is shown to always take on a special structure , that is , a full - rank matrix minus a rank - one matrix .
note that this observation plays a key role in proving the rank - one property of @xmath391 , which is also identified in ( * ? ? ?
* appendix c ) .      finally ,
multiplying both sides of by @xmath391 , as per , we obtain the following equation : @xmath393 which further implies that @xmath394 .
in addition , since the optimality of @xmath177 suggests that @xmath395 , @xmath396 is thus proved .
as @xmath391 can be decomposed as @xmath397 by evd , results in @xmath398 , which further implies that @xmath399 therefore , admits a unique solution @xmath182 up to a scaling factor , which is given by @xmath400 consequently , we have @xmath401 , where @xmath402 . on the other hand , by plugging @xmath401 into the equality constraint of @xmath177 , we have @xmath403 , which yields @xmath404    at last , we show @xmath405 of proposition  [ prop : rank - one solution with cj for sps ] . for the case of @xmath406 ,
it is obvious that @xmath407 . for the case of @xmath408 ,
only a sketch of the proof is provided here due to the length constraint . according to ,
first it is provable that @xmath409 is a full - rank matrix when @xmath177 obtains its optimum value ; next , observing that @xmath410 , it follows that @xmath411 as a result of @xmath412 ; then according to , @xmath413 is thus obtained .          for convenience of the proof ,
the optimum value for @xmath237 and @xmath415 are denoted by @xmath416 and @xmath417 , respectively . assuming that @xmath418 is the optimal solution to @xmath237 , it is easily verified to be feasible for @xmath415 as well , which implies that @xmath419 . on the other hand , if @xmath415 returns an optimal solution of @xmath420 , by defining @xmath421 and @xmath422 , @xmath24 , we show that @xmath423 is also feasible for @xmath237 as follows . as for , @xmath424 where @xmath163 is due to , and
@xmath164 comes from @xmath425 and @xmath426 .
similarly , @xmath427 + @xmath428 can be proved to satisfy as well .
in addition , @xmath429 , @xmath24 , i.e. , holds true .
these feasibility implies that @xmath430 . by combining the above two facts
, we have @xmath431 and complete the proof .    then
, we apply the charnes - cooper transformation again to @xmath415 , the result of which is denoted by @xmath432 , to study the property of @xmath433 and @xmath250 .
it is noteworthy that the charnes - cooper transformed constraint of admits the form given by @xmath434 , @xmath24 , where @xmath435 is the column vector inside @xmath7 of the lhs of and @xmath436 indicates the rhs expression . since it is easily checked that @xmath437 as a result of feasibility , we have @xmath438 , which implies that @xmath439 according to schur complement .
it thus follows that holds true , @xmath440 : @xmath434 . as such
, we can show that can be recast into a constraint without @xmath441 , i.e. , not related to @xmath442 , @xmath443 , following the same procedure as ( * ? ? ?
* appendix iii ) by exploiting ( * ? ? ?
* _ lemma 2 _ ) .
similarly , the charnes - cooper transformed constraint of can also be rewritten without @xmath442 , @xmath443 .
hence , the partial lagrangian for @xmath432 in terms of @xmath442 and @xmath443 can be expressed as ( [ eq : lagrangian of ( p2.1-sdp - eqv ) ] ) ( see top of next page ) ,      where @xmath445 denotes a tuple comprising all the associated primal and dual variables : @xmath344 , @xmath345 , and @xmath349 are lagrangian multipliers associated with @xmath442 , @xmath443 , and , @xmath93 , respectively ; @xmath346 is the dual variable associated with the only equality constraint ; @xmath446_{i=1}^n)$ ] and @xmath447_{i=1}^n)$ ] denote those associated with @xmath448 and @xmath449 , @xmath24 , respectively ; finally , the diagonal entry of @xmath450_{i=1}^n)$ ] denotes the dual variable associated with @xmath451 , @xmath24 .
the kkt conditions related to are accordingly given by        1 .
* case i * : @xmath456 , @xmath133 . in this case
, since @xmath457 ( c.f .  )
due to the strong duality , it is easily verified that @xmath458 and therefore @xmath459 .
* case ii * : @xmath460 such that @xmath461 .
according to , @xmath462 , which is a positive semidefinite diagonal matrix . according to the similar argument made in * case iii * of appendix  [ appendix : proof of rank - one solution with cj for sps ] , it is thus shown that @xmath249 has maximum one zero diagonal entry and the positive definiteness of @xmath463 can be proved by definition without difficulty .    as @xmath344 ( c.f .  )
again complies with the difference between a positive definite matrix and a rank one matrix , i.e. , @xmath464 , it turns out that @xmath465 according to .
then , following the same procedure as that in appendix  [ appendix : proof of rank - one solution with cj for sps ] , @xmath466 of proposition  [ prop : rank - one solution with cj for dps ] can be proved ( details omitted for brevity ) .
finally , it is verified that @xmath238 is related to @xmath250 merely with its diagonal entries , viz @xmath467_{i , i}$ ] , @xmath166 , ( c.f .  , ) .
furthermore , denoting @xmath468_{i , i}^{1/2}]_{i=1}^n$ ] by @xmath251 , it is easily checked that the diagonal entries remain the same after we replace @xmath250 with @xmath252 . combining the above two facts , we arrive at the conclusion that such modification does not affect its optimality while returning a rank - one @xmath250 for @xmath238 , which completes the proof for @xmath405 .
e. tekin , s. member , and a. yener , `` the general gaussian multiple - access and two - way wiretap channels : achievable rates and cooperative jamming , '' _ ieee trans .
inf . theory _
54 , no . 6 , pp
. 27352751 , jun .
2008 .
j. zhang , m. c. gursoy , `` relay beamforming strategies for physical - layer security , '' in _ proc .
ieee annual conference on information sciences and systems ( ciss10 ) _ , princeton , usa , mar .
2010 , pp.16 .
y. yang , q. li , w .- k .
ma , j. ge , and p. c. ching , `` cooperative secure beamforming for af relay networks with multiple eavesdroppers , '' _ ieee signal process .
1 , pp . 3538 , jan .
2013 .
h. xing , k .-
k wong , z. chu , and a. nallanathan , `` to harvest and jam : a paradigm of self - sustaining friendly jammers for secure af relaying , '' _ to appear in ieee trans . signal process .
_ , available on - line at arxiv:1502.07066v3 .
z. ding , k. k. leung , d. l. goeckel , and d. towsley , `` opportunistic relaying for secrecy communications : cooperative jamming vs. relay chatting , '' _ ieee trans .
10 , no . 6 , pp .
17251729 , jun . 2011 .
y. yang , q. li , w .- k .
ma , j. ge , and m. lin , `` optimal joint cooperative beamforming and artificial noise design for secrecy rate maximization in af relay networks , '' in _ proc .
ieee workshop on signal processing advances in wireless communications ( spawc ) _
, darmstadt , de , jun .
2013 , pp.360364 .
q. li , y. yang , w .- k .
ma , m. lin , j. ge , and j. lin ,  robust cooperative beamforming and artificial noise design for physical - layer secrecy in af multi - antenna multi - relay networks , " _ ieee trans . signal process .
1 , pp . 206220 , jan .
k. ishibashi , `` dynamic harvest - and - forward : new cooperative diversity with rf energy harvesting , '' in _ proc .
ieee international conference on wireless communications and signal processing ( wcsp ) _ , hefei , cn , oct .
2014 , pp . 15 .
a. nasir , x. zhou , s. durrani , and r. kennedy , `` relaying protocols for wireless energy harvesting and information processing , '' _ ieee trans.wireless commun .
36223636 , jul .
z. ding , s. perlaza , i. esnaola , and h. poor , `` power allocation strategies in energy harvesting wireless cooperative networks , '' _ ieee trans .
wireless commun .
_ , vol . 13 , no .
2 , pp . 846860 , feb . 2014 .
s. timotheou , i. krikidis , g. zheng , and b. ottersten , `` beamforming for miso interference channels with qos and rf energy transfer , '' _ ieee trans .
wireless commun .
_ , vol.13 , no.5 , pp .
26462658 , may 2014 .
m. zhao , x. wang , and s. feng , `` joint power splitting and secure beamforming design in the multiple non - regenerative wireless - powered relay networks , '' _ ieee commun .
9 , pp . 15401543 , sept . 2015 .
d. ng , e. lo , and r. schober , `` robust beamforming for secure communication in systems with wireless information and power transfer , '' _ ieee trans .
wireless commun .
_ , vol . 13 , no . 8 , pp
45994615 , aug . 2014 .
q. li , q. zhang , and j. qin , `` secure relay beamforming for simultaneous wireless information and power transfer in non - regenerative relay networks , '' _ ieee trans .
5 , pp . 24622467 , jun . 2014 .
q. li and w .- k .
ma , `` spatially selective artificial - noise aided transmit optimization for miso multi - eves secrecy rate maximization , '' _ ieee trans . signal process .
27042717 , may 2013 .