 $\rho_0$. (a)-(f) correspond to $\rho_0=0.06$, $0.08$, $0.10$, $0.12$, $0.14$ and $0.16$. []{data-label="fig_PFC_vac"}](fig3f "fig:"){width="0.33\columnwidth"}
  -------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------

With the vacancy term, Eq. (\[eqn\_vac\_f\]), we can numerically verify the analytical calculation for the coexistence between the periodic phase and vacancies. Fig. (\[fig\_PFC\_vac\]) provides results from simulations with $r=-0.9$ and different values of $\rho_0$, showing clearly that the number of atoms increases with $\rho_0$. In addition, Fig. (\[fig\_pfc\_density\]) shows that the PFC atomic density (*i.e.,* the number of atoms per unit area.) indeed increases linearly with $\rho_0$, as expected. The curve starts to saturate at around $\rho_0=0.15$, as opposed to the prediction from Eq (\[eqn\_vac\_phase\]), $\rho_{min} \approx 0.134$. This discrepancy is not surprising for several reasons: In the analytical calculation, we consider only the one-mode approximation in the ansatz; we did not account for the surface energy between the triangular phase and the vacancies; and we did not account for thermal fluctuations, introduced in the simulation to help the system equilibrate faster.

![ The PFC atomic density increases linearly with the order parameter, $\rho_0$, when the vacancy term is added to the model. $r=-0.9$ is used. The curve starts to saturate at around $\rho_0=0.15$, as opposed to the analytical prediction $\rho_0=\sqrt{(-48-56r)/133}=0.134$. []{data-label="fig_pfc_density"}](fig4){width="0.75\columnwidth"}

With the modifications described above,the PFC simulation operates very much like a molecular dynamics simulation, but on diffusive time scales many orders of magnitude faster than pure molecular dynamics[@Elder04]. We can control the number of atoms and the temperature in the system by adjusting $\rho_0$ and the magnitude of thermal noise, $\eta$, respectively. The interaction potential between individual PFC atoms is specified by the PFC free energy
$ \rho_0$. (a)-(f) correspond to $ \rho_0=0.06 $, $ 0.08 $, $ 0.10 $, $ 0.12 $, $ 0.14 $ and $ 0.16$. [ ] { data - label="fig_PFC_vac"}](fig3f " fig:"){width="0.33\columnwidth " } 
   ------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------- 

 With the vacancy term, Eq.   (\[eqn\_vac\_f\ ]), we can numerically control the analytic calculation for the coexistence between the periodic phase and vacancies. Fig.   (\[fig\_PFC\_vac\ ]) provides results from simulation with $ r=-0.9 $ and different values of $ \rho_0 $, show distinctly that the number of atom increases with $ \rho_0$. In addition, Fig.   (\[fig\_pfc\_density\ ]) show that the PFC atomic density (* i.e., * the phone number of atom per unit area .) indeed increases linearly with $ \rho_0 $, as expected. The curve start to saturate at around $ \rho_0=0.15 $, as opposed to the prediction from Eq   (\[eqn\_vac\_phase\ ]), $ \rho_{min } \approx 0.134$. This discrepancy is not surprising for several rationality: In the analytical calculation, we consider only the one - mode estimate in the ansatz; we did not account for the airfoil energy between the triangular phase and the vacancies; and we did not account for thermal fluctuation, introduced in the simulation to help the system equilibrate faster. 

! [ The PFC atomic density increases linearly with the order parameter, $ \rho_0 $, when the vacancy term is added to the model. $ r=-0.9 $ is used. The curve start to saturate at around $ \rho_0=0.15 $, as opposed to the analytic prediction $ \rho_0=\sqrt{(-48 - 56r)/133}=0.134$. [ ] { datum - label="fig_pfc_density"}](fig4){width="0.75\columnwidth " } 

 With the change described above, the PFC simulation operate on very much like a molecular dynamics simulation, but on diffusive meter scales many orders of magnitude faster than pure molecular dynamics[@Elder04 ]. We can control the phone number of atoms and the temperature in the system by adjusting $ \rho_0 $ and the magnitude of thermal randomness, $ \eta$, respectively. The interaction potential between individual PFC atom is assign by the PFC free energy
 $\rhl_0$. (a)-(f) correspond to $\rho_0=0.06$, $0.08$, $0.10$, $0.12$, $0.14$ and $0.16$. []{data-labgl="dig_PFC_tac"}](fig3f "fig:"){widtf="0.33\columnwidth"}
  -------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------

With the vacanry twrm, Ew. (\[eqn\_vac\_f\]), we can numerizally verpfy the abalyuical calculation for the coexiabence uetween the perlodic phase and vacancies. Fkg. (\[yig\_PFC\_vac\]) provides results from simujations wlth $r=-0.9$ and diffgrent daluss of $\rho_0$, showing clearly that the number of atoms incteases with $\rho_0$. In additioj, Fih. (\[fig\_pfc\_density\]) shoas that the PFC qtomic densigy (*i.e.,* the number of atkms per unit area.) indeed increares lnnearly wity $\eho_0$, ds expected. The burve starts bp satusate at around $\rho_0=0.15$, as op'osee to the prediction fcom Eq (\[eqn\_vac\_phase\]), $\rhj_{min} \apprmx 0.134$. This discrepanct us noj sur[riskbg wor sxvedal rewsois: In the ahalytical cqlculation, we consicew only the one-jode a[pwoximation in the ansatz; we did not accmunf for the surface energt between the triangupar phase and the vacancies; and we did not account for thesmal hljctbqtionr, ijtroduced in the simulation to help the systei eauplibrate faster.

![ Tme PFC atomic densotj ogcreases linexrly wnfh the order parametfr, $\rho_0$, rhen rhe vacansy trrm is added to the model. $r=-0.9$ is used. The xurve starts to sacurate at aruund $\rho_0=0.15$, ss opposed to the analycical lrediction $\gho_0=\sqrt{(-48-56r)/133}=0.134$. []{sxta-label="fig_pfc_devsiny"}](fic4){width="0.75\columnwidth"}

With the iodificatmons bescribea abpve,the PFC simulwtion operates very much like a molecular dynamics simulation, but on difhnsive time scslas kany ordzrs of magnitude fafter than pure molecujar dhnamics[@Eldvr04]. We can control thq number of adlms and the vemperatuwe ib thw systeo by adjusting $\tho_0$ and tkt magnitudw of thermal noise, $\etx$, respectively. Tkt ibteraction potemtixl fenwexn inqhvidual PFC dtomr ir speckfied by thc PWC ftee energy
$\rho_0$. (a)-(f) correspond to $\rho_0=0.06$, $0.08$, $0.10$, and []{data-label="fig_PFC_vac"}](fig3f "fig:"){width="0.33\columnwidth"} With the vacancy numerically the analytical calculation the coexistence between periodic phase and vacancies. Fig. (\[fig\_PFC\_vac\]) results from simulations with $r=-0.9$ and different values of $\rho_0$, showing clearly that number of atoms increases with $\rho_0$. In addition, Fig. (\[fig\_pfc\_density\]) shows that the atomic (*i.e.,* number atoms per unit area.) indeed increases linearly with $\rho_0$, as expected. The curve starts to saturate around $\rho_0=0.15$, as opposed to the prediction from (\[eqn\_vac\_phase\]), $\rho_{min} \approx 0.134$. discrepancy is not surprising for reasons: the analytical we only one-mode approximation in ansatz; we did not account for the surface energy between the triangular phase and the vacancies; and did not thermal fluctuations, in simulation help the system ![ The PFC atomic density increases order parameter, $\rho_0$, when the vacancy term is to the $r=-0.9$ is used. The curve starts saturate at around $\rho_0=0.15$, as opposed to the prediction $\rho_0=\sqrt{(-48-56r)/133}=0.134$. []{data-label="fig_pfc_density"}](fig4){width="0.75\columnwidth"} With the modifications described above,the PFC simulation operates very much like a simulation, but on diffusive scales many orders magnitude than molecular We can the number of atoms and the temperature in the system by $\rho_0$ and the magnitude of thermal noise, $\eta$, respectively. The between PFC atoms is by the PFC free
 $\rho_0$. (a)-(f) correspond to $\rho_0=0.06$, $0.08$, $0.10$, $0.12$, $0.14$ and $0.16$. []{daTa-label="fig_pFC_vaC"}](fiG3f "fIg:"){WidtH="0.33\colUmnwidth"}
  -------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------

With thE VacaNcy term, Eq. (\[eqn\_vac\_f\]), we can nUmeriCaLLy veRIfY the aNalyticAL cALCulAtIoN foR tHE cOexisTenCe betweEn the perioDic PhAse and vacancIEs. fig. (\[fig\_PFC\_vAc\]) pRovides resulTs fRom simUlAtiONs witH $r=-0.9$ aNd difFerent VAlues oF $\rho_0$, showiNg CLearly THat the nUMBeR of aToms increases with $\RHo_0$. iN addition, Fig. (\[fiG\_pfc\_deNsITy\]) SHOws ThaT the PFC atoMiC densITy (*i.e.,* the NUmBER Of aTOms per unit areA.) indeed incrEAseS lineaRlY wiTH $\rho_0$, as ExpecTeD. the Curve starts To saTurate at aRound $\rHO_0=0.15$, as oppoSEd to the PredicTioN frOm Eq (\[EQn\_VaC\_phAsE\]), $\Rho_{MIn} \AppROx 0.134$. THis discrEpAnCy is nOt suRPRISing For SeveRal reAsons: In the anaLytIcal CAlcUlatiOn, we cOnsiDeR only The one-Mode aPpRoximation in the AnsaTz; we did noT acCoUnt FoR the sURface eNerGy bEtween tHe trianGUlaR pHASE aNd the vacancies; and wE dID NoT account For theRMaL fLUctuatioNs, IntRoduCED in thE simULaTion to heLp the sYStEm EquilibRaTe fastEr.

![ the pFC AtomiC DensIty incReases liNearlY With the order paRAmeter, $\rho_0$, when THe VACaNCy teRm iS added to the ModeL. $R=-0.9$ is uSed. THE cUrvE StartS to saTuRAtE At around $\rho_0=0.15$, as opposeD tO the anAlytiCal prediction $\Rho_0=\sqrt{(-48-56r)/133}=0.134$. []{daTA-LAbel="fig_pFc_deNSiTY"}](fig4){width="0.75\columNwidtH"}

With the moDIficatioNs desCribed abOve,the PFC SIMulation OpeRatEs vEry MUCh Like a moleculaR DYnamIcS simulaTioN, but on dIffUsiVe tIme ScAles many oRders of mAgNiTuDe FasTer thAN pure molEcUlaR dYnaMics[@ELDer04]. We cAn conTrol ThE nUMbeR of atomS AnD THe teMpErAturE in ThE systEm by ADjuSting $\rhO_0$ and the maGniTUde oF tHeRmal noiSe, $\eta$, respectiVeLy. The interAcTioN potenTIAl betweeN individual PFC atoms is spECified bY thE PFC fRee eNergy
 $\rho_0$. (a)-(f) corresp ond to $\r ho_0= 0.0 6$,  $ 0.08 $, $ 0.10$, $0.12$, $0.1 4$ and $0.16$. []{data -labe l= " fig_ P FC _vac" }](fig3 f  " f i g:" ){ wi dth =" 0 .3 3\col umn width"} 
  ------- --- -- ------------ - -- ---------- --- ------------ --- ------ -- --- - ----- --- ----- ------ - ------ --------- -- - ------ - ------- - - -- ---- ----------------- - -- - -------------- ------ -- - -- - - --- --- ---------- -- ----- - ------- - -- - - - --- - ------------- ----------- - --- ------ -- --- - ------ ----- -- - --- ----------- ---- --------- ------ - ------- - ------- ------ --- --- ---- - -- -- --- -- - --- - -- --- - --- -------- -- -- ----- ---- - - - - ---- --- ---- ----- ------------- --- ---- - --- ----- ----- ---- -- ----- ------ ----- -- --------------- ---- --------- --- -- --- -- ----- - ------ --- --- ------- ------- - --- -- - - - -- ------------------ -- - - -- -------- ------ - -- -- - -------- -- --- ---- - - ----- ---- - -- -------- ------ - -- -- ------- -- ------ -- --- --- ----- - ---- ------ -------- ----- - -------------- - ------------- - -- - - -- - ---- --- ----------- ---- - ---- ---- - -- --- - ----- ----- -- - -- - -----------

With t he  vacan cy te rm, Eq. (\[eq n\_vac\_f\ ] ) ,  we cannume r ic a lly verify the  anal ytical cal c ulationfor t he coexi stence be t w een theper iod icpha s e  a nd vacancies. F ig.(\ [fig\_P FC\ _vac\])  pr ovi des  re su lts fromsimulati on swi th  $r =-0.9 $  and dif fe ren tval ues o f  $\rho _0$,show in gc lea rly tha t  t h e  num be rof a tom sincre ases wit h $\rho _0$. In a ddi t ion,  F ig . (\[fi g\_pfc\_densi ty \]) showsth atthe PF C atomic d ensity (*i.e.,* the num b er of a tom s per  uni t area.)ind eed in cre a ses li nearly  with  $ \rh o _ 0$, a s ex pec te d. The cur v e  st artsto  sat urate a t around $\rho_0=0 . 15$ , as opposedtothep r ed ict i on fro mE q ( \ [ eqn\_vac\_phase \]), $\rho _{ m in } \approx0 .13 4$ . Thisdiscrep ancyi s not s urprising  for seve ra l re a s ons : In the a nalytica l calcula t ion,w econsi der  onlyth e o ne-mo de app r oxi matio n in t he  ansat z; we  d id not a ccount for the surfaceenergy  betw een  the tria ngu l arphase and  the  vacancies ; a ndwe di d n o t acc ount fo r t h ermal  flu c tuations, in tro d u ce d in the si m u l ati on to  he l p thesyst em equilibrate fa s ter.

![ The P FC a t o mic  de n sity  i ncreases linea rly  w i t h the or de r parameter , $\rho_ 0$ ,  when  the v acancy  term i s ad d ed tothemod el. $r=-0 .9$  i s  used.Th ec urve s tart sto sat uratea t ar o u nd $\rho_0=0.15$ , aso p posed tothe a na lytical pred iction $\r ho_0=\sqrt{ (-48-5 6r)/ 133}= 0.134$.  [ ]{data -la be l="fig_pfc _ density"} ](fig 4){widt h= "0.7 5\c olumnw idth " } 

Wit h th emod ification s de s cr ib e d a bove ,thePF C si mulationo peratesver y  much l ik e a m olecul a rd y namics sim ula tion, b ut on diff u sive ti m e sca les ma ny ord ers ofm agn it ude fas ter t han puremolecular dyn am ics[ @Elder04 ].  Weca n c on trol the number  of atom s and th e temper atu r e in the sy st em by adjust ing $ \r h o_0$ an d the  magnit u deofthe rmal  no ise, $\eta$, respectively. Th e in terac t ion  pote ntial  b etween i n divi dual PFC atoms is spec ified by the PF C  fr e e en erg y 
 $\rho_0$._(a)-(f) correspond_to $\rho_0=0.06$, $0.08$, $0.10$,_$0.12$, $0.14$_and_$0.16$. []{data-label="fig_PFC_vac"}](fig3f_"fig:"){width="0.33\columnwidth"}
_ -------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------

With the_vacancy term, Eq. (\[eqn\_vac\_f\]),_we can numerically verify_the analytical calculation_for_the coexistence between the periodic phase and vacancies. Fig. (\[fig\_PFC\_vac\]) provides results from simulations with_$r=-0.9$_and different_values_of_$\rho_0$, showing clearly that the_number of atoms increases with_$\rho_0$. In_addition, Fig. (\[fig\_pfc\_density\]) shows that the PFC atomic density_(*i.e.,*_the number of_atoms per unit area.) indeed increases linearly with $\rho_0$,_as expected. The curve starts to_saturate at around_$\rho_0=0.15$,_as_opposed to the prediction_from Eq (\[eqn\_vac\_phase\]), $\rho_{min} \approx 0.134$. This_discrepancy is not surprising for several_reasons: In the analytical calculation, we consider_only the one-mode approximation in the_ansatz; we did not account_for the_surface energy between the triangular_phase and the_vacancies; and_we did not_account for thermal fluctuations, introduced in_the simulation to_help the system equilibrate faster.

![ The_PFC_atomic density increases_linearly_with_the order_parameter, $\rho_0$, when_the_vacancy term_is_added to the model. $r=-0.9$ is_used._The curve starts to saturate at around_$\rho_0=0.15$, as opposed to_the_analytical prediction $\rho_0=\sqrt{(-48-56r)/133}=0.134$. []{data-label="fig_pfc_density"}](fig4){width="0.75\columnwidth"}

With_the modifications described above,the PFC_simulation operates very much like a_molecular dynamics_simulation, but_on diffusive time scales many orders of magnitude faster than pure_molecular dynamics[@Elder04]. We can control the_number of atoms and_the temperature_in_the system by_adjusting_$\rho_0$ and_the magnitude of thermal noise, $\eta$, respectively._The interaction_potential between individual PFC atoms is_specified by the PFC_free_energy