materials with correlated electrons exhibit some of the most intriguing quantum states in condensed matter physics @xcite . since the number of electric charge carriers essentially determines such states , external electric or light field can be applied for controllable manipulations .
stability of electric field or light induced states has been demonstrated in some novel systems , where switching occurs between neighboring equilibrium thermodynamic states @xcite .
this powerful characteristic can be applied in electric devices such as transistors and memories , which are of great importance to not only the fundamental physics research but also the information processing technology @xcite . correlated materials with a rich set of quantum states delicately balanced on a similar energy scale
will be promising platforms to realize such devices .
( a ) crystal structure of @xmath0-tas@xmath1 .
( b ) top view of ta - ta plane of @xmath0-tas@xmath1 .
the unit cells of conventional @xmath0 structure and low temperature ccdw structure are denoted with green solid and dashed lines , respectively .
the black hexagrams are the so - called `` david star '' pattern in ccdw phase .
( c ) schematic diagrams of `` david star '' ( left ) , ccdw ( middle ) , and nccdw ( right ) in @xmath0-tas@xmath1 . ]    the transition metal dichalcogenide ( tmd ) @xmath0-tas@xmath1 is one of the promising candidate due to the multiple competing ground states in it @xcite . as shown in fig .
[ fig_1_structure ] ( a ) , @xmath0-tas@xmath1 shows a cdi@xmath1-type layered crystal structure with ta atoms octahedrally coordinated by s atoms .
a unit layer consists of one ta layer sandwiched between two s layers . at low temperatures ,
strong @xmath4-dependent electron - phonon coupling induced periodic lattice distortion makes a @xmath5 superlattice @xcite , in which ta atoms displace to make `` david star '' clusters ( fig .
[ fig_1_structure ] ( b ) ) .
the outer twelve atoms within each star move slightly towards the atom at the center , leading to the commensurate charge - density - wave ( ccdw ) ground state .
in particular , in ccdw state , the correlation effect of @xmath6 electrons of ta atoms turns the system into mott insulating state @xcite . upon heating to 225 k , it undergoes a sequence of first order phase transition to a nearly commensurate ( nc ) cdw . the nccdw phase is composed of metallic incommensurate ( ic ) network and mott insulating ccdw domains .
the ccdw domains shrink upon heating and finally disappear at 355 k , while the system transforms to iccdw state .
the standard metallic @xmath0 structure appears above 535 k. moreover , when ccdw state is suppressed , superconductivity emerges in this system @xcite .
one can expect that the controllable switching between those states will be helpful for figuring out the mechanism of cdw and superconductivity , and realizing the high performance memory and transistor in future technology . to meet this goal , many groups performed investigations on this novel material .
yu et al . reported a gate - controlled li ion intercalation can suppress ccdw and introduce superconducivity @xcite .
tsen et al .
@xcite , hollander et al .
@xcite , yoshida et al .
@xcite , and mihailovic et al .
@xcite reported the in - plane electric field induced transition from ccdw state to nccdw or some metastable hidden state .
cho et al . applied perpendicular electric pulse on @xmath0-tas@xmath1 and found positive electric pulse can introduce ic network to suppress mott state at low temperature @xcite .
moreover , zhang et al .
@xcite , mihailovic et al .
@xcite , and han et al .
@xcite suggested that light can introduce a transition from ccdw state to nccdw or hidden state as well .
in all those manipulations , the transitions are considered to be determined by the charge carrier doping .
however , the mechanisms of such transitions are not fully clear yet .    in this work , we simulated by first - principles calculations the charge carrier doping effect on cdw in @xmath0-tas@xmath1 .
we found that ccdw is stable upon electron doping , while hole doping significantly suppresses ccdw instability , implying different mechanisms of recently reported electric and photoelectric manipulations of cdw in @xmath0-tas@xmath1 .
we figured out such mechanism by analysis of carrier doping effects .
furthermore , we show that superconductivity with @xmath2 about @xmath3 k can be introduced by hole doping in the system .
the density functional theory ( dft ) calculations were carried out using quantum espresso package @xcite with ultrasoft pseudopotentials .
the exchange - correlation interaction was treated with the generalized gradient approximation ( gga ) with pw91 parametrization @xcite .
the energy cutoff for the plane - wave basis set was 35 ry .
the marzari - vanderbilt fermi smearing method @xcite with a smearing parameter of @xmath7 ry was used for the calculations of the total energy and electron charge density . to simulate the monolayer ,
a vacuum layer more than 15 @xmath8 was introduced . for the bulk sample
, brillouin zone sampling is performed on the monkhorst - pack ( mp ) mesh @xcite of @xmath9 , while a denser @xmath10 grid was used in the electron phonon coupling calculations .
phonon dispersions were calculated using density functional perturbation theory ( dfpt ) @xcite with an @xmath11 mesh of @xmath12-points . for the monolayer sample , @xmath13-points grids of @xmath14 and @xmath15 and @xmath12-points grids of @xmath16 are used . since the frequency and stability of soft phonon modes were found to be very sensitive to the plane - wave cutoff for the charge density , a large cutoff of 1500 ry was used for the charge density in the phonon calculations @xcite .
all the calculation parameters were well tested . in the investigation of carrier doping effects , the electron / hole doping was simulated by increasing / decreasing the total electron numbers of the system , together with a compensating uniform positive / negative background to maintain the charge neutrality .
the crystal structure of each doped sample was optimized with respect to lattice parameters and atomic positions . for comparisons , we also calculated some doped samples by using the lattice parameters of pristine tas@xmath1 and relaxing the atomic positions , as described later .
the low symmetry cdw structure is usually considered as the high symmetry phase with distortion introduced by some instability . therefore , we firstly investigated carrier doping effects in the bulk @xmath0-tas@xmath1 .
the electronic structures and phonon properties for the samples with different doping level were calculated . as examples , figure [ fig_2_bulk ]
shows the calculated band structures and fermi surfaces of the pristine bulk @xmath0-tas@xmath1 , the doped bulk @xmath0-tas@xmath1 with the doping level of @xmath17 electrons / f.u . , and the doped bulk @xmath0-tas@xmath1 with @xmath17 holes / f.u . as expected , electron doping increases the lattice parameters and hole doping shrinks the lattice ( figs .
[ fig_2_bulk ] ( j ) and ( k ) ) . for the pristine tas@xmath1 ,
our results are of good agreement with the previous calculations @xcite .
there is a gap of @xmath18 ev below fermi energy ( @xmath19 ) ( fig .
[ fig_2_bulk ] ( b ) ) . in the @xmath20-@xmath21 direction ,
the two bands around the gap are nearly flat due to the quasi-2d nature of the layered structure .
the band above the gap crosses fermi energy ( @xmath19 ) , forming a 2d electron pocket around @xmath22 point ( figs .
[ fig_2_bulk ] ( b ) and ( e ) ) .
such gap increases upon electron doping and decreases upon hole doping ( figs .
[ fig_2_bulk ] ( a ) and ( c ) ) .
electron doping increases @xmath19 . as shown in figs .
[ fig_2_bulk ] ( a ) and ( d ) , for the doped bulk @xmath0-tas@xmath1 with @xmath17 electrons / f.u . the rise of @xmath19 expands and opens up the original electron pockets around @xmath22 , leaving hole pockets centered at @xmath23 point . besides the original band crossing @xmath19 , a band with higher energy starts to cross @xmath19 , forming a 2d cylinder - like electron pocket around the zone center . on the contrary
, the hole doping enhances the dispersion in the @xmath20-@xmath21 direction and weakens the quasi-2d nature by shrinking the lattice and decreasing the interlayer distance . as shown in figs .
[ fig_2_bulk ] ( c ) and ( f ) , for the doped bulk @xmath0-tas@xmath1 with @xmath17 holes / f.u .
, the bands close to @xmath19 are not flat anymore .
hole doping reduces @xmath19 .
the decrease of @xmath19 shrinks the original electron pockets around m point .
moreover , a lower energy band starts to cross @xmath19 , forming a 3d hole pocket around @xmath20 point .
the phonon instability of the high symmetry structure is considered to be directly related to the cdw distortion : at high temperatures , the phonon of the high symmetry structure softens at cdw vector ( @xmath24 ) . above the transition temperature the phonon frequency near @xmath24 drops but does not go to zero . just below the transition temperature
the phonon frequency near @xmath24 is imaginary , meaning there is a restructuring of lattice with a superlattice vector of @xmath24 @xcite .
therefore , phonon dispersion without imaginary frequency implies that the structure is stable compared to cdw structure .
the phonon calculation is proved to be an effective method to simulate the cdw instability : the calculated phonon dispersions show instability just locating at @xmath24 of some tmds @xcite .
more specially , for the present @xmath0-tas@xmath1 , experimental report showed cdw can be suppressed by pressure @xcite , which is correctly simulated by liu s phonon calculation @xcite . in this work
we also performed the phonon calculations on each sample .
for the pristine @xmath0-tas@xmath1 , our calculation is of good agreement with liu s calculation ( fig .
[ fig_2_bulk ] ( h ) ) .
the phonon dispersion show instability very close to the ccdw vector @xmath25 .
this instability persists at all values of @xmath26 , as shown in inset of fig .
[ fig_2_bulk ] ( h ) . for the electron doped sample ,
the acoustic branches become more unstable . as shown in fig .
[ fig_2_bulk ] ( g ) , the unstable modes in the doped bulk @xmath0-tas@xmath1 with with @xmath17 electrons / f.u .
appear in @xmath23-@xmath22 and @xmath23-@xmath20 directions , indicating the area of instability are largely expanded .
on the contrary , the hole doping significantly stabilizes lattice . as shown in fig .
[ fig_2_bulk ] ( i ) , no unstable mode can be found in the phonon dispersion of the doped bulk @xmath0-tas@xmath1 with @xmath17 holes / f.u .
since the lowest mode in the pristine tas@xmath1 locates near @xmath27 , we used such mode as an indicator of doping effect on ccdw in tas@xmath1 .
the frequency variation of such mode under doping is shown in fig .
[ fig_2_bulk ] ( l ) .
one can note that upon electron doping the mode is always unstable in bulk @xmath0-tas@xmath1 , while hole doping significantly suppresses the instability . according to our calculation ,
the lattice becomes completely stable when the doping level is higher than @xmath28 electrons / f.u .
since hole doping shrinks the lattice as the high pressure does , we should figure out the suppression of ccdw under hole doping is majorly attributed to whether the lattice compression or to the pure doping effect .
we calculated a hypothetic doped @xmath0-tas@xmath1 with n=0.3 holes / f.u . ,
in which the lattice parameters are fixed to those of the undoped pristine @xmath0-tas@xmath1 . as shown in fig .
[ fig_2_bulk ] ( l ) , ccdw can be suppressed by just doping holes into @xmath0-tas@xmath1 without changing its lattice volume .
that demonstrates the suppression of ccdw in the present case is predominately by the hole doping effect .
furthermore , we investigated the doping induced energy and the distortion variations in the ccdw state for the pristine tas@xmath1 , the doped tas@xmath1 with @xmath17 electrons / f.u . , and the doped tas@xmath1 with @xmath17 holes / f.u .
we defined the ccdw formation energy @xmath29 as @xmath30 where @xmath31 and @xmath32 are the total energies of the relaxed ccdw structure and the @xmath0 structure .
the distortion rations @xmath33 can be expressed as    @xmath34    and @xmath35 where @xmath36 is the in - plane lattice parameter of the undistorted @xmath0 structure , @xmath37 and @xmath37 are radius of the inner and outer circles of `` david star '' , as shown in fig . 3 ( b ) . for @xmath0-tas@xmath1 with ccdw structure ,
the layers stacking order is not clear yet .
therefore , we simply construct the @xmath38 supercell with 39 atoms to simulate the ccdw structure .
we relaxed the atomic positions using the lattice parameters @xmath39 and @xmath40 , where @xmath36 and @xmath40 are from figs .
[ fig_2_bulk ] ( j ) and ( k ) , respectively . in the relaxed ccdw structures
, we found that for the pristine tas@xmath1 and the doped tas@xmath1 with @xmath17 electrons / f.u .
[ fig_3_charge ] ( b ) and ( a ) ) , the distortion ratios are large , while the formation energies of those are negative . on the other hand , the distortion ratios of the doped tas@xmath1 with @xmath17 holes / f.u .
are nearly zero ( fig .
[ fig_3_charge ] ( c ) ) , and the energy difference between the relaxed supercell and the @xmath0 unit cell is very small ( @xmath41 mev ) , i.e. the ccdw structure relaxes to high symmetry @xmath0 structure .
we also relaxed the atomic positions with @xmath36 and @xmath40 fixed to those in the pristine tas@xmath42 to prevent the volume variation . with the lattice parameters of the pristine tas@xmath42 , the doped tas@xmath42 with @xmath17 electrons / f.u
can be seen as under pressure . however , the ccdw is still stable . on the other hand , without doping introduced volume shrinking , in the doped tas@xmath42 with @xmath17 holes / f.u .
the ccdw is still not favorable in energy .
such results of the ccdw supercell coincide well with our deductions based on the analysis of the phonon instabilities of the @xmath0 unit cell .
if we simply consider that the roles of electron / hole doping on the electronic structures are to increase / decrease @xmath19 , we can estimate the population of the added / removed electrons when electron / hole doping and understand the calculated results more clearly .
figure [ fig_3_charge ] ( d ) shows the density of states ( dos ) near @xmath19 of the pristine tas@xmath42 in the ccdw structure calculated directly by dft .
based on the integrations of the dos upwards / downwards from @xmath19 , we can see that if one electron is doped into a `` david star '' ( the doping level is equal to @xmath43 electron / f.u . ) , it will be mainly added into the center of the `` david star '' .
obviously , it will enhance the clustering of the charge density .
therefore , ccdw is stable upon electron doping . on the contrary ,
doping one hole into a `` david star '' will notably decrease the charge density at the center and inner atoms of the `` david star '' , which will weaken the charge density clustering . therefore , ccdw is suppressed upon hole doping .
[ cols="^,^,^,^,^,^,^",options="header " , ]     the doping effects in the bulk tas@xmath1 investigated using the low symmetry ccdw structure .
( a ) , ( b ) , and ( c ) are the charge density in ta - ta plane of the doped tas@xmath1 with @xmath17 electrons / f.u . , pristine tas@xmath1 , and the doped tas@xmath1 with @xmath17 holes / f.u .
, respectively .
( d ) dos of the pristine tas@xmath1 directly calculated in ccdw structure .
the colored areas below and above @xmath19 denote the electron density which can be integrated to one electron .
the related electron charge densities are plotted in the insets of ( d ) . ]
the doping effects in monolayer @xmath0-tas@xmath1 . for the doped
tas@xmath1 with @xmath17 electrons / f.u . , pristine tas@xmath1 , and the doped tas@xmath1 with @xmath17 holes / f.u . , ( a ) , ( b ) , and ( c ) are the band dispersions of those samples , ( d ) , ( e ) , and ( f ) are the fermi surfaces , and ( g ) , ( h ) , and ( i ) are the phonon dispersions , respectively .
the red dashed lines show the phonon dispersion of undoped monolayer @xmath0-tas@xmath1 with tensile strain of 5% . ]
schematic picture of suppression of ccdw by hole doping .
( a ) the situation of tas@xmath1 in the ccdw state when holes are doped locally .
the solid orange `` david stars '' are the area under ccdw distortion .
the hollow orange `` david stars '' are the area in which only one hole per a `` david star '' is doped , which are still under ccdw distortion .
the dashed light orange `` david stars '' are the area with more doped holes , in which the ccdw distortion is fully suppressed .
the black arrows denote the diffusion of the doped holes .
( b ) the schematic band structure of tas@xmath1 with ccdw distortion .
( c ) the schematic band structure of tas@xmath1 with ccdw distortion when one hole per `` david star '' is doped .
( d ) the schematic band structure of tas@xmath1 in @xmath0 structure . ]        we also investigate the doping effect in monolayer @xmath0-tas@xmath1 . as examples , fig .
[ fig_4_monolayer ] shows the calculated electronic structures and phonon properties of the monolayer pristine tas@xmath1 , the doped tas@xmath1 with @xmath17 electrons / f.u . ,
and the doped tas@xmath1 with @xmath17 holes / f.u in @xmath0 structure . for the undoped pristine monolayer @xmath0-tas@xmath1 ,
the band structure is similar to that of pristine bulk sample .
one band crosses @xmath19 , forming a 2d electron pocket around @xmath22 point ( figs .
[ fig_4_monolayer ] ( b ) and ( e ) ) .
there is a gap of  0.8 ev below the band crossing @xmath19 . in the electrons and holes doped monolayer samples , the gap slightly changes , which is different from the case of bulk samples .
for the doped monolayer @xmath0-tas@xmath1 with @xmath17 electrons / f.u . , the band structure is slightly different from that of the doped bulk @xmath0-tas@xmath1 with @xmath17 electrons / f.u .
[ fig_4_monolayer ] ( a ) ) .
although @xmath19 changes with electron doping , there is still only one band crossing @xmath19 .
the 2d cylinder - like electron pocket in the doped bulk @xmath0-tas@xmath1 with @xmath17 electrons / f.u .
does not appear in monolayer sample ( fig .
[ fig_4_monolayer ] ( d ) ) .
for the doped monolayer @xmath0-tas@xmath1 with @xmath17 holes / f.u . , there is only one band crossing @xmath19 ( fig .
[ fig_4_monolayer ] ( c ) ) .
similar to the case of bulk samples , the phonon calculation shows that the cdw instability in the monolayer @xmath0-tas@xmath1 can not be suppressed by electron doping ( fig .
[ fig_4_monolayer ] ( g ) ) , while it can be suppressed under hole doping ( fig .
[ fig_4_monolayer ] ( i ) ) .
the result indicates the suppression of ccdw in @xmath0-tas@xmath1 is not due to the hole doping enhanced band dispersion along @xmath20-@xmath21 direction ( figs .
[ fig_2_bulk ] ( c ) and ( f ) ) .
the suppression should be attributed to the weakening of the electron - phonon coupling at @xmath24 upon hole doping .    besides
, one may note in the monolayer sample , a small instability near @xmath20 can be found in the phonon dispersion ( fig .
[ fig_4_monolayer ] ( h ) ) . in a very recent calculation of phonon dispersion of monolayer @xmath0-tas@xmath1 by zhang et al .
@xcite , there is a similar instability near @xmath20 .
such instability is consistent with the instability against long - wavelength transversal waves @xcite .
this instability is suggest to be fixed by defects , such as ripples or grain boundaries , which do not allow these waves by limiting the size @xcite .
furthermore , we found tensile strain can suppress such instability near @xmath20 , which can be easily applied to monolayer materials @xcite .
as an example , in fig .
[ fig_4_monolayer ] ( h ) , the red dashed lines show the phonon dispersion of @xmath0-tas@xmath1 under a tensile strain of 5% , in which the instability near @xmath20 is suppressed .
the observation indicates that the application of tensile strain is helpful for stabilizing the experimentally exfoliated monolayer or fewlayer @xmath0-tas@xmath42 .
based on our calculations of phonon properties in @xmath0 structure , ccdw formation energies and distortion ratios in ccdw structure , we can conclude that when the doping level is above @xmath28 holes / f.u .
( 2.6 holes per a `` david stars '' ) , ccdw can be completely suppressed in @xmath0-tas@xmath1 for bulk and monolayer , i.e. all the `` david stars '' should be melted .
experimentally , holes can directly be introduced by a positive electric field perpendicular to the sample @xcite , or by light @xcite .
one can expect that in reality the doped holes will firstly distribute in a local area and gradually diffuse out , i.e. even the doping level might be much lower than 0.2 holes / f.u .
, in the local area , the doping level could be high enough to melt the `` david stars '' in this area and destroy the long - range ccdw , as shown in fig .
[ fig_5_hole_mechanism ] ( a ) .
we describe the possible picture of such process here : in @xmath0 structure , the electronic structure near @xmath19 are formed predominantly from a single ta @xmath44 band , which splits into subbands by the formation of ccdw state .
six of these subbands are fulfilled with 12 electrons per new ccdw unit cell , forming a manifold of occupied states .
the 13th leftover electron is localized on the central ta atom of the `` david stars '' , forming a half - filled subband at @xmath19 .
this half - filled subband further splits by the coulomb interaction into upper and lower hubbard bands ( uhb and lhb ) , as shown in fig .
[ fig_5_hole_mechanism ] ( b ) .
light or positive perpendicular electric filed can firstly excite one electron from lhb to uhb , and create one hole in lhb ( fig .
[ fig_5_hole_mechanism ] ( c ) ) . in real space ,
the hole doped by light or positive perpendicular electric filed locates at the center of the `` david stars '' , leaving a polaron with an excess charge .
the other 12 electrons in this polaron are still star - shaped around the central , thus screening the excess charge .
when more holes are doped into the `` david stars '' , the ccdw distortion is suppressed locally , the structure transforms into @xmath0 in the doping area .
the splitting subbands merge to a single band again ( fig .
[ fig_5_hole_mechanism ] ( d ) ) . in this case
, the `` david stars '' shaped clusters are annihilated and can not screen the holes any more .
so the holes diffuse into neighboring `` david stars '' .
therefore , upon hole doping in @xmath0-tas@xmath42 in ccdw phase , ccdw should firstly transform to a nccdw or metastable phase composed of ccdw domains and iccdw network .
for example , cho et al .
@xcite applied a very small positive perpendicular voltage pulse on @xmath0-tas@xmath1 single crystal sample within a typical scanning tunneling microscope ( stm ) set - up .
since local hole concentration under the stm tip is largely enhanced , a pulse creates a textured cdw domain of a few tens of nanometers with an irregular domain wall network inside .
they considered that such network is consistent with those in the thermally excited nccdw phase @xcite .
besides , di / dv measured implies the weakening and broadening of the hubbard states together with the reduction of the mott gap inside the textured cdw domain induced by positive perpendicular electric filed @xcite .
that is corresponding to our picture : holes diffuse from the domain wall into the neighboring area inside the domain .
although the concentration of the diffused holes is not high enough to suppress ccdw distortion inside the domain , it can create the hole - electron pair in lhb and uhb to reduce the mott gap .
moreover , the photoexcitation of electrons can be considered as directly doping holes into system .
mihailovic s group found a laser pulse can introduce a transition from ccdw to a hidden state in @xmath0-tas@xmath1 thin flake @xcite .
the raman spectrum of such hidden state is completely different from that in ccdw and nccdw phases @xcite .
such hidden state is demonstrated to be a metastable phase , which is different from nccdw , but should be composed of ccdw domains and iccdw network as well @xcite .
han et al .
observed the electron diffraction spots varies from ccdw case to iccdw case when applying laser to tas@xmath1 in ccdw phase @xcite .
very recently , zhang et al . showed laser can introduce a ccdw - nccdw transition in bulk @xmath0-tas@xmath1 single crystal @xcite . on the other hand , in experiments upon hole doping in @xmath0-tas@xmath1 in high temperature phase , the ccdw transition temperature is significantly lowered @xcite .
both the two kinds of structural evidences can be well explained by our calculation . for the sake of illustration
we draw , we drew a schematic diagram to describe the mechanism of the hole doping induced ccdw - nccdw / metastable phase transition , as shown in figs .
[ fig_6_mechanism ] ( a ) , ( b ) and ( c ) . in the pristine @xmath0-tas@xmath1 ,
ccdw phase has the lowest energy ( fig .
[ fig_6_mechanism ] ( a ) ) .
the hole doping by perpendicular field or by light can stabilize lattice and largely increase the energy of ccdw phase . in this case
, the system transforms to nccdw or other metastable phase ( figs .
[ fig_6_mechanism ] ( b ) and ( c ) ) .    on the other hand , according to our calculation , ccdw lattice distortion is stable upon electron doping . according to the report by cho et al .
, negative perpendicular voltage pulses could not make change in ccdw state @xcite , which is consistent with our estimation . however , some recent experimental works suggest the opposite results .
doped electron into @xmath0-tas@xmath1 by intercalation of li ions and found ccdw can be suppressed @xcite .
the intercalated ions can not only carry electrons , but also strongly influence the lattice structures and induce disorder , which can suppress cdw in tmds @xcite .
therefore , one can not simply attribute the suppression to the electron doping . to demonstrate the pure electron doping effect on cdw ,
more direct doping experiments by negative perpendicular electric field or liquid - gated method should perform .
some works reported the suppression of ccdw by in - plane electric field @xcite . by measuring the temperature dependence of resistivity ( @xmath45 ) , yoshida et al . found the in - plane field can not affect the nccdw / iccdw transition in @xmath45 curve , but can introduce a metastable state with very low resistivity at low temperatures @xcite . if we consider that the @xmath45 in ccdw , nccdw , and iccdw reflects the related structural characteristics , one can infer that the in - plane field induced the metastable state with different structural characteristics , i.e. the in - plane field can suppress ccdw distortion as well .
the suppressions of ccdw by the in - plane field are usually explained as electron injection effect @xcite .
however , according to our calculation , ccdw could be stable upon electron doping . therefore , the effect of in - plane field should be more complex .
we consider some potential mechanisms of the suppressions by in - plane field : besides the pure charge carrier injection , one can expect that the in - plane electric field can also force the electrons in the `` david star '' to be delocalized .
moreover , the in - plane electric field might depinning the cdw , which is suggested in @xmath0-tas@xmath1 bulk single crystals @xcite .
the suppression of ccdw by in - plane electric field might also be due to the thermal activation by local joule heating as current flows through the material .
the real mechanism of the suppression of ccdw by in - plane field needs to be figured out by more experimental and theoretical works .
here we only offer a general description of such process : although ccdw state might still have the lower energy , the in - plane electric field can drive the system to cross over the energy barrier from ccdw to nccdw or other metastable phase , as shown in the schematic diagram in fig .
[ fig_6_mechanism ] ( d ) .
eliashberg function ( left ) and the integrated electron - phonon coupling strength ( right ) for ( a ) the doped bulk @xmath0-tas@xmath1 with @xmath46 holes / f.u . and ( b ) the doped monolayer @xmath0-tas@xmath1 with @xmath47 holes / f.u . under a tensile strain of 2.5% , respectively . ]    as described above , once long range ccdw is suppressed , the transition between mott insulating state to metallic state can be observed : when ccdw is suppressed , domain walls show up , which can be seen as conducting channels to induce metallic state @xcite . in the ccdw area near the domain walls , the hubbard states are weakening and broadening and the mott gap is reduced .
@xcite . furthermore
, superconductivity can emerge in percolated metallic ic network @xcite .
while we are not able to directly model the textured ic phase , we can still estimate the superconductivity using @xmath0 structure qualitatively . approximately using the high symmetry doped bulk @xmath0-tas@xmath1 with @xmath46 holes / f.u . and the doped monolayer @xmath0-tas@xmath1 with @xmath47 holes / f.u . under a tensile strain of 2.5% , for which all the phonon instabilities are just suppressed , we estimated the potential carrier doping induced superconductivity by the electron - phonon coupling calculation .
figure [ fig_7_a2f ] shows the calculated eliashberg spectral function @xmath48 where @xmath49 is the density of states at @xmath19 , @xmath50 is phonon frequency , @xmath51 is electronic energy , and @xmath52 is electron - phonon coupling matrix element .
the total electron - phonon coupling strength is then @xmath53 the calculated @xmath54 for the two samples are 0.81 and 0.96 , respectively .
we estimated @xmath2 based on the modified mcmillan formula @xcite : @xmath55 where the coulomb pseudopotential @xmath56 is set to a typical value of @xmath57 .
the logarithmically averaged characteristic phonon frequency @xmath58 is defined as @xmath59 the calculated @xmath58 for the two samples are 136.7 and 98.4 k , respectively . using those parameters , we can estimate that the superconductivity with @xmath2 of @xmath3 k can emerge in the hole doped @xmath0-tas@xmath1 . to observe the superconductivity we predicted ,
the doping induced metallic network between ccdw domain in nccdw or other metastable phase must be percolated , which requires experimental devices with ability to dope more holes .
recently , suda et al . suggested a high level hole doping technique using photoactive electric double layer @xcite , which might be applied to verify our prediction .
in conclusion , based on the first - principles calculation , we simulated the carrier doping effect in the bulk and monolayer @xmath0-tas@xmath1 .
we found that ccdw is stable upon electron doping , while hole doping can notably suppress the ccdw . according to our analysis ,
the different mechanisms of the reported electric and photoelectric manipulations of ccdw in @xmath0-tas@xmath1 are figured out : light or positive perpendicular electric field induced hole doping significantly increases the energy of ccdw , so that the system transforms to nccdw or similar metastable state . although the ccdw distortion is more stable upon the in - plane electric field induced electron injection , some accompanied effects can drive the system to cross over the energy barrier from ccdw to nccdw or similar metastable state .
we also estimated that a potential superconductivity with @xmath2 of @xmath3 k can be introduced by hole doping in the system . by tuning the carrier density , controllable switching of different states such as ccdw / mott insulating state , metallic state , and even the superconducting state
can be realized in @xmath0-tas@xmath1 , which makes the novel material have a very promising application in the future electronic devices .
we thank prof .
z. g. sheng and dr .
y. z. zhang for helpful discussion .
this work was supported by the national key research and development program under contract no .
2016yfa0300404 , the national nature science foundation of china ( grant no .
11674326 , 11304320 , 11274311 , 11404340 , 1408085ma11 , and u1232139 ) .
10 m. imada , a. fujimori , and y. tokura , rev .
70 , 1039 ( 1998 ) .
e. dagotto , rev .
66 , 763 ( 1994 ) .
p. a. lee , n. nagaosa , and x. -g .
wen , rev .
78 , 17 ( 2006 ) .
h. ohno , d. chiba , f. matsukura , t. omiya , e. abe , t. dietl , y. ohno , and k. ohtani , nature 408 , 944 ( 2000 ) .
y. yamada , k. ueno , t. fukumura , h. t. yuan , h. shimotani , y. iwasa , l. gu , s. tsukimoto , y. ikuhara , and m. kawasaki , science 332 , 1065 ( 2011 ) .
r. e. glover , iii and m. d. sherrill , phys .
lett . 5 , 248 ( 1960 ) .
n. takubo , y. ogimoto , m. nakamura , h. tamaru , m. izumi , and k. miyano , phys .
95 , 017404 ( 2005 ) .
a. zakery , s. r. elliott , optical nonlinearities in chalcogenide glasses and their applications ( springer , new york , 2007 ) . m. nakano , k. shibuya , d. okuyama , t. hatano , s. ono , m. kawasaki , y. iwasa , and y. tokura , nature 487 , 459 ( 2012 ) .
j. a. wilson , f. j. di salvo , and s. mahajan , adv .
24 , 117 ( 1975 ) .
a. y. liu , phy .
b 79 , 220515(r ) ( 2009 ) .
y. ge and amy y. liu , phy . rev .
155133 ( 2010 ) .
y. liu , d. f. shao , l. j. li , w. j. lu , x. d. zhu , p. tong , r. c. xiao , l. s. ling , c. y. xi , l. pi , h. f. tian , h. x. yang , j. q. li , w. h. song , x. b. zhu , y. p. sun , phys .
b 94 , 045131 ( 2016 ) .
p. darancet , a. j. millis , and c.a .
marianetti , phys .
b 90 , 045134 ( 2014 ) .
r. ang , y. tanaka , e. ieki , k. nakayama , t. sato , l. j. li , w. j. lu , y. p. sun , and t. takahashi , phys .
109 , 176403 ( 2012 ) .
r. ang , y. miyata , e. ieki , k. nakayama , t. sato , y. liu , w. j. lu , y. p. sun , and t. takahashi , phys .
b 88 , 115145 ( 2013 ) .
b. sipos , a. f. kusmartseva , a. akrap , h. berger , l. forr , and e. tuti , nature mater . 7 , 960 ( 2008 ) .
y. liu , r. ang , w. j. lu , w. h. song , l. j. li , and y. p. sun , appl .
102 , 192602 ( 2013 ) .
l. j. li , w. j. lu , x. d. zhu , l. s. ling , z. qu , and y. p. sun , europhys .
97 , 67005 ( 2012 ) .
y. yu , f. yang , x. f. lu , y. j. yan , y. -h .
cho , l. ma , x. niu , s. kim , y. -w . son , d. feng , s. li , s. -w .
cheong , x. h. chen , and y. zhang , nat .
10 , 270 ( 2015 ) .
a. w. tsen , r. hovden , d. wang , y. d. kim , j. okamoto , k. a. spoth , y. liu , w. j. lu , y. p. sun , j. hone , l. f. kourkoutis , p. kim , and a. n. pasupathy , pnas 8 , 15054(2015 )
. m. j. hollander , y. liu , w. j. lu , l. j. li , y. p. sun , j. a. robinson , and s. datta , nano lett .
15 , 1861 ( 2015 ) .
masaro yoshida , ryuji suzuki , yijin zhang , masaki nakano , and yoshihiro iwasa , sci . adv . 1 , e1500606 ( 2015 ) . i. vaskivskyi , i. a. mihailovic , s. brazovskii , j. gospodaric , t. mertelj , d. svetin , p. sutar and d. mihailovic , arxiv:1409,3794 ( 2014 ) . i. vaskivskyi , j. gospodaric , s. brazovskii , d. svetin , p , sutar , e. goreshnik , i. a. mihailovic , t. mertelj , and d. mihailovic , sci . adv . 1 , e1500168 ( 2015 ) .
d. cho , s. cheon , k. -s .
kim , s. -h .
lee , y. -h .
cho , s. -w .
cheong , and h. w. yeom , nature comm . 7 , 10453 ( 2016 ) .
y. z. zhang , z. g. sheng , et al . , private communication .
l. stojchevska , i. vaskivskyi , t. mertelj , p. kusar , d. svetin , s. brazovskii , and d. mihailovic , science 344 , 177 ( 2014 ) . t. -r .
t. han , f.zhou , c. d. malliakas , p. m. duxbury , s. d. mahanti , m. g. kanatzidis , and c. -y .
ruan , sci . adv . 1 , e1400173 ( 2015 ) .
p. giannozzi , et al .
, j. phys . : condens .
matter 21 , 395502 ( 2009 ) .
j. p. perdew , j. a. chevary , s. h. vosko , k. a. jackson , m. r. pederson , d. j. singh , and c. fiolhais , phys . rev .
b 46 , 6671 ( 1992 ) .
n. marzari , d. vanderbilt , a. de vita , and m. c. payne , phys .
82 , 3296 ( 1999 ) . h. j. monkhorst and j. d. pack , phys .
b 13 , 5188 ( 1976 ) .
s. baroni , s. de gironcoli , a. dal corso , and p. giannozi , rev .
73 , 515 ( 2001 ) .
m. d. johannes , i. i. mazin , and c. a. howells , phys . rev .
b 73 , 205102 ( 2006 ) . c. battaglia , h. cercellier , f. clerc , l. despont , m. g. garnier , c. koitzsch , p. aebi , h. berger , l. forr , and c. ambrosch - draxl , phys .
b 72 , 195114 ( 2005 ) .
r. bianco , m. calandra and f. mauri , phys .
b 92 , 094107 ( 2015 ) .
q. zhang , l. -y .
gan , y. cheng , and u. schwingenschlcgl , phys .
b 90 , 081103(r ) ( 2014 ) .
s. cahangirov , m. topsakal , e. aktrk , h. ahin , and s. ciraci , phys .
102 , 236804 ( 2009 ) .
h.ahin , s. cahangirov , m. topsakal , e. bekaroglu , e. akturk , r. t. senger , and s. ciraci , phys .
b 80 , 155453 ( 2009 ) .
a. j. mannix , x. -f .
zhou , b. kiraly , j. d. wood , d. alducin , b. d. myers , x. liu , b. l. fisher , u. santiago , j. r. guest , m. jose yacaman , a. ponce , a. r. oganov , m. c. hersam , and n. p. guisinger , science 350 , 1513 ( 2015 ) .
m. a. bissett , s. konabe , s. okada , m. tsuji , and h. ago , acs nano 7 , 10335 ( 2013 ) .
s. -i . park , j. -h .
ahn , x. feng , s. wang , y. huang , and j. a. rogers , adv .
18 , 2673 ( 2008 ) .
e. morosan , h. w. zandbergen , b. s. dennis , j. w. g. bos , y. onose , t. klimczuk , a. p. ramirez , n. p. ong and r. j. cava , nature phys . 2 , 544 ( 2006 )
. k. e. wagner , e. morosan , y. s. hor , j. tao , y. zhu , t. sanders , t. m. mcqueen , h. w. zandbergen , a. j. williams , d. v. west , and r. j. cava , phys .
b 78 , 104520 ( 2008 ) .


the in - medium properties of the vector mesons ( @xmath1 and @xmath0 ) in hot and dense matter have an important role to play in the low mass dilepton production resulting from heavy ion collision experiments .
this has , hence , been a topic of great interest in the recent past , both experimentally @xcite and theoretically @xcite .
the experimental observation of enhanced dilepton production @xcite in the low invariant mass regime could be due to a reduction in the vector meson masses in the medium .
brown and rho suggested the hypothesis that the vector meson masses drop in the medium according to a simple ( br ) scaling law @xcite , given as @xmath2 , where @xmath3 is the pion decay constant and asterisk refers to in - medium quantities . in the framework of quantum hadrodynamics ( qhd ) as a description of the hot hadronic matter , it is seen that the dropping of the vector meson masses has its dominant contribution arising from the vacuum polarisation effects from the baryon sector @xcite , which is not observed in the mean field approximation .
this means that the quantum effects do play an important role in the medium modification of vector meson properties .
there have also been approaches based on qcd sum rules @xcite which confirm a scaling hypothesis @xcite as suggested by brown and rho , with a saturation scheme that leads to a delta function at the vector meson pole and a continuum for higher energies for the hadronic spectral function .
it is , however , seen that such a simple saturation scheme for finite densities does not work and a more realistic description for the hadronic spectral function is called for @xcite . using an effective lagrangian model to calculate the hadronic spectral function
, it is seen that such a universal scaling law is not suggested for in - medium vector meson masses @xcite .
the medium modifications of the vector mesons have been thus a subject of several investigations .
it has been emphasized in the literature @xcite that the properties of the hadrons are modified due to their interactions with the thermal bath , and such modifications are reflected in the dilepton and photon spectra emitted from a hot and dense matter @xcite .
dileptons are interesting probes for the study of evolving matter arising from relativistic heavy ion collisions since they do not interact strongly and escape unthermalized from the hot and dense matter at all stages of the evolution .
the temperature @xcite and density @xcite modifications of the dileptons from hot hadronic matter as well as from a quark gluon plasma ( qgp ) resulting from heavy ion collisions have been a subject of extensive investigations in the recent past .
broadly , two types of modifications for the hadrons in the medium are expected : shift in the pole position for the hadron propagator , giving rise to a modification of the mass , and broadening of the spectral function . in the present work , we shall attempt to study the medium modification of vector meson masses and decay widths in the hot hadronic matter within the framework of quantum hadrodynamics ( qhd ) taking the vacuum polarisation effects into account , and their subsequent effect on the dilepton spectra .    in the conventional hadronic models
@xcite , the masses of the vector mesons stay constant or increase slightly , in the mean field approximation , i.e. , when the polarization from the dirac sea is neglected @xcite . with the inclusion of quantum corrections from the baryonic sector ,
however , in the walecka model one observes a drop in the vector meson masses in the medium @xcite .
this medium modification of the vector meson masses plays an important role in the enhanced dilepton yield @xcite for masses below the @xmath1 resonance in the heavy ion collision experiments @xcite .
it has been emphasized recently that the dirac sea contribution ( taken into account through summing over baryonic tadpole diagrams in the relativistic hartree approximation ( rha ) ) dominates over the fermi sea contribution and @xmath4 is caused by a large dressing of @xmath5 cloud in the medium @xcite .
further quantum effects arising from the scalar mesons have been considered @xcite along with the rha for the baryon sector in the context of vector meson mass modifications in strange hadronic matter @xcite .
the present method is a step in that direction of studying medium modifications of the vector meson properties including such quantum effects from the scalar mesons in a self  consistent manner , along with those arising from the baryon sector .
it was earlier demonstrated in a nonperturbative formalism that a realignment of the ground state with baryon - antibaryon condensates is equivalent to the relativistic hartree approximation @xcite .
the ground state for the nuclear matter was extended to include sigma condensates to take into account the quantum correction effects from the scalar meson sector @xcite .
such a formalism includes multiloop effects and is self consistent @xcite .
the methodology was then generalized to consider hot nuclear matter @xcite as well as to the situation of hyperon - rich dense matter @xcite relevant for neutron stars .
the effect of vacuum polarisations on the vector meson masses has also been recently studied @xcite . in the present work ,
we study the effect of such quantum corrections on the in - medium vector meson properties .
this apart , for studying @xmath1 meson properties , we take into account the @xmath1-@xmath6-@xmath6 interactions , and their effects on the dilepton production in the low invariant mass regime .
the pionic interactions are seen to modify the @xmath1 meson mass only marginally , as already stated in the literature @xcite . due to scalar meson contributions ,
there is broadening in the @xmath0 and @xmath1 peaks , as compared to rha .
this is seen to be quite pronounced for the @xmath0 meson , which leads to smearing and ultimate disappearance of the @xmath0 peak at finite densities @xcite .
we organize the paper as follows .
we first briefly recapitulate the nonperturbative framework used for studying the quantum correction effects in hot nuclear matter in section ii .
the medium modification of the @xmath0 and @xmath1 meson masses and decay widths are considered in section iii .
we discuss the effect of these medium modifications on the low mass dilepton spectra in section iv . in section
v , we summarize the results of the present work and give an outlook .
we briefly recapitulate here the vacuum polarisation effects arising from the nucleon and scalar meson fields in hot nuclear matter in a nonperturbative variational framework @xcite .
the method of thermofield dynamics ( tfd ) @xcite is used here to study the  ground state " ( the state with minimum thermodynamic potential ) at finite temperature and density within the walecka model with a quartic scalar self interaction .
the temperature and density dependent baryon and sigma masses are also calculated in a self - consistent manner in the formalism .
the ansatz functions involved in such an approach are determined through functional minimisation of the thermodynamic potential .
the lagrangian density in the walecka model is given as @xmath7 in the above , @xmath8 , @xmath9 , and @xmath10 are the fields for the nucleon , @xmath9 , and @xmath0 mesons with masses @xmath11 , @xmath12 , and @xmath13 respectively
. since we are interested in symmetric nuclear matter , the isovector rho meson does not contribute .
the quartic coupling term in @xmath9 is necessary for the sigma condensates to exist , through a vacuum realignment @xcite .
our calculations thus include the quantum effects arising from the sigma meson in addition to the mean field contribution from the the quartic self interaction of the scalar meson .
we retain the quantum nature of both the nucleon and the scalar meson fields , whereas the vector @xmath0 meson is treated as a classical field , using the mean field approximation for the @xmath0 meson , given as @xmath14 .
the reason is that without higher  order term for the @xmath0 meson , the quantum effects generated due to the @xmath0 meson through the present variational ansatz turn out to be zero .
we then write down the expressions for the thermodynamic quantities including the quantum effects .
the details regarding the formalism can be found in earlier references @xcite .
the energy density , after carrying out the renormalisation procedures for the baryonic and scalar meson sectors @xcite , is obtained as @xmath15 with    @xmath16    @xmath17    @xmath18    @xmath19    as the mean field result , contribution from the dirac sea , and contributions from the @xmath0 and @xmath9 mesons , respectively .
we might note here that the quantum effects arising here from the scalar meson sector through @xmath9 meson condensates amounts to a sum over a class of multiloop diagrams and , does not correspond to the one meson loop approximation for scalar meson quantum effects considered earlier @xcite . in the above , @xmath20 and @xmath21 are the thermal distribution functions for the baryons and antibaryons , given in terms of the effective energy , @xmath22 , and the effective chemical potential , @xmath23 .
@xmath24 is the effective nucleon mass and @xmath25 is the baryon number density after subtracting out the vacuum contribution .
the spin - isospin degeneracy factor is @xmath26 for symmetric nuclear matter . in the expression for the energy density arising from the scalar meson sector ,
the field dependent effective sigma mass , @xmath27 , satisfies the gap equation in terms of the renormalised parameters as @xmath28 where @xmath29 where @xmath30 .
the effective meson mass @xmath31 is determined through a self consistent solution of eq .
( [ mm2 ] ) . in eq .
( [ vph0 ] ) , @xmath32 and @xmath33 are the expressions as given by eqs .
( [ mm2 ] ) and ( [ if ] ) with @xmath34 .
we might note here that the gap equation given by eq .
( [ mm2 ] ) is identical to that obtained through resumming the daisy and superdaisy graphs @xcite and hence includes higher  order corrections from the scalar meson field in a self - consistent manner .
extremisation of the thermodynamic potential , @xmath35 , with respect to the meson fields @xmath36 and @xmath37 yields the self
consistency conditions for @xmath36 ( and hence for the effective nucleon mass @xmath38 ) and for the vector meson field @xmath37 .
we now examine the medium modification to the masses and decay widths of the @xmath0 and @xmath1 mesons including the quantum correction effects in hot nuclear matter in the relativistic random phase approximation .
the interaction vertices for these mesons with nucleons are given as    @xmath39    where @xmath40 or @xmath41 , @xmath42 is the free nucleon mass , @xmath8 is the nucleon field and @xmath43 or @xmath44 , @xmath45 being the pauli matrices . @xmath46 and
@xmath47 correspond to the couplings due to the vector and tensor interactions for the corresponding vector mesons to the nucleon fields .
the vector meson self energy is expressed in terms of the nucleon propagator , @xmath48 modified by the quantum effects .
this is given as @xmath49,\ ] ] where @xmath50 is the isospin degeneracy factor for nuclear matter , and @xmath51 represents the meson - nucleon vertex function .
the vector meson self energy can be written as the sum of two parts @xmath52 where @xmath53 is the contribution arising from the vacuum polarisation effects , described by the coupling to the @xmath54 excitations and @xmath55 is the density dependent contribution to the vector self energy . for the @xmath0 meson ,
the tensor coupling is generally small as compared to the vector coupling to the nucleons @xcite , and is neglected in the present work .
the feynman part of the self energy , @xmath56 is divergent and needs renormalization .
we use dimensional regularization to separate the divergent parts and the renormalization condition @xmath57 . because of the tensor interaction , the vacuum self energy for the @xmath1 meson is not renormalizable .
we employ a phenomenological subtraction procedure @xcite to extract the finite part using the renormalization condition @xmath58 , with ( @xmath59 ) .
the effective mass of the vector meson is obtained by solving the equation , with @xmath60 , @xmath61      the decay width for the process @xmath62 is calculated from the imaginary part of the self energy using the cutkosky rule , and in the rest frame of the @xmath1 meson is given by @xmath63 \label { gmrho}\ ] ] where , @xmath64^{-1}$ ] is the bose - einstein distribution function .
the first and the second terms in the equation ( [ gmrho ] ) represent the decay and the formation of the resonance , @xmath1 .
the medium effects have been shown to play a very important role for the @xmath1 meson decay width . in the calculation for the @xmath1 decay width
, the pion has been treated as free , and any modification of the pion propagator due to effects such as delta - nucleon hole excitation @xcite to yield a finite decay width for the pion , have not been taken into account .
the coupling @xmath65 is fixed from the decay width of @xmath1 meson in vacuum ( @xmath66=151 mev ) decaying to two pions .    for the nucleon - rho couplings , we use the vector and tensor couplings as obtained from the n - n forward dispersion relation @xcite . with the couplings as described above , we consider the temperature and density dependence of the @xmath0 and @xmath1 meson properties in the hot nuclear matter due to quantum correction effects .
to calculate the decay width for the @xmath0 meson , we write down the effective lagrangian for the @xmath0 meson as @xcite @xmath67 in the above , the first term refers to the direct coupling of the vector meson @xmath0 to the photon , and hence to the dilepton pairs , as given by the vector dominance model .
the decay width of the @xmath0 meson in vacuum is dominated by the channel @xmath68 . in the medium , the decay width for @xmath69 is given as @xmath70,\end{aligned}\ ] ] where @xmath71 , @xmath72 and @xmath73 s are 4-momenta and energies for the pions , and @xmath74 s are their thermal distributions .
the matrix element @xmath75 has contributions from the channels @xmath76 ( described by the second term in ( [ lomg ] ) ) and the direct decay @xmath69 resulting from the contact interaction ( last term in ( [ lomg ] ) ) @xcite . for the @xmath77 coupling
we take the value @xmath78 = 2 according to ref .
@xcite , which is similar in value to the coupling taken in @xcite .
we fix the point interaction coupling @xmath79 by fitting the partial decay width @xmath68 in vacuum ( 7.49 mev ) to be 0.24 .
the contribution arising from the direct decay is upto around 15 % which is comparable to the results of @xcite .    with the modifications of the vector meson masses in the hot and dense medium ,
a new channel becomes accessible , i.e. , the decay mode @xmath80 for @xmath81 , which is not possible in the vacuum ( since @xmath82 ) .
the in - medium decay width for this process is found from the cutting rules at finite temperature to be @xmath83\ ] ] where @xmath84 arises from phase space considerations , and @xmath85 s are the thermal distribution functions for the pions and @xmath1 mesons .
the first and second terms refer to the emission and absorption of pion and @xmath1 meson .
the pion annihilation ( @xmath86 ) dominates dilepton production in the low invariant mass region .
the dropping of the vector meson masses in the hot and dense matter also gives rise to substantial contribution to the dilepton yield from the direct decays of the vector mesons into dileptons ( @xmath87 , @xmath88 ) .
the thermal production rate per unit four volume for lepton pairs is related to the imaginary part of the one  particle irreducible photon self energy by @xcite    @xmath89    where @xmath90 is the 4-momentum of the virtual photon , and @xmath91 is the fine structure constant .      the invariant mass distribution for the lepton pairs from pion annihilation is given by @xmath92 where , @xmath93 , @xmath94 . in the above , @xmath95 is the cross section for the pion annihilation , and is given by    @xmath96    where @xmath97 is the electromagnetic form factor of the pion .
the photon can also couple to the hadrons through vector mesons , which is the basis of the vector dominance model @xcite .
the photon can thus interact with the pions through a vector meson , which must be an isovector and is taken to be the @xmath1 meson .
the effective lagrangian for the photon - pion - rho interaction can then be written as @xcite @xmath98 where @xmath99 is the field tensor for the @xmath1 field and @xmath100 .
the coupling @xmath65 in the above can be determined from the decay @xmath62 , and the direct photon - rho coupling @xmath101 is determined from the process @xmath102 @xcite .
these , including a finite decay width for the @xmath1 meson , then lead to the @xmath1-dominated pion form factor as @xmath103 in the rest frame of the @xmath1 meson , the invariant mass distribution of the lepton pair of invariant mass @xmath11 from pion annihilation is given as @xmath104 with the cross section @xmath105 for the pion annihilation , given in terms of the pion form factor at @xmath106 .
the summation over @xmath107 arises from the bose - einstein distribution function for the @xmath1 meson .
since we are mostly in the regime @xmath108 , the first term in the series , corresponding to the boltzmannian distribution function is dominant , with quantum corrections being of the order of upto 10 % .
a vector meson @xmath109 in a thermal bath can decay into leptons and photons , described by the final state @xmath110 say , which escape the bath without thermalization . with @xmath111 as the four momentum of the nonhadronic state @xmath110 , a resonance peak should appear in the invariant mass ( @xmath11 ) plot for the number of lepton pair events , versus @xmath106 at @xmath112 , where @xmath113 is the mass of the resonance in the thermal bath .
hence , the medium modifications of the vector mesons get reflected in the dilepton spectra through the shift and broadening of the peaks .    the invariant mass distribution of lepton pairs from the decays of vector mesons is given by a generalization of the breit wigner formula as @xcite @xmath114 in the above
, @xmath115 is the number of lepton pair events , say in the space time and four momentum element @xmath116 , @xmath117 being the spin of the resonance , @xmath118 is the formation width , @xmath119 is the total width and @xmath120 is the partial decay width for the off - shell @xmath109 ( i.e. of mass @xmath11 ) to go into the nonhadronic state @xmath110 .
the total width @xmath121 is the rate at which the particles equilibriate and relax to chemical equilibrium .
the invariant mass distribution of the lepton pairs from decays of vector mesons propagating with a momentum @xmath122 is obtained from equation ( [ gbw ] ) as    @xmath123 q_0 \gamma _ { v \rightarrow e^+ e^-}^{\rm vac}\ ] ]    in the rest frame of the vector meson , this reduces to @xmath124 in the above , @xmath125 is the partial width for the leptonic decay for the off - shell vector meson @xmath109 with mass @xmath11 , and is given as @xmath126 the coupling @xmath46 refers to the direct coupling of the vector meson @xmath109 to the photon as given in the equations ( [ lomg ] ) and ( [ lrho ] ) for the @xmath0 and @xmath1 mesons respectively .
the coupling of @xmath0 to the photon taken to be one third of the coupling of @xmath1 to the photon gives reasonable agreement with the ratio of the observed partial widths @xmath127 @xcite .
in the following we present numerical results on the influence of vacuum effects on hadronic masses , decay widths of vector mesons , and dilepton emission from hot hadronic matter . as described in the previous sections , the study is based entirely on a hadronic scenario and makes use of the original walecka model , extended by a nonlinear sigma self coupling .    in this work we concentrate on one particular aspect , namely on the role of scalar vacuum effects . for
each quantity studied , we compare the results of the relativistic hartree approximation ( rha ) , which was already treated in ref .
@xcite , with those additional effects arising from a scalar vacuum condensate .
the quartic self coupling @xmath128 is taken as a free parameter ; we have used @xmath129 and @xmath130 as representative values .
the values of the incompressibility corresponding to these self couplings are 401 mev ( for @xmath128=1.8 ) and 329 mev ( for @xmath128=5 ) .
these values are smaller than the mean field result of 545 mev as well as the relativistic hartree result ( @xmath128=0 ) of 450 mev .
this is because an increase of @xmath128 leads to weaker potentials of the @xmath9 and @xmath0 mesons and to a softening of the equation of state .
the parameters of the model are adjusted to fit nuclear matter , corresponding to a saturation density of @xmath131 = 0.193/@xmath132 .
the fitted values for the scaled coupling constants are @xmath133 183.3 , 167.5 and 137.9 , and @xmath134 114.7 , 96.45 , and 63.7 for @xmath135 ( i.e. , rha ) , 1.8 , and 5 , respectively .
when discussing dilepton production , we do not attempt to give a realistic description of specific experimental results . for this , the dynamic space - time evolution for the hot hadronic matter in the course of a nuclear collision would have to be addressed and additional effects contributing to the dilepton spectrum would have to be included . however , our study should help to elucidate , in which way in - medium vacuum effects can have an influence on experimentally observable quantities .
we first discuss the effect of quantum corrections on the in - medium nucleonic properties for hot nuclear matter @xcite . in fig .
1a , we plot the temperature dependence of the nucleon mass for the baryon density @xmath136 and the nuclear matter density @xmath137 .
the quantum corrections arising from the sigma field have the effect of softening the equation of state , and produce a higher value for the effective nucleon mass as compared to the rha .
the effective nucleon mass was seen to increase with the vacuum polarisation effects arising from the baryons as compared to the mft calculations @xcite .
though for @xmath136 , the effective nucleon mass decreases monotonically with temperature , for @xmath138 , the nucleon mass first increases and then decreases as a function of temperature , which has also been observed earlier to be the trend for higher densities @xcite .
the variation in the nucleon mass is very slow up to a temperature @xmath139 150 mev , beyond which there is a fast decrease .
though qualitative features are similar for the in - medium mass with inclusion of quantum effects from the scalar meson sector , it is observed that such effects lead to a higher value for the nucleon mass .
we next study the temperature and density dependence of the vector mesons in hot nuclear matter .
the @xmath0 meson mass as a function of the temperature for @xmath140=0 and @xmath131 , is plotted in fig .
1b . in the walecka model ,
the dirac sea has been shown to have a significant contribution over the fermi sea , leading to a substantial drop in the omega meson mass in nuclear matter @xcite .
specifically , at saturation density and at zero temperature , the decrease in the omega meson mass from the vacuum value , is around 150 mev in rha , whereas with sigma quantum effects , the drop in the mass reduces to around 117 mev or 61 mev for the two values of @xmath128 studied .
a similar reduction in the mass drop , as compared to the rha was also observed when scalar field quantum corrections were included within a one  loop approximation @xcite , unlike the approximation adopted here , which is self - consistent and includes multi - loop effects .
this indicates that the quantum corrections do play an important role in the medium modification of the vector meson masses . in the present work
, it is seen that the quantum corrections from the scalar mesons , over those arising from the baryonic sector , have the effect of giving rise to a higher value for the @xmath0 meson mass .
this again is related to the fact that the effective nucleon mass is higher when such quantum effects are taken into account .    in figure [ figf2 m ]
, we illustrate the medium modification for the @xmath1 meson mass with the vector and tensor couplings to the nucleons being fixed from the nn forward dispersion relation @xcite . the values for these couplings
are given as @xmath141=0.55 and @xmath142=6.1 .
we notice that the decrease in the @xmath1 meson mass with increase in temperature is much sharper than that of the @xmath0 meson .
such a behaviour of the @xmath1 meson undergoing a much larger medium modification was also observed earlier @xcite within the relativistic hartree approximation for the nucleons .
this indicates that the tensor coupling , which is negligible for the @xmath0 meson , plays a significant role for the @xmath1 meson . with inclusion of quantum corrections from the scalar meson
, the qualitative comparison between the @xmath0 and @xmath1 vector mesons remains the same , though these effects are generally seen to lead to larger values for the vector meson masses .
we have also calculated the contribution to the @xmath1 meson mass arising from pion loops through the @xmath143 coupling , but this was found to lead to a marginal increase of at most 1% , in agreement with @xcite .    in figure [ figf3 m ] , the decay widths for the @xmath1 meson with and without the scalar quantum corrections , over the relativistic hartree approximation , are considered .
such quantum effects lead to an appreciable increase in the decay widths , which , as we shall see , is observed as a broadening of the @xmath1-peak in the dilepton spectra .    in figure [ figf4 m ] , the decay width for the @xmath0 meson is seen to have a positive contribution arising from the @xmath9-meson quantum effect as compared to rha .
the decay width is calculated as arising from the processes @xmath69 as well as @xmath80 , a channel which becomes possible with the mass modifications of the vector mesons in the medium .
it is seen that the channel @xmath144 opens for temperature greater than around 180 mev for @xmath136 , whereas it is observed even at zero temperature at nuclear matter saturation density .
we might note that the collisional modification to the @xmath0 width , arising from the process @xmath145 @xcite , has not been taken into account in the present calculations .
the dilepton spectra for the low invariant mass regime , at t=200 mev are shown in figures [ figf5 m ] and [ figf6 m ] calculated for densities @xmath136 and @xmath138 .
the total dilepton spectra along with contributions arising from pion annihilation and direct decays of the @xmath1 and @xmath0 mesons , are displayed in figures 5a and 6a .
the effect of sigma condensates shifts the peak position for the vector mesons to higher values as compared to rha , and also results in a broadening of the peaks in the dilepton spectra as can be seen from figures 5b and 6b .
the @xmath1 and @xmath0 peaks are well separated due to the fact that the @xmath1 meson mass gets much larger medium modification as compared to the @xmath0 meson , which would not be the case in the absence of medium effects .
the omega peak is seen to be smeared due to scalar meson quantum corrections .
the broadening of the @xmath0 peak due to @xmath9-quantum effects is more clearly visible for @xmath138 than for zero density , in the dilepton spectra .
the dominant contributions to the dilepton enhancement are from pion annihilation and @xmath1 meson decay , which are plotted in figures 5c and 6c in comparison with the rha results .
the effect of @xmath9 condensates leads to a higher value for the peak position as well as broadening of the @xmath1-peak .
the dilepton production rates due to the @xmath0 meson decay are displayed in figures 5d and 6d .
it is observed that the broadening of the @xmath0 peak due to scalar quantum corrections is larger for @xmath131 than for zero density .
this is a reflection of the fact that there is a larger modification of the omega width due to scalar quantum effects for higher densities .
there have been various mechanisms under discussion which can lead to an increase in the @xmath0 width .
it has been argued that the @xmath0 decay width can increase by an order of magnitude due to the reduction in the pion decay constant related with chiral symmetry restoration @xcite . the medium modification due to scalar meson quantum effects giving rise to higher values for the @xmath1 meson mass and decay width are reflected in the dilepton spectra .
the medium modification for the vector meson masses depend more sensitively on the density than on temperature
. this can give rise to an enhancement of the dilepton production rate when matter is highly compressed at the initial stage of a heavy ion collision .
in the dls experiments for a - a collisions at 1 agev , there have been models that focus on density induced change for the @xmath1 meson mass to describe the dilepton pair production for invariant mass of dilepton pair , @xmath146gev @xcite . even though there is dilepton enhancement arising from the density  modified @xmath1 mass , the results still underestimate the dls data by a factor of 2 - 3 @xcite .    the mass modifications for the vector meson masses with density are plotted in figure [ figf7 ] , for temperature t= 100 mev . the trend is a sharp decrease of the masses with density up to @xmath147 , and then a rise for higher values of @xmath140 . the observed trend is due to the competing effects of vacuum polarisation which leads to a reduction in the mass , and the density dependent dressing for the vector self energy causing an increase in the meson mass .
such a behaviour of the vector meson mass for the @xmath0 meson was also observed earlier @xcite
. a higher value of the scalar self interaction coupling , @xmath128 , leads to a higher value for the vector meson masses .
the change in the @xmath1 meson mass for higher densities is seen to be rather slow which means that there is a delicate cancellation between the feynman and the density dependent parts of the vector self energy . as depicted in figure [ figf7 ]
, the @xmath0 meson mass seems to be more dominated by the density dependent contributions for higher densities .
the decay widths are plotted in figure [ figf8 ] . the slow change in @xmath1 meson mass for higher densities
is reflected in the in - medium @xmath1-decay width , while the @xmath0 decay width @xmath148 grows strongly with density and gets strongly enhanced by quantum effects from the @xmath9 meson .
the dilepton spectra for @xmath149 mev and @xmath150 are plotted in figures [ figf91 ] and [ figf92 ] .
these values should be typical for the compression zone of a nuclear collision at energies of a few gev per nucleon .
for example , for inclusive ca+ca collisions at 1 agev studied by the dls collaboration @xcite , maximum temperatures and densities in this range were extracted from a transport model @xcite . in figure 9a ,
the total dilepton spectrum is plotted , for @xmath128=1.8 , displaying the individual contributions .
figure 9b shows that the sigma quantum corrections lead to a broadening of the peaks in the dilepton spectra , which is seen to be more appreciable for the higher value of the scalar self coupling , @xmath128 .
the individual contributions to the dilepton spectra , plotted in figure [ figf92 ] , illustrate the effect of @xmath9 quantum corrections .
finally , the density dependence for the dilepton spectra is shown for temperature @xmath149 mev in figures [ figf101 ] and [ figf102 ] for couplings @xmath128 as 1.8 and 5 , respectively , in comparison with the rha results .
the vector meson masses decreasing with density up to around twice the nuclear matter density and then slowly increasing or staying constant , are reflected in the dilepton spectra through the peak positions .
one might notice that the broadening of the peaks due to @xmath9- quantum effects is larger for higher densities and is more pronounced for the @xmath0 peak .
the curves for densities twice and four times nuclear matter density in figure [ figf102 ] are coincident up to a value of the invariant mass of about 0.6 gev .
this is due to the fact that the rho meson properties whose effects are relevant for that mass region , remain almost constant for @xmath128=5 . beyond that however , they are different around the invariant mass region for the omega meson . due to substantial increase in the decay width of the @xmath0 meson with density the @xmath0 peak
is almost smeared out even at nuclear matter density and the omega loses its quasi particle nature .
the corrections from the scalar meson quantum effects are thus seen to be appreciable .
to summarize , we have considered in the present paper , the modifications of the vector meson ( @xmath1 and @xmath0 ) properties due to quantum correction effects in hot and dense nuclear matter and their effect on dilepton spectra .
the baryonic properties , as modified due to such effects , subsequently determine the vector meson masses and decay widths in the hot and dense hadronic matter .
it has been recently emphasized that the dirac sea contribution dominates over the fermi sea , thus implying that the vacuum polarisation effects arising from the baryonic sector do play an important role in the vector meson properties in a medium .
we consider the hadronic properties as modified arising due to the quantum effects from the baryon and scalar meson sectors and investigate their effect on the dilepton production rate .
the effect of the scalar meson quantum corrections leads to an increase in the vector meson masses and decay widths .
the effect on the omega meson decay width is particularly more pronounced at higher densities , than in the rha .
the value for @xmath149 mev is modified from 67 mev ( for rha ) to 125 mev ( for @xmath128=1.8 ) and 251 mev ( for @xmath128=5 ) for @xmath151=2 , and for @xmath151=4 the rha value of 101 mev is modified to 175 mev and 370 mev , for @xmath128=1.8 and 5 , respectively .    in the present work ,
we have studied the effect of quantum corrections arising from the scalar meson over those resulting from the baryon sector ( given through the summing over baryonic tadpole diagrams in the relativistic hartree approximation ) .
the quantum corrections from the @xmath9-meson give rise to a higher value for the vector meson masses as compared to keeping only quantum effects from the baryons .
there is also an appreciable increase in the decay widths in the presence of these quantum corrections , and this is reflected in the invariant mass plot for dilepton spectra through the broadening of the peaks . for the @xmath0 meson , an appreciable broadening arising due to sigma meson quantum effects leads to the disappearance of the peak from the spectra for high densities when @xmath0 ceases to be exist as a quasi particle .
the density dependence of in - medium properties of the vector mesons and their subsequent effects on the dilepton spectra have been studied , which will be relevant for highly compressed matter at the early stage of a heavy ion collision .
one of the authors ( am ) acknowledges financial support from alexander von humboldt foundation when the work was initiated .
am would like to thank s. scherer , h. weber , d. dietrich for discussions , and would like to acknowledge the institut fr theoretische physik for warm hospitality and dfg ( sto 161/8 - 1 ) for financial support .
99 n. masera for the helios-3 collaboration , nucl
. phys . * a 590 * , 93c ( 1995 ) .
g. agakichiev et al ( ceres collaboration ) , phys .
rev . lett . * 75 * , 1272 ( 1995 ) ; g. agakichiev et al ( ceres collaboration ) , phys . lett .
* b 422 * , 405 ( 1998 ) ; g. agakichiev et al ( ceres collaboration ) , nucl . phys . * a 661 * , 23c ( 1999 ) .
r. j. porter et al ( dls collaboration ) , phys .
. lett . * 79 * , 1229 ( 1997 ) ; w. k. wilson et al ( dls collaboration ) , phys .
c 57 , 1865 ( 1998 ) .
d. p. morrison ( phenix collaboration ) , nucl .
phys . * a 638 * , 565c ( 1998 ) .
j. stroth ( hades collaboration ) , advances nuclear dynamics * 5 * , 311 ( 1999 ) .
g. e. brown and m. rho , phys .
lett . * 66 * , 2720 ( 1991 ) . t. hatsuda and su h. lee , phys . rev .
* c 46 * , r34 ( 1992 ) ; t. hatsuda , s. h. lee and h. shiomi , phys . rev .
* c 52 * , 3364 ( 1995 ) .
x. jin and d. b. leinweber , phys .
* c 52 * , 3344 ( 1995 ) ; t. d. cohen , r. d. furnstahl , d. k. griegel and x. jin , prog .
* 35 * , 221 ( 1995 ) ; r. hofmann , th .
gutsche , a. faessler , eur .
j. * c 17 * , 651 ( 2000 ) ; s. mallik and k. mukherjee , phys . rev . *
d 58 * , 096011 ( 1998 ) .
s. mallik and a. nyffeler , phys . rev .
* c 63 * , 065204 ( 2001 ) .
f. klingl , n. kaiser , w. weise , nucl . phys . * a 624 * , 527 ( 1997 ) . c. ernst , s. a. bass , m. belkacem , h. stcker and w. greiner , phys
* c 58 * , 447 ( 1998 ) .
h. shiomi and t. hatsuda , phys .
b 334 * , 281 ( 1994 ) . t. hatsuda , h. shiomi and h. kuwabara , prog .
. phys . * 95 * , 1009 ( 1996 ) .
jeans , j. piekarewicz and a. g. williams , phys . rev .
* c 49 * , 1981 ( 1994 ) ; k. saito , k. tsushima , a. w. thomas , a. g. williams , phys .
* b 433 * , 243 ( 1998 ) .
jan - e alam , s. sarkar , p. roy , b. dutta - roy and b. sinha , phys . rev .
* c 59 * , 905 ( 1999 ) .
r. rapp and j. wambach , adv .
phys . * 25 * , 1 ( 2000 ) . c. m. ko , v. koch and g. q. li , ann .
* 47 * , 505 ( 1997 ) .
jan - e alam , p. roy , s. sarkar and b. sinha , nucl - th/0106038 ; jan - e alam , s. sarkar , p. roy , t. hatsuda and b. sinha , ann .
phys . * 286 * , 159 ( 2000 ) ; f. karsch , e. laermann , p. petreczky , s. stickan , i. wetzorke , hep - lat/0110208 . k. redlich , j. cleymans , v. v. goloviznin , in proceedings on nato advanced workshop on hot hadronic matter : theory and experiment , n.y . plenum press , 1995 , 562p ; a. dumitru , d. h. rischke , th .
schnfeld , l. winckelmann , h. stcker and w. greiner , phys .
lett . * 70 * , 2860 ( 1993 ) ; p. jaikumar , r. rapp and i. zahed , hep - ph/0112308 ; song gao , ru - keng su , xue- qian li , comm .
* 28 * , 207 ( 1997 ) ; d. dutta , k. kumar , a. k. mohanty , r. k. choudhury , phys . rev . * c 60 * , 014905 , 1999 . b. d. serot and j. d. walecka , adv .
nucl . phys .
* 16 * , 1 ( 1986 ) .
s. a. chin , ann . phys . * 108 * , 301 ( 1977 ) ; m. asakawa , c. m. ko , p. levai and x. j. qiu , phys . rev .
* c 46 * , r1159 ( 1992 ) .
s. pal , song gao , h. stcker and w. greiner , phys
* b 465 * , 282 ( 1999 ) .
a. mishra , p. k. panda , s. schramm , j. reinhardt and w. greiner , phys .
* c 56 * , 1380 ( 1997 ) .
a. mishra and h. mishra , j. phys .
* g 23 * , 143 ( 1997 ) .
s.y . pi and m. samiullah , phys . rev . *
d 36 * , 3121 ( 1987 ) ; g. a. camelia and s.y .
pi , phys . rev . * d 47 * , 2356 ( 1993 )
. a. mishra , p. k. panda and w. greiner , jour .
g 27 * , 1561 ( 2001 ) .
a. mishra , p. k. panda and w. greiner , jour .
* g 28 * , 67 ( 2002 ) .
a. mishra , j. c. parikh and w. greiner , jour .
* g 28 * , 151 ( 2002 ) . c. gale and j. i. kapusta , nucl .
* b 357 * , 65 ( 1991 ) .
r. a. schneider , t. renk and w. weise , nucl .
* a 698 * , 428 ( 2002 ) ; t. renk , r. a. schneider and w. weise , phys . rev . * c 66 * , 014902 ( 2002 ) . h. umezawa , h. matsumoto and m. tachiki , _ thermofield dynamics and condensed states _ ( north - holland , amsterdam , 1982 ) .
m. asakawa , c. m. ko , p. levai and x. j. qiu , phys . rev .
* c 46 * , r1159 ( 1992 ) ; m. herrmann , b. l. friman and w. nrenberg , nucl .
phys . * a 560 * , 411 ( 1993 ) ; g. chanfray and p. shuck , nucl . phys . * a 545 * , 271c ( 1992 ) .
w. grein , nucl .
* b 131 * , 255 ( 1977 ) ; w. grein and p. kroll , nucl .
phys . * a 338 * , 332 ( 1980 ) . j. j. sakurai , _ currents and mesons _
( the university od chicago press , chicago , 1969 ) .
m. gell - mann , d. sharp , and w. d. wagner , phys .
rev . lett .
* 8 * , 261 ( 1962 ) .
b. a. li , phys . rev . * d 52 * , 5165 ( 1995 ) . f. klingl , n. kaiser and w. weise , z. phys . * a 356 * , 193 ( 1996 ) .
. kaymakcalan , s. rajeev and i. schechter , phys . rev .
* d 30 * , 594 ( 1984 ) . m. le bellac , _ thermal field theory _
( cambridge university press , new york , 1996 ) .
h. b. connell , b. c. pearce , a. w. thomas and a. g. williams , prog .
. phys . * 39 * , 201 ( 1997 ) .
h. a. weldon , ann .
( n. y. ) * 228 * , 43 ( 1993 ) .
g. q. li , c. m. ko , g. e. brown , nucl . phys . * a 606 * , 568 ( 1996 ) .
s. sarkar , jan - e alam , p. roy , a.k .
dutt - mazumder , b. dutta - roy and b. sinha , nucl . phys . * a 634 * , 206 ( 1998 ) .
r. d. pisarski , hep - ph/9503330 .
g. chanfrey and p. schuck , nucl .
phys . * a 555 * , 329 ( 1993 ) ; g. q. li and c. m. ko , nucl . phys . * a 582 * , 731 ( 1995 ) .
e. l. bratkovskaya , w. cassing , r. rapp and j. wambach , nucl . phys . * a 634 * , 168 ( 1998 ) .
j. zhang and c. gale , phys . rev .
* c 50 * , 1617 ( 1994 ) .