Medidas da Amostra de DPPH⁽¹⁾

Para expressar o espectro de EPR em termos do fator g (gFactor) efetivo é preciso realizar a seguinte conversão para cada ponto: \[ gFactor=\frac{714.55*v(GHz)}{B(G)} \] As medidas da amostra de 5 mg de DPPH sem a utilização do Teslameter foram realizadas nos dias 11 de fevereiro de 2014 e 11 de novembro de 2016. Foram utilizados os mesmos parâmetros para a configuração do equipamento. A amostra foi posicionada no centro da cavidade.

(1) a amostra de DPPH será caracterizada por outras técnicas.

Gráficos:

Espectros das duas medidas:

#echo=FALSE para nao publicar o cod
# primeira medida
dsc <- readLines("dpph5mg.DSC",encoding="UTF-8")
info<-grep("XPTS",dsc,value=T);campos<-strsplit(info,'\t')[[1]];XPTS<-as.numeric(campos[2])
info<-grep("XMIN",dsc,value=T);campos<-strsplit(info,'\t')[[1]];XMIN<-as.numeric(campos[2])
info<-grep("XWID",dsc,value=T);campos<-strsplit(info,'\t')[[1]];XWID<-as.numeric(campos[2])
info<-grep("MWFQ",dsc,value=T);campos<-strsplit(info,'    ')[[1]];MF<-as.numeric(campos[2])/1e9
info<-grep("DATE", dsc,value=T);campos<-strsplit(info,'    ')[[1]];data1<-(campos[2])
passo<-round(XWID/XPTS,1)
b<-seq(1:XPTS)
gg<-seq(1:XPTS)
for(i in 1:XPTS){
  b[i]<-XMIN+(passo*(i-1))
}
b<-data.frame(b)
g<-714.55*MF/b
s<-readBin("dpph5mg.DTA", double(), n = XPTS, size = 8, endian = "big");
s<-data.frame(s);
dados1<-cbind(b,g,s);colnames(dados1)[2] <- "g"
###
# segunda medida
dsc <- readLines("dpph5mg_11nov2016.DSC",encoding="UTF-8")
info<-grep("XPTS",dsc,value=T);campos<-strsplit(info,'\t')[[1]];XPTS<-as.numeric(campos[2])
info<-grep("XMIN",dsc,value=T);campos<-strsplit(info,'\t')[[1]];XMIN<-as.numeric(campos[2])
info<-grep("XWID",dsc,value=T);campos<-strsplit(info,'\t')[[1]];XWID<-as.numeric(campos[2])
info<-grep("MWFQ",dsc,value=T);campos<-strsplit(info,'    ')[[1]];MF<-as.numeric(campos[2])/1e9
info<-grep("DATE", dsc,value=T);campos<-strsplit(info,'    ')[[1]];data2<-(campos[2])
passo<-round(XWID/XPTS,1)
b<-seq(1:XPTS)
g<-seq(1:XPTS)
for(i in 1:XPTS){
  b[i]<-XMIN+(passo*(i-1))
}
b<-data.frame(b)
g<-714.55*MF/b
s<-readBin("dpph5mg_11nov2016.DTA", double(), n = XPTS, size = 8, endian = "big");
s<-data.frame(s);
dados2<-cbind(b,g,s);colnames(dados2)[2] <- "g"
###
#plots
campo1<-dados1[which.max(dados1[,2]),1]
campo2<-dados2[which.max(dados2[,2]),1]
dif<-campo1-campo2
#
plot(dados1$b,dados1$s,type='l',xlim=c(3450,3550),ylim=c(-3.5e6,3.5e6),xlab='',ylab='',axes=FALSE,col='blue')
par(new=T)

plot(dados2$b,dados2$s,type='l',xlim=c(3450,3550),ylim=c(-3.5e6,3.5e6),xlab="Magnetic Field (G)",ylab="EPR Signal (a.u.)",col="red")
par(new=F)

grid (NULL,NULL, lty = 4, col = "black") 
legend(3530, 3e6, legend=c(data1, data2), col=c("blue", "red"), lty=1:1, cex=0.75)

Em termos do fator g:

gmax1<-dados1[which.max(dados1[,3]),2]
gmax2<-dados2[which.max(dados2[,3]),2]
gdif<-round(100*(gmax2-gmax1)/gmax1,3)
plot(dados1$g,dados1$s,type='l',xlim=c(2.04,1.96),ylim=c(-3.5e6,3.5e6),xlab='',ylab='',axes=FALSE,col="blue")
par(new=T)

plot(dados2$g,dados2$s,type='l',xlim=c(2.04,1.96),ylim=c(-3.5e6,3.5e6),xlab='gFactor',ylab='EPR Signal (a.u.)',col="red")
par(new=F)

legend(1.977,2.6e6, legend=c(data1, data2), col=c("blue","red"), lty=1:1,cex=0.75)
text(2.015, 2.6e6, bquote(paste(Delta,"g = ", .(gdif), "%")))

hpp1<-max(dados[,3])-min(dados[,3])
hpp2<-max(dados2[,3])-min(dados2[,3])
hdif<-round(100*(hpp2-hpp1)/hpp1,3)
text(2.015, -1.6e6, bquote(paste(Delta,"h = ", .(hdif), "%")))

Fator g do primeiro pico:

• 02/11/14 = 2.0032281
  
• 11/11/16 = 2.0034934