program compound implicit none C comment: C C 2024.07.10 Y.Satou C C Calculate the mass fraction of each element in G10-FR4 : LBL C Calculate I_eff etc. of G10-FR4 : LBL C C 1 : E_Glass, mass fraction = 0.60 C 2 : Epotek301-1, mass fraction = 0.40 C cf) http://personalpages.to.infn.it/~tosello/EngMeet/ITSmat/SDD/SDD_G10FR4.html C C In E_Glass, TiO2, Na2O, K2O, Fe2O3, F2 components are ignored, C as in the LBL PDG Web page. C C constant: real*8 A_H,A_C,A_O,A_N ! Atomic weight parameter(A_H=1.008D0,A_C=12.0107D0,A_O=15.9990D0,A_N=14.0070D0) !LBL real*8 A_Si,A_Ca,A_Al,A_Mg,A_B,A_Ti,A_Na,A_K,A_Fe,A_F parameter(A_Si=28.0855D0,A_Ca=40.078D0,A_Al=26.9815D0, & A_Mg=24.3050D0,A_B=10.8100D0,A_Ti=47.867D0, & A_Na=22.9898D0,A_K=39.0983D0,A_Fe=55.8450D0, & A_F=18.9984D0) real*8 Z_H,Z_C,Z_O,Z_N ! Atomic charge parameter(Z_H=1.D0,Z_C=6.D0,Z_O=8.D0,Z_N=7.D0) real*8 Z_Si,Z_Ca,Z_Al,Z_Mg,Z_B,Z_Ti,Z_Na,Z_K,Z_Fe,Z_F parameter(Z_Si=14.D0,Z_Ca=20.D0,Z_Al=13.D0, & Z_Mg=12.D0,Z_B=5.D0, Z_Ti=22.D0, & Z_Na=11.D0,Z_K=19.D0, Z_Fe=26.D0, & Z_F=9.D0) real*8 w_1,w_2 ! Mass fraction parameter(w_1=0.60D0,w_2=0.40D0) C real*8 I_H,I_C,I_O,I_N ! H:liquid,C:amorphous parameter(I_H=21.8D0,I_C=81.0D0,I_O=95.0D0,I_N=82.0D0) real*8 I_Si,I_Ca,I_Al,I_Mg,I_B,I_Ti,I_Na,I_K,I_Fe,I_F parameter(I_Si=173.0D0,I_Ca=191.0D0,I_Al=166.0D0, & I_Mg=156.0D0,I_B=76.0D0,I_Ti=233.0D0, & I_Na=149.0D0,I_K=190.0D0,I_Fe=286.0D0, & I_F=115.0D0) real*8 L_H,L_C,L_O,L_N parameter(L_H=63.05D0,L_C=42.70D0,L_O=34.24D0,L_N=37.99D0) real*8 L_Si,L_Ca,L_Al,L_Mg,L_B,L_Ti,L_Na,L_K,L_Fe,L_F parameter(L_Si=21.82D0,L_Ca=16.14D0,L_Al=24.01D0, & L_Mg=25.03D0,L_B=52.68D0,L_Ti=16.16D0, & L_Na=27.74D0,L_K=17.32D0,L_Fe=13.84D0, & L_F=32.93D0) C real*8 I_E_Glass,I_Epotek301_1 ! Mean ionization potential CC parameter(I_E_Glass=143.86D0,I_Epotek301_1=67.94D0) CC parameter(I_E_Glass=143.4D0, I_Epotek301_1=76.7D0) !LBL CC parameter(I_E_Glass=143.90D0,I_Epotek301_1=76.7D0) !LBL(Epoteck) C C variables: real*8 N_1,N_2 ! Atomic fraction real*8 A_eff,Z_eff real*8 N_H,N_C,N_O,N_N real*8 N_Si,N_Ca,N_Al,N_Mg,N_B real*8 w_H,w_C,w_O,w_N real*8 w_Si,w_Ca,w_Al,w_Mg,w_B C real*8 A1_eff,Z1_eff real*8 N1_Si,N1_Ca,N1_Al,N1_Mg,N1_B, & N1_O real*8 A2_eff,Z2_eff real*8 N2_H,N2_C,N2_O,N2_N C real*8 A_eff_d,Z_eff_d,I_eff_alt,ln_I_eff_alt C real*8 I_eff,ln_I_eff real*8 L_rad,recip_L_rad real*8 A_avr,A_avr_d,Z_avr C begin: call E_Glass(A1_eff,Z1_eff,N1_Si,N1_Ca,N1_Al,N1_Mg,N1_B,N1_O, & I_E_Glass) call Epotek301_1(A2_eff,Z2_eff,N2_H,N2_C,N2_O,N2_N,I_Epotek301_1) C I_Epotek301_1=76.7D0 C write(*,'(A)')'***********************************************'// &'*********************************' write(*,'(A)')'G10-FR4:' write(*,'(A)')'Calculate various parameters from '// & 'mass fractions of elements (LBL-like way)' write(*,'(A,F6.2,A)')' with I_E_Glass= ',I_E_Glass, ' (ev)' write(*,'(A,F6.2,2A)') & ' with I_Epotek301_1=',I_Epotek301_1,' (ev)', & ' (<-- From LBL Epotk301-1)' write(*,'(A)')'TiO2, Na2O, K2O, Fe2O3, F2 components'// & ' are ignored.' write(*,'(A)')'***********************************************'// &'*********************************' write(*,'(A,F7.3)')'Fraction by weight ([E_Glass])= ',w_1 write(*,'(A,F7.3)')'Fraction by weight ([Epotek301-1])= ',w_2 C write(*,'(A,F10.6)')'Mass of E_Glass (arb-unit))= ', & A1_eff write(*,'(A,F10.6)')'Mass of Epotek301-1 (arb-unit))= ', & A2_eff C N_2=1.D0 A_eff_d = A2_eff/w_2 N_1 = w_1*A_eff_d/A1_eff C write(*,'(A,F10.6)')'Molecular fraction ([E_Glass])= ',N_1 write(*,'(A,F10.6)')'Molecular fraction ([Epotek301-1])= ',N_2 C Z_eff_d = N_1*Z1_eff + N_2*Z2_eff ln_I_eff_alt = (N_1 *Z1_eff *log(I_E_Glass) & + N_2 *Z2_eff *log(I_Epotek301_1))/Z_eff_d I_eff_alt = exp(ln_I_eff_alt) write(*,'(A18,F10.5,A5)')'I_eff_alt = ',I_eff_alt,' (eV)' write(*,'(A18,F10.5)')'A_eff_d = ',A_eff_d write(*,'(A18,F10.5)')'Z_eff_d = ',Z_eff_d write(*,'(A18,F10.5)')'Z_eff_d/A_eff_d = ',Z_eff_d/A_eff_d C N_Si = N_1*N1_Si N_Ca = N_1*N1_Ca N_Al = N_1*N1_Al N_Mg = N_1*N1_Mg N_B = N_1*N1_B C N_H = N_2*N2_H N_C = N_2*N2_C N_O = N_1*N1_O + N_2*N2_O N_N = N_2*N2_N C write(*,'(A,F10.6)')'Atomic fraction (H)= ',N_H write(*,'(A,F10.6)')'Atomic fraction (C)= ',N_C write(*,'(A,F10.6)')'Atomic fraction (O)= ',N_O write(*,'(A,F10.6)')'Atomic fraction (N)= ',N_N C write(*,'(A,F10.6)')'Atomic fraction (Si)= ',N_Si write(*,'(A,F10.6)')'Atomic fraction (Ca)= ',N_Ca write(*,'(A,F10.6)')'Atomic fraction (Al)= ',N_Al write(*,'(A,F10.6)')'Atomic fraction (Mg)= ',N_Mg write(*,'(A,F10.6)')'Atomic fraction (B)= ',N_B C A_eff = N_O *A_O & + N_Si*A_Si + N_Ca*A_Ca + N_Al*A_Al & + N_Mg*A_Mg + N_B *A_B & + N_H *A_H + N_C *A_C + N_N *A_N Z_eff = N_O *Z_O & + N_Si*Z_Si + N_Ca*Z_Ca + N_Al*Z_Al & + N_Mg*Z_Mg + N_B *Z_B & + N_H *Z_H + N_C *Z_C + N_N *Z_N ln_I_eff = (N_O *Z_O *log(I_O) + N_Si*Z_Si*log(I_Si) & + N_Ca*Z_Ca*log(I_Ca) + N_Al*Z_Al*log(I_Al) & + N_Mg*Z_Mg*log(I_Mg) + N_B *Z_B *log(I_B) & + N_H *Z_H *log(I_H) + N_C *Z_C *log(I_C) & + N_N *Z_N *log(I_N) & )/Z_eff I_eff = exp(ln_I_eff) write(*,'(A18,F10.5,A5)')'I_eff = ',I_eff,' (eV)' write(*,'(A18,F10.5)')'A_eff = ',A_eff write(*,'(A18,F10.5)')'Z_eff = ',Z_eff write(*,'(A18,F10.5)')'Z_eff/A_eff = ',Z_eff/A_eff C w_H = N_H *A_H /A_eff w_C = N_C *A_C /A_eff w_O = N_O *A_O /A_eff w_N = N_N *A_N /A_eff w_Si = N_Si*A_Si/A_eff w_Ca = N_Ca*A_Ca/A_eff w_Al = N_Al*A_Al/A_eff w_Mg = N_Mg*A_Mg/A_eff w_B = N_B *A_B /A_eff C write(*,'(A,F10.6)')'Mass fraction (H)= ',w_H write(*,'(A,F10.6)')'Mass fraction (C)= ',w_C write(*,'(A,F10.6)')'Mass fraction (O)= ',w_O write(*,'(A,F10.6)')'Mass fraction (N)= ',w_N C write(*,'(A,F10.6)')'Mass fraction (Si)= ',w_Si write(*,'(A,F10.6)')'Mass fraction (Ca)= ',w_Ca write(*,'(A,F10.6)')'Mass fraction (Al)= ',w_Al write(*,'(A,F10.6)')'Mass fraction (Mg)= ',w_Mg write(*,'(A,F10.6)')'Mass fraction (B)= ',w_B C recip_L_rad=w_O /L_O + w_Si/L_Si + w_Ca/L_Ca + w_Al/L_Al & + w_Mg/L_Mg + w_B /L_B & + w_H /L_H + w_C /L_C + w_N /L_N L_rad=1.D0/recip_L_rad write(*,'(A18,F10.5,A9)')'L_rad = ',L_rad,' (g/cm^2)' C Z_avr = Z_O *w_O + Z_Si*w_Si + Z_Ca*w_Ca + Z_Al*w_Al & + Z_Mg*w_Mg + Z_B *w_B & + Z_H *w_H + Z_C *w_C + Z_N *w_N A_avr = A_O *w_O + A_Si*w_Si + A_Ca*w_Ca + A_Al*w_Al & + A_Mg*w_Mg + A_B *w_B & + A_H *w_H + A_C *w_C + A_N *w_N A_avr_d = Z_avr/(Z_eff/A_eff) write(*,'(A18,F10.5)')' = ',A_avr write(*,'(A18,F10.5,A21)')' = ',A_avr_d, & ' <- /(Z_eff/A_eff)' write(*,'(A18,F10.5)')' = ',Z_avr end