  We denote the abstract simplicial complex containing all matchings in the
complete graph of order $n$ by $M_n$. Bouc first studied the topological
properties of $M_n$ in connection with the Quillen complexes. Later Bj\"orner
et al. showed that $M_n$ is homotopically $(\nu_n-1)$-connected, where
$\nu_n=\lfloor{\frac{n+1}{3}}\rfloor-1$, and conjectured that this connectivity
bound is sharp. Shareshian and Wachs settled the conjecture by inductively
showing that the $\nu_n$-dimensional homology group of $M_n$ is nontrivial,
with Bouc's "hand" calculation of $H_1(M_7)$ serving as the pivotal base step.
In general, the topology of $M_n$ is not very well-understood, even for a small
$n$. In the present article, we look into the topology of $M_n$, and $M_7$ in
particular, in the light of discrete Morse theory as developed by Forman. This
theory gives a combinatorial generalization of the classical notion of a smooth
gradient vector field, called (discrete) gradient vector field. Similar to the
smooth case, the critical simplices capture the topology of the complex. We
first construct a gradient vector field on $M_n$ that doesn't admit any
critical simplices up to dimension $\nu_n-1$, except one $0$-simplex, which
also leads to the aforementioned $(\nu_n-1)$-connectedness of $M_n$ in a purely
combinatorial way. However, for an efficient homology computation, we need to
work with a gradient vector field that admits a low number of critical
simplices. An optimal gradient vector field is one with the least number of
critical simplices, but the problem of finding one is an NP-hard problem. We
improve the gradient vector field constructed on $M_7$ in particular to a
near-optimal one, and then with the help of this improved gradient vector
field, compute the homology groups of $M_7$ in an efficient and algorithmic
manner. We also augment this near-optimal gradient vector field to an optimal
one.
