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())); }
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())); }