This comment is in response to Nieves and Pal [1] who dispute our claim [2] that a classical gravitational field may possibly distinguish between Dirac and Majorana neutrinos, described in terms of Gaussian wave packets propagating in a Lense-Thirring background, where the distinction is manifested in spin-gravity corrections to the neutrino oscillation length. They contend that our model for Majorana neutrinos is incorrect and that any conclusions relating to this neutrino type are not reliable. Furthermore, they suggest that any distinction between the two neutrino types will be suppressed by factors of m/E where m is the neutrino mass and E is its mean energy of propagation, and claim that this distinction is unobservable when m/E ≪ 1. We beg to disagree.As noted above, our model for the Dirac and Majorana neutrinos is motivated by a wave packet approach in quantum mechanics, as opposed to a plane-wave expansion in quantum field theory. This is an essential detail which Nieves and Pal fail to acknowledge. In addition, the gravitational field [2] is incorporated in terms of a gravitational phase Φ G , giving rise to an interaction Hamiltonian H ΦG with spin-dependent features, to be evaluated in terms of time-independent perturbation theory. These two details are important for framing the context underpinning our reply. Regarding the technical concerns, we agree that the Majorana condition they note in their eq. (1) is certainly true for a fermion field operator. However, we again emphasize that our perspective is quantum mechanical, so our treatment of the Majorana condition must be described in terms of wave functions. Adopting their notation, our approach is to identify [3,4] Unlike what Nieves and Pal claim, it indeed follows [3] that W 1,2 is a solution of the free particle equation / k W 1,2 = ± m W 1,2 [5] for wave functions. However, a more precise implementation of (1) with our notation leads to a Majorana wave packet model with the formwhereand ξ(k) is the Gaussian function [2]. Clearly, it follows from (2) that |ψ c 1(2) Maj. = ± |ψ 1(2) Maj. . This leads to a modification of our eq. (13) for the Majorana matrix element [2], which isOur plots [2], as applied to the SN 1987A data, are completely unaffected by the adjustments because the corrections only alter the contributions coupled to M ΩR 2 /r 2 , which are all exponentially damped compared to the M/r contributions. As for the m/E suppression issue raised by Nieves and Pal, this is of no relevance because all such terms are automatically excluded within the construction of |W 1,2 (k) Maj. , so everything presented in [2] is of leading order. As shown in (5), and also present in our eq. (13)