The incorporation of radionuclides into low-temperature mineral hosts may strongly influence the concentration and migration of radioactive contaminants in the subsurface. One difficulty in evaluating the thermodynamics of incorporation is that experiments are often performed at high supersaturations and typically do not reach equilibrium. An alternative way to obtain the equilibrium thermodynamics is the quantum-mechanical analysis of the mineral host and the incorporated species before and after incorporation. In this contribution, density functional theory is used to calculate the energetics, resulting structures, and electronic configuration of uranyl (UO 2 2+) and neptunyl (NpO 2 +) incorporation into sulfate and carbonate minerals. In each host mineral, gypsum (CaSO 4 2H 2 O), anhydrite (CaSO 4), anglesite (PbSO 4), celestine (SrSO 4), barite (BaSO 4), calcite (CaCO 3), aragonite (CaCO 3), cerussite (PbCO 3), strontianite (SrCO 3), and witherite (BaCO 3), a divalent cation is replaced with either UO 2 2+ or NpO 2 + (in the case of neptunyl, charge balance is maintained with an additional hydrogen ion). The source of the actinyl ion and the sink for the host cation are modeled as both solid and aqueous phases, the latter of which requires an expansion of previous descriptions of incorporation. By combining periodic and cluster computational methods, this newly-developed approach enables the quantum-mechanical simulation of reactions between charged, aqueous molecular species and solid mineral phases. Among the host minerals considered, gypsum and aragonite are the most favorable hosts for both uranyl and neptunyl uptake (∆E ୷୮,ୟ୯