the presence of a tight correlation between the rest - frame optical color and luminosity of early - type galaxies , the so - called color - magnitude relation ( cmr ) or red - sequence , is well established in the literature .
the relation implies that the properties of the stellar populations in early - type galaxies are strongly dependent on the total stellar mass of the system . the tight correlation observed along this relationship ( e.g. * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ) requires a strong coherence between both the age and metallicity of stellar populations present in early - type galaxies .
measurements of the slope , zero - point , and dispersion around the color - magnitude relation probe the star formation histories of early - type galaxies and can provide physical insight into galaxy formation and evolution .
one possible explanation for the slope of the cmr is that massive , and thus more luminous , galaxies retain more metals than less massive ones
. the deeper potential wells of massive galaxies may allow fewer metals to be expelled in supernova driven winds @xcite .
alternatively , @xcite have shown that the color - magnitude relation can be reproduced in hierarchical merging models if the epoch of mergers occurs at very early times . since the age of a stellar population is degenerate with its metallicity on the cmr , however
, metallicity may not be the sole parameter behind the relationship @xcite . while both age and metallicity effects
may create the cmr , studies of the evolution of the relationship toward high - redshift argue that the relation is driven primarily by metallicity and can not be generated by age effects alone @xcite .
the color - magnitude relation of galaxies in clusters has been studied extensively . using high precision @xmath2 and @xmath3 photometry , @xcite determined the dispersion around the cmr for early - type galaxies in coma and virgo to be small , @xmath4 milli - mag ( mmag ) .
this small scatter implies early - type galaxies in clusters must be uniformly old with formation epochs earlier than a @xmath5 .
@xcite found @xmath6 mmag for clusters at @xmath7 ; this scatter , combined with the early age of the universe at the observed epoch , indicates that elliptical galaxies must have stopped forming stars in these clusters by a @xmath8 for a simple approximation to the color evolution of stellar populations .
the cluster ms 1054 - 03 at a @xmath9 has a 29 mmag dispersion in the @xmath10 color @xcite , while @xcite and @xcite showed that the scatter in single clusters is only 40 and 30 mmag at @xmath11 .
these small dispersions , seemingly independent of redshift , indicate that the coordination of star formation between galaxies in clusters is quite strong . using a sample of 158 clusters with @xmath12 drawn from the early data release of the sloan digital sky survey ( sdss )
, @xcite found that the cmr between clusters is very homogeneous .
locally , the cmr is universal between clusters of various masses as well @xcite .
while the mean color of red sequence galaxies in over - dense regions is only slightly redder than for similar galaxies in less dense regions and the slope of the cmr is independent of environment @xcite , several studies show that the scatter around the cmr is dependent on galaxy environment .
a re - analysis of the data presented in @xcite indicates the dispersion around the cmr in clusters is 2-@xmath13 smaller than that for galaxies in groups or in the field @xcite . the scatter in galaxy colors in the core of the cluster cl 1358 + 62 at @xmath14
have @xmath15 mmag while early - types at large radii ( @xmath16 mpc ) have nearly double that value , indicating that the dispersion of galaxy colors in the field may be larger than that measured in clusters @xcite .
early - type galaxies in the hubble deep field have a scatter of 120 @xmath17 60 mmag in rest - frame @xmath18 with an average redshift of 0.9 @xcite . while this measurement is rather uncertain , it is larger than that measured in clusters , further indicating that field galaxies may have less coordinated star formation than galaxies in clusters .
a majority of the work on the cmr has focused on early - type galaxies in clustered environments with luminosities near @xmath19 .
this is natural , as statistically significant samples are more easily selected in these over - dense regions and @xmath19 galaxies are more common than their more massive analogs . in this work ,
we extend the analysis of the color - magnitude relation to a sample of nearly 20,000 very luminous red field galaxies from the sloan digital sky survey .
these galaxies are not selected to reside in clustered environments and are chosen to be quite luminous with @xmath20 .
we create two galaxy samples from the sdss spectroscopic data . at low - redshifts ( @xmath21 )
we use galaxies from the entire survey region while at higher redshifts ( @xmath22 ) , where the galaxies are apparently fainter , we select galaxies from a 270 deg@xmath23 region which has been imaged an average of 10 times and up to 29 times
. the repeated imaging of our moderate - redshift galaxies makes comparisons between the galaxy samples possible with similar fidelity .
our sample allows us to probe a new area in parameter space  luminous galaxies in the field  and compare our results to those found for galaxies in clustered environments .
the paper is organized as follows :  2 describes the sloan digital sky survey and the selection for the galaxies used here . in  3 , we present our measurements of the slope and scatter of the color - magnitude relations .
we discuss simple star formation history models in  4 and use our measurements of the scatter around the cmr as a constraint on the evolutionary history of massive early - type galaxies in the field . in  5 , we discuss the role environment plays on the cmr of galaxies .
we construct the composite spectra of galaxies across the red sequence in  6 before closing in  7 . throughout this work ,
we use @xmath24 = @xmath25 km s@xmath26 mpc@xmath26 and ( @xmath27 , @xmath28 ) = ( 0.3,0.7 ) to calculate luminosities and distances .
when calculating ages of stellar populations , we set @xmath29 . all magnitudes used here
have been corrected for reddening using the @xcite dust maps .
the sloan digital sky survey @xcite is imaging @xmath30 steradians of the northern sky through 5 passbands @xcite .
the imaging is conducted with a ccd mosaic in drift - scanning mode @xcite on a dedicated 2.5 m telescope located at apache point observatory in new mexico .
after the images are processed @xcite and calibrated @xcite , targets are selected for spectroscopy with two double - spectrographs mounted on the same telescope using an automated spectroscopic fiber assignment algorithm @xcite .
two spectroscopic galaxy samples are created using the sdss imaging .
the main galaxy sample @xcite is a complete , flux limited , sample of galaxies with @xmath31 .
this cut is nearly five magnitudes brighter than the sdss detection limit of @xmath32 22.5 , and thus the photometric quantities for these galaxies are measured with signal - to - noise of a few hundred .
the luminous red galaxy ( lrg ) sample @xcite selects luminous early - type galaxies out to @xmath33 with @xmath34 19.5 using several color - magnitude cuts in @xmath35 , @xmath36 , and @xmath37 .
the combination of these two samples allows us to study the broadband colors of luminous field galaxies at moderate - redshifts with statistically significant samples .
while the focal point for the sdss is a contiguous survey of the northern galactic cap , the sdss also conducts a deep imaging survey , the sdss southern survey , by repeatedly imaging an area on the celestial equator in the southern galactic cap .
currently , objects in the 270 deg@xmath23 region have been imaged an average of 10 times and up to 29 times , resulting in improved photometric precision for faint galaxies .
objects detected in each observational epoch were matched using a tolerance of 0.5 arcseconds to create the final coadded catalog used here .
the photometric measurements from each epoch were combined by converting the reported asinh magnitudes @xcite into flux and then calculating the mean value .
errors on each parameter are simply the standard deviation of the flux measurements .
the southern survey is an ideal region to compare the properties of faint galaxies at moderate - redshift with brighter low - redshift analogs drawn from the entire sdss survey area with photometry of similar fidelity .
several methods are used to measure galaxy fluxes in sdss .
below , we briefly describe the two flux measurements used throughout this paper . the petrosian ratio , @xmath38 , the ratio of the local surface brightness at radius @xmath39 to the average surface brightness at that radius , is given by @xmath40}{\left .
\int_0^\theta dr ' 2\pi \theta ' i(\theta ' ) \right/ [ \pi \theta^2]},\ ] ] where @xmath41 is the azimuthally averaged surface brightness profile of a galaxy and @xmath42 and @xmath43 are chosen to be 0.85 and 1.25 for sdss .
the petrosian flux is given by the flux within a circular aperture of @xmath44 , where @xmath45 is the radius at which @xmath38 falls below 0.2 .
in sdss , @xmath45 is determined in the @xmath36 band then subsequently used in each of the other bands .
this flux measurement contains a constant fraction of a galaxy s light , independent of its size or distance , in the absence of seeing .
more details of sdss petrosian magnitudes can be found in @xcite , and @xcite .    for each galaxy imaged by sdss ,
two seeing - convolved models , a pure @xcite profile and a pure exponential profile , are fit to the galaxy image .
the best - fit model in the @xmath36 band is used to measure the flux of a galaxy in each of the other bands .
these model magnitudes are unbiased in the absence of color gradients and provide higher signal - to - noise ratio colors than petrosian colors .
a more complete description of model magnitudes is given in @xcite . throughout this paper
, we use petrosian magnitudes when calculating galaxy luminosities while model magnitudes are used to determine galaxy colors .
we calculate photometric @xmath46-corrections using the method of @xcite ( ` kcorrect v3_2 ` ) .
this program constructs a linear combination of carefully chosen spectral templates in order to best match the observed photometry at the measured redshift of the galaxy .
the rest - frame colors of these best - fit spectra are used to derive corrections to the observed galaxy colors . with this method
, it becomes simple to transform between bandpasses and , if necessary , photometric systems given a detailed understanding of the filter characteristics .
after calculating @xmath46-corrections with ` kcorrect ` , we remove any mean @xmath46-corrected color versus redshift trends with a low order polynomial to account for the passive evolution of stellar populations between observed epochs .
this evolutionary correction is normalized to a redshift of 0.3 , near the median of the lrg redshift distribution .
the color evolution correction is small ( less than 2% for all bands in both samples ) and does not affect any of the results presented here .
we further assume a galaxy dims by one magnitude per @xmath47 change in redshift in the @xmath35-band due to the passive evolution of its stellar populations .
again , this correction is normalized to @xmath48 . in the remainder of this paper , all of the color and luminosity measurements refer to quantities which have been adjusted for the @xmath46-correction and evolution .
cc @xmath49 & 4026 + @xmath50 & 5307 + @xmath51 & 6441 + @xmath52 & 7687 + @xmath53 & 3408 + @xmath54 & 4494 + @xmath55 & 5454 + @xmath56 & 6509 +    throughout this paper ,
we work in the @xmath49 , @xmath50 , @xmath51 , @xmath52 and @xmath53 , @xmath54 , @xmath55 , @xmath56 systems for the low- and moderate - redshift galaxy samples , respectively . here
, the notation @xmath57 indicates the reported ab magnitudes are derived through a standard @xmath35 filter that has been blueshifted by @xmath58 . for a filter system which has been shifted by @xmath58 ,
the @xmath46-correction for galaxies at that redshift is trivial , @xmath59 , and is independent of the galaxy spectral energy distribution . in choosing shifts which nearly match the median redshift of each sample , we minimize the uncertainty introduced through the @xmath46-corrections @xcite .
table 1 lists the effective wavelengths for each of the bandpasses used here .      in order to investigate the color - magnitude relation of massive galaxies in the range @xmath60 , we create two samples from the sdss spectroscopic database to be analyzed separately ; figure [ sample ]
shows the selection cuts used to define each of these samples .
from the main sample of galaxies , we extracted 16622 galaxies in the redshift range @xmath21 satisfying the rest - frame color cut @xmath61 which approximates a cut of constant stellar mass @xcite .
this selection criterion results in a cut diagonally across the red sequence , which could bias our measurements of the scatter in red sequence colors and will greatly affect any slope measurement of the color - magnitude relation . in order to avoid biased slope measurements and test the consistency of our scatter measurements , we also construct a sample of 5100 galaxies in the luminosity range @xmath62 with @xmath63 for comparison ( shown by the dashed region in figure [ sample ] ) .
we further select 2782 lrg galaxies from the sdss southern survey in the redshift range @xmath22 which satisfy @xmath64 and have been imaged 6 or more times .
equation 3 is the same rest - frame color cut applied to the low - redshift galaxies ( equation 2 ) shifted to the @xmath65 color  the effective wavelengths of the filters that define each of these colors are similar , and thus these colors probe similar regions in the rest - frame galaxy spectrum ( see table 1 ) . of these 2782 sample galaxies
, 62% were imaged more than 9 times while 20% of the galaxies were imaged more than 12 times in the sdss southern survey . as in the low - redshift sample
, we adopt a simple constant color cut , @xmath66 and @xmath62 , to define a subsample of 1380 high - redshift galaxies to test the consistency of our scatter measurements and allow us to determine the slope of the color - magnitude relation without a strong cut across the red sequence .
we find no differences within our quoted error between the dispersion around the cmr measured from galaxies defined using the stellar mass cut and the constant color cut ; throughout the remainder of this paper , scatter measurements will be based upon galaxies drawn in the first manner while all slope measurements are referenced to samples created using the latter method .
figure [ zhist ] shows the redshift distribution of the galaxies used to measure the scatter around the color - magnitude relation in this paper .    in detail ,
the red sequence is populated by elliptical and s0 galaxies as well as late - type galaxies reddened by dust .
the arrow in the lower right of figure [ sample ] shows the reddening vector given by the @xcite extinction curve . in order for a dusty late - type galaxy to scatter into our sample
, it must have a large unobscured luminosity as we only consider the most luminous red - sequence galaxies in this work .
there is a clear paucity of very luminous blue galaxies in figure [ sample ] and thus the contamination from dusty spiral galaxies is likely small in our sample .
furthermore , @xcite constructed the average spectrum of massive galaxies , such as those used in our analysis , and found that only 5% of the galaxy spectra in the @xmath67 luminosity range were contaminated by emission lines uncharacteristic of early - type galaxies . in the remainder of this paper
, we will use the red - sequence and early - type classifications synonymously .    throughout this paper
, we will refer to the galaxies considered here as field galaxies as they are not specifically selected to reside in clustered environments .
it is important to note , however , that galaxies which are of similar color and luminosity as those in our sample tend to lie in over - dense environments @xcite .
@xcite found the clustering strengths and mean separations of lrgs are comparable to those of poor clusters or rich groups . for comparison , of the lrgs in the volume surveyed by @xcite to find galaxy clusters with
13 or more detected cluster members , 16% are projected less than 500 kpc from the inferred location of a cluster and within 0.05 of the photometric redshift for that cluster .
a detailed examination of the density dependence of the scatter around the cmr is conducted in  5 .
when determining the slope of the color - magnitude relation , we only consider the galaxies selected with our simple luminosity and constant color cuts as any slope measurements based on samples with a strong cut across the red sequence will be quite biased .
the galaxy samples considered in this work only span a small range of luminosities and thus are not ideally suited for measurements of the slope of the cmr  we simply use our measurement of this quantity to remove the mean relationship before determining the scatter around it .
the slopes and zero - points for the cmr for several colors are listed in table 2 .
figure [ fig : slopeplot ] shows the measured slopes of the cmr as a function of the effective wavelengths of the two bandpasses used to define the color .
the slope in bluer colors is more pronounced than that for redder bands which tend to show small or negligible slopes .
this trend is consistent with a metallicity sequence as the primary driver along the red - sequence . as noted by @xcite using models from @xcite , the flattening of the color - magnitude relation toward redder colors can also be reproduced by an age trend along the red sequence but only if a very specific relationship between galaxy mass and star formation time holds with very little scatter . in the bluest colors , which have been used historically in the literature , the slopes measured here are in general agreement with past work ; for comparison , @xcite report a slope of @xmath68 mag mag@xmath26 in the @xmath69 color for a large sample of galaxies at @xmath70 .
rcrrcrrcr[hb ] @xmath71 & 1.162 & @xmath72 & & 1.162 & @xmath73 & & 1.164 & @xmath74 + @xmath75 & 0.449 & @xmath76 & & 0.449 & @xmath77 & & 0.451 & @xmath78 + @xmath79 & 0.329 & @xmath80 & & 0.329 & @xmath81 & & 0.329 & @xmath82 + @xmath83 & 1.611 & @xmath84 & & 1.610 & @xmath85 & & 1.612 & @xmath86 + @xmath87 & 1.750 & @xmath88 & & 1.750 & @xmath89 & & 1.754 & @xmath90 + @xmath91 & 0.613 & @xmath92 & & 0.612 & @xmath93 & & 0.614 & @xmath94 + @xmath95 & 0.379 & @xmath96 & & 0.377 & @xmath97 & & 0.381 & @xmath98 + @xmath65 & 0.991 & @xmath99 & & 0.980 & @xmath100 & & 0.983 & @xmath101 +    figure [ fig : sigmaresid ] illustrates the typical scatter about the mean color - magnitude relation for the @xmath75 and @xmath91 colors ; the measured dispersion about this relationship is quite small .
the second column of table 3 lists the measured dispersion around the cmr for each of the colors used in this study . for
each color studied , we fit a gaussian to the observed color distribution using poisson errors on each color bin and define the dispersion of the best - fit gaussian to be the scatter in galaxy colors on the red sequence . by fitting the histograms rather than calculating the standard deviation directly , we place more weight on the core of the distribution compared to the wings .
thus , our scatter measurements will not be heavily affected by any non - red - sequence interlopers in the wings of the color distributions . as this method
could potentially depend on the binning used , we are careful to vary the bin size for each fit and ensure the measured scatter is not dependent on our choice of bin size .
it should be noted that past work has used the standard deviation to quantify the scatter in galaxy colors . while the standard deviation and gaussian widths may give systematically different dispersions in the case of non - gaussian distributions , figure [ fig : sigmaresid ]
shows the color distribution on the red sequence is well described by a gaussian and thus the two measurements are comparable .
bootstrap resampling of the data ( e.g. * ? ? ?
* ) was employed to determine the uncertainty on our scatter measurements .
it is important to note that the scatter measurements in the second column of table 3 represent upper limits to the intrinsic scatter in galaxy colors as these values are not corrected for systematic or experimental effects , which could be significant .    in order to investigate the level of skewness of the observed color distribution , we define @xmath102 to be the absolute value of the color difference between the @xmath103 quantile and the median of the distribution
the level of asymmetry in the distribution is then estimated using @xmath104 for this statistic , a maximally skewed distribution would have @xmath105 such that a tail of galaxies to the red would result in @xmath106 . a symmetric distribution would have @xmath107 .
we adopt this approach rather than a traditional third moment calculation as the calculation of higher moments of a distribution is strongly dependent on the wings of the distribution where contamination may be important . for the histograms shown in figure [ fig : sigmaresid ]
, we find that @xmath108 for the @xmath75 color and @xmath109 for the @xmath91 color .
both of these measurements reflect the slight over - abundance of red galaxies compared to a gaussian in both panels of figure [ fig : sigmaresid ] .
this red tail is due to the diagonal cut across the red sequence imposed by our stellar mass selction criterion ; red galaxies are selected to lower luminosities than blue galaxies .
since the luminosity function of galaxies is steeply rising toward lower luminosities in this range , more red galaxies are selected than blue galaxies . for comparison ,
we measure @xmath110 for the @xmath75 color and @xmath111 for the @xmath91 color when the simple luminosity cuts are used to select our sample galaxies .    since our intermediate redshift sample is drawn from regions that have been repeatedly imaged by the sloan digital sky survey , we can assess the level to which the dispersion we observe in the color - magnitude relation is dependent on measurement noise .
if the scatter is contaminated by measurement noise , we would expect the squared scatter to decrease inversely with the number of observations : @xmath112 figure [ fig : richardson ] shows the results of this test ; the methods used for the ensemble sample were used in the same manner to calculate the scatter and error for each bin . in each bin of @xmath113-measurements , we include galaxies which were imaged a total of @xmath113 times , and thus each bin contains a unique subsample of galaxies .
for the moderate - redshift sample , the third and fourth columns in table 3 lists the instrumental contribution and intrinsic galaxy scatter predicted by extrapolating the fits to equation 2 to infinite observations for each of the three colors .
cccccc[hb ] @xmath71 & @xmath114 & @xmath115 & @xmath116 & @xmath117 & @xmath118 + @xmath75 & @xmath119 & @xmath120 & @xmath121 & @xmath122 & @xmath123 + @xmath79 & @xmath124 & @xmath125 & @xmath125 & @xmath126 & @xmath127 + @xmath83 & @xmath128 & @xmath129 & @xmath130 & @xmath131 & @xmath132 + @xmath87 & @xmath133 & @xmath134 & @xmath135 & @xmath136 & @xmath137 + @xmath91 & @xmath138 & @xmath139 & @xmath140 & @xmath141 & @xmath142 + @xmath95 & @xmath143 & @xmath144 & @xmath145 & @xmath146 & @xmath147 + @xmath65 & @xmath148 & @xmath149 & @xmath150 & @xmath151 & @xmath152 +    while the intrinsic dispersion around the color - magnitude relation at moderate redshift can be extracted from the measured scatter based on subsamples of galaxies with different number of measurements , the low - redshift galaxies do not allow for this test as most galaxies in this sample have been imaged only once . in order to estimate the intrinsic scatter in galaxy colors for this sample , we identified 5839 of our low - redshift galaxies that have been imaged more than once under photometric conditions .
we adopt the average root - mean squared variation of each photometric measurement as the instrumental error for each color of interest .
we find the mean measurement uncertainty to be 17.6 , 15.1 , and 20.9 mmag in the @xmath71 , @xmath75 , and @xmath79 bands , respectively .
the internal dispersion in each color is the difference , in quadrature , between the measured scatter and inferred mean uncertainty in our single pass photometry .
for the low - redshift galaxies , the third and fourth columns of table 3 lists the observational error and intrinsic scatter determined by this correction .
while our estimate of the instrumental contribution to the observed color dispersion corrects for uncorrelated errors on each independent measurement , we do not correct for systematic errors common to all observations which could cause an overestimation of the intrinsic dispersion in galaxy colors .
figure [ fig : sigmaplot ] shows our scatter measurements as a function of the effective wavelength of the filters used to construct each color .
the scatter about the cmr is clearly wavelength dependent ; the bluest colors , which are most affected by small changes in the age , metallicity , or dust attenuation of the system , have systematically larger dispersions than the reddest colors where changes in metallicity and age have less impact .
here , we consider three simple star formation histories and derive the range of permitted parameters for these histories following a method similar to that of @xcite . in these toy models , we will assume that only age variations between red - sequence galaxies creates the scatter in early - type galaxy colors .
in reality , the scatter across the cmr is likely driven by some combination of age , metallicity , dust , and possibly @xmath153-enhancement variations across the red - sequence .
these models will thus illustrate the broadest range of ages allowed by our measurements ; the inclusion of any variation in the other galaxy properties will only constrain the age more strongly .    in this section
, we will primarily consider the high - redshift sample , as these are the intrinsically youngest galaxies in our sample and hence should provide the tightest constraints on the past star formation history of massive early - type galaxies .
an identical analysis of the low - redshift sample provides results consistent with those presented below . in all three model cases ,
the low - redshift galaxies provide looser constraints than the high - redshift sample ; star formation is allowed to continue to slightly later times .    first , we consider a star formation history composed of a single delta - function burst .
the probability of a burst occurring is distributed uniformly between the earliest allowed time for galaxy formation , @xmath154 , set to @xmath155 , and @xmath156 , the time at which the youngest early - type galaxies form . using a grid of monte - carlo simulations with 5000 galaxies per realization ,
we create samples of galaxy spectra synthesized using the method of @xcite with solar metallicity , a @xcite initial mass function , and varying values of @xmath156 .
we find that realizations allowing galaxies younger than 8.8 gyr ( @xmath157 ) generate color distributions with standard deviations larger than that measured .
we perform a similar test using a scheme in which 80% of the stars in a galaxy are formed in a single burst at @xmath155 and 20% are formed in a second burst which occurs at a random time between @xmath155 and a minimum age , as done by @xcite . for this scenario , the secondary bursts must occur between @xmath158 to reconstruct the observed scatter in the color - magnitude relation . for bursts that generate less than 5% of a galaxy s stars ,
the time scale for the secondary burst is unconstrained - small amounts of star formation can occur at late epochs without violating the constraints placed by our measurements .
finally , we create a grid of models in which star formation begins at @xmath155 and continues at a constant rate until it is truncated at @xmath159 . for each realization , we allow @xmath159 to range from @xmath155 to a critical cutoff time at which all galaxies have stopped forming stars . star formation in this manner must have ended by @xmath160 in order to reproduce the observed scatter in the color - magnitude relation at @xmath161 .
the allowed range of parameters and the distribution of galaxy colors found in each model are shown graphically in figure [ fig : models ] . to allow comparison with figure [ fig : sigmaresid ] , typical observational errors have been included in the color distributions in the lower panels of the figure .
interestingly , all of these models predict similar scatter trends to those seen in figure [ fig : sigmaplot ] . in all three cases ,
the scatter in the @xmath87 color is nearly double that predicted for the @xmath91 and @xmath95 colors .
while each of the three models succeeds in reproducing the observed scatter in early - type galaxy colors , none of the models adequately recreates the color distribution measured from real galaxies , a sign that the models chosen do not perfectly track the true evolution of the star formation history of early - type galaxies .
the @xmath162 parameter , as defined in equation 4 , for each of the model distributions is shown in figure [ fig : models ] .
the three models produce color distributions with @xmath163 compared to the value of @xmath164 calculated for the observed sample of galaxies at the same redshift .
the difference between the observed skewness and that determined from our simple models reiterates the mismatch between the observed and predicted color distributions . in our simulations , we assume all galaxies have the same metallicity and no dust attenuation . also , the onset of star formation is uniform in the second two models . in reality , these parameters are not constant , and thus the sharp truncations seen on the red edges of model color distributions will be blurred by the addition of other complications to the models .
however , any of these corrections will increase the observed scatter in galaxy colors thus driving the last epoch of star formation needed to satisfy our constraints to earlier times .
a majority of the work on the color - magnitude relation of early - type galaxies has focused on cluster galaxies at various redshifts . while some work ( e.g. * ? ? ? * ; * ? ? ?
* ; * ? ? ?
* ) tackles the problem in the field , these studies suffer from small numbers of galaxies .
the measurement of the scatter around the cmr presented here , on the other hand , focuses on the most massive field galaxies at moderate redshifts using a statistically significant sample .
we can now , for the first time , statistically compare the dispersion around the cmr in both clustered and field environments .    in order to explore the environmental dependence of the color - magnitude relation within our data set , we use sdss imaging of 5305 deg@xmath23 to identify a comparison sample of 16 million normal galaxies down to @xmath165 . at the location of each lrg galaxy in our sample , we count the number of neighboring galaxies within a @xmath166 mpc ( proper ) aperture which have colors similar to those expected for an @xmath19 galaxy on the red sequence , and luminosities in the @xmath167 to @xmath168 range ; the details of the luminosity and color cuts are given in @xcite .    using the number of neighbor galaxies as a proxy for the environment of each lrg in our sample , we divide each of our galaxy sets into a clustered subsample and a low density subsample .
we define lrg galaxies with more than 5 ( 7 ) neighboring galaxies in the low- ( moderate- ) redshift samples to be in clustered environments while a field subsample is comprised of galaxies with fewer neighbors .
this threshold was chosen such that the clustered sample in both redshift bins composes 30% of the total galaxy sample at that redshift . to place these values in physical context , at @xmath169 ( @xmath161 ) there are 25.3 ( 170 ) tracer galaxies per square degree . given the angular diameter distance of 398 ( 712 ) @xmath170 mpc to that redshift , this yields a total of 1.65 ( 3.46 ) background galaxies per aperture on average .
the average number of neighbors in the field sample is 2.3 ( 3.9 ) galaxies while the clustered galaxies have a mean of 9.1 ( 11.4 ) neighbors .
thus our clustered samples are over - dense by a factor of 4 ( 3 ) compared to the field galaxies .
@xcite found that the slope and zeropoint of the color - magnitude relation do not depend strongly on the local galaxy environment .
we find the same result ; table 2 lists the measurements of the slope and zeropoint of the relationship for the ensemble sample as well as for the field and clustered subsamples .
the variation between the two samples is small . in the following calculations ,
we use the ensemble slope and zeropoint measurements as these are measured with the highest signal - to - noise .
the results are unaffected if the individual measurements of the slope and zeropoint are used for each subsample .
the fifth and sixth columns of table 3 list the measured scatter in galaxy colors for the field and clustered galaxies respectively . in all of the bands ,
the scatter for the field subsample is larger than that in clustered environments .
if we assume the scatter introduced by instrumental effects is the same for both the field and clustered subsamples , the same correction determined for the ensemble galaxy sample can be applied to these subsamples , as well . with this correction , we find that the scatter in galaxy colors in dense environments is @xmath171 smaller than for that of isolated systems in both redshift intervals studied here .
our clustered galaxies reside in a range of densities , not simply rich clusters , as have often been used for similar studies , and thus we likely underestimate the significance of the difference in scatter between the field and strongly clustered environments .
furthermore , our field sample is not completely composed of isolated galaxies as the average number of neighbor galaxies in this sample is larger than the average density of background galaxies .
again , this leads to an underestimation of the true difference in the dispersion around the color - magnitude relation for field and clustered galaxies .
the larger scatter in field galaxy colors supports previous claims of an environmental dependence for the scatter around the color - magnitude relation for early - type galaxies ( e.g. * ? ? ?
* ; * ? ? ?
* ) , though the effect measured here is weak .
comparisons of the scatter in galaxy colors for clustered galaxies often draw on galaxies from a single cluster whereas our clustered sample is created by combining galaxies from a number of clustered environments .
it is possible that the star formation histories of galaxies in a single cluster are more highly coordinated than the star formation histories between different clusters , thus the scatter measured using galaxies from a single cluster would be smaller than when several clusters are combined . as the fraction of blue galaxies is larger in the field compared to clustered environments , the field subsample may have a larger contamination fraction compared the the clustered subsample .
this would also cause an increased scatter in field galaxies colors compared to the galaxies in more dense environments .
having measured the scatter in galaxy colors , we next investigate the spectral differences between galaxies by examining the average spectrum of massive early - type galaxies across the red sequence .
the composite early - type galaxy spectrum has been investigated as a function of environment , luminosity , redshift , size , and velocity dispersions @xcite ; here , we construct the composite spectrum for red sequence galaxies of different mean colors .
since the galaxies in our low - redshift sample are brighter , and thus have higher quality spectra , and are more numerous than the higher redshift galaxies , we only consider these galaxies in our analysis .    to construct the average spectrum of the galaxies in each color bin , we first shift the spectra according to the observed redshift and then normalize the observed sdss spectrum by the flux at rest - frame 5000 @xmath172 .
we augment the sdss pixel mask to include all pixels within 10 @xmath172 of the bright sky lines at 5577 , 5890 , 5896 , 6300 , 6364 , 7245 , and 7283 @xmath172 .
since the galaxies of interest are in the redshift range @xmath21 , the masked pixels occur at different rest - frame wavelengths for each galaxy , leaving the mean spectrum unaffected .
we then find the number - weighted mean of the spectra , ignoring all masked pixels in the mean .
galaxies are divided into two distinct color bins  red galaxies having @xmath71 colors between 1-@xmath13 and 3-@xmath13 redward of the mean and a blue sample in the same range blueward of the mean .
the red galaxy bin contains 1321 galaxies with a mean @xmath71 color of 1.20 and @xmath173 while the blue sample contains 1161 galaxies with mean color and luminosity of 1.09 and @xmath174 .
figure [ fig : coadd](a ) shows the coadded spectra of the red and blue galaxies .
the gross properties of each composite spectrum are similar .
figure [ fig : coadd](b ) shows the difference between the blue and red spectra .
this difference spectrum shows clear signs of enhanced balmer absorption in the blue galaxies compared to the red composite , indicating that , on average , the blue galaxies contain more young stars , in agreement with past studies ( e.g. * ? ? ?
* ; * ? ? ?
* ) . while the red composite shows stronger [ n  @xmath175ii ] emission lines relative to the blue composite , no other prominent emission lines , such as [ o  @xmath175ii ] , [ o  @xmath175iii ] , and h@xmath153 are present ; this emission line difference is likely due to a larger fraction of agn in the red galaxies compared to the blue
we consider two possible evolutionary scenarios to explain the spectral differences between the blue and red galaxies on the red sequence : pure age ( i.e. passive ) evolution of a single age stellar population and a recent epoch of star formation mixed with an existing old population . in the first scenario
, we fit each composite spectrum with a solar metallicity , single age , stellar synthesis model created using the method of @xcite .
the red and blue spectra are best fit by 6.3 gyr and 3.7 gyr populations , respectively . in the second scenario
, we find the blue composite spectrum is well fit by a 6.3 gyr old stellar population mixed with a 1 gyr population containing 3.5% of the stellar mass and we model the red composite with a simple 6.3 gyr population
. the difference spectra for both scenarios are shown in figure [ fig : coadd](b ) .
while each of the two model hypotheses would suggest quite different evolutionary histories for massive early - type galaxies , the predicted difference spectra are remarkably similar .
the model consisting of a recent burst of star formation slightly over - estimates the continuum in the difference spectrum near 3900 @xmath172 and thus provides a marginally poorer fit to the data than the passive evolution model .
the residuals between the data and the model difference spectra are shown in figure [ fig : coadd](c ) .
both of the models fail to match the strengths of the @xmath176 and mg @xmath177 features near 4140 and 5160 @xmath172 .
@xcite interpret these differences between observations and synthetic spectra as a enhanced @xmath153-to - iron ratio present in the galaxies but not included in the synthetic models .
also , several weak iron features are evident in the two model spectra near 4300 and 5300 @xmath172 which are missing in the observed difference spectrum . in a pure - age model for the evolution across the red sequence ,
the enhanced blue continuum provided by the hotter stars in the young population dilutes the metal lines seen in the older population .
this dilution leads to differences in the strengths of the metal lines of the two population models , thus causing the residual iron features in both model difference spectra .
the absence of these residuals in our observed difference spectrum indicates that a pure - age model for the spectral differences between red and blue galaxies on the red sequence is an incomplete model . in order to match the variations of the spectral features across the red sequence ,
a more detailed model is required .
this model should include both age and metallicity variation and , more importantly , non - solar abundance ratios .
until spectral synthesis models with realistic @xmath153-enhancement prescriptions are available , more quantitative analysis of our observed difference spectrum is difficult as the interplay between age , metallicity , and @xmath153-to - iron ratio variations must be included for a full understanding of the physical differences between these galaxies .
another clear difference between the models and the data is the strong na d feature near 5892 @xmath172 . here
, the red composite spectrum has considerably more na d absorption than predicted by the stellar synthesis models . as the model spectra do not contain absorption due to the interstellar medium
, this discrepancy may indicate that the red galaxies contain more neutral gas than their blue counterparts .
we present new multi - color observations of a large sample of galaxies in the range @xmath1 gathered from the sloan digital sky survey .
these galaxies represent the most massive galaxies in the universe and are not selected to be in clustered environments .
past studies of the color - magnitude relationship have focused on @xmath178 galaxies in clusters , thus our work focuses on an area of parameter space previously unexplored .
using the robust photometry of bright galaxies provided by the sdss as well as multiple observations of moderate - redshift galaxies , we are able to construct the color - magnitude relations for massive field galaxies in several colors . while our data only span one magnitude in luminosity , and thus are not ideal for in - depth analysis of the slope and zeropoint of the cmr , our large sample of galaxies allows us to measure the dispersion about this relationship quite well .
the scatter around the red sequence is 20 mmag in the reddest bands while the bluest colors have nearly double this value .
we construct several toy model star formation histories in order to place constraints on the evolutionary processes at work in massive galaxies .
we find that the majority of star formation in these early - type galaxies must have occurred before @xmath179 , in agreement with other authors . for secondary bursts of star formation with strengths less than 5% of the total stellar mass of the system , our data are consistent with bursts of star formation at any epoch .
however , one should note that none of the toy models reproduce the observed distribution of galaxy colors .
furthermore , we make comparisons between our new measurement of the scatter around the color - magnitude relation for massive field galaxies and early - type galaxies in clusters .
an internal comparison shows that the scatter around the cmr in dense environments is @xmath180 smaller than that measured using more isolated galaxies .
the composite spectra of lrg galaxies redward and blueward of the mean galaxy color are physically distinct .
the blue galaxies show stronger balmer absorption than the red galaxies , consistent with the presence of a younger stellar population .
simple models slightly favor passive evolution of stellar populations as the cause of the observed differences in spectra across the red sequence compared to a more recent injection of young stars onto an old stellar population , though the difference between the predicted spectra from these models is small .
we would like to thank michael blanton , scott burles , douglas finkbeiner , david hogg , david schlegel , john moustakas , amelia stutz , and jane rigby for software and helpful comments during the preparation of this paper .
we further thank the anonymous referee for a thorough and critical reading of our paper which lead to substantial improvements to the text .
this research made use of the nasa astrophysics data system .
rjc is funded through a national science foundation graduate research fellowship .
dje is supported by nsf grant ast 339020 and an alfred p. sloan research fellowship .
funding for the creation and distribution of the sdss archive has been provided by the alfred p. sloan foundation , the participating institutions , the national aeronautics and space administration , the national science foundation , the u.s .
department of energy , the japanese monbukagakusho , and the max planck society .
the sdss web site is http://www.sdss.org/.    the sdss is managed by the astrophysical research consortium ( arc ) for the participating institutions .
the participating institutions are the university of chicago , fermilab , the institute for advanced study , the japan participation group , the johns hopkins university , the korean scientist group , los alamos national laboratory , the max - planck - institute for astronomy ( mpia ) , the max - planck - institute for astrophysics ( mpa ) , new mexico state university , university of pittsburgh , university of portsmouth , princeton university , the united states naval observatory , and the university of washington .
abazajian , k. , et al .
2003 , , 126 , 2081 abazajian , k. , et al .
2004 , , 128 , 502 abazajian , k. , et al .
2005 , , 129 , 1755 andreon , s.  2003 , , 409 , 37 arimoto , n.  & yoshii , y.  1987 , , 173 , 23 bahcall , n.  a. , et al .
2003 , , 148 , 243 baldry , i.  k. , glazebrook , k. , brinkmann , j. , ivezi ' c ,  . ,
lupton , r.  h. , nichol , r.  c. , & szalay , a.  s.  2004 , , 600 , 681 bell , e.  f. , & de jong , r.  s.  2000 , , 312 , 497 bernardi , m. , et al .
2003 , , 125 , 1882 bernardi , m. , sheth , r .
k. , nichol , r.  c. , schneider , d.  p. , & brinkmann , j.  2005 , , 129 , 61 blakeslee , j.  p. , et al .
2003 , , 596 , l143 blanton , m.  r. , et al .
2001 , , 121 , 2358 blanton , m.  r. , et al .
2003 , , 594 , 186 blanton , m.  r. , lin , h. , lupton , r.  h. , maley , f.  m. , young , n. , zehavi , i. , & loveday , j.  2003 , , 125 , 2276 blanton , m.  r. , et al .
2003 , , 125 , 2348 bower , r.  g. , lucey j.  r. , & ellis , r.  s.  1992 , , 254 , 601 bressan , a. , chiosi , c. , & tantalo , r.  1996 , , 311 , 425 bruzual , g.  & charlot , s.  2003 , , 344 , 1000 de lucia , g. , springel , v. , white , s.  d.  m. , croton , d. , & kauffmann , g.  2005 , arxiv astrophysics e - prints , arxiv : astro - ph/0509725 de vaucouleurs , g.   1948 , annales dastrophysique , 11 , 247 eisenstein , d.  j. , blanton , m. , zehavi , i. , bahcall , n. , brinkmann , j. , loveday , j. , meiksin , a. , & schneider , d.  2005 , , 619 , 178 eisenstein , d.  j. , et al .
2003 , , 585 , 694 eisenstein , d.  j. , et al .
2001 , , 122 , 2267 ellis , r.  s. , smail , i. , dressler , a. , couch , w.  j. , oemler , a.  j. , butcher , h. , & sharples , r.  m.  1997 , , 483 , 582 faber , s.  m.  1972 , , 20 , 361 faber , s.  m.  1973 , , 179 , 731 fukugita , m. , ichikawa , t. , gunn , j.  e. , doi , m. , shimasaku , k. , & schneider , d.  p.  1996
, , 111 , 1748 gladders , m.  d. , lopez - cruz , o. , yee , h.  k.  c. , & kodama , t.  1998 , , 501 , 571 gunn , j.  e. , et al .  1998 , , 116 , 3040 hogg , d.  w. , finkbeiner , d.  p. , schlegel , d.  j. , & gunn , j.  e.  2001 , , 122 , 2129 hogg , d.  w. , et al .  2003 , , 585 , l5 hogg , d.  w. , et al .
2004 , , 601 , l29 ivezi ' c ,  . , et al .
2004 , astronomische nachrichten , 325 , 583 kauffmann , g.  & charlot , s.  1998 , , 294 , 705 kodama , t. 1997 , ph.d .
thesis , univ .
tokyo kodama , t.  & arimoto , n.  1997 , , 320 , 41 kodama , t. , bower , r.  g. , & bell , e.  f.  1999 , , 306 , 561 larson , r.  b.  1974 , , 166 , 585 larson , r .
b. , tinsley , b.  m. , & caldwell , c.  n.  1980 , , 237 , 692 lugger , p.  m.  1984 , , 278 , 51 l ' opez - cruz , o. , barkhouse , w.  a. , & yee , h.  k.  c.  2004 , , 614 , 679 lupton , r.  h. , gunn , j.  e. , ivezi ' c , z. , knapp , g.  r. , kent , s. , & yasuda , n.  2001 , asp conf .
238:astronomical data analysis software and systems x , 10 , 269 lupton , r.  h. , gunn , j.  e. , & szalay , a.  s.  1999 , , 118 , 1406 matteucci , f.  & tornambe , a.  1987 , , 185 , 51 mcintosh , d.  h. , zabludoff , a.  i. , rix , h. , & caldwell , n.  2005 , , 619 , 193 oconnell , r.  w.  1980 , , 236 , 430 odonnell , j.  e.  1994 , , 422 , 158 pier , j.  r. , munn , j.  a. , hindsley , r.  b. , hennessy , g.  s. , kent , s.  m. , lupton , r.  h. , & ivezi ' c ,  .
2003 , , 125 , 1559 press , w.  h.  2002 , numerical recipes in c++ : the art of scientific computing by william h.  press .
xxviii , 1,002 p.  : ill .  ; isbn : 0521750334 , rose , j.  a.  1985 , , 90 , 1927 salpeter , e.  e.  1955 , , 121 , 161 sandage , a.  & visvanathan , n.  1978 , , 225 , 742 schlegel , d.  j. , finkbeiner , d.  p. , & davis , m.  1998 , , 500 , 525 smith , j.  a. , et al .
2002 , , 123 , 2121 stoughton , c. , et al .
2002 , , 123 , 485 stoughton , c. , et al .
2002 , , 4836 , 339 strauss , m.  a. , et al .
2002 , , 124 , 1810 terlevich , a.  i. , kuntschner , h. , bower , r.  g. , caldwell , n. , & sharples , r.  m.  1999 , , 310 , 445 terlevich , a.  i. , caldwell , n. , & bower , r.  g.  2001 , , 326 , 1547 trager , s.  c. , faber , s .
m. , worthey , g. , & gonz ' alez , j.  j.  2000 , , 119 , 1645 van dokkum , p.  g.,franx , m. , kelson , d.  d. , illingworth , g.  d. , fisher , d. , & fabricant , d.  1998 , , 500 , 714 van dokkum , p.  g. , franx , m. , fabricant , d. , illingworth , g.  d. , & kelson , d.  d.  2000 , , 541 , 95 van dokkum , p.  g. , stanford , s.  a. , holden , b.  p. , eisenhardt , p.  r. , dickinson , m. , & elston , r.  2001 , , 552 , l101 visvanathan , n.  & sandage , a.  1977 , , 216 , 214 worthey , g. , trager , s.  c. , & faber , s.  m.  1995 , asp conf .
ser .   86 : fresh views of elliptical galaxies , 86 , 203 york , d.  g. , et al .
2000 , , 120 , 1579 zehavi , i. , et al .  2005 , , 621 , 22 zepf , s.  e. , whitmore , b.  c. , & levison , h.  f.  1991 , , 383 , 524


discovery of the higgs boson as predicted by the standard model ( sm ) by the atlas and the cms experiments became the moment of triumph for particle physics @xcite .
such a historic discovery together with decades of eletroweak precision data have well established the validity of sm up to accessable energies .
however , there is no verified explanation of the origin of the small neutrino masses and no viable candidate for the dark matter in the sm . due to these unwavering issues
, various extensions of sm have been proposed .
the secret of neutrino masses may lie in some form of seesaw mechanism , where a sm singlet right - handed neutrinos with large majorana masses cause the light neutrino masses ( type - i seesaw ) @xcite or a sm weak scalar triplet with a tiny induced vacuum expectation value ( vev ) generates the small neutrino masses ( type - ii seesaw ) @xcite .
if neutrino masses are generated as loop corrections , the masses will naturally be suppressed and such extensions of sm are both theoretically well motivated and phenomenologically viable @xcite .
searches for a stable dark matter candidate have also been motivation for various extensions of the sm .
some form of symmetry usually stabilizes the dark matter .
simple discrete symmetries such as r - parity in supersymmetric models @xcite can perform an excellent job of preventing the particle from decaying .
kaluza - klein parity @xcite in universal extra dimension models and t - parity in the littlest higgs models @xcite can stabilize the lightest particle ; turning them into promising dark matter candidates .
a similar role is played by a @xmath5 symmetry for the case of inert doublet models @xcite or scotogenic models @xcite .
sm extended by a scalar singlet carrying a discrete @xmath5 parity is yet another example for a simple dark matter model . instead of being an adhoc symmetry , this @xmath5 symmetry can be a remnant of the @xmath6 generator of @xmath7 grand unified theories ( guts ) @xcite .
@xmath7 guts provide one of the most lucrative frameworks , where one can incorporate many of the aforementioned extensions of sm along with a beautiful unified picture of sm gauge couplings . among the classes of @xmath7 guts , supersymmetric versions have multiple features such as successful unification of gauge couplings and natural dark matter candidate owing to an automatic r - parity , while it solves the gauge hierarchy problem based on symmetry principle .
in addition , supersymmetric models offer a mechanism for triggering electroweak symmetry breaking ( ewsb ) via radiative effects @xcite . in this scenario ,
the positive mass parameter of the higgs boson at high energy becomes negative at low energy due to the renormalization group flow which dictates how the parameters evolve with scale .
the purpose of this paper is to explore the possibility of radiative ewsb breaking in non - supersymmetric ( non - susy ) extensions of the sm .
since this is an attractive mechanism to trigger ewsb , checking its viability in non - susy models is of great interest .
as we argue below , such a radiative ewsb may be necessary in certain unified theories which have two stages of symmetry breaking . in a general context a positive mass parameter for the higgs boson turning negative
also enhances the available parameter space .    in some extensions of sm ,
the higgs boson is a part of the larger multiplet , which breaks some higher symmetry .
this occurs in trinification models based on @xmath8 gauge symmetry @xcite broken by two @xmath9 higgs multiplets .
these multiplets contain sm singlet components which acquire large vevs breaking the gauge symmetry down to sm .
the same @xmath9 multiplets also contain the higgs boson of the sm which should develop a negative squared mass to trigger ewsb .
consistency of the high scale symmetry breaking , however , would demand that all physical higgs bosons , including the sm higgs , have positive squared masses at a high energy .
one could introduce new higgs fields to break the electroweak symmetry in which case the model looses its minimality and predictivity .
for such class of models one might employ radiative loop corrections to turn higgs mass parameter negative at low energy from positive value it obtained at high energy and thus cause eletroweak symmetry breaking .
a second example is provided by a class of @xmath7 models with the symmetry breaking sectors containing @xmath10 along with either a @xmath11 or a @xmath12 @xcite where flavor mixing is induced by vector - like fermions in the @xmath13 representation . in such models ,
the sm singlet from @xmath10 acquires a gut scale vev , breaking @xmath7 down to @xmath14 .
the @xmath10 also contains a sm higgs doublet which must have positive squared mass at the gut scale .
this positive mass term can turn negative at low energy due to renormalization group flow .
similar arguments can be applied for the case where a sm singlet of @xmath15-representation breaks @xmath7 down to sm @xcite .
the higgs doublet is also part of @xmath15 , which should have a positive squared mass at the gut scale .    in this paper
we explore the possibility of radiative electroweak symmetry breaking in several popular extensions of the sm .
the mechanism fails in the sm , as reviewed in sec .
[ sec : sm and type - i ] . in type - i seesaw models which includes right - handed neutrinos to the spectrum of sm
, radiative ewsb is not achieved - in fact the effect of @xmath16 fields is to provide _ positive _ corrections to the higgs mass parameter in evolving from high to low energies .
the situation is different in type - ii seesaw models which contain a weak scalar triplet , if the mass of the triplet is around the tev scale .
the key difference is the cross coupling between the sm higgs boson and an additional scalar field in the higgs potential .
the need for this additional scalar field to be at the tev scale arises from the needed magnitude of the @xmath17 parameter : @xmath18 . as the correction to @xmath17 from the scalar cross coupling grows as @xmath19 of the new scalar
, this scalar should not be much heavier than a tev , assuming that the quartic cross couplings are not extremely weak .
radiative mass generation is a popular mechanism for neutrino masses where one assumes new scalars at the tev scale for lepton number violation .
such models are amenable to radiative ewsb .
dark matter models employing a scalar singlet or an inert doublet also exhibit radiative ewsb .
finally , we propose and analyze a universal seesaw model wherein a new @xmath0 symmetry is broken at tev scale , which also shows radiative ewsb .    in non - supersymmtric extensions of the sm such as the ones studied here , the gauge hierarchy problem has to be somehow solved .
here we simply assume that this is done by fine - tuning .
the dimensionful parameters of the sm extensions are given by @xmath20 where the @xmath21 donate mass parameters for additional scalar fields that may be present and @xmath22 is the cosmological constant .
these dimensionful parameters may take , for reasons not understood , special values , rather than their  natural values " which are of order the planck scale .
once the scalar masses are set at these special ( or fine - tuned ) values , we assume that the corrections to @xmath17 arising from other particles present in the model do not exceed the physical mass of @xmath23 .
a positive mass parameter turning negative via rge flow leads to dimensional transmutation as can be seen in a coleman - weinberg @xcite analysis of the effective potential .
the rge evolution that we employ is in one to one correspondence with the effective potential , where the minimization is performed at a momentum scale close to the mass of the higgs scalar .
the outline of the paper is as follows . in sec .
[ sec : sm and type - i ] , we discuss the absence of such radiative ewsb in sm and type - i seesaw model . even though such a mechanism fails for type - i seesaw models , in sec .
[ sec : type - ii seesaw ] we show that the presence of a tev scale weak triplet makes radiative ewsb a success for the case of type - ii seesaw models . in sec .
[ sec : two - loop neutrino mass model ] , we show that for a two - loop neutrino mass model positive higgs mass parameter at a high energy scale turns negative in the renormalization group flow to low energy . in sec .
[ sec : idm ] and sec .
[ sec : scalar singlet dm model ] , we show that simple dark matter models such as inert doublet model and scalar singlet dark matter model also exhibit radiative ewsb .
when the models have tev scale scalar(s ) coupled to the sm higgs boson . in sec .
[ sec : universal seesaw ] , we study the radiative ewsb for a universal seesaw model . finally in sec .
[ sec : conclusion ] we conclude .
the mechanism of radiative ewsb occurs when the renormalization group flow of higgs mass parameter @xmath24 receives enough negative contribution from various parameters of the model turning the positive quantity into a negative one while evolving from high to low energies .
unfortunately such can not be the case in the sm .
the higgs potential of the sm is given by : @xmath25 and the renormalization group equation ( rge ) of mass parameter @xmath24 is given by @xcite : @xmath26 the evolution of the higgs mass parameter @xmath24 is dominated by the gauge couplings and the top quark yukawa couplings .
however , these corrections are proportional to @xmath17 itself , which implies that a positive @xmath17 can not turn into negative @xmath17 in rge evolution making radiative ewsb an impossibility within sm .    for type - i seesaw models ,
the additional lagrangian is given by : @xmath27 this part of the lagrangian manages to contribute in the renormalization group flow as it adds a new term to the rge of higgs mass parameter : @xmath28 unfortunately , the contributions coming from the @xmath16 fields make the situation worse as they only strengthen the positivity of the mass parameter as it evolves from high to low energies .
one should also notice that if we want to use the criterion of naturalness , ie .
the correction to the higgs mass parameter @xmath29 the scale of the majorana mass of the right - handed neutrino should not exceed @xmath30 @xcite .    in supersymmetry , the radiative ewsb becomes successful @xcite as the scalar partners of the fermions contribute significantly in the renormalization group flow of the higgs mass parameter in the right direction , making the positive term negative as it evolves from high to low energies .
similar incidents occur in other extensions of sm such as type - ii seesaw models where tev scale particle(s ) manages to dominate the renormalization group flow and turn the positive value of the higgs mass parameter into negative value triggering radiative ewsb .
a common feature of these extensions is the presence of new scalar(s ) at the tev scale as we show in the next five sections .
while type - i seesaw needs right - handed ( rh ) neutrinos which are neutral under the standard model gauge group with large majorana masses , the minimal type - ii seesaw mechanism requires the existence of a weak scalar triplet .
the most natural source for such triplets is provided by the left - right symmetric theories which can be realized either at low energy or can be embedded in grand unified theories ( guts ) such as @xmath7 or @xmath31 . here
we study the potential radiative ewsb scenario of the type - ii seesaw extension of the sm .
we consider the possibility that the weak scalar triplet @xmath32 is the only low - energy remnant of the new physics beyond the sm and the neutral component @xmath33 acquires a very small induced vev at low energy .
the sm electroweak doublet @xmath34 and the electroweak triplet @xmath35 are denoted by : @xmath36 where @xmath37 s are the pauli matrices .
the most general renormalizable tree - level scalar potential is @xmath38\\   & + \lambda_4 \ ; \phi^\dagger \phi \ ; \tr ( \delta^\dagger \delta)+ \lambda_5 \ ; \phi^\dagger [ \delta^\dagger , \delta ] \phi + \ { \dfrac{\mu}{\sqrt{2}}\ ; \phi^t\ ; i \sigma_2 \ ; \delta^\dagger \phi + h.c.\}. \end{split}\ ] ] the weak triplet also generates a majorana mass term for the neutrinos through the yukawa lagrangian : @xmath39 with the vev of the electroweak doublet @xmath40 , an effective dimension 5 operator generates the neutrino masses through the small but non - zero induced vev , @xmath41 and /or @xmath42 is so small that @xmath43 as the electroweak triplet decouples .
the neutrino mass matrix is given by : @xmath44 here one assumes @xmath45 so that @xmath46 is induced only after @xmath47 is generated .
one also needs to realize the fact that integrating out the heavy scalar triplet in the tree level approximation will also have an effect on the sm higgs quartic coupling . the effective quartic coupling below the scale @xmath48 is given by : @xmath49 this is the connecting formula for the standard model quartic coupling @xmath50 at the scale @xmath48 .
one needs to be careful while studying the solution to the set of rges of the parameters of a model .
the parameters in the scalar potential have to satisfy certain conditions at all energy scales which ensures that the potential is bounded from below . for
that purpose one must identify the necessary and sufficient conditions for boundedness of the potential . for type - ii seesaw models ,
the stability conditions for the potential to be bounded from below can be derived to be :    1 .    2 .    3 .
@xmath51    this is compact set of constraints that is necessary and sufficient as we show in the appendix [ appendix : stability cond ] . for a less compact set of constraints
the couplings of the lagrangian have to maintain these stability conditions upto the energy scale of new physics such as guts .
using vertex corrections and the wave function renormalization factors , we can calculate the complete set of @xmath52-functions and renormalization group equations ( rges ) .
we have also determined the rges for the mass parameters of the model which are related to the anomalous dimensions @xmath53 of the scalar masses by @xmath54 where @xmath55 and @xmath56 is the running scale .
the set of rges for the mass parameters is given by : and @xmath57 in the rge of the mass parameter @xmath42 given in eq .
( 17 ) of ref .
we also disagree with the coefficient in front of the mass parameter @xmath58 in the rge of @xmath59 given in eq .
( 18 ) of ref
] @xmath60 \mu^2_\phi + 6 \lambda_4 \mu_\delta^2 + 6|\mu|^2;\\ 16\pi^2 \dfrac{d \mu_\delta^2 } { d t }   = & \left[\l(-\frac{18}{5}g_1 ^ 2 - 12g_2 ^ 2\r)+8\lambda_2 + 2\lambda_3 + 2\tr\left(\mathbf{y_\delta^\dagger y_\delta}\right)\right]\mu_\delta^2 + 4\lambda_4 \mu_\phi^2 + 2|\mu|^2;\\ 16\pi^2 \dfrac{d \mu } { d t } = & \left[\lambda_1 + 4\lambda_4 - 8\lambda_5-\frac{27}{10}g_1 ^ 2-\frac{21}{2}g_2 ^ 2    +
2 t +   \tr\left(\mathbf{y_\delta^\dagger y_\delta}\right)\right]\mu;\\ \end{split}\ ] ] where @xmath61 .\ ] ]    a complete set of rges for all the couplings of the lagrangian is given in the appendix  [ appendix : rges for type - ii ] . with this set of rges
we proceed towards its numerical solution .
we can already see that the contribution from the cubic coupling @xmath42 and quartic couplings such as @xmath62 has the potential to drive the positive mass parameter @xmath17 at high energy scale towards a negative value at lower energy scale to trigger radiative ewsb . for a tev scale
scalar triplet mass one realizes that to get the correct order of neutrino mass the cubic parameter @xmath42 needs to be very small ( @xmath63  gev ) , which makes the contribution of @xmath42 term in the rges of mass parameters irrelevant .
this choice is natural , since quantum corrections to @xmath42 are proportional to @xmath42 itself - see eq .
( [ mass rges type - ii ] ) .      to analyze the evolution of the mass parameters , one needs to solve the set of rges which in turn requires one to define the relevant couplings at some energy scale @xcite . in this case , all the sm gauge couplings , yukawa couplings and higgs quartic couplings were evaluated at two - loop level upto the energy scale corresponding to the scalar triplet mass . at energies above the triplet mass the gauge couplings
were evolved continuously but with the updated rges given in the set of rges in eq .
( [ rges type - ii gauge ] ) .
the quartic coupling of the sm electroweak doublet higgs has a discontinuity at the triplet mass scale due to the matching condition of the parameter given in eq .
( [ quartic matching condition ] ) . above the energy scale @xmath64
, the full set of rges was used to evolve all the parameters of the model .
.quartic couplings and mass parameter values for the sample point used for the type - ii seesaw model in fig .
[ fig : sample rge type - ii seesaw ] . [ cols="^,^,^,^",options="header " , ]     .5 . in fig .
[ fig : sub1vf ] the evolution of the absolute value of the sm higgs mass parameter @xmath65 and @xmath66 along with the masses parameters @xmath67 and @xmath42 have been plotted . the point at which @xmath68 touches the horizontal line
corresponding to mass = @xmath69 , is the energy scale where radiative ewsb is triggered as the sign of the sm higgs mass - squared parameter @xmath70 switches from positive to negative while evolving from high to low energies . for the sample point , this happens around tev energy scale .
note that the mass parameter of the scalar @xmath71 also turns negative around the renormalization energy scale @xmath72 indicating that the @xmath73 symmetry also gets broken radiatively . in fig .
[ fig : sub2vf ] shows the evolution of all the quartic couplings of the model from planck scale down to weak scale .
note that the model remains perturbative all the way for the selected sample point.,title="fig:",width=302 ]    .5 . in fig .
[ fig : sub1vf ] the evolution of the absolute value of the sm higgs mass parameter @xmath65 and @xmath66 along with the masses parameters @xmath67 and @xmath42 have been plotted .
the point at which @xmath68 touches the horizontal line corresponding to mass = @xmath69 , is the energy scale where radiative ewsb is triggered as the sign of the sm higgs mass - squared parameter @xmath70 switches from positive to negative while evolving from high to low energies .
for the sample point , this happens around tev energy scale .
note that the mass parameter of the scalar @xmath71 also turns negative around the renormalization energy scale @xmath72 indicating that the @xmath73 symmetry also gets broken radiatively . in fig .
[ fig : sub2vf ] shows the evolution of all the quartic couplings of the model from planck scale down to weak scale .
note that the model remains perturbative all the way for the selected sample point.,title="fig:",width=321 ]      to find the solution of the set of rges , we took a more simplified case where we kept all the yukawa coupling @xmath74 to be small and negligible and yukawa coupling @xmath75 .
the numerical solution was hunted for the case where one of the eigenvalue of the scalar mass matrix @xmath76 corresponds to the sm higgs boson and another one corresponds to scalar boson of mass @xmath77 . here
, we used @xmath78 gev just as an example as similar universal seesaw model @xcite was used to explain apparent diphoton excess @xcite which eventually became statistically insignificant @xcite .
the vector - like fermion mass was kept around tev scale where the mass is approximated by @xmath79 .
one such sample point is given by table [ table : sample point vector - like fermion model ] .
as the first new particle in this model is at @xmath78 gev , the sm rges were evolved at two - loop level upto the scale @xmath80 and then the new set of rges was deployed to do the evolution of the couplings and mass parameters .
[ fig : sample rge vector - like model ] shows that the both the vevs ( electroweak vev and vev for the single @xmath71 ) can be generated by radiative correction .
we have presented in this paper various extensions of the sm where electroweak symmetry breaking is triggered by the renormalization group flow . even though such symmetry breaking fails to occur in sm ,
the scalar extensions are able to incorporate this attractive mechanism .
extensions like type - ii seesaw models , loop induced neutrino mass models and scalar dark matter models all have this built - in feature . a common shared feature of all models where radiative ewsb is realized is the presence of new scalars at the tev scale . these tev scale scalars may be detected in the large hadron collider ( lhc ) in the near future .
ksb would like to thank bartol research institute for hospitality during a visit where the work was initiated .
ksb and sk would like to thank the organizers of cetup * 2016 for hospitality and the participants for helpful discussions and comments . the work of ksb and sk is supported in part by the u. s. department of energy grant no .
de - sc0016013 and the work of ig is supported in part by bartol research institute .
the gauge transformation of the fields defined in eqs .
( [ eq : type - ii scalar potential])-([eq : type ii yukawa ] ) can be written as : @xmath81 where @xmath82 is a unitary matrix . under this gauge transformation
the yukawa term remains invariant : @xmath83    let us define @xmath84 : @xmath85 the gauge transformation of the field @xmath86 can be written as @xmath87 now , one can diagonalize @xmath86 by gauge transformation and as @xmath88 , we can write @xmath86 in such a basis as @xmath89    let us define the higgs doublet as @xmath90 in terms of these real fields the quartic part of the scalar potential ( given in eq .
( [ eq : type - ii scalar potential]))becomes @xmath91 where @xmath92 the quartic couplings of the potential form a vector space spanned by the real valued vector , @xmath93 and the quartic couplings of the scalar potential can be written as @xmath94 where @xmath95 we use known results from the copositivity conditions of real symmetric matrices @xcite to determine the boundedness conditions as    1 .    2 .    3 .      where @xmath96 as this set of boundedness conditions need to be satisfied for all values of @xmath52 , this set of conditions reduces to the set of inequalities ( [ cond1 ] , [ cond2 ] and [ cond3 ] )
note that , all the conditions mentioned in inequalities ( [ cond1 ] ) and ( [ cond2 ] ) and at least one of the conditions in inequality ( [ cond3 ] ) need to be satisfied for the potential to be bounded from below .
the rges for the yukawa couplings are given by @xcite : @xmath97 + \mathbf{y_d}\left [ t- \frac{1}{4 } g_1 ^ 2 - \frac{9}{4 } g_2 ^ 2 - 8\,g_3 ^ 2\right ] ; \\
16\pi^2 \dfrac{d\mathbf{y}_u}{d t }   = & \mathbf{y_u } \left [ \frac{3}{2 } \mathbf{y_u^\dagger y_u } - \frac{3}{2}\ , \mathbf{y_d^\dagger y_d}\right]+\mathbf{y_u}\left [ t-\frac{17}{20 } g_1 ^ 2 - \frac{9}{4 } g_2 ^ 2 - 8\,g_3 ^ 2\right ] ; \\ 16\pi^2 \dfrac{d \mathbf{y}_e}{d t}=&\mathbf{y_e}\left[\frac{3}{2}\mathbf{y_e^\dagger y_e}+\frac{3}{2}\mathbf{y_\delta^\dagger y_\delta}\right ] + \mathbf{y_e}\left[t-\frac{9}{4}g_1 ^ 2-\frac{9}{4}g_2 ^ 2\right];\\ 16\pi^2 \dfrac{d \mathbf{y}_\delta}{d t}=&\left[\frac{1}{2}\mathbf{y_e^\dagger y_e}+\frac{3}{2}\mathbf{y_\delta^\dagger y_\delta}\right]^t \mathbf{y_\delta } + \mathbf{y_\delta}\left[\frac{1}{2}\mathbf{y_e^\dagger   y_e}+\frac{3}{2}\mathbf{y_\delta^\dagger y_\delta}\right ] \\&+\left[-\frac{3}{2}\left(\frac{3}{5}g_1 ^ 2 + 3g_2 ^ 2\right)+\tr\left(\mathbf{y_\delta^\dagger   y_\delta}\right)\right]\mathbf{y_\delta}.\\ \end{split}\ ] ] the rges for the quartic couplings of the lagrangians are given by @xcite : @xmath98\lambda_4\nonumber + 8\lambda_5 ^ 2;\\ 16\pi^2 \dfrac{d \lambda_5}{d t}=&-\frac{9}{2}g_1 ^ 2\lambda_5-\frac{33}{2}g_2 ^ 2\lambda_5-\frac{18}{5}g_1 ^ 2g_2 ^ 2+\[2\lambda_2 - 2\lambda_3 + 2\lambda_1 + 8\lambda_4 + 2 t + 2 \tr\left(\mathbf{y_\delta^\dagger y_\delta}\right)\]\lambda_5.\nonumber\\ \end{split}\ ] ] @xmath99 ; \\ h & = \tr \[\mathbf{y_e^\dagger y_e y_e^\dagger y_e }   + 3 \mathbf{y_d^\dagger y_d y_d^\dagger y_d } + 3 \mathbf{y_u^\dagger y_u y_u^\dagger y_u } \ ] .
\end{split}\ ] ] the rges for the mass parameters are given by eq .
( [ mass rges type - ii ] ) .
+ beyond the energy scale corresponding to the mass of the triplet @xmath100 the sm gauge coupling evolution also needs to be recalculated due to the triplet s contribution .
while the weak triplet does not effect the evolution of the @xmath101 gauge coupling evolution , it does change the rges of the other gauge couplings .
the rges for the gauge couplings are given by : @xmath102 where @xmath103 are the three gauge couplings with the one loop @xmath52-function coefficient @xmath104 .
for the two - loop neutrino mass model , among the gauge couplings only the hypercharge gauge coupling is modified due the additional scalar particles .
so , the rges for the sm gauge couplings are given by : @xmath105 where @xmath103 are the three gauge couplings with the one loop @xmath52-function coefficient @xmath106 + the rges for the yukawa couplings are given by @xcite : @xmath107 the rges for the quartic scalar couplings are given by @xcite : @xmath108 - 8 \tr \l(\mathbf{f^\dagger f y^\dagger_\ell y_\ell}\r ) ;   \end{split}\ ] ] @xmath109 - 8 \tr \l(\mathbf{y^\dagger_\ell y_\ell h^\dagger h } \r ) ; \\ 16 \pi^2 \dfrac{d\lambda_6}{d t}=&\ ; 4 \lambda_2 \lambda_6 + 4
\lambda_3 \lambda_6 + 4 \lambda_4 \lambda_5 + 4 \lambda_6 ^ 2-\dfrac{90}{5 } \lambda_6 g_1 ^ 2 + \dfrac{432}{25 } g_1 ^ 4 + 8 \lambda_6 \ [ \tr ( \mathbf{f^\dagger f})+\tr ( \mathbf{h^\dagger h})\].\\ \end{split}\ ] ] the rges for the mass parameters are given by eq .
( [ mass rges two - loop neutrino ] ) .
the one loop rges for the inert doublet model have already been computed .
the sm gauge coupling rges are given by : @xmath110 where @xmath111 are the @xmath52-coefficients of the sm gauge couplings updated with the added particles .
+ the quark sector of the model remains unchanged , while the leptonic sector needs to be revisited .
the rges for the leptonic yukawa couplings are @xcite : @xmath112    for the quartic scalar coupling we find the following set of rges @xcite :
@xmath113 @xmath114;\\ \end{split}\ ] ] @xmath115 ; \hspace{.5 cm } t_{4\nu } \equiv   \tr \ [ \mathbf{h^\dagger h h^\dagger h}\];\hspace{.5 cm } t_{\nu e } \equiv   \tr \ [ \mathbf{h^\dagger h y_e^\dagger y_e}\].\\ \end{split}\ ] ]    the rges for the mass parameters are given by eq .
( [ mass rges for idm ] ) .
one notices from the set of rges that the evolution of majorana mass @xmath116 , the new yukawa coupling @xmath117 and the scalar quartic coupling @xmath57 are proportional to the respective quantities themselves .
the upshot of this setting is that these parameters remain small if they are small at any energy scale .
this feature of the model becomes self - explanatory upon realization that if any of these parameters becomes zero , the neutrino becomes massless and global @xmath0 symmetry conserving the lepton number is restored .
while the rges for the mass parameters are given by eq .
( [ mass rges sm + real singlet scalar ] ) , the rges for the quartic couplings are given by : @xmath118 rges for the yukawa couplings and the gauge couplings remain the same as sm .
the rges for the mass parameters are given by eq .
( [ mass rges for vector - like fermion model ] ) .
the rges for the gauge couplings are given by : @xmath119 where @xmath120 .
the set of rges for all the yukawa couplings is given by : @xmath121;\\ 16 \pi^2 \dfrac{d\mathbf{{y}_d}}{d t}=&\ ; \mathbf{y_d } \ [   \dfrac{3}{2 } \l ( \mathbf{y_d^\dagger y_d } - \mathbf{y_u^\dagger y_u } \r)+\dfrac{1}{2}\mathbf{g_d^\dagger g_d}+ t- \dfrac{1}{4}g_1 ^ 2-\dfrac{9}{4}g_2 ^ 2 - 8g_3 ^ 2 - 3g_4 ^ 2\];\\ 16 \pi^2 \dfrac{d\mathbf{{y}_e}}{d t}=&\ ; \mathbf{y_e } \ [   \dfrac{3}{2 }   \mathbf{y_e^\dagger y_e } + \dfrac{1}{2}\mathbf{g_e^\dagger g_e}+ t- \dfrac{9}{4}g_1 ^ 2-\dfrac{9}{4}g_2 ^ 2 - 3g_4 ^ 2\];\\ 16 \pi^2\dfrac{d\mathbf{{f}_u}}{d t}= & \ ; \mathbf{f_u } \[\mathbf{f_u^\dagger f_u}+t_f-\dfrac{8}{5}g_1 ^ 2 - 8g_3 ^ 2 - 15 g_4 ^ 2 \]+\dfrac{1}{2}\mathbf{g_u g_u^\dagger f_u};\\ 16 \pi^2 \dfrac{d\mathbf{{f}_d}}{d t}= & \ ; \mathbf{f_d } \[\mathbf{f_d^\dagger f_d}+t_f-\dfrac{2}{5}g_1 ^ 2 - 8g_3 ^ 2 - 15 g_4 ^ 2 \]+\dfrac{1}{2}\mathbf{g_d g_d^\dagger f_d};\\ 16 \pi^2 \dfrac{d\mathbf{{f}_e}}{d t}= & \ ; \mathbf{f_e } \[\mathbf{f_e^\dagger f_e}+t_f-\dfrac{18}{5}g_1 ^ 2 - 15 g_4
^ 2 \]+\dfrac{1}{2}\mathbf{g_e g_e^\dagger f_e};\\ 16 \pi^2 \dfrac{d\mathbf{{g}_u}}{d t}= & \ ; \mathbf{g_u } \[\mathbf{g_u^\dagger g_u}+\mathbf{y_u^\dagger y_u}+t_g-\dfrac{8}{5}g_1 ^ 2 - 8 g_3 ^ 2 - 6 g_4 ^ 2 \]+\dfrac{1}{2}\mathbf{f_u f_u^\dagger g_u};\\ 16 \pi^2 \dfrac{d\mathbf{{g}_d}}{d t}= & \ ; \mathbf{g_d } \[\mathbf{g_d^\dagger g_d}+\mathbf{y_d^\dagger y_d}+t_g-\dfrac{2}{5}g_1 ^ 2 - 8 g_3 ^ 2 - 6 g_4 ^ 2 \]+\dfrac{1}{2}\mathbf{f_d f_d^\dagger g_d};\\ 16 \pi^2\dfrac{d\mathbf{{g}_e}}{d t}= & \ ; \mathbf{g_e } \[\mathbf{g_e^\dagger g_e}+\mathbf{y_e^\dagger y_e}+t_g-\dfrac{18}{5}g_1 ^ 2 - 6 g_4 ^ 2 \]+\dfrac{1}{2}\mathbf{f_e f_e^\dagger g_e}.\\ \end{split}\ ] ] the rges for scalar quartic couplings are given by : @xmath122 @xmath123 where @xmath124 ; \\ t_f & = \tr \[\mathbf{f_e^\dagger f_e }   + 3 \mathbf{f_d^\dagger f_d } + 3 \mathbf{f_u^\dagger f_u } \ ] ; \\ t_g & = \tr \[\mathbf{g_e^\dagger g_e }   + 3 \mathbf{g_d^\dagger g_d } + 3 \mathbf{g_u^\dagger g_u } \ ] ; \\ h & = \tr \[\mathbf{y_e^\dagger y_e y_e^\dagger y_e }   + 3 \mathbf{y_d^\dagger y_d y_d^\dagger y_d } + 3 \mathbf{y_u^\dagger y_u y_u^\dagger y_u } \];\\
h_f & = \tr \[\mathbf{f_e^\dagger f_e f_e^\dagger f_e }   + 3 \mathbf{f_d^\dagger f_d f_d^\dagger f_d } + 3 \mathbf{f_u^\dagger f_u f_u^\dagger f_u } \];\\
h_g & = \tr \[\mathbf{g_e^\dagger g_e g_e^\dagger g_e }   + 3 \mathbf{g_d^\dagger g_d g_d^\dagger g_d } + 3 \mathbf{g_u^\dagger g_u g_u^\dagger g_u } \];\\ h_{yg } & = \tr \[\mathbf{y_e^\dagger y_e g_e^\dagger g_e }   + 3 \mathbf{y_d^\dagger y_d g_d^\dagger g_d } + 3 \mathbf{y_u^\dagger y_u g_u^\dagger g_u } \];\\ h_{fg } & = \tr \[\mathbf{f_e^\dagger f_e g_e^\dagger g_e }   + 3 \mathbf{f_d^\dagger f_d g_d^\dagger g_d } + 3 \mathbf{f_u^\dagger f_u g_u^\dagger g_u } \ ] .
\end{split}\ ] ]    99 s.  chatrchyan _ et al . _
[ cms collaboration ] , phys .
b * 716 * , 30 ( 2012 ) [ arxiv:1207.7235 [ hep - ex ] ] ; g.  aad _ et al .
_ [ atlas collaboration ] , phys .
b * 716 * , 1 ( 2012 ) [ arxiv:1207.7214 [ hep - ex ] ] .
g.  aad _ et al . _ [ atlas and cms collaborations ] , phys .
lett .   * 114 * , 191803 ( 2015 ) [ arxiv:1503.07589 [ hep - ex ] ] .
p.  minkowski , phys .
b * 67 * , 421 ( 1977 ) ; t.  yanagida , proceedings of the _ workshop on unified theories and baryon number in the universe _ , tsukuba , 1979 , eds .
a.  sawada , a.  sugamoto ; s.  glashow , in cargese 1979 , _ proceedings , quarks and leptons _
( 1979 ) ; m.  gell - mann , p.  ramond and r.  slansky , conf .
c * 790927 * , 315 ( 1979 ) [ arxiv:1306.4669 [ hep - th ] ] ; r.  n.  mohapatra and g.  senjanovic , phys .
lett .   * 44 * , 912 ( 1980 ) .
j.  schechter and j.  w.  f.  valle , phys .
d * 22 * , 2227 ( 1980 ) ; g.  lazarides , q.  shafi and c.  wetterich , nucl .  phys .
b * 181 * , 287 ( 1981 ) ; r.  n.  mohapatra and g.  senjanovic , phys .
d * 23 * , 165 ( 1981 ) ; c.  wetterich , nucl .
b * 187 * , 343 ( 1981 ) ; j.  schechter and j.  w.  f.  valle , phys .
d * 25 * , 774 ( 1982 ) .
a.  zee , phys .
b * 93 * , 389 ( 1980 ) erratum : [ phys .  lett .
b * 95 * , 461 ( 1980 ) ] .
k.  s.  babu , phys .
b * 203 * ( 1988 ) 132 .
g.  r.  farrar and p.  fayet , phys .
b * 76 * , 575 ( 1978 ) .
h.  c.  cheng , k.  t.  matchev and m.  schmaltz , phys .
d * 66 * , 036005 ( 2002 ) [ hep - ph/0204342 ] . h.  c.  cheng and i.  low , jhep * 0309 * , 051 ( 2003 ) [ hep - ph/0308199 ] ; h.  c.  cheng and i.  low , jhep * 0408 * , 061 ( 2004 ) [ hep - ph/0405243 ] ; i.  low , jhep * 0410 * , 067 ( 2004 ) [ hep - ph/0409025 ] . n.  g.  deshpande and e.  ma , phys .  rev .
d * 18 * , 2574 ( 1978 ) .
e.  ma , phys .
d * 73 * , 077301 ( 2006 ) [ hep - ph/0601225 ] .
t.  w.  b.  kibble , g.  lazarides and q.  shafi , phys .
b * 113 * , 237 ( 1982 ) ; r.  n.  mohapatra , phys .  rev .
d * 34 * , 3457 ( 1986 ) ; a.  font , l.  e.  ibanez and f.  quevedo , phys .
b * 228 * , 79 ( 1989 ) ; l.  m.  krauss and f.  wilczek , phys .
.  lett .   * 62 * , 1221 ( 1989 ) ; l.  e.  ibanez and g.  g.  ross , phys .
b * 260 * , 291 ( 1991 ) ; l.  e.  ibanez and g.  g.  ross , nucl .
b * 368 * , 3 ( 1992 ) ; s.  p.  martin , phys .
d * 46 * , 2769 ( 1992 ) [ hep - ph/9207218 ] .
l.  e.  ibanez and g.  g.  ross , phys .
b * 110 * , 215 ( 1982 ) ; k.  inoue , a.  kakuto , h.  komatsu and s.  takeshita , prog .
phys .   *
68 * , 927 ( 1982 ) erratum : [ prog .
phys .   * 70 * , 330 ( 1983 ) ] ; l.  e.  ibanez , phys .
b * 118 * , 73 ( 1982 ) ; j.  r.  ellis , d.  v.  nanopoulos and k.  tamvakis , phys .
b * 121 * , 123 ( 1983 ) ; j.  r.  ellis , j.  s.  hagelin , d.  v.  nanopoulos and k.  tamvakis , phys .  lett .
b * 125 * , 275 ( 1983 ) ; l.  alvarez - gaume , j.  polchinski and m.  b.  wise , nucl .
b * 221 * , 495 ( 1983 ) .
a.  de  rujula , h.  georgi , and s.  l.  glashow , in _ fifth workshop on grand unification _ , edited by k.  kang , h.  fried , and p.   frampton ( world scientific , singapore , 1984 ) , p. 88 .
for an earlier discussion , see y.  achiman and b.  stech , in _ new phenomena in lepton - hadron physics _ , edited by d. e. c. fries and j. wess ( plenum , new york , 1979 ) , p. 303 .
k.  s.  babu , x.  g.  he and s.  pakvasa , phys .
d * 33 * , 763 ( 1986 ) ; g.  m.  pelaggi , a.  strumia and s.  vignali , jhep * 1508 * , 130 ( 2015 ) [ arxiv:1507.06848 [ hep - ph ] ] . k.  s.  babu , b.  bajc and s.  saad , arxiv:1605.05116 [ hep - ph ] . k.  s.  babu , i.  gogoladze , p.  nath and r.  m.  syed , phys .
d * 72 * , 095011 ( 2005 ) [ hep - ph/0506312 ] .
s.  r.  coleman and e.  j.  weinberg , phys .
d * 7 * , 1888 ( 1973 ) .
h.  arason , d.  j.  castano , b.  keszthelyi , s.  mikaelian , e.  j.  piard , p.  ramond and b.  d.  wright , phys .
d * 46 * , 3945 ( 1992 ) ; m.  x.  luo and y.  xiao , phys .
lett .   * 90 * , 011601 ( 2003 ) [ hep - ph/0207271 ] .
f.  vissani , phys .
d * 57 * , 7027 ( 1998 ) [ hep - ph/9709409 ] .
a.  arhrib , r.  benbrik , m.  chabab , g.  moultaka , m.  c.  peyranere , l.  rahili and j.  ramadan , phys .
d * 84 * , 095005 ( 2011 ) [ arxiv:1105.1925 [ hep - ph ] ] .
m.  a.  schmidt , phys .
d * 76 * , 073010 ( 2007 ) erratum : [ phys .
d * 85 * , 099903 ( 2012 ) ] [ arxiv:0705.3841 [ hep - ph ] ] .
i.  gogoladze , n.  okada and q.  shafi , phys .
d * 78 * , 085005 ( 2008 ) [ arxiv:0802.3257 [ hep - ph ] ] ; a.  kobakhidze and a.  spencer - smith , jhep * 1308 * , 036 ( 2013 ) [ arxiv:1305.7283 [ hep - ph ] ] .
k.  s.  babu and j.  julio , aip conf .
proc .   * 1604 * , 134 ( 2014 ) .
k.  s.  babu and c.  macesanu , phys .
d * 67 * , 073010 ( 2003 ) [ hep - ph/0212058 ] .
m.  nebot , j.  f.  oliver , d.  palao and a.  santamaria , phys .
d * 77 * , 093013 ( 2008 ) [ arxiv:0711.0483 [ hep - ph ] ] ; j.  herrero - garcia , m.  nebot , n.  rius and a.  santamaria , nucl .
b * 885 * , 542 ( 2014 ) [ arxiv:1402.4491 [ hep - ph ] ] ; d.  schmidt , t.  schwetz and h.  zhang , nucl .
b * 885 * , 524 ( 2014 ) [ arxiv:1402.2251 [ hep - ph ] ] ; j.  herrero - garcia , m.  nebot , n.  rius and a.  santamaria , arxiv:1410.2299 [ hep - ph ] .
r.  barbieri , l.  j.  hall and v.  s.  rychkov , phys .
d * 74 * , 015007 ( 2006 ) [ hep - ph/0603188 ] ; j.  a.  casas , j.  r.  espinosa and i.  hidalgo , nucl .
b * 777 * , 226 ( 2007 ) [ hep - ph/0607279 ] .
e.  ma , mod .
a * 21 * , 1777 ( 2006 ) [ hep - ph/0605180 ] .
t.  hambye and m.  h.  g.  tytgat , phys .
b * 659 * , 651 ( 2008 ) [ arxiv:0707.0633 [ hep - ph ] ] .
a.  merle and m.  platscher , phys .
d * 92 * , no .
9 , 095002 ( 2015 ) [ arxiv:1502.03098 [ hep - ph ] ] .
v.  silveira and a.  zee , phys .
* 161b * , 136 ( 1985 ) . j.  mcdonald , phys .  rev .
d * 50 * , 3637 ( 1994 ) [ hep - ph/0702143 [ hep - ph ] ] .
m.  kadastik , k.  kannike and m.  raidal , phys .
d * 80 * , 085020 ( 2009 ) erratum : [ phys .
d * 81 * , 029903 ( 2010 ) ] [ arxiv:0907.1894 [ hep - ph ] ] ; y.  mambrini , n.  nagata , k.  a.  olive and j.  zheng , phys .
d * 93 * , no .
11 , 111703 ( 2016 ) [ arxiv:1602.05583 [ hep - ph ] ] .
s.  m.  barr , phys .
lett .   * 53 * , 329 ( 1984 ) ; a.  e.  nelson , phys .
* 136b * , 387 ( 1984 ) .
k.  s.  babu and r.  n.  mohapatra , phys .
lett .   * 62 * , 1079 ( 1989 ) .
k.  s.  babu and r.  n.  mohapatra , phys .
d * 41 * , 1286 ( 1990 ) .
r.  n.  mohapatra and y.  zhang , jhep * 1406 * , 072 ( 2014 ) [ arxiv:1401.6701 [ hep - ph ] ] .
s.  b.  dev , r.  n.  mohapatra and y.  zhang , jhep * 1602 * , 186 ( 2016 ) [ arxiv:1512.08507 [ hep - ph ] ] .
cms collaboration [ cms collaboration ] , collisions at 13tev ,  cms - pas - exo-15 - 004 ; the atlas collaboration , atlas - conf-2015 - 081 ; m.  aaboud _ et al .
_ [ atlas collaboration ] , arxiv:1606.03833 [ hep - ex ] ; v.  khachatryan _ et al . _ [ cms collaboration ] , arxiv:1606.04093 [ hep - ex ] .
cms collaboration [ cms collaboration ] , cms - pas - exo-16 - 027 ; the atlas collaboration [ atlas collaboration ] , atlas - conf-2016 - 059 . k.  p.
hadeler , linear  algebra  appl .
,   * 49 * , 79 ( 1983 ) .
k.  g.  klimenko , theor .
phys . ,
* 62 * , 58 ( 1985 ) .