The frozen natural orbital ͑FNO͒ coupled-cluster method increases the speed of coupled-cluster ͑CC͒ calculations by an order of magnitude with no consequential error along a potential energy surface. This method allows the virtual space of a correlated calculation to be reduced by about half, significantly reducing the time spent performing the coupled-cluster ͑CC͒ calculation. This paper reports the derivation and implementation of analytical gradients for FNO-CC, including all orbital relaxation for both noncanonical and semicanonical perturbed orbitals. These derivatives introduce several new orbital relaxation contributions to the CC density matrices. FNO-CCSD͑T͒ and FNO-⌳CCSD͑T͒ are applied to a test set of equilibrium structures, verifying that these methods are capable of reproducing geometries and vibrational frequencies accurately, as well as energies. Several decomposition pathways of nitroethane are investigated using CCSD͑T͒ and ⌳CCSD͑T͒ with 60% of the FNO virtual orbitals in a cc-pVTZ basis, and find differences on the order of 5 kcal/ mol with reordering of the transition state energies when compared to B3LYP 6-311 +G͑3df ,2p͒.