Based on the Nambu–Jona–Lasinio (NJL) model, we develop a framework for calculating the spin alignment of vector mesons and applied it to study ϕ\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\phi $$\\end{document} mesons in a magnetic field. We calculate mass spectra for ϕ\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\phi $$\\end{document} mesons and observe mass splitting between the longitudinally polarized state and transversely polarized states. The ϕ\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\phi $$\\end{document} meson in a thermal equilibrium system is preferred to occupy the state with spin λ=0\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\lambda =0$$\\end{document} than those with spin λ=±1\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\lambda =\\pm 1$$\\end{document}, because the former state has a smaller energy. As a consequence, we conclude that the spin alignment will be larger than 1/3 if one measures along the direction of the magnetic field, which is qualitatively consistent with the recent STAR data. Around the critical temperature TC=150MeV\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$T_{C}=150~\\hbox {MeV}$$\\end{document}, the positive deviation from 1/3 is proportional to the square of the magnetic field strength, which agrees with the result from the non-relativistic coalescence model. Including the anomalous magnetic moments for quarks will modify the dynamical masses of quarks and thus affect the mass spectra and spin alignment of ϕ\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\phi $$\\end{document} mesons. The discussion of spin alignment in the NJL model may help us better understand the formation of hadron’s spin structure during the chiral phase transition.