Esempio n. 1
0
void SemiconductorPhase::getChemPotentials(doublereal* mu) const
{
    getActivityConcentrations(DATA_PTR(m_work));
    doublereal r = m_work[0]/nc();
    mu[0] = ec() + GasConstant*temperature()*(JoyceDixon(r));
    mu[1] = ev() + GasConstant*temperature()*(log(m_work[1]/nv()));
}
Esempio n. 2
0
void SemiconductorPhase::getChemPotentials(doublereal* mu) const
{
    getActivityConcentrations(m_work.data());
    mu[0] = ec() + RT()*(JoyceDixon(m_work[0]/nc()));
    mu[1] = ev() + RT()*(log(m_work[1]/nv()));
}