) laser pulses [@Wood:2003]. Later, investigations of laser-induced damage in the picosecond (ps) and fs regime gained significance when the advent of the chirped-pulse amplification (CPA) technique [@Strickland:1985] made it feasible to develop fs and petawatt-class laser systems [@Sung:2010; @Yu:2012]. In particular, laser ablation occurring on the fs time scale became critical because a laser pulse duration of few tens of fs is much shorter than the time scale for electron energy transfer to the lattice and subsequent lattice heating. In 1995, Stuart $\mathit{et \, al.}$, investigated the laser-induced damage threshold at $1053 \, \mathrm{nm}$ and $526 \, \mathrm{nm}$ for pulse durations ranging from $270 \, \mathrm{fs}$ to $1 \, \mathrm{ns}$, through a theoretical model based on CB electron production via multiphoton ionization, Joule heating and collisional ionization [@Stuart:1995]. Subsequent studies by other groups were conducted for a more accurate analysis of the damage and ablation threshold by including energy dependence of the CB electrons [@Rethfeld:2002] and nonlinear pulse propagation effect in a medium [@Penano:2005; @Petrov:2008; @Gulley:2010; @Apalkov:2012] in the fs regime. However, all these studies were based on theoretical models that used experimental and/or empirical values of the material parameters such as ionization rate, refractive index, relaxation rate, and band structure. Hence, the need for developing a method that uses non-empirical values of the material parameters grew continuously in the search for a comprehensive and reliable method of investigating laser–matter interactions in the intermediate laser intensity range.

In this paper, we employ an alternative method to compute the optical response and the ablation threshold of an optical medium. In contrast to the previous studies, our method is based on first-principles simulations computed from fundamental equations. A multiscale approach using the wave equation and the time-dependent density functional theory (TDDFT) is applied to calculate directly the density of the CB electrons generated in the optical medium. The first report on the use of such a multiscale approach for investigation of the interaction between a laser pulse and an optical medium was made for crystalline silicon, where it was said to yield reliable results [@Yabana:2012]. In this approach, no empirical parameters and approximations were used except for information on the crystal structure and on
) laser pulses [ @Wood:2003 ]. Later, investigations of laser - induced price in the picosecond (ps) and fs government gained significance when the advent of the tweedle - pulse amplification (CPA) proficiency [ @Strickland:1985 ] made it feasible to develop fs and petawatt - class laser organization [ @Sung:2010; @Yu:2012 ]. In particular, laser ablation occurring on the fs prison term scale became critical because a laser pulsation duration of few tens of fs is much shorter than the time scale for electron energy transfer to the wicket and subsequent lattice heating. In 1995, Stuart $ \mathit{et \, al.}$, investigated the laser - induce damage threshold at $ 1053 \, \mathrm{nm}$ and $ 526 \, \mathrm{nm}$ for pulsation durations crop from $ 270 \, \mathrm{fs}$ to $ 1 \, \mathrm{ns}$, through a theoretical model base on CB electron production via multiphoton ionization, Joule heating and collisional ionization [ @Stuart:1995 ]. Subsequent studies by early groups were conducted for a more accurate analysis of the damage and ablation threshold by including energy dependence of the CB electrons [ @Rethfeld:2002 ] and nonlinear pulse propagation effect in a medium [ @Penano:2005; @Petrov:2008; @Gulley:2010; @Apalkov:2012 ] in the fs regime. However, all these studies were based on theoretical models that used experimental and/or empiric value of the material parameter such as ionization pace, deflective index, relaxation rate, and band structure. Hence, the need for developing a method acting that uses non - empirical values of the material parameters grew continuously in the search for a comprehensive and reliable method of investigating laser – matter interactions in the average laser intensity stove. 

 In this newspaper, we employ an alternative method to compute the ocular response and the ablation threshold of an optical culture medium. In contrast to the previous studies, our method is base on first - principles simulation computed from cardinal equations. A multiscale approach using the wave equation and the prison term - dependent density functional theory (TDDFT) is applied to account directly the concentration of the CB electrons generated in the optical medium. The beginning report card on the use of such a multiscale approach for probe of the interaction between a laser pulse and an optical medium was made for crystalline silicon, where it was said to yield reliable resultant role [ @Yabana:2012 ]. In this approach, no empirical parameters and approximations were used except for information on the quartz glass structure and on
) lader pulses [@Wood:2003]. Later, inyestigations of laser-invuced dzmage in the picosecond (ps) and fs rejime gaintb significance when ghe advenn of the xhirked-pulse amplificefion (CPA) techhlque [@Vvrickland:1985] made lt feasible to develop fs avd petawatt-class laser systems [@Sung:2010; @Yt:2012]. In pattlcular, laser aflatpog ocdlrving on the fs time scale became criticel because a laxer pulse duration of few hens of fs is much shogter than tye tyne scale for electron tnzrgy transfgr to the lattice and subsequent lattnce heating. Ib 1995, Djuart $\mathit{xt \, al.}$, investigated the laver-induved damage thrcsholv at $1053 \, \mathrm{nm}$ and $526 \, \matirm{nm}$ for pulse durajions ranghny from $270 \, \mathrm{fs}$ to $1 \, \nathrk{ns}$, dhrojth x tgeprstical movel based oh CB electrin production via mileplhoton ionizztion, Totle heating and collisional ionization [@Vtuzrt:1995]. Subsequent studies vy other groups were fonducted for a more accurate analysis of the damage and atlatikv tkvcshoue hy including energy dependence of the CB elecedoms [@Rethfeld:2002] and konlinear pulse prppwgsjion effect in a mednhm [@Penano:2005; @Petrov:2008; @Gulpey:2010; @Apajkov:2012] un the fs regome. However, all these studiws were baseb ob theoretical modeps that useb expetimentsl and/or empirical valuzs of fhe materiap parametsfs such as ionizxtipn rate, refractive index, relwxation rete, aud band rtruvture. Rence, the jeed nmr developing a mehhod jhat uves non-emplrical values of the material pecameters grew cmntpnuously nn the search for a comprehensive and rejiabld method or invesvigating lasqr–matter inteswctions in tie intermqdiare lqser ingdnsity range.

In this papvr, we wmploy an alternatlve mgtgod to compute cke optical responxe xnd tje afnation thresvold of sn opgical mtbiun. In vontrast to the prevhous studies, our methoc ls based in first-[rinciples sikulations computed from fuidamenyal equations. A multiscale approadh using hhe wave equatijn akd tre time-depzndent density functional theory (TDDFT) iw applied to calcolate directly the deuslty of the CU elecerons genarated in the opticao medium. The firsu report on the use of auch a multlscale approach for investigation of the interaction between a laser pulsw and en optical mesium was kabe yor cryftalpmne silicon, where it was said to yield reliable cesults [@Yataua:2012]. In this approach, no empirivau parameters xnd approximations were hsed excrpt for information on the crysyal structure and on
) laser pulses [@Wood:2003]. Later, investigations of in picosecond (ps) fs regime gained the amplification (CPA) technique made it feasible develop fs and petawatt-class laser systems @Yu:2012]. In particular, laser ablation occurring on the fs time scale became critical a laser pulse duration of few tens of fs is much shorter than time for energy to the lattice and subsequent lattice heating. In 1995, Stuart $\mathit{et \, al.}$, investigated the laser-induced threshold at $1053 \, \mathrm{nm}$ and $526 \, for pulse durations ranging $270 \, \mathrm{fs}$ to $1 \mathrm{ns}$, a theoretical based CB production via multiphoton Joule heating and collisional ionization [@Stuart:1995]. Subsequent studies by other groups were conducted for a more accurate of the ablation threshold including dependence the CB electrons nonlinear pulse propagation effect in a @Gulley:2010; @Apalkov:2012] in the fs regime. However, all studies were on theoretical models that used experimental empirical values of the material parameters such as rate, refractive index, relaxation rate, and band structure. Hence, the need for developing a method non-empirical values of the parameters grew continuously the for comprehensive reliable method investigating laser–matter interactions in the intermediate laser intensity range. In this we employ an alternative method to compute the optical response ablation of an optical In contrast to the studies, method is based on computed fundamental approach the equation and the time-dependent functional theory (TDDFT) is applied calculate directly the density in the optical medium. The first report on use of such a multiscale approach for of the interaction between a laser pulse and an optical medium was for crystalline it was said to yield reliable results [@Yabana:2012]. this approach, no empirical and approximations were used except for information on the structure on
) laser pulses [@Wood:2003]. Later, invesTigations oF laseR-inDucEd DamaGe in The picosecond (pS) And fS regime gained significaNce whEn THe adVEnT of thE chirpeD-PuLSE amPlIfIcaTiON (CpA) tecHniQue [@StriCkland:1985] made It fEaSible to develOP fS and petawaTt-cLass laser sysTemS [@Sung:2010; @YU:2012]. IN paRTiculAr, lAser aBlatioN OccurrIng on the fS tIMe scalE Became cRITiCal bEcause a laser pulse DUrATion of few tens oF fs is mUcH ShORTer ThaN the time scAlE for eLEctron eNErGY TRanSFer to the lattiCe and subseqUEnt LatticE hEatINg. In 1995, StUart $\mAtHIt{eT \, al.}$, investigAted The laser-iNduced DAmage thREshold aT $1053 \, \mathrM{nm}$ And $526 \, \MathRM{nM}$ fOr pUlSE duRAtIonS RanGing from $270 \, \MaThRm{fs}$ tO $1 \, \matHRM{NS}$, thrOugH a thEoretIcal model baseD on cB elECtrOn proDuctiOn viA mUltipHoton iOnizaTiOn, Joule heating aNd coLlisional IonIzAtiOn [@stuarT:1995]. subseqUenT stUdies by Other grOUps WeRE COnDucted for a more accuRaTE AnAlysis of The damAGe AnD Ablation ThResHold BY IncluDing ENeRgy depenDence oF ThE Cb electrOnS [@RethfElD:2002] anD noNlineAR pulSe propAgation eFfect IN a medium [@Penano:2005; @pEtrov:2008; @Gulley:2010; @ApALkOV:2012] In THe fs RegIme. However, aLl thESe stUdieS WeRe bASed on TheorEtICaL Models that used experImEntal aNd/or eMpirical valueS of the mateRIAL parametErs sUCh AS ionization ratE, refrActive indeX, RelaxatiOn ratE, and band Structure. hENce, the neEd fOr dEveLopING a Method that useS NOn-emPiRical vaLueS of the mAteRiaL paRamEtErs grew coNtinuousLy In ThE sEarCh for A CompreheNsIve AnD reLiablE Method Of invEstiGaTiNG laSer–mattER iNTEracTiOnS in tHe iNtErmedIate LAseR intensIty range.

IN thIS papEr, We Employ aN alternative mEtHod to compuTe The OpticaL REsponse aNd the ablation threshold oF An opticAl mEdium. in coNtrast to tHe pReviouS stUDies, ouR methoD is baSeD on FIRst-prINCiPleS sImulations COMpuTed frOm FundAmental Equations. A multiscaLE apProach using thE waVe eqUATiOn aND tHE tiMe-DEpeNDEnt density functIonal theorY (TddFt) is applied TO caLcUlate diRectly tHe denSIty of thE CB electrOns generaTeD in tHE OptIcal medium. the first Report on tHE use oF SuCh a muLtiScale aPpRoaCh for InvestIGatIon of The intErAction BetweEn A laser puLse and an optical medium waS made fOr cryStaLline siliCon, WHerE it was saiD to yIeld reliabLe rEsuLts [@YaBanA:2012]. in thiS appROaCh, nO EmpirIcal PArameters ANd AppROXiMations were USED exCept fOr iNFormatIon oN the crystal structURe and on
) laser pulses [@Wood:2003 ]. Later,inves tig ati on s of  las er-induced dam a ge i n the picosecond (ps)and f sr egim e  g ained  signif i ca n c e w he nthe  a d ve nt of  th e chirp ed-pulse a mpl if ication (CPA )  t echnique [ @St rickland:198 5]made i tfea s ibletodevel op fsa nd pet awatt-cla ss lasers ystems[ @ Su ng:2 010; @Yu:2012]. I n  p a rticular, lase r abla ti o no c cur rin g on the f stimes cale be c am e c rit i cal because a  laser puls e  du ration  o f f e w tens  of f si s m uch shorter  tha n the tim e scal e  for el e ctron e nergytra nsf er t o  t he  la tt i cea nd  su b seq uent lat ti ce  heat ing. I n 1995 , S tuar t $\m athit{et \, a l.} $, i n ves tigat ed th e la se r-ind uced d amage  t hreshold at $10 53 \ , \mathrm {nm }$  an d$526\ , \mat hrm {nm }$ forpulse d u rat io n s ra nging from $270 \,  \ m a th rm{fs}$to $1\ ,\m a thrm{ns} $,  th roug h a the oret i ca l modelbasedo nCB  electr on  produ ct ion  vi a mul t ipho ton io nization , Jou l e heating andc ollisional io n iz a t io n  [@S tua rt:1995]. S ubse q uent  stu d ie s b y  othe r gro up s  w e re conducted for amo re acc urate  analysis ofthe damage a n d ablati on t h re s hold by includ ing e nergy depe n dence of  theCB elect rons [@Re t h feld:200 2]and  no nli n e ar  pulse propag a t ionef fect in  amedium[@P ena no: 200 5;  @Petrov: 2008; @G ul le y: 20 10;  @Apa l kov:2012 ]inth e f s reg i me. Ho wever , al lth e sestudies we r e  bas ed  o n th eor et icalmode l s t hat use d experim ent a l an d/ or  empiri cal values of  t he materia lpar ameter s such asionization rate, refrac t ive ind ex,  rela xati on rate,and  bandstr u cture.  Hence , the  n eed f or de v e lo pin ga method t h a t u ses n on -emp iricalvalues of the mate r ial  parameters g rew  con t i nu ous l yi n t he sea r c h for a compreh ensive and  r e li able metho d  of  i nvestig ating l aser– m atter i nteractio ns in the  i nter m e dia te laser i ntensity  range.
I n thi s  p aper,  we  emplo yanalter native met hod t o comp ut e theoptic al  respons e and the ablation thre sholdof an  op tical med ium .  In  contrast  tothe previo usstu dies,  ou r  meth od i s  b ase d  on f irst - principle s  s imu l a ti ons compute d f rom  fund ame n tal eq uati ons. A multiscale approach using  the w ave  eq u atio nand the time-d epe nd e n t densit yfunctionaltheory ( TD D FT) i s appl ied to  calcul a t ed irectl y th e d ensity of  th eC B elect ro ns genera tedin  the o ptical medi u m . The first repo rt on t he us e  of  such  a  multis c aleapproach f or investig ationof t he in teracti on  betwe enalaser puls e  and an o ptica l mediu mwasmad e forcrys t a lline  sil ic on,  where it w as sa id toyiel d rel ia bleresults [ @ Yabana:2 012 ] . In th is  ap p r oach,n oe m pirical pa ram eters a nd approxi m atio n sw ere u sed ex cept f or info r mat io n on th e c r y stal stru cture and on
) laser_pulses [@Wood:2003]._Later, investigations of laser-induced_damage in_the_picosecond (ps)_and_fs regime gained_significance when the_advent of the chirped-pulse_amplification (CPA) technique_[@Strickland:1985]_made it feasible to develop fs and petawatt-class laser systems [@Sung:2010; @Yu:2012]. In particular,_laser_ablation occurring_on_the_fs time scale became critical_because a laser pulse duration_of few_tens of fs is much shorter than the_time_scale for electron_energy transfer to the lattice and subsequent lattice heating._In 1995, Stuart $\mathit{et \, al.}$,_investigated the laser-induced_damage_threshold_at $1053 \, \mathrm{nm}$_and $526 \, \mathrm{nm}$ for pulse_durations ranging from $270 \, \mathrm{fs}$_to $1 \, \mathrm{ns}$, through a theoretical_model based on CB electron production_via multiphoton ionization, Joule heating_and collisional_ionization [@Stuart:1995]. Subsequent studies by_other groups were_conducted for_a more accurate_analysis of the damage and ablation_threshold by including_energy dependence of the CB electrons_[@Rethfeld:2002]_and nonlinear pulse_propagation_effect_in a_medium [@Penano:2005; @Petrov:2008;_@Gulley:2010;_@Apalkov:2012] in_the_fs regime. However, all these studies_were_based on theoretical models that used experimental_and/or empirical values of_the_material parameters such as_ionization rate, refractive index, relaxation_rate, and band structure. Hence, the_need for_developing a_method that uses non-empirical values of the material parameters grew continuously_in the search for a comprehensive_and reliable method of_investigating laser–matter_interactions_in the intermediate_laser_intensity range.

In_this paper, we employ an alternative method_to compute_the optical response and the ablation_threshold of an optical_medium._In contrast to the previous studies,_our method is based on first-principles_simulations computed from fundamental equations._A_multiscale_approach using the wave equation_and the time-dependent density functional theory_(TDDFT) is applied_to calculate directly the density of the_CB_electrons generated in the optical medium._The_first report on the use of_such_a_multiscale approach for investigation of_the interaction between a laser pulse_and an optical medium was made for crystalline silicon,_where it was_said to yield reliable results_[@Yabana:2012]._In_this approach, no empirical parameters and approximations were used except_for information_on the crystal_structure and on