HCN is a key ingredient for synthesizing biomolecules such as nucleobases and amino acids. We calculate 42 reaction rate coefficients directly involved with or in competition with the production of HCN in the early Earth or Titan atmospheres. These reactions are driven by methane and nitrogen radicals produced via UV photodissociation or lightning. For every reaction in this network, we calculate rate coefficients at 298 K using canonical variational transition state theory (CVT) paired with computational quantum chemistry simulations at the BHandHLYP/augcc-pVDZ level of theory. We also calculate the temperature dependence of the rate coefficients for the reactions that have barriers from 50-400 K. We present 15 new reaction rate coefficients with no previous known value. 93% of our calculated coefficients are within an order of magnitude of the nearest experimental or recommended values. Above 320 K, the rate coefficient for the new reaction H2CN − −→ HCN + H dominates. Contrary to experiments, we find the HCN reaction pathway, N + CH3 − −→ HCN + H2, to be inefficient, and suggest the experimental rate coefficient actually corresponds to an indirect pathway, through the H2CN intermediate. We present CVT using energies computed with density functional theory as a feasible and accurate method for calculating a large network of rate coefficients of small-molecule reactions. * Corresponding author:pearcbe@mcmaster.ca these sources has errors associated with it, and there are often a range of experimental and theoretical values to choose from for a single reaction. As a result, atmospheric HCN compositions can vary by orders of magnitude from one simulation to the next. Therefore, it is perhaps unsurprising that, as of yet, no simulation has matched the HCN profile of Titan completely.There are also several reactions without past experimental, theoretical, or suggested values that are missing in these networks that may play important roles in HCN formation (e.g. 1 CH 2 + 2 N − −→ H 2 CN and H 2 CN − −→ HCN + H).The focus of this work is to create a theoretical reduced HCN chemical network, where all the rate coefficients are consistently calculated with the same theoretical and computational method. Using this strategy, all reactions can be theoretically validated before being employed in a chemical network, and key reaction pathways with previously unknown rate coefficients can be included. Furthermore, by constructing a model chemistry 19,20 the errors for consistently calculated rate coefficients are expected to be similar, thus employing such a network has a chance to improve accuracy.The limitation of calculating a consistent theoretical network is that one cannot feasibly include a large number of molecular species. For every additional species, there is a potential additional reaction with all the existing species in the network. Therefore in this work, we focus only on the small set of reactions involved in the production of HCN from methane and nitrogen dissociation radicals, as well as the direct competing ...