Kris's Research Notes

July 6, 2011

Droplet Statistics and Fast Adatom Diffusion

Filed under: GaAs Simulations — Kris Reyes @ 4:34 pm

In this note, we examine the effect of Ga-Ga bonds on droplet statistics and describe changes to the code that speed up adatom diffusion.

Droplet Statistics

During the June 23 meeting, we discussed the results of a large run that looped over several values for the five energy parameters \gamma_{GG}, \gamma_{GA}, \gamma_{HH}, \lambda_{D}, \mu_{AS}, which affect G-G bonds, G-A bonds, atom-atom exchanges, and diffusion through the droplet and As desorption, respectively. The results of the run can be found here. Observe that the size of the droplets smaller than what is expected in experiments. For example, here is a crystallized droplet formed by depositing Ga at a rate of 0.1 monolayers/second for 35 seconds then exposed to an As flux of 1 monolayer/second at a temperature of 473K with parameter values \gamma_{GG} = 0.3 eV, \gamma_{GA} = 0.5 eV, \gamma_{HH} = -0.6 eV, \lambda_{D} = 0.7 eV and \mu_{AS} = 0.4 eV:

The domain is 512 atoms wide, and the droplet is at most 120 atoms wide. However, in experiments (see Z. Gong et al., Appl. Phys. Lett. 87, 093116 (2005) ) droplets/nanorings of diameter 150 nm (about 600 atoms) form at 423K, much larger than the above example.

To address this, we consider the effect of \gamma_{GG} and \gamma^\prime_{GG} on the size of the droplets. The idea is that on a Ga-terminated surface, by lowering these energies, diffusing Ga adatoms are more likely to find each other over a larger domain. We deposit 3.5 monolayers of Ga on a domain of width 4096 atoms at a temperature of 423K. Using the values in Gong et al. (namely the average diameter is 360 atoms and average height 24 atoms) is and assuming the droplets are elliptical in shape, we expect each droplet to have a mass of about 6786 atoms, half of which are Ga atoms. Since we deposit 14336 Ga atoms total, we want between 4 and 5 droplets to form in our domain, resulting in a droplet density between 3.9 \times 10^4 cm^{-1} and 4.8\times 10^4 cm^{-1} (compare this to the experimental density of 2.5 \times 10^9 cm^{-2}).

We vary \gamma_{GG} \in \left\{ 0.20, 0.22, \hdots, 0.30\right\} eV and \gamma^\prime_{GG} \in \left\{ 0.1, 0.2, 0.3\right\} eV. For each pair of parameters, Ga is deposited at a rate of 0.1 monolayers/second for 35 seconds with substrate temperature of 423K. Statistics are collected over an ensemble of 16 identical runs. Widths are determined using the height profile auto-correlation function. (All results here)

First is a plot droplet width (which we want to be 360/2 = 180 atoms) vs. \gamma_{GG}. The different lines correspond to different values of \gamma^\prime_{GG}:

Second is a plot of number of droplets (which we want around 4 or 5) vs. \gamma_{GG}.

These two plots suggest parameter values of \gamma_{GG}, \gamma^\prime_{GG} on the low side of values tested. Specifically,

  • if \gamma^\prime_{GG} = 0.1 eV (red), then 0.2 \leq \gamma_{GG} \leq 0.28 eV;
  • if \gamma^\prime_{GG} = 0.2 eV (blue), then 0.2 \leq \gamma_{GG} \leq 0.24 eV;
  • if \gamma^\prime_{GG} = 0.3 eV (green), then 0.2 \leq \gamma_{GG} \leq 0.22 eV.

Here is a histogram of how much Ga (monolayers) was actually deposited in the 35 seconds:

Note the standard deviation of 0.06 monolayers Ga. We can repeat the same experiment, this time fixing the amount of Ga deposited to exactly 3.5 monolayers (results here). Here is droplet width vs \gamma_{GG} in that case:

Number of droplets vs \gamma_{GG}:

The behavior is similar to the previous case.

Fast Adatom Diffusion

In order for the above simulations to run in a reasonable amount of time, we changed the code to allow for fast adatom diffusion. This means that when an adatom hops, it may hop several one direction before any bookkeeping is performed. An adatom hops in one direction until its local environment is different from the one it started out in, or until 10 hops are performed. To correct for this, we simply sample adatom hops 10 times less than before. This has a significant impact on the speed of the simulation, but also affects the diffusivity of adatoms in a non-trivial way.

To illustrate this, fix \gamma_{GG} = 0.30, \gamma^\prime_{GG} = 0.10, \gamma_{GA} = 0.6, \gamma_{HH} = -1.0, \mu_{AS} = 0.4, \lambda_{D} = 0.6 eV and consider the result of depositing 1 monolayer of Ga at a rate of 0.1 monolayer/second. Here is the result with the above fast-adatom diffusion mechanism turned on. Note that there are nine Ga droplets formed. The simulation performed 194,749,977 updates. Here is the result with the fast diffusion turned off. In this case, thirteen droplets have begun to form. Here, the simulation performed 624,706,417 updates, about three times more than when fast adatom diffusion was turned on. We conclude that the Ga adatoms diffuse more rapidly with fast adatom diffusion turned onn even when the rates are scaled by a factor of 0.1.


Leave a Comment »

No comments yet.

RSS feed for comments on this post. TrackBack URI

Leave a Reply

Fill in your details below or click an icon to log in: Logo

You are commenting using your account. Log Out /  Change )

Google+ photo

You are commenting using your Google+ account. Log Out /  Change )

Twitter picture

You are commenting using your Twitter account. Log Out /  Change )

Facebook photo

You are commenting using your Facebook account. Log Out /  Change )


Connecting to %s

%d bloggers like this: