原子質量単位のメモ

原子質量単位の値メモ

なにかに使えれば。 使ったことのある値だけ書いてます。 追加情報・間違い等あればコメントください。

PROGRAM AMU
IMPLICIT NONE

REAL(8), PARAMETER :: c         = 299792458d0
REAL(8), PARAMETER :: e         = 1.602176634d-19
REAL(8), PARAMETER :: AMU_U     = 1.66053906660d-27
REAL(8), PARAMETER :: AMU_e     = 0.00054858d0
REAL(8), PARAMETER :: AMU_p     = 1.0073d0
REAL(8), PARAMETER :: AMU_n     = 1.008665d0
REAL(8), PARAMETER :: AMU_1H    = 1.007825d0
REAL(8), PARAMETER :: AMU_D     = 2.014103d0
REAL(8), PARAMETER :: AMU_T     = 3.0160492d0
REAL(8), PARAMETER :: AMU_3He   = 3.01603
REAL(8), PARAMETER :: AMU_4He   = 4.0026022d0
REAL(8), PARAMETER :: AMU_238U  = 238.0508d0
REAL(8), PARAMETER :: AMU_235U  = 235.04392d0
REAL(8), PARAMETER :: AMU_95Mo  = 94.945d0
REAL(8), PARAMETER :: AMU_139La = 138.995d0
REAL(8) :: M1, M2, SUM_M, SUM_MM, D_M, 
REAL(8) :: MM1, MM2, MM3, E1, E2, dE, E3

STOP
CONTAINS
! =============================================== !
SUBROUTINE CALC_E_22
    IMPLICIT NONE
    SUM_M = M1 + M2
    SUM_MM = MM1 + MM2
    D_M = SUM_M - SUM_MM
    dE = d_M * amu_U * c**2
    E1 = dE * MM2 / (MM1 + MM2)
    E2 = dE * MM1 / (MM1 + MM2)
END SUBROUTINE CALC_E_22
! =============================================== !
SUBROUTINE CALC_E_23
    IMPLICIT NONE
    SUM_M = M1 + M2             ; PRINT*, SUM_M
    SUM_MM = MM1 + MM2 + MM3    ; PRINT*, SUM_MM
    D_M = SUM_M - SUM_MM        ; PRINT*, D_M
    dE = d_M * amu_U * c**2             ; PRINT*, dE
    E1 = dE * (MM2+MM3) / (MM1 + MM2 + MM3)
    E2 = dE * MM1 / (MM1 + MM2 + MM3)
END SUBROUTINE CALC_E_23
! =============================================== !